The R2R3-MYB transcription factor ZeMYB32 negatively regulates anthocyanin biosynthesis in Zinnia elegans

KEY MESSAGE

This study identified an R2R3-MYB from Zinnia elegans, ZeMYB32, which negatively regulates anthocyanin biosynthesis.

Molecular Mechanism of Resistance to Alternaria alternata Apple Pathotype in Apple by Alternative Splicing of Transcription Factor MdMYB6-like

As a fruit tree with great economic value, apple is widely cultivated in China. However, apple leaf spot disease causes significant damage to apple quality and economic value. In our study, we found that MdMYB6-like is a transcription factor without auto-activation activity and with three alternative spliced variants. Among them, MdMYB6-like-β responded positively to the pathogen infection. Overexpression of MdMYB6-like-β increased the lignin content of leaves and improved the pathogenic resistance of apple flesh callus. In addition, all three alternative spliced variants of MdMYB6-like could bind to the promoter of MdBGLU H. Therefore, we believe that MdMYB6-like plays an important role in the infection process of the pathogen and lays a solid foundation for breeding disease-resistant cultivars of apple in the future.

Identification and Expression Analysis of R2R3-MYB Transcription Factors Associated with Flavonoid Biosynthesis in Panax quinquefolius

Panax quinquefolius L. is an important medicinal plant, and flavonoids are among its main secondary metabolites. The R2R3-MYB transcription factor plays an irreplaceable role in plant growth, development, and secondary metabolism. In our study, we identified 159 R2R3-MYBs and analyzed their physical and chemical properties in P. quinquefolius. The protein length of 159 PqMYBs varied from 107 to 1050 amino acids. The molecular weight ranged from 12.21 to 116.44 kDa. The isoelectric point was between 4.57 and 10.34. We constructed a phylogenetic tree of P. quinquefolius and Arabidopsis thaliana R2R3-MYB family members, and PqMYB members were divided into 33 subgroups. Transcriptome data analysis showed that the expression patterns of PqMYBs in root, leaf, and flower were significantly different. Following the MeJA treatment of seedlings, five candidate PqMYB genes demonstrated a response. A correlation analysis of PqMYBs and candidate flavonoid pathway genes showed that PqMYB2, PqMYB46, and PqMYB72 had correlation coefficients that were higher than 0.8 with PqCHS, PqANS4, and PqCCoAMT10, respectively. Furthermore, a transient expression assay confirmed that the three PqMYBs were localized in the nucleus. We speculated that these three PqMYBs were related to flavonoid biosynthesis in P. quinquefolius. These results provided a theoretical basis and a new perspective for further understanding the R2R3-MYB gene family and the biosynthesis mechanism of secondary metabolites in P. quinquefolius.

Absence of long-term balancing selection on variation in EuMYB3, an R2R3-MYB gene responsible for the anther-color polymorphism in Erythronium umbilicatum

Balancing selection has been shown to be common in plants for several different types of traits, such as self-incompatibility and heterostyly. Generally, for these traits balancing selection is generated by interactions among individuals or between individuals and other species (e.g., pathogens or pollinators). However, there are phenotypic polymorphisms in plants that do not obviously involve types of interactions that generate balancing selection. Little is known about the extent to which balancing selection also acts to preserve these polymorphisms. Here we ask whether balancing selection preserves an anther-color polymorphism in Erythronium umbilicatum (Liliaceae). We identified a major gene underlying this polymorphism. We then attempted to detect signatures of balancing selection on that gene by developing a new coalescence test for balancing selection. We found that variation in anther color is in large part caused by variation in a paralog of EuMYB3, an anthocyanin-regulating R2R3-MYB transcription factor. However, we found little evidence for balancing selection having acted historically on EuMYB3. Our results thus suggest that plant polymorphisms, especially those not involved in interactions that are likely to generate negative frequency-dependent selection, may reflect a transient state in which one morph will eventually be fixed by either genetic drift or directional selection. Our results also suggest that regulation of the anthocyanin pathway is more evolutionarily labile than is generally believed.

Transcriptome-Based Identification of the SaR2R3-MYB Gene Family in Sophora alopecuroides and Function Analysis of SaR2R3-MYB15 in Salt Stress Tolerance

As one of the most prominent gene families, R2R3-MYB transcription factors significantly regulate biochemical and physiological processes under salt stress. However, in Sophora alopecuroides, a perennial herb known for its exceptional saline alkali resistance, the comprehensive identification and characterization of SaR2R3-MYB genes and their potential functions in response to salt stress have yet to be determined. We investigated the expression profiles and biological functions of SaR2R3-MYB transcription factors in response to salt stress, utilizing a transcriptome-wide mining method. Our analysis identified 28 SaR2R3-MYB transcription factors, all sharing a highly conserved R2R3 domain, which were further divided into 28 subgroups through phylogenetic analysis. Some SaR2R3-MYB transcription factors showed induction under salt stress, with SaR2R3-MYB15 emerging as a potential regulator based on analysis of the protein-protein interaction network. Validation revealed the transcriptional activity and nuclear localization of SaR2R3-MYB15. Remarkably, overexpression of SaR2R3-MYB15 in transgenic plants could increase the activity of antioxidant enzymes and the accumulation of proline but decrease the content of malondialdehyde (MDA), compared with wild-type plants. Moreover, several salt stress-related genes showed higher expression levels in transgenic plants, implying their potential to enhance salt tolerance. Our findings shed light on the role of SaR2R3-MYB genes in salt tolerance in S. alopecuroides.

Comparative transcriptome analysis and identification of candidate R2R3-MYB genes involved in anthraquinone biosynthesis in Rheum palmatum L

BACKGROUND

Rheum palmatum L. has important medicinal value because it contains biologically active anthraquinones. However, the key genes and TFs involved in anthraquinone biosynthesis and regulation in R. palmatum remain unclear.

METHODS

Based on full length transcriptome data, in this study, we screened the differentially expressed genes in the anthraquinone biosynthesis pathway. The R2R3-MYB family genes of R. palmatum were systematically identified based on full-length transcriptome sequencing followed by bioinformatics analyses. The correlation analysis was carried out by using co-expression analysis, protein interaction analysis, and real-time fluorescence quantitative analysis after MeJA treatment. The RpMYB81 and RpMYB98 genes were amplified by RT-PCR, and their subcellular localization and self-activation characteristics were analyzed.

RESULTS

Comparative transcriptome analysis results revealed a total of 3525 upregulated and 6043 downregulated DEGs in the CK versus MeJA group; 28 DEGs were involved in the anthraquinone pathway. Eleven CHS genes that belonged to the PKS family were differentially expressed and involved in anthraquinone biosynthesis. Twelve differentially expressed MYBs genes were found to be co-expressed and interact with CHS genes. Furthermore, 52 MYB genes were identified as positive regulators of anthraquinone biosynthesis and were further characterized. Three MYB genes including RpMYB81, RpMYB98, and RpMYB100 responded to MeJA treatment in R. palmatum, and the levels of these genes were verified by qRT-PCR. RpMYB81 was related to anthraquinone biosynthesis. RpMYB98 had an interaction with genes in the anthraquinone biosynthesis pathway. RpMYB81 and RpMYB98 were mainly localized in the nucleus. RpMYB81 had self-activation activity, while RpMYB98 had no self-activation activity.

CONCLUSION

RpMYB81, RpMYB98, and RpMYB100 were significantly induced by MeJA treatment. RpMYB81 and RpMYB98 are located in the nucleus, and RpMYB81 has transcriptional activity, suggesting that it might be involved in the transcriptional regulation of anthraquinone biosynthesis in R. palmatum.

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

Genetic basis of nectar guide trichome variation between bumblebee- and self-pollinated monkeyflowers (Mimulus): role of the MIXTA-like gene GUIDELESS

Nectar guide trichomes play crucial ecological roles in bee-pollinated flowers, as they serve as footholds and guides for foraging bees to access the floral rewards. However, the genetic basis of natural variation in nectar guide trichomes among species remains poorly understood. In this study, we performed genetic analysis of nectar guide trichome variation between two closely related monkeyflower (Mimulus) species, the bumblebee-pollinated Mimulus lewisii and self-pollinated M. parishii. We demonstrate that a MIXTA-like R2R3-MYB gene, GUIDELESS, is a major contributor to the nectar guide trichome length variation between the two species. The short-haired M. parishii carries a recessive allele due to non-synonymous substitutions in a highly conserved motif among MIXTA-like MYB proteins. Furthermore, our results suggest that besides GUIDELESS, additional loci encoding repressors of trichome elongation also contribute to the transition from bumblebee-pollination to selfing. Taken together, these results suggest that during a pollination syndrome switch, changes in seemingly complex traits such as nectar guide trichomes could have a relatively simple genetic basis, involving just a few genes of large effects.

The R2R3-MYB Transcriptional Repressor TgMYB4 Negatively Regulates Anthocyanin Biosynthesis in Tulips (Tulipa gesneriana L.)

Anthocyanins play a paramount role in color variation and significantly contribute to the economic value of ornamental plants. The conserved activation complex MYB-bHLH-WD40 (MBW; MYB: v-myb avian myeloblastosis viral oncogene homolog; bHLH: basic helix-loop-helix protein; WD40:WD-repeat protein) involved in anthocyanin biosynthesis has been thoroughly researched, but there have been limited investigations into the function of repressor factors. In this study, we characterized TgMYB4, an R2R3-MYB transcriptional repressor which is highly expressed during petal coloration in red petal cultivars. TgMYB4-overexpressing tobaccos exhibited white or light pink petals with less anthocyanin accumulation compared to control plants. TgMYB4 was found to inhibit the transcription of ANTHOCYANIDIN SYNTHASE (TfANS1) and DIHYDRO-FLAVONOL-4-REDUCTASE (AtDFR), although it did not bind to their promoters. Moreover, the TgMYB4 protein was able to compete with the MYB activator to bind to the :bHLHprotein, thereby suppressing the function of the activator MBW complex. These findings demonstrate that TgMYB4 plays a suppressive role in the regulation of anthocyanin synthesis during flower pigmentation.

Sugar Signals and R2R3-MYBs Participate in Potassium-Repressed Anthocyanin Accumulation in Radish

Anthocyanin biosynthesis in plants is influenced by a wide range of environmental factors, such as light, temperature and nutrient availability. In this study, we revealed that the potassium-repressed anthocyanin accumulation in radish hypocotyls was associated with altered sugar distribution and sugar signaling pathways rather than changes in oxidative stress status. Sugar-feeding experiments suggested a hexokinase-independent glucose signal acted as a major contributor in regulating anthocyanin biosynthesis, transport and regulatory genes at the transcriptional level. Several R2R3-MYBs were identified as anthocyanin-related MYBs. Phylogenetic and protein sequence analyses suggested that RsMYB75 met the criteria of subgroup 6 MYB activator, while RsMYB39 and RsMYB82 seemed to be a non-canonical MYB anthocyanin activator and repressor, respectively. Through yeast-one-hybrid, dual-luciferase and transient expression assays, we confirmed that RsMYB39 strongly induced the promoter activity of anthocyanin transport-related gene RsGSTF12, while RsMYB82 significantly reduced anthocyanin biosynthesis gene RsANS1 expression. Molecular models are proposed in the discussion, allowing speculation on how these novel RsMYBs may regulate the expression levels of anthocyanin-related structural genes. Together, our data evidenced the strong impacts of potassium on sugar metabolism and signaling and its regulation of anthocyanin accumulation through different sugar signals and R2R3-MYBs in a hierarchical regulatory system.

The Characterization of R2R3-MYB Genes in Ammopiptanthus nanus Uncovers That the miR858-AnaMYB87 Module Mediates the Accumulation of Anthocyanin under Osmotic Stress

R2R3-MYB transcription factors (TFs) participate in the modulation of plant development, secondary metabolism, and responses to environmental stresses. Ammopiptanthus nanus, a leguminous dryland shrub, tolerates a high degree of environmental stress, including drought and low-temperature stress. The systematic identification, structural analysis, evolutionary analysis, and gene profiling of R2R3-MYB TFs under cold and osmotic stress in A. nanus were performed. Up to 137 R2R3-MYB TFs were identified and clustered into nine clades, with most A. nanus R2R3-MYB members belonging to clade VIII. Tandem and segmental duplication events drove the expansion of the A. nanus R2R3-MYB family. Expression profiling revealed that multiple R2R3-MYB genes significantly changed under osmotic and cold stress conditions. MiR858 and miR159 targeted 88 R2R3-MYB genes. AnaMYB87, an miR858-targeted clade VIII R2R3-MYB TF, was up-regulated under both osmotic and cold stress. A transient expression assay in apples showed that the overexpression of AnaMYB87 promoted anthocyanin accumulation. A luciferase reporter assay in tobacco demonstrated that AnaMYB87 positively affected the transactivation of the dihydroflavonol reductase gene, indicating that the miR858-MYB87 module mediates anthocyanin accumulation under osmotic stress by regulating the dihydroflavonol reductase gene in A. nanus. This study provides new data to understand the roles of R2R3-MYB in plant stress responses.

Identification of HpMYB1 inducing anthocyanin accumulation in Hippeastrum Hybridum tepals by RNA-seq

BACKGROUND

Cultivated Hippeastrum × hybridum is a popular ornamental plant with large and colorful flowers, long flowering duration, and high commercial value. As its main ornamental feature, its flower color is related to the anthocyanin content in the tepals. However, the molecular regulatory mechanisms of anthocyanin biosynthesis in H. × hybridum have not yet been elucidated.

RESULTS

In the present study, 12 cDNA libraries of four stages of H.× hybridum 'Royal Velvet' tepal development were used for RNA-seq, obtaining 79.83 gigabases (GB) of clean data. The data were assembled into 148,453 unigenes, and 11,262 differentially expressed genes were identified. Forty key enzymes participating in anthocyanin biosynthesis were investigated, and the results showed that most of the anthocyanin structural genes were expressed at low levels in S1 and were markedly upregulated in S2 and S3. The expression profiles of 12 selected genes were verified by qRT-PCR. Furthermore, the R2R3-MYB transcription factor (TF), HpMYB1, involved in the regulation of anthocyanin biosynthesis was identified by sequence, expression pattern, and subcellular localization analyses. Its overexpression in tobacco significantly increased the anthocyanin levels in various tissues and activated anthocyanin-related genes.

CONCLUSIONS

Using RNA-seq technology, we successfully identified a potential R2R3-MYB gene, HpMYB1, that regulates anthocyanin biosynthesis in H.× hybridum 'Royal Velvet'. Our findings provide basic transcript information and valuable transcriptome data for further identification of key genes involved in anthocyanin biosynthesis and can be applied in the artificial breeding of new H. × hybridum cultivars with enhanced ornamental value.

Ethylene inhibits ABA-induced stomatal closure via regulating NtMYB184-mediated flavonol biosynthesis in tobacco

Stomatal movement can be regulated by ABA signaling through synthesis of reactive oxygen species (ROS) in guard cells. By contrast, ethylene triggers the biosynthesis of antioxidant flavonols to suppress ROS accumulation and prevent ABA-induced stomatal closure; however, the underlying mechanism remains largely unknown. In this study, we isolated and characterized the tobacco (Nicotiana tabacum) R2R3-MYB transcription factor NtMYB184, which belongs to the flavonol-specific SG7 subgroup. RNAi suppression and CRISPR/Cas9 mutation (myb184) of NtMYB184 in tobacco caused down-regulation of flavonol biosynthetic genes and decreased the concentration of flavonols in the leaves. Yeast one-hybrid assays, transactivation assays, EMSAs, and ChIP-qPCR demonstrated that NtMYB184 specifically binds to the promoters of flavonol biosynthetic genes via MYBPLANT motifs. NtMYB184 regulated flavonol biosynthesis in guard cells to modulate ROS homeostasis and stomatal aperture. ABA-induced ROS production was accompanied by the suppression of NtMYB184 and flavonol biosynthesis, which may accelerate ABA-induced stomatal closure. Furthermore, ethylene stimulated NtMYB184 expression and flavonol biosynthesis to suppress ROS accumulation and curb ABA-induced stomatal closure. In myb184, however, neither the flavonol and ROS concentrations nor the stomatal aperture varied between the ABA and ABA+ethylene treatments, indicating that NtMYB184 was indispensable for the antagonism between ethylene and ABA via regulating flavonol and ROS concentrations in the guard cells.

Association Analysis of Transcriptome and Targeted Metabolites Identifies Key Genes Involved in Iris germanica Anthocyanin Biosynthesis

The anthocyanin biosynthetic pathway is the main pathway regulating floral coloration in Iris germanica, a well-known ornamental plant. We investigated the transcriptome profiles and targeted metabolites to elucidate the relationship between genes and metabolites in anthocyanin biosynthesis in the bitone flower cultivar 'Clarence', which has a deep blue outer perianth and nearly white inner perianth. In this study, delphinidin-, pelargonidin-, and cyanidin-based anthocyanins were detected in the flowers. The content of delphinidin-based anthocyanins increased with the development of the flower. At full bloom (stage 3), delphinidin-based anthocyanins accounted for most of the total anthocyanin metabolites, whereas the content of pelargonidin- and cyanidin-based anthocyanins was relatively low. Based on functional annotations, a number of novel genes in the anthocyanin pathway were identified, which included early biosynthetic genes IgCHS, IgCHI, and IgF3H and late biosynthetic genes Ig F3'5'H, IgANS, and IgDFR. The expression of key structural genes encoding enzymes, such as IgF3H, Ig F3'5'H, IgANS, and IgDFR, was significantly upregulated in the outer perianth compared to the inner perianth. In addition, most structural genes exhibited their highest expression at the half-color stage rather than at the full-bloom stage, which indicates that these genes function ahead of anthocyanins synthesis. Moreover, transcription factors (TFs) of plant R2R3-myeloblastosis (R2R3-MYB) related to the regulation of anthocyanin biosynthesis were identified. Among 56 R2R3-MYB genes, 2 members belonged to subgroup 4, with them regulating the expression of late biosynthetic genes in the anthocyanin biosynthetic pathway, and 4 members belonged to subgroup 7, with them regulating the expression of early biosynthetic genes in the anthocyanin biosynthetic pathway. Quantitative real-time PCR (qRT-PCR) analysis was used to validate the data of RNA sequencing (RNA-Seq). The relative expression profiles of most candidate genes were consistent with the FPKM of RNA-seq. This study identified the key structural genes encoding enzymes and TFs that affect anthocyanin biosynthesis, which provides a basis and reference for the regulation of plant anthocyanin biosynthesis in I. germanica.

Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

Light-mediated anthocyanin biosynthesis in rose petals involves a balanced regulatory module comprising transcription factors RhHY5, RhMYB114a, and RhMYB3b

Roses are significant botanical species with both ornamental and economic value, displaying diverse floral traits, particularly an extensive array of petal colors. The red pigmentation of rose petals is predominantly attributed to anthocyanin accumulation. However, the underlying regulatory mechanism of anthocyanin biosynthesis in roses remains elusive. This study presents a novel light-responsive regulatory module governing anthocyanin biosynthesis in rose petals, which involves the transcription factors RhHY5, RhMYB114a, and RhMYB3b. Under light conditions (1000-1500 μmol m-2 s-1), RhHY5 represses RhMYB3b expression and induces RhMYB114a expression, positively regulating anthocyanin biosynthesis in rose petals. Notably, activation of anthocyanin structural genes probably involves an interaction and synergy between RhHY5 and the MYB114a-bHLH3-WD40 complex. Additionally, RhMYB3b is activated by RhMYB114a to prevent excessive accumulation of anthocyanin. Conversely, under low light conditions (<10 μmol m-2 s-1), the degradation of RhHY5 leads to down-regulation of RhMYB114a and up-regulation of RhMYB3b, which in turn inhibits the expression of both RhMYB114a and anthocyanin structural genes. Additionally, RhMYB3b competes with RhMYB114a for binding to RhbHLH3 and the promoters of anthocyanin-related structural genes. Overall, our study uncovers a complex light-mediated regulatory network that governs anthocyanin biosynthesis in rose petals, providing new insights into the molecular mechanisms underlying petal color formation in rose.

Genome-wide identification of the longan R2R3-MYB gene family and its role in primary and lateral root

R2R3-MYB is an important transcription factor family that regulates plant growth and development. Root development directly affects the absorption of water and nutrients by plants. Therefore, to understand the regulatory role of R2R3-MYB transcription factor family in root development of longan, this study identified the R2R3-MYB gene family members at the genome-wide level, and analyzed their phylogenetic characteristics, physical and chemical properties, gene structure, chromosome location and tissue expression. The analysis identified 124 R2R3-MYB family members in the longan genome. Phylogenetic analysis divided these members into 22 subfamilies, and the members of the unified subfamily had similar motifs and gene structures. The result of qRT-PCR showed that expression levels of DlMYB33, DlMYB34, DlMYB59, and DlMYB77 were significantly higher in main roots than in lateral as opposed to those of DlMYB35, DlMYB69, DlMYB70, and DlMYB83, which were significantly lower. SapBase database prediction and miRNAs sequencing results showed that 34 longan miRNAs could cleave R2R3-MYB, including 17 novel miRNAs unique to longan. The qRT-PCR and subcellular localization experiments of DlMYB92 and DlMYB98 showed that DlMYB92 is a key factor that regulates transcription in the nucleus and participates in the regulation of longan lateral root development. Longan also has a conserved miRNA-MYB-lateral root development regulation mechanism. This study provides a reference for further research on the transcriptional regulation of the miRNA-R2R3-MYB module in the root development of longan.

Light-Induced TaHY5-7A and TaBBX-3B Physically Interact to Promote PURPLE PERICARP-MYB 1 Expression in Purple-Grained Wheat

Purple-grained wheat (Triticum aestivum L.) is an important germplasm source in crop breeding. Anthocyanin biosynthesis in the pericarps of purple-grained wheat is largely light-dependent; however, the regulatory mechanisms underlying light-induced anthocyanin accumulation in the wheat pericarp remain unknown. Here we determined that anthocyanins rapidly accumulate in the pericarps of the purple-grained wheat cultivar Heixiaomai 76 (H76) at 16 days after pollination under light treatment. Using transcriptome sequencing, differential gene expression analysis, and phylogenetic analysis, we identified two key genes involved in light signaling in wheat: ELONGATED HYPOCOTYL 5-7A (TaHY5-7A) and B-BOX-3B (TaBBX-3B). TaHY5-7A and TaBBX-3B were highly expressed in purple-grained wheat pericarps. The heterologous expression of TaHY5-7A partially restored the phenotype of the Arabidopsis (Arabidopsis thaliana) hy5 mutant, resulting in increased anthocyanin accumulation and a shortened hypocotyl. The heterologous expression of TaBBX-3B in wild-type Arabidopsis had similar effects. TaHY5-7A and TaBBX-3B were nucleus-localized, consistent with a function in transcription regulation. However, TaHY5-7A, which lacks a transactivation domain, was not sufficient to activate the expression of PURPLE PERICARP-MYB 1 (TaPpm1), the key anthocyanin biosynthesis regulator in purple pericarps of wheat. TaHY5-7A physically interacted with TaBBX-3B in yeast two-hybrid and bimolecular fluorescence complementation assays. Additionally, TaHY5-7A, together with TaBBX-3B, greatly enhanced the promoter activity of TaPpm1 in a dual luciferase assay. Overall, our results suggest that TaHY5-7A and TaBBX-3B collaboratively activate TaPpm1 expression to promote light-induced anthocyanin biosynthesis in purple-pericarp wheat.

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis

BACKGROUND

The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis.

RESULTS

In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR.

CONCLUSIONS

This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

Function of the R2R3-MYB Transcription Factors in Dalbergia odorifera and Their Relationship with Heartwood Formation

R2R3-MYB transcription factors (TFs) form one of the most important TF families involved in regulating various physiological functions in plants. The heartwood of Dalbergia odorifera is a kind of high-grade mahogany and valuable herbal medicine with wide application. However, the role of R2R3-MYB genes in the growth and development of D. odorifera, especially their relevance to heartwood formation, has not been revealed. A total of 126 R2R3-MYBs were screened from the D. odorifera genome and named DodMYB1-126 based on their location on 10 chromosomes. The collinearity results showed that purification selection was the main driving force for the evolution of the R2R3-MYB TFs family, and whole genome/fragment replication event was the main form for expanding the R2R3-MYB family, generating a divergence of gene structure and function. Comparative phylogenetic analysis classified the R2R3-MYB TFs into 33 subfamilies. S3-7,10,12-13,21 and N4-7 were extensively involved in the metabolic process; S9,13,16-19,24-25 and N1-3,8 were associated with the growth and development of D. odorifera. Based on the differential transcriptional expression levels of R2R3-MYBs in different tissues, DodMYB32, DodMYB55, and DodMYB89 were tentatively screened for involvement in the regulatory process of heartwood. Further studies have shown that the DodMYB89, localized in the nucleus, has transcriptional activation activity and is involved in regulating the biosynthesis of the secondary metabolites of heartwood by activating the promoters of the structural genes DodI2'H and DodCOMT. This study aimed to comprehensively analyze the functions of the R2R3-MYB TFs and screen for candidate genes that might be involved in heartwood formation of D. odorifera.

Two R2R3-MYB transcription factors from Chinese cedar (Cryptomeria fortunei Hooibrenk) are involved in the regulation of secondary cell wall formation

As the most abundant renewable energy source, wood comprises the secondary cell wall (SCW). SCW biosynthesis involves lignin and cellulose deposition. Increasing studies have illustrated that R2R3-MYB transcription factors (TFs) play pivotal roles in affecting lignin accumulation and SCW formation. Nevertheless, the regulatory roles of R2R3-MYBs are still unresolved in Cryptomeria fortunei Hooibrenk cambium and wood formation. To dissect the potentials of CfMYBs, we successfully cloned and intensively studied the functions of CfMYB4 and CfMYB5 in SCW formation and abiotic stress response. They both contained the conserved MYB domain capable of forming a special structure that could bind to the core motifs of downstream genes. The phylogenetic tree implied that two CfMYBs clustered into different evolutionary branches. They were predominantly expressed in the stem and were localized to the nucleus. Furthermore, CfMYB4 functioned as an activator to enhance lignin and cellulose accumulation, and increase the SCW thickness by elevating the expression levels of SCW-related genes. By contrast, CfMYB5 negatively regulated lignin and cellulose biosynthesis, and decreased SCW formation by reducing the expression of SCW biosynthetic genes. Our data not only highlight the regulatory functions of CfMYBs in lignin deposition but also provide critical insights into the development of strategies for the genetic improvement of Cryptomeria fortunei wood biomass.

Three R2R3-MYB transcription factors from banana (Musa acuminata) activate structural anthocyanin biosynthesis genes as part of an MBW complex

OBJECTIVE

Bananas are one of the most popular fruits in the world, providing food security and employment opportunities in several developing countries. Increasing the anthocyanin content of banana fruit could improve the health-promoting properties. Anthocyanin biosynthesis is largely regulated at the transcriptional level. However, relatively little is known about the transcriptional activation of anthocyanin biosynthesis in banana.

RESULTS

We analysed the regulatory activity of three Musa acuminata MYBs that were predicted by bioinformatic analysis to transcriptionally regulate anthocyanin biosynthesis in banana. MaMYBA1, MaMYBA2 and MaMYBPA2 did not complement the anthocyanin-deficient phenotype of the Arabidopsis thaliana pap1/pap2 mutant. However, co-transfection experiments in A. thaliana protoplasts showed that MaMYBA1, MaMYBA2 and MaMYBPA2 function as components of a transcription factor complex with a bHLH and WD40 protein, the so called MBW complex, resulting in the activation of the A. thaliana ANTHOCYANIDIN SYNTHASE and DIHYDROFLAVONOL 4-REDUCTASE promoters. The activation potential of MaMYBA1, MaMYBA2 and MaMYBPA2 was increased when combined with the monocot Zea mays bHLH ZmR instead of the dicot AtEGL3. This work paves the path towards decoding the MBW complex-mediated transcriptional activation of anthocyanin biosynthesis in banana. It will also facilitate research towards increased anthocyanin content in banana and other monocot crops.

Global dissection of R2R3-MYB in Pogostemon cablin uncovers a species-specific R2R3-MYB clade

MYB family is one of the largest transcription factor families in plants and plays a crucial role in regulating plant biochemical and physiological processes. However, R2R3-MYBs in patchouli have not been systematically investigated. Here, based on the gene annotation of patchouli genome sequence, 484 R2R3-MYB transcripts were detected. Further in-depth analysis of the gene structure and expression of R2R3-MYBs supported the tetraploid hybrid origin of patchouli. When combined with R2R3-MYBs from Arabidopsis, a phylogenetic tree of patchouli R2R3-MYBs was constructed and divided into 31 clades. Interestingly, a patchouli-specific R2R3-MYB clade was found and confirmed by homologous from other Lamiaceae species. The syntenic analysis demonstrated that tandem duplication contributed to its evolution. This study systematically analysed the R2R3-MYB family in patchouli, providing information on its gene characterization, functional prediction, and species evolution.

An R2R3-MYB network modulates stem strength by regulating lignin biosynthesis and secondary cell wall thickening in herbaceous peony

Stem strength is an important agronomic trait affecting plant lodging, and plays an essential role in the quality and yield of plants. Thickened secondary cell walls in stems provide mechanical strength that allows plants to stand upright, but the regulatory mechanism of secondary cell wall thickening and stem strength in cut flowers remains unclear. In this study, first, a total of 11 non-redundant Paeonia lactiflora R2R3-MYBs related to stem strength were identified and isolated from cut-flower herbaceous peony, among which PlMYB43, PlMYB83 and PlMYB103 were the most upregulated differentially expressed genes. Then, the expression characteristics revealed that these three R2R3-MYBs were specifically expressed in stems and acted as transcriptional activators. Next, biological function verification showed that these P. lactiflora R2R3-MYBs positively regulated stem strength, secondary cell wall thickness and lignin deposition. Furthermore, yeast-one-hybrid and dual luciferase reporter assays demonstrated that they could bind to the promoter of caffeic acid O-methyltransferase gene (PlCOMT2) and/or laccase gene (PlLAC4), two key genes involved in lignin biosynthesis. In addition, the function of PlLAC4 in increasing lignin deposition was confirmed by virus-induced gene silencing and overexpression. Moreover, PlMYB83 could also act as a transcriptional activator of PlMYB43. The findings of the study propose a regulatory network of R2R3-MYBs modulating lignin biosynthesis and secondary cell wall thickening for improving stem lodging resistance, and provide a resource for molecular genetic engineering breeding of cut flowers.

CgMYB1, an R2R3-MYB transcription factor, can alleviate abiotic stress in an annual halophyte Chenopodium glaucum

MYB transcription factors (TFs) are important regulators of the stress response in plants. In the present study, we characterized the CgMYB1 gene in Chenopodium glaucum, a member of the R2R3-MYB TF family. CgMYB1 was located in the nucleus with an activating domain at the C terminus. The CgMYB1 gene could be induced by salt and cold stress in C. glaucum. Overexpressing CgMYB1 in Arabidopsis significantly enhanced salt and cold tolerance, probably by improving physiological performance and stress-related gene expression. Further analysis suggests that the positive response of CgMYB1 to abiotic stress may partially be attributed to the interaction between CgMYB1 and the CgbHLH001 promoter followed by activation of downstream stress-responsive genes, which mediates stress tolerance. Our findings should contribute to further understanding of the function of R2R3 MYB TF in response to abiotic stress.

Mechanisms underlying the formation of complex color patterns on Nigella orientalis (Ranunculaceae) petals

Complex color patterns on petals are widespread in flowering plants, yet the mechanisms underlying their formation remain largely unclear. Here, by conducting detailed morphological, anatomical, biochemical, optical, transcriptomic, and functional studies, we investigated the cellular bases, chromogenic substances, reflectance spectra, developmental processes, and underlying mechanisms of complex color pattern formation on Nigella orientalis petals. We found that the complexity of the N. orientalis petals in color pattern is reflected at multiple levels, with the amount and arrangement of different pigmented cells being the key. We also found that biosynthesis of the chromogenic substances of different colors is sequential, so that one color/pattern is superimposed on another. Expression and functional studies further revealed that a pair of R2R3-MYB genes function cooperatively to specify the formation of the eyebrow-like horizontal stripe and the Mohawk haircut-like splatters. Specifically, while NiorMYB113-1 functions to draw a large splatter region, NiorMYB113-2 functions to suppress the production of anthocyanins from the region where a gap will form, thereby forming the highly specialized pattern. Our results provide a detailed portrait for the spatiotemporal dynamics of the coloration of N. orientalis petals and help better understand the mechanisms underlying complex color pattern formation in plants.

bracteatus

The R2R3-MYB proteins comprise the largest class of MYB transcription factors, which play an essential role in regulating anthocyanin synthesis in various plant species. Ananas comosus var. bracteatus is an important colorful anthocyanins-rich garden plant. The spatio-temporal accumulation of anthocyanins in chimeric leaves, bracts, flowers, and peels makes it an important plant with a long ornamental period and highly improves its commercial value. We conducted a comprehensive bioinformatic analysis of the R2R3-MYB gene family based on genome data from A. comosus var. bracteatus. Phylogenetic analysis, gene structure and motif analysis, gene duplication, collinearity, and promoter analysis were used to analyze the characteristics of this gene family. In this work, a total of 99 R2R3-MYB genes were identified and classified into 33 subfamilies according to phylogenetic analysis, and most of them were localized in the nucleus. We found these genes were mapped to 25 chromosomes. Gene structure and protein motifs were conserved among AbR2R3-MYB genes, especially within the same subfamily. Collinearity analysis revealed four pairs of tandem duplicated genes and 32 segmental duplicates in AbR2R3-MYB genes, indicating that segmental duplication contributed to the amplification of the AbR2R3-MYB gene family. A total of 273 ABRE responsiveness, 66 TCA elements, 97 CGTCA motifs, and TGACG motifs were the main cis elements in the promoter region under response to ABA, SA, and MEJA. These results revealed the potential function of AbR2R3-MYB genes in response to hormone stress. Ten R2R3-MYBs were found to have high homology to MYB proteins reported to be involved in anthocyanin biosynthesis from other plants. RT-qPCR results revealed the 10 AbR2R3-MYB genes showed tissue-specific expression patterns, six of them expressed the highest in the flower, two genes in the bract, and two genes in the leaf. These results suggested that these genes may be the candidates that regulate anthocyanin biosynthesis of A. comosus var. bracteatus in the flower, leaf, and bract, respectively. In addition, the expressions of these 10 AbR2R3-MYB genes were differentially induced by ABA, MEJA, and SA, implying that these genes may play crucial roles in hormone-induced anthocyanin biosynthesis. Our study provided a comprehensive and systematic analysis of AbR2R3-MYB genes and identified the AbR2R3-MYB genes regulating the spatial-temporal anthocyanin biosynthesis in A. comosus var. bracteatus, which would be valuable for further study on the anthocyanin regulation mechanism of A. comosus var. bracteatus.

Transcriptomic analysis reveals biosynthesis genes and transcription factors related to leaf anthocyanin biosynthesis in Aglaonema commutatum

BACKGROUND

Aglaonema commutatum 'Red Valentine', as a foliage ornamental plant, is widely used for interior and exterior decoration because of its easy cultivation and management. However, reduced proportion of red foliage during large-scale production of A. commutatum seedlings is a frequent occurrence, which has considerable implications on the plant's ornamental and market value. However, the molecular mechanisms underlying this phenomenon remain unclear.

RESULTS

To explore the molecular basis of the variation in leaf color of A. commutatum Red Valentine, we performed transcriptome sequencing with the Illumina platform using two different varieties of A. commutatum, namely Red Valentine and a green mutant, at three different stages of leaf development. We annotated 63,621 unigenes and 14,186 differentially expressed genes by pairwise comparison. Furthermore, we identified 26 anthocyanin biosynthesis structural genes. The transcript per million (TPM) values were significantly higher for Red Valentine than for the green mutant in all three developmental stages, consistent with the high anthocyanin content of Red Valentine leaves. We detected positive transcription factors that may be involved in the regulation of anthocyanin biosynthesis using BLAST and through correlation analysis. Downregulation of these transcription factors may downregulate the expression of anthocyanin genes. We obtained full-length cDNA of the anthocyanin biosynthesis and regulatory genes and constructed phylogenetic trees to ensure accuracy of the analysis.

CONCLUSIONS

Our study provides insights into the molecular mechanisms underlying leaf variation in A. commutatum Red Valentine and may be used to facilitate the breeding of ornamental cultivars with high anthocyanin levels.

An R2R3-MYB FtMYB11 from Tartary buckwheat has contrasting effects on abiotic tolerance in Arabidopsis

R2R3-MYB transcription factors play important roles in response to abiotic stresses in planta, such as salt, drought, and osmotic stress. However, the role of FtMYB11 in Tartary buckwheat (Fagopyrum tataricum) in drought and osmotic tolerance has not yet been elucidated. In this study, we found that FtMYB11 was markedly induced by exogenous abscisic acid (ABA), salinity, and mannitol. Further, FtMYB11-overexpressing Arabidopsis showed hypersensitivity to ABA-mediated seed germination and seedling establishment through regulating transcripts of AtCBF1, AtDREB2A, and AtRD20, compared with wild type, indicating that FtMYB11 plays a positive role in ABA signaling. In contrast, transgenic lines overexpressing FtMYB11 were sensitive to mannitol and NaCl treatments, suggesting that FtMYB11 plays a negative role in osmotic tolerance. Intriguingly, the transcripts of ABA biosynthetic enzyme genes were significantly elevated in plants overexpressing FtMYB11 after exposure to osmotic stresses, such as AtABA3 and AtNCED3. In addition, flavonoid biosynthesis genes were also upregulated in transgenic Arabidopsis under ABA, salt, and drought treatments, including AtC4H, AtF3H, AtANS, AtFLS, and At4CL. The drought tolerance assay showed that plants overexpressing FtMYB11 displayed greater tolerance to water deficit through regulating MDA and proline content. Taken together, FtMYB11 has opposite roles in response to abiotic stresses, but it may mediate flavonoid biosynthesis through regulation of related enzyme genes.

To stripe or not to stripe: the origin of a novel foliar pigmentation pattern in monkeyflowers (Mimulus)

The origin of phenotypic novelty is one of the most challenging problems in evolutionary biology. Although genetic regulatory network rewiring or co-option has been widely recognised as a major contributor, in most cases how such genetic rewiring/co-option happens is completely unknown. We have studied a novel foliar pigmentation pattern that evolved recently in the monkeyflower species Mimulus verbenaceus. Through genome-wide association tests using wild-collected samples, experimental crosses of laboratory inbred lines, gene expression analyses, and functional assays, we identified an anthocyanin-activating R2R3-MYB gene, STRIPY, as the causal gene triggering the emergence of the discrete, mediolateral anthocyanin stripe in the M. verbenaceus leaf. Chemical mutagenesis revealed the existence of upstream activators and repressors that form a 'hidden' prepattern along the leaf proximodistal axis, potentiating the unique expression pattern of STRIPY. Population genomics analyses did not reveal signatures of positive selection, indicating that nonadaptive processes may be responsible for the establishment of this novel trait in the wild. This study demonstrates that the origin of phenotypic novelty requires at least two separate phases, potentiation and actualisation. The foliar stripe pattern of M. verbenaceus provides an excellent platform to probe the molecular details of both processes in future studies.

'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

Interplay between R2R3 MYB-type activators and repressors regulates proanthocyanidin biosynthesis in banana (Musa acuminata)

Proanthocyanidins are oligomeric flavonoids that promote plant disease resistance and benefit human health. Banana is one of the world's most extensively farmed crops and its fruit pulp contain proanthocyanidins. However, the transcriptional regulatory network that fine tunes proanthocyanidin biosynthesis in banana remains poorly understood. We characterised two proanthocyanidin-specific R2R3 MYB activators (MaMYBPA1-MaMYBPA2) and four repressors (MaMYBPR1-MaMYBPR4) to elucidate the mechanisms underlying the transcriptional regulation of proanthocyanidin biosynthesis in banana. Heterologous expression of MaMYBPA1 and MaMYBPA2 partially complemented the Arabidopsis thaliana proanthocyanidin-deficient transparent testa2 mutant. MaMYBPA1 and MaMYBPA2 interacted physically with MaMYCs to transactivate anthocyanin synthase, leucoanthocyanidin reductase, and anthocyanidin reductase genes in vitro and form functional MYB-bHLH-WD Repeat (MBW) complexes with MaTTG1 to transactivate these promoters in vivo. Overexpression of MaMYBPAs alone or with MaMYC in banana fruits induced proanthocyanidin accumulation and transcription of proanthocyanidin biosynthesis-related genes. MaMYBPR repressors are also shown to interact with MaMYCs forming repressing MBW complexes, and diminished proanthocyanidin accumulation. Interestingly overexpression of MaMYBPA induces the expression of MaMYBPR, indicating an agile regulation of proanthocyanidin biosynthesis through the formation of competitive MBW complexes. Our results reveal regulatory modules of R2R3 MYB- that fine tune proanthocyanidin biosynthesis and offer possible targets for genetic manipulation for nutritional improvement of banana.

Comparative analysis of R2R3-MYB transcription factors in the flower of Iris laevigata identifies a novel gene regulating tobacco cold tolerance

Breeding for flower cold resistance is a priority for flower breeding research in northern China. The identification of cold resistance genes will not only provide genetic resources for cold resistance breeding, but also form a basis for the study of plant cold resistance mechanisms. Based on the flower transcriptome of Iris laevigata, 20 R2R3-MYBs were identified and comprehensive analysis, including conservative domain, phylogenetic analyses and functional distribution, were performed for R2R3-MYBs. Expression patterns of the abiotic stress genes under cold stress were detected, the upregulated gene was genetically transformed into tobacco, and the related physiological indicators of the transgenic tobacco were measured. A novel cold resistance gene, IlMYB306, was obtained. qRT-PCR indicated that IlMYB306 was dramatically induced by cold stress and was significantly upregulated in roots. The free proline content, MDA, SOD and POD activity of the transgenic tobacco improved after cold stress, and the chlorophyll content decreased slowly. In addition, overexpression of IlMYB306 improved cold resistance of the seeds. SEM results showed leaves of transgenic tobacco had obvious folds, more grooves and bulges on the lower leaf surface. Overall, we report a novel cold resistance R2R3-MYB gene, IlMYB306, in the flower of I. laevigata, which could improve tobacco cold stress tolerance by thickening the waxy layer, increasing antioxidant activity and the content of proline.

Seeing the unseen: a trifoliate (MYB117) mutant allele fortifies folate and carotenoids in tomato fruits

In tomato (Solanum lycopersicum), mutations in the gene encoding the R2R3-MYB117 transcription factor elicit trifoliate leaves and initiate the formation of axillary meristems; however, their effects on fruit ripening remain unexplored. The fruits of a new trifoliate (tf) mutant (tf-5) were firmer and had higher °Brix values and higher folate and carotenoid contents. The transcriptome, proteome, and metabolome profiling of tf-5 reflected a broad-spectrum change in cellular homeostasis. The tf-5 allele enhanced the fruit firmness by suppressing cell wall softening-related proteins. tf-5 fruit displayed a substantial increase in amino acids, particularly γ-aminobutyric acid, with a parallel reduction in aminoacyl-tRNA synthases. The increased lipoxygenase protein and transcript levels seemingly elevated jasmonic acid levels. In addition, increased abscisic acid hydrolase transcript levels coupled with reduced precursor supply lowered abscisic acid levels. The upregulation of carotenoids was mediated by modulation of methylerythreitol and plastoquinone pathways and increased the levels of carotenoid isomerization proteins. The upregulation of folate in tf-5 was connoted by the increase in the precursor p-aminobenzoic acid and transcript levels of several folate biosynthesis genes. The reduction in pterin-6-carboxylate levels and γ-glutamyl hydrolase activity indicated that reduced folate degradation in tf-5 increased folate levels. Our study delineates that in addition to leaf development, MYB117 also influences fruit metabolism. The tf-5 allele can be used to increase γ-aminobutyric acid, carotenoid, and folate levels in tomato.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT-qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.

DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium hybrid

BACKGROUND

Flower pigment and shape are determined by the coordinated expression of a set of structural genes during flower development. R2R3-MYB transcription factors are known regulators of structural gene expression. The current study focused on two members of this large family of transcription factors that were predicted to have roles in pigment biosynthesis and organ shape development in orchids.

METHODS

Phylogenetic analysis was used to identify candidate Dendrobium catenatum R2R3-MYB (DcaMYB) sequences associated with pigment and cell shape development. Gene silencing of candidate DhMYBs in Dendrobium hybrid by direct application of dsRNA to developing flowers was followed by observation of gene expression level and flower phenotypes. Silencing of the structural gene chalcone synthase was used as a comparative control.

KEY RESULTS

Ten candidate flower-associated DcaMYBs were identified. Flowers treated with dsRNA of DhMYB22 and DhMYB60 sequences were less pigmented and had relatively low expression of anthocyanin biosynthetic genes (F3'H and DFR), lower total anthocyanin concentration and markedly lower levels of cyanidin-3-glucoside and cyanidin-3-rutinoside. Petals of DhMYB22-treated flowers and sepals of DhMYB60-treated flowers showed the greatest colour difference relative to the same organs in untreated flowers. DhMYB22-treated flowers had relatively narrow and constricted lips, while DhMYB60-treated flowers had narrow and constricted sepals. No significant difference in shape was observed for DhCHS-treated or untreated flowers.

CONCLUSIONS

Our results demonstrate that DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium. This is a first report of MYB regulation of floral organ shape in orchids.

The R2R3-MYB gene family in Cicer arietinum: genome-wide identification and expression analysis leads to functional characterization of proanthocyanidin biosynthesis regulators in the seed coat

We identified 119 typical CaMYB encoding genes and reveal the major components of the proanthocyanidin regulatory network. CaPARs emerged as promising targets for genetic engineering toward improved agronomic traits in C. arietinum. Chickpea (Cicer arietinum) is among the eight oldest crops and has two main types, i.e., desi and kabuli, whose most obvious difference is the color of their seeds. We show that this color difference is due to differences in proanthocyanidin content of seed coats. Using a targeted approach, we performed in silico analysis, metabolite profiling, molecular, genetic, and biochemical studies to decipher the transcriptional regulatory network involved in proanthocyanidin biosynthesis in the seed coat of C. arietinum. Based on the annotated C. arietinum reference genome sequence, we identified 119 typical CaMYB encoding genes, grouped in 32 distinct clades. Two CaR2R3-MYB transcription factors, named CaPAR1 and CaPAR2, clustering with known proanthocyanidin regulators (PARs) were identified and further analyzed. The expression of CaPAR genes correlated well with the expression of the key structural proanthocyanidin biosynthesis genes CaANR and CaLAR and with proanthocyanidin levels. Protein-protein interaction studies suggest the in vivo interaction of CaPAR1 and CaPAR2 with the bHLH-type transcription factor CaTT8. Co-transfection analyses using Arabidopsis thaliana protoplasts showed that the CaPAR proteins form a MBW complex with CaTT8 and CaTTG1, able to activate the promoters of CaANR and CaLAR in planta. Finally, transgenic expression of CaPARs in the proanthocyanidin-deficient A. thaliana mutant tt2-1 leads to complementation of the transparent testa phenotype. Taken together, our results reveal main components of the proanthocyanidin regulatory network in C. arietinum and suggest that CaPARs are relevant targets of genetic engineering toward improved agronomic traits.

Systematic Analysis and Functional Characterization of R2R3-MYB Genes in Scutellaria baicalensis Georgi

R2R3-MYB transcription factors participate in multiple critical biological processes, particularly as relates to the regulation of secondary metabolites. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine and possesses various bioactive attributes including anti-inflammation, anti-HIV, and anti-COVID-19 properties due to its flavonoids. In the current study, a total of 95 R2R3-MYB genes were identified in S. baicalensis and classified into 34 subgroups, as supported by similar exon–intron structures and conserved motifs. Among them, 93 R2R3-SbMYBs were mapped onto nine chromosomes. Collinear analysis revealed that segmental duplications were primarily responsible for driving the evolution and expansion of the R2R3-SbMYB gene family. Synteny analyses showed that the ortholog numbers of the R2R3-MYB genes between S. baicalensis and other dicotyledons had a higher proportion compared to that which is found from the monocotyledons. RNA-seq data indicated that the expression patterns of R2R3-SbMYBs in different tissues were different. Quantitative reverse transcriptase-PCR (qRT-PCR) analysis showed that 36 R2R3-SbMYBs from different subgroups exhibited specific expression profiles under various conditions, including hormone stimuli treatments (methyl jasmonate and abscisic acid) and abiotic stresses (drought and cold shock treatments). Further investigation revealed that SbMYB18/32/46/60/70/74 localized in the nucleus, and SbMYB18/32/60/70 possessed transcriptional activation activity, implying their potential roles in the regulatory mechanisms of various biological processes. This study provides a comprehensive understanding of the R2R3-SbMYBs gene family and lays the foundation for further investigation of their biological function.

The CaMYB340 transcription factor induces chilling injury in post-harvest bell pepper by inhibiting fatty acid desaturation

Bell pepper (Capsicum annuum L.) is a tradable and desirable crop; however, its perishable nature requires low-temperature handling. Paradoxically, cold causes chilling injury (CI) and post-harvest waste. Current knowledge about CI in pepper is limited. The mechanism of CI is multi-faceted; therefore, we focused on fatty acid (FA) desaturation. We identified an upstream nuclear transcription factor (TF), CaMYB340, belonging to the R2R3 MYB subfamily, that negatively regulates FA desaturation and CaCBF3 expression and whose gene and protein expression is induced by low temperature (4°C). Specifically, McrBC treatment and bisulfite sequencing PCR indicate that exposure to cold triggers DNA methylation on one of the CHH sites in the CaMYB340 promoter. This epigenetic event at least partly contributes to the upregulation of CaMYB340 transcript levels. Increased expression of CaMYB340 results in the formation of protein complexes with CabHLH93 and CaMYB1R1, which in turn downregulate the expression of downstream genes. For peppers held at low temperature, transient overexpression of CaMYB340 reduced unsaturated FA content and membrane fluidity, resulting in cold-induced poor peel texture. Transient CaMYB340 silencing increased FA desaturation and lowered electrolyte leakage, enhancing cold tolerance in CaMYB340 knockdown fruits. Overall, these results underscore the intricacy of transcriptional networks in plants and highlight the role of CaMYB340 in CI occurrence in pepper fruits.

AcMYB1 Interacts With AcbHLH1 to Regulate Anthocyanin Biosynthesis in Aglaonema commutatum

Aglaonema commutatum is one of the most popular foliage plants with abundant leaf phenotypes; therefore, anthocyanin coloration is a vital economic trait in A. commutatum. However, the molecular mechanisms underlying anthocyanin biosynthesis and its regulation remain unclear. In this study, AcMYB1 and AcbHLH1, transcription factor genes related to an R2R3-myeloblast (MYB) and a basic helix-loop-helix (bHLH), respectively, were isolated from A. commutatum "Red Valentine" and functionally characterized. AcMYB1 and AcbHLH1 were found to interact by Y2H and BiFC assay. AcMYB1 was grouped into the AN2 subgroup and shared high homology with the known regulators of anthocyanin biosynthesis. Gene expression analysis showed that both AcMYB1 and AcbHLH1 have similar expression patterns to anthocyanin structural genes and correlate with anthocyanin distribution in different tissues of A. commutatum. Light strongly promoted anthocyanin accumulation by upregulating the expression of anthocyanin-related genes in A. commutatum leaves. Ectopic expression of AcMYB1 in tobacco remarkably increased anthocyanin accumulation in both vegetative and reproductive tissues at various developmental stages. These results provide insights into the regulation of anthocyanin biosynthesis in A. commutatum and are useful for breeding new A. commutatum cultivars with enhanced ornamental value.

R2R3-MYB gene family: Genome-wide identification provides insight to improve the content of proanthocyanidins in Trifolium repens

The R2R3-MYB family is one of largest transcription factor families in plants playing significant roles in regulating anthocyanin and proanthocyanidin biosynthesis. Proanthocyanidins are one of major objectives to improve the quality of white clover (Trifolium repens L.), which have a beneficial effect on ruminant to prevent the lethal pasture bloat. A total of 133 TrR2R3-MYB genes were identified and distributed on all 16 chromosomes based on the whole genome information of white clover. Also, by exploring the gene structure, motifs and duplication events of TrR2R3-MYBs, as well as the evolutionary relationship with TrR2R3-MYB genes of other species, 10 TrR2R3-MYB genes with the potential to regulate the anthocyanins and proanthocyanidins biosynthesis were screened. These TrR2R3-MYB genes responded significantly to low temperature in white clover. In addition, they have different expression patterns in leaves, petioles and inflorescences of white clover. Importantly, TrMYB116 and TrMYB118 may positively regulate anthocyanin accumulation and low temperature response in white clover. TrMYB118 may also be associated with anthocyanin pigmentation pattern in Purple leaves. This study provides a basis for verifying the function of TrR2R3-MYB and breeding white clover cultivars with high proanthocyanidins.

TcMYB29a, an ABA-Responsive R2R3-MYB Transcriptional Factor, Upregulates Taxol Biosynthesis in Taxus chinensis

Paclitaxel (Taxol), a highly modified diterpene agent mainly obtained from Taxus species, is the most widely used anticancer drug. Abscisic acid (ABA) is a well-known stress hormone that plays important roles in the secondary metabolism of plants, and it can also induce the accumulation of taxol in Taxus cell suspension cultures. However, the mechanism behind the regulation of taxol biosynthesis by ABA remains largely unknown. In previous research, a R2R3 MYB transcription factor (TF) TcMYB29a was observed to show a significant correlation with taxol biosynthesis, indicative of its potential role in the taxol biosynthesis. In this study, the TcMYB29a encoded by its gene was further characterized. An expression pattern analysis revealed that TcMYB29a was highly expressed in the needles and roots. Overexpression of TcMYB29a in Taxus chinensis cell suspension cultures led to an increased accumulation of taxol, and upregulated expression of taxol-biosynthesis-related genes, including the taxadiene synthase (TS) gene, the taxane 5α-hydroxylase (T5OH) gene, and the 3'-N-debenzoyl-2'-deoxytaxol-N-benzoyltransferase (DBTNBT) gene as compared to the controls. Chromatin immunoprecipitation (ChIP) assays, yeast one-hybrid (Y1H) assays, electrophoretic mobility shift assays (EMSAs), and dual-luciferase reporter assays verified that TcMYB29a could bind and activate the promoter of TcT5OH. Promoter sequence analysis of TcMYB29a revealed that its promoter containing an AERB site from -313 to -319 was a crucial ABA-responsive element. Subsequently, the ABA treatment assay showed that TcMYB29a was strongly upregulated at 6 h after ABA pretreatment. Furthermore, TcMYB29a was strongly suppressed at 3 h after the methyl jasmonate (MeJA) treatment and was depressed to the platform at 12 h. Taken together, these results reveal that TcMYB29a is an activator that improves the accumulation of taxol in Taxus chinensis cells through an ABA-medicated signaling pathway which is different from JA-medicated signaling pathways for the accumulation of taxol. These findings provide new insights into the potential regulatory roles of MYBs on the expression of taxol biosynthetic genes in Taxus.

Genome-Wide Identification and Expression Analysis of the R2R3-MYB Transcription Factor Family Revealed Their Potential Roles in the Flowering Process in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is a productive fruit crop with high nutritional and medical value in tropical and subtropical regions. The MYB gene family is one of the most widespread plant transcription factor (TF) families participating in the flowering regulation. However, little is known about the MYB TFs involved in the flowering process in longan and its regulatory network. In this study, a total of 119 DlR2R3-MYB genes were identified in the longan genome and were phylogenetically grouped into 28 subgroups. The groupings were supported by highly conserved gene structures and motif composition of DlR2R3-MYB genes in each subgroup. Collinearity analysis demonstrated that segmental replications played a more crucial role in the expansion of the DlR2R3-MYB gene family compared to tandem duplications, and all tandem/segmental duplication gene pairs have evolved under purifying selection. Interspecies synteny analysis among longan and five representative species implied the occurrence of gene duplication events was one of the reasons contributing to functional differentiation among species. RNA-seq data from various tissues showed DlR2R3-MYB genes displayed tissue-preferential expression patterns. The pathway of flower development was enriched with six DlR2R3-MYB genes. Cis-acting element prediction revealed the putative functions of DlR2R3-MYB genes were related to the plant development, phytohormones, and environmental stresses. Notably, the orthologous counterparts between Arabidopsis and longan R2R3-MYB members tended to play conserved roles in the flowering regulation and stress responses. Transcriptome profiling on off-season flower induction (FI) by KClO₃ indicated two up-regulated and four down-regulated DlR2R3-MYB genes involved in the response to KClO₃ treatment compared with control groups. Additionally, qRT-PCR confirmed certain genes exhibited high expression in flowers/flower buds. Subcellular localization experiments revealed that three predicted flowering-associated MYB proteins were localized in the nucleus. Future functional studies on these potential candidate genes involved in the flowering development could further the understanding of the flowering regulation mechanism.

[Identification and expression profiling of R2R3-MYB transcription factors in Erigeron breviscapus]

R2 R3-MYB transcription factors are ubiquitous in plants, playing a role in the regulation of plant growth, development, and secondary metabolism. In this paper, the R2 R3-MYB transcription factors were identified by bioinformatics analysis of the genomic data of Erigeron breviscapus, and their gene sequences, structures, physical and chemical properties were analyzed. The functions of R2 R3-MYB transcription factors were predicted by cluster analysis. Meanwhile, the expression patterns of R2 R3-MYB transcription factors in response to hormone treatments were analyzed. A total of 108 R2 R3-MYB transcription factors, named EbMYB1-EbMYB108, were identified from the genome of E. breviscapus. Most of the R2 R3-MYB genes carried 2-4 exons. The phylogenetic tree of MYBs in E. breviscapus and Arabidopsis thaliala was constructed, which classified 234 MYBs into 30 subfamilies. The MYBs in the five MYB subfamilies of A.thaliala were clustered into independent clades, and those in E. breviscapus were clustered into four clades. The transcriptome data showed that MYB genes were differentially expressed in different tissues of E. breviscapus and in response to the treatments with exogenous hormones such as ABA, SA, and GA for different time. The transcription of 13 R2 R3-MYB genes did not change significantly, and the expression patterns of some genes were up-regulated or down-regulated with the extension of hormone treatment time. This study provides a theoretical basis for revealing the mechanisms of R2 R3-MYB transcription factors in regulating the growth and development, stress(hormone) response, and active ingredient accumulation in E. breviscapus.

Functional Characterization of Hedychium coronarium J. Koenig MYB132 Confers the Potential Role in Floral Aroma Synthesis

The R2R3-MYB transcription factors (TFs) play several key roles in numerous plant biological processes. Hedychium coronarium is an important ornamental plant well-known for its elegant flower shape and abundant aroma type. The floral aroma of H. coronarium is due to the presence of a large amount of terpenes and benzenoids. However, less is known about the role of R2R3-MYB TFs in the regulatory mechanism of floral aroma production in this breed. Herein, we isolate and functionally characterize the R2R3-MYB TF HcMYB132, which is potentially involved in regulating floral aroma synthesis. Sequence alignment analysis revealed that it includes a nuclear localization signal NLS(s) and a 2R, 3R motif signature in the sequences. A subcellular localization assay revealed that HcMYB132 protein localizes to the nucleus. Real-time qPCR assays showed that HcMYB132 is specifically expressed in flowers and its expression pattern correlates with the emission of floral volatile compounds. In HcMYB132-silenced flowers, the levels of floral volatile compounds were significantly reduced, and the expression of key structural volatile synthesis genes was downregulated compared to control. Collectively, these results suggest that HcMYB132 might play a significant role in the regulation of terpenoid biosynthesis in H. coronarium.

A Novel R2R3-MYB Gene LoMYB33 From Lily Is Specifically Expressed in Anthers and Plays a Role in Pollen Development

Lily (Lilium spp.) is an important commercial flower crop, but its market popularity and applications are adversely affected by severe pollen pollution. Many studies have examined pollen development in model plants, but few studies have been conducted on flower crops such as lily. GAMYBs are a class of R2R3-MYB transcription factors and play important roles in plant development and biotic resistance; their functions vary in different pathways, and many of them are involved in anther development. However, their function and regulatory role in lily remain unclear. Here, the GAMYB homolog LoMYB33 was isolated and identified from lily. The open reading frame of LoMYB33 was 1620 bp and encoded a protein with 539 amino acids localized in the nucleus and cytoplasm. Protein sequence alignment showed that LoMYB33 contained a conserved R2R3 domain and three BOX motifs (BOX1, BOX2, and BOX3), which were unique to the GAMYB family. LoMYB33 had transcriptional activation activity, and its transactivation domain was located within 90 amino acids of the C-terminal. LoMYB33 was highly expressed during the late stages of anther development, especially in pollen. Analysis of the promoter activity of LoMYB33 in transgenic Arabidopsis revealed that the LoMYB33 promoter was highly activated in the pollen of stage 12 to 13 flowers. Overexpression of LoMYB33 in Arabidopsis significantly retarded growth; the excess accumulation of LoMYB33 also negatively affected normal anther development, which generated fewer pollen grains and resulted in partial male sterility in transgenic plants. Silencing of LoMYB33 in lily also greatly decreased the amount of pollen. Overall, our results suggested that LoMYB33 might play an important role in the anther development and pollen formation of lily.

Genome-wide analysis and characterization of R2R3-MYB family in pigeon pea (Cajanus cajan) and their functional identification in phenylpropanoids biosynthesis

Thirty CcMYB were identified to involve in flavonoid and lignin biosynthesis in pigeon pea genome. A comprehensive analysis of gene structure, phylogenetic relationships, distribution on chromosomes, gene duplication, and expression patterns was performed. MYB transcription factor is one of the largest gene families in plants and plays critical roles in plant growth and development, as well as resistance to biotic and abiotic stress. However, the function of MYB genes in pigeon pea (Cajanus cajan) remains largely unknown. Here, 30 R2R3-MYB which involved flavonoid and lignin biosynthesis were identified in the pigeon pea genome and were classified into five groups based on phylogenetic analysis. Simultaneously, another 122 key enzyme genes from biosynthetic pathways of flavonoid and lignin were identified and all of them were mapped on 11 chromosomes with the co-linearity relationship. Among these genes, the intron/exon organization and motif compositions were conserved and they have undergone a strong purifying selection and tandem duplications during evolution. Expression profile analysis demonstrated most of these genes were expressed in different tissues and responded significantly to MeJA, RNA-seq analysis revealed clear details of genes varied with time of induction. Ten key genes from the phenylpropanoid pathway were selected to further verify whether they responded to induction under different abiotic stress conditions (UV-B, cold, heat, salt, drought, and GA₃). This study elaborates on potential regulatory relationships between R2R3-MYB genes and some key genes involved in flavonoid and lignin biosynthesis under MeJA treatment, as well as adding to the understanding of improving abiotic stress tolerance and regulating the secondary metabolism in woody crops. A simplified discussion model for the different regulation networks involved with flavonoid and lignin biosynthesis in pigeon pea is proposed.

Anthocyanin synthesis in orange carrot cv. Danvers is activated by transgene expression of the transcription factors DcMYB113_NB and DcEGL1_NB from black carrot cv. Nightbird

Black carrots are potent sources of anthocyanin for the natural food color industry as their anthocyanins contain very high percentages of acylated anthocyanins which are much more stable than non-acylated anthocyanins. Anthocyanins are synthesized by a specific branch of the phenylpropanoid pathway activated by a triad of R2R3-MYB, bHLH and WD40 transcription factors (TFs). Recent studies in black carrots have elucidated major anthocyanin related structural genes and also regulatory TFs. However, the active TFs responsible for anthocyanin production in black carrots differ between cultivars. We have previously shown by RNAseq that DcMYB113 (LOC108213488), a R2R3-MYB TF, was up-regulated in colored as compared to non-colored tissues of the black carrots 'Superblack' and 'CH05544' and that this upregulation was positively correlated with anthocyanin content. However, this gene showed no upregulation in the black carrot 'Nightbird' also included in that study. In the present study, we present a novel R2R3-MYB DcMYB113_NB (LOC108212072) and a complementary bHLH DcEGL1_NB (LOC108210744) isolated from the RNA of 'Nightbird'. Their functionality as anthocyanin regulators was confirmed by their simultaneous expression under the control of a constitutive promoter in the background of the orange carrot 'Danvers 126'. Transformants showed activation of the structural anthocyanin genes and accumulation of anthocyanins across leaves, stems and taproots. Interestingly, the anthocyanin profile of the transformants showed increases of 20 to 30% in acylated anthocyanins as compared to 'Nightbird' resulting in transformants with almost 100% acylated anthocyanins.

R2R3-MYB transcription factor family in tea plant (Camellia sinensis): Genome-wide characterization, phylogeny, chromosome location, structure and expression patterns

MYB transcription factors play essential roles in many biological processes and environmental stimuli. However, the functions of the MYB transcription factor family in tea plants have not been elucidated. Here, a total of 122 CsR2R3-MYB genes were identified from the chromosome level genome of tea plant (Camellia sinensis). The CsR2R3-MYB genes were phylogenetically classified into 25 groups. Results from the structure analysis of the gene, conserved motifs, and chromosomal distribution supported the relative conservation of the R2R3-MYB genes family in the tea plant. Synteny analysis indicated that 122, 34, and 112 CsR2R3-MYB genes were orthologous to Arabidopsis thaliana, Oryza sativa and C. sinensis var. 'huangdan' (HD), respectively. Tissue-specific expression showed that all CsR2R3-MYB genes had different expression patterns in the tea plant tissues, indicating that these genes may perform diverse functions. The expression patterns of representative R2R3-MYB genes and the regulatory network of the main anthocyanin components were analyzed, which suggested that CsMYB17 may played a key role in the regulation of cya-3-O-gal, del-3-O-gal, cya-3-O-glu and pel-3-O-glu. Results from the qRT-PCR validation of selected genes suggested that CsR2R3-MYB genes were induced in response to drought, cold, GA, and ABA treatments. Overall, this study provides comprehensive and systematic information for research on the function of R2R3-MYB genes in tea plants.