Betalain biosynthesis in red pulp pitaya is regulated via HuMYB132: a R-R type MYB transcription factor
BMC Plant Biology volume 23, Article number: 28 (2023)
Multiple MYB transcription factors (TFs) are involved in the regulation of plant coloring. Betalain is a kind of natural plant pigment and its biosynthesis is regulated by a number of enzymes. Despite this, little is known about the molecular properties and roles of MYB TFs in pitaya betalain biosynthesis.
In the present study, we identified a 1R-MYB gene, HuMYB132, which is preferentially expressed in red-pulp pitaya at the mature stage. It was clustered with Arabidopsis R-R-type genes and had two DNA-binding domains and a histidine-rich region. The expression assays in N. benthamiana and yeast indicated that HuMYB132 is a nucleus-localized protein with transcriptional activation activity. Dual luciferase reporter assay and electrophoretic mobility shift assays (EMSA) demonstrated that HuMYB132 could promote the transcriptional activities of HuADH1, HuCYP76AD1–1, and HuDODA1 by binding to their promoters. Silencing HuMYB132 reduced betalain accumulation and the expression levels of betalain biosynthetic genes in pitaya pulps.
According to our findings, HuMYB132, a R-R type member of 1R-MYB TF subfamily, positively regulates pitaya betalain biosynthesis by regulating the expression of HuADH1, HuCYP76AD1–1, and HuDODA1. The present study provides a new theoretical reference for the management of pitaya betalain biosynthesis and also provides an essential basis for future regulation of betalain biosynthesis in Hylocereus.
Transcription factors (TFs) play important roles in plant growth and development. Based on the pitaya genome and transcriptome data, a comprehensive regulatory network of betalain biosynthesis in pitaya fruit was constructed to provide multiple potential candidate TFs in the pitaya betalain biosynthesis . MYB TFs are among the most abundant plant-specific TFs which are classified into four subfamilies according to the presence of one to four highly conserved MYB repeats: MYB-related (1R-MYB), R2R3-MYB (2R-MYB), R1R2R3-MYB (3R-MYB), and 4R-MYB . Surprisingly, the pitaya genome contains 75 1R-MYB, 105 R2R3-MYB, four R1R2R3-MYB, and one 4R-MYB TFs. HuMYB1 is an R2R3-MYB negative regulator of pitaya betalain biosynthesis that inhibits the transcription activities of HuADH1, HuCYP76AD1–1, and HuDODA1 . Furthermore, an inventory of the Arabidopsis genome reveals that it has 126 R2R3-MYB, five R1R2R3-MYB, 64 MYB-like, and three atypical MYB genes . Many 3R-MYB and 2R-MYB genes have been extensively studied in A. thaliana, giving us a better understanding of MYB gene functions in other plant species . For instance, a homologous R2R3-MYB TF of AtPAP1/2 in beet, BvMYB1, enhanced betalain biosynthesis by activating the expression levels of BvCYP76AD1 and BvDODA1 . However, the identification and role of other MYB TFs regulating betalain biosynthesis remains largely unknown.
Pitaya is a popular tropical fruit belonging to the Hylocereus and Seleniereus in the Cactaceae family of the Caryophyllales order. According to peel and pulp color, pitayas are mainly categorized into three species: H. undatus (red peel with white pulp), H. monacanthus or H. polyrhizus (red peel with red pulp), and H. megalanthus or S. megalanthus (yellow peel with white pulp) . Pitaya is the only commercial cultivation of fruit containing abundant betalains for consumers. Betalain biosynthesis and accumulation are responsible for different peel and pulp colors in pitaya. Interestingly, betalains are restricted to the angiosperm order of core Caryophyllales, with the exception of Caryophyllaceae, Molluginaceae, Kewaceae, Limeaceae, Macarthuriaceae, and Simmondsiaceae families, which have been reported to produce anthocyanins . Betalains and anthocyanins have never been found in the same plant species and share a mutual exclusion relationship . Interestingly, the precursors of betalains and anthocyanins are Tyr and Phe, respectively, and both of them are synthesized from arogenate [10, 11]. Moreover, a recent study inferred that anthocyanins in most plants of Caryophyllales were replaced by betalains, probably due to betalain biosynthetic genes, i.e., ADH, CYP76AD, and DODA happening in several evolutionary events . However, the mechanism of the mutual exclusion between betalains and anthocyanins still remains unknown.
Betalains are vacuole-localized, water-soluble, nitrogen-containing, and tyrosine-derived pigments which are divided into purple-red betacyanins and yellow betaxanthins . Three key enzymatic reactions have been characterized, including the tyrosine biosynthesis by arogenate dehydrogenase (ADH) , the betalamic acid biosynthesis by 4,5-DOPA extradiol dioxygenase (DODA)  and the cyclo-DOPA synthesis by a cytochrome P450 enzyme (CYP76AD1) , which were critical for betalain biosynthesis. Betacyanins are synthesized via glucosyltransferase (GT) and acyltransferase (ATs), while betaxanthins are produced through spontaneous reactions [15,16,17,18]. So far, candidate genes, proteins and metabolites of betalain biosynthesis were isolated in green-, red-, and yellow-peeled pitayas based on microRNAs, transcriptomics, proteomics, and metabolomics analyses [19,20,21,22,23].
Our previous studies suggested that betalain-related structural genes, i.e. HuADH1, HuCYP76AD1–1, and HuDODA1 were gathered with other Caryophyllales ADHα, CYP76AD1α, and DODAα genes, respectively, on the same chromosome in pitaya genome . Moreover, a predicted R-R-type member of 1R-MYB TF, HuMYB132 (HU04G01397.1), showed a co-expression pattern relationship with betalain biosynthesis genes by weighted correlation network analysis (WGCNA) [1, 3]. As a result, this study was designed for the first time to elucidate the molecular mechanism of HuMYB132 TF which was further explored to be responsible for regulating the expression of betalain biosynthesis genes of HuADH1, HuCYP76AD1–1, and HuDODA1. Our findings could shed light on the regulatory mechanism of 1R-MYB in pitaya betalain biosynthesis.
Identification of HuADH1, HuCYP76AD1–1, and HuDODA1 involved in pitaya Betalain biosynthesis
Our previous studies showed that HuADH1, HuCYP76AD1–1, and HuDODA1 were respectively gathered with ADHα, CYP76AD1α, and DODAα genes and their expression levels in the pulp of H. monacanthus (red peel with red pulp) were significantly higher than those of H. undatus (red peel with white pulp) . HuADH1, HuCYP76AD1–1, and HuDODA1 are co-localized on chromosome 3 according to the H. undatus genome data. A transient expression assay was performed in N. benthamiana leaves to elucidate the roles of HuADH1, HuCYP76AD1–1, and HuDODA1 in betalain biosynthesis. As shown in Fig. 1 (Additional file 1: Fig. S1), no betalain accumulation was detected in N. benthamiana leaves when HuADH1, HuCYP76AD1–1, and HuDODA1 were individually overexpressed, which is consistent with the phenomenon of the empty vector. Co-infiltration of A. tumefaciens containing HuCYP76AD1–1 and HuDODA1 resulted in the production of betalains compared to no visible pigmentation for co-expression of HuADH1 and HuCYP76AD1–1. Betalain contents in N. benthamiana leaves of co-expression of HuADH1, HuCYP76AD1–1, and HuDODA1 were two times higher than that of co-expression of HuCYP76AD1–1 and HuDODA1. These results suggested that HuCYP76AD1–1 and HuDODA1 are key structural genes for betalain biosynthesis in pitaya, while HuADH1 is a key structural gene for promoting betalain accumulation.
HuMYB132 exhibits a preferentially expression level in red-pulp pitaya at mature stage
Based on the WGCNA of pitaya transcriptome data, a 1R-MYB gene named HuMYB132 was obtained, and it was predominantly expressed in the mature fruit of red-pulp pitaya [1, 3]. The highest expression level of HuMYB132 was detected in mature fruit (32 d), which is consistent with RNA-Seq data during pitaya fruit development (Fig. 2A). HuMYB132 belonged to the R-R-type MYB TFs and was phylogenetically related to AT1G49010.1 (AtMYBL) involving in abiotic stress (Fig. 2B) . It encoded 296 amino acid residues with a theoretical pI of 9 and a molecular weight of 32 kDa, as well as two DNA binding domains (DBD) (named DBD I and DBD II) at the N-terminal and a histidine-rich region at the C-terminal (Fig. 2C). Based on these results, it seems that HuMYB132 is a candidate 1R-MYB TF responsible for betalain biosynthesis in pitaya.
HuMYB132 is a nuclear protein with transcriptional activation activity
To investigate the subcellular localization of HuMYB132, its encoding region was fused with the GFP gene under the control of the CaMV promoter. The GFP-HuMYB132 fluorescence was exclusively observed in the nucleus of N. benthamiana protoplasts while the positive control fluorescence was dispersed throughout the cytosol and nucleus (Fig. 3A). These results suggested that HuMYB132 is a nuclear-localized protein.
To explore whether HuMYB132 has transcriptional activation activity, the full-length cDNA was fused to the GAL4BD domain of the BD-62SK vector as an effector and transiently expressed in N. benthamiana leaves with reporter (Fig. 3B). The BD-62SK-HuMYB132 and positive control (BD-62SK-VP16) showed a significantly higher value of LUC/REN ratio than the negative control (BD-62SK) (Fig. 3C). Besides, full-length cDNA of HuMYB132 was fused to the GAL4BD domain of the pGBKT7 vector and introduced into yeast cells. Compared with negative control (pGBKT7), the transformed yeast cells of pGBKT7-HuMYB132 and positive control (pGBKT7-P53 + pGADT7-T) grew well on SD/−Trp-His-Ade medium and showed α-Gal activity (Fig. 3D). These results suggested that HuMYB132 has transcriptional activation activity.
HuMYB132 directly binds to the promoters of Betalain biosynthesis genes and activates their transcription
According to PlantCARE database (http://bioinformatics.Psb.ugent.be/webtools/ plantcare/html/), all promoters of HuADH1, HuCYP76AD1–1, and HuDODA1 have MYB cis-elements, including MRE (CAACCA and ACCTAA) and MYB-binding sites (MBS, T/CAACTG). To elucidate whether HuMYB132 directly targets the promoters of HuADH1, HuCYP76AD1–1, and HuDODA1 genes, their promoter regions were cloned into a pABAi vector and introduced into yeast cells. As shown in Fig. 4A, the bait yeast cells respectively harboring HuADH1, HuCYP76AD1–1, and HuDODA1 promoter plasmids could not grow on the SD/−Ura medium supplemented with 200 ng/mL AbA. Subsequently, the full-length of HuMYB132 was inserted into the pGADT7 vector and transformed into the bait yeast cells. Yeast cells co-transformed with HuMYB132 and HuADH1 or HuCYP76AD1–1 or HuDODA1 grew well on the SD/−Leu medium with 200 ng/mL AbA, indicating that HuMYB132 could directly bind to the promoters of HuADH1, HuCYP76AD1–1 and HuDODA1 (Fig. 4B). Besides, an in vitro EMSA assay was performed to verify the binding ability of HuMYB132 to HuADH1, HuCYP76AD1–1, and HuDODA1 promoters. Clear mobility shifts in bands were detected when DNA fragments from HuADH1, HuCYP76AD1–1, and HuDODA1 promoters containing MRE and MBS elements were biotin-labelled and incubated with purified recombinant HuMYB132 protein (Fig. 4C-E; Additional file 1: Fig. S4-S6). In addition, the shift bands became weaker with the addition of the unlabeled probes compared to no significant change happened with the presence of the mutant probes. These results indicated that HuMYB132 can directly bind to promoters of HuADH1, HuCYP76AD1–1, and HuDODA1 via the MRE and MBS elements.
To verify the effects of HuMYB132 on the transcription of betalain biosynthetic genes, the coding region of HuMYB132 and the promoter regions of HuADH1, HuCYP76AD1–1, and HuDODA1 were inserted into the dual luciferase assay system (Fig. 5A). Compared with the empty control, the infiltration of HuMYB132 could activate the HuADH1, HuCYP76AD1–1, and HuDODA1 promoters and their LUC/REN ratio values were significantly increased (Fig. 5B). These results indicated that HuMYB132 could activate the transcriptions of HuADH1, HuCYP76AD1–1, and HuDODA1 responsible for pitaya betalain biosynthesis.
Silencing HuMYB132 results in reduction of Betalain production
A gene silencing assay was performed to further explore the function of HuMYB132 in the pulp of ‘Guanhuahong’ pitaya. Compared with the empty vector, silencing HuMYB132 repressed the process of pulp coloration due to the reduction of betalains, especially for betacyanins (Fig. 6A-C). Results from RT-qPCR showed that silencing HuMYB132 resulted in significant decreases in expression levels of HuADH1, HuCYP76AD1–1, and HuDODA1 (Fig. 6D). These results suggested that HuMYB132 acts as a positive regulator involved in the betalain biosynthesis of pitaya.
ADH, an enzyme responsible for tyrosine synthesis, underwent a gene duplication event that resulted in the deregulated ADHα clade and the canonical tyrosine-sensitive ADHβ clade . Due to gene duplication events, CYP76AD1, a novel cytochrome P450 (Cyt P450), is divided into betacyanin-functional CYP76AD1-α, betaxanthin-functional CYP76AD1-β, and unknown-functional CYP76AD1-γ . DODA is an enzyme catalyzing betalamic acid synthesis, which is classified into functional DODAα and nonfunctional DODAβ based on gene duplication events . In B. vulgaris, C. quinoa, Amaranthus tricolor, Mirabilis jalapa, Stegnosperma halimifolium, Carnegiea gigantea and Mesembryanthemum crystallinum, the functions of CYP76AD1-α and DODAα clade homologous genes have been figured out [17, 28,29,30,31]. However, the involvement of ADHα in betalain biosynthesis has only been investigated in beet . Our previous study suggested that HuADH1, HuCYP76AD1–1, and HuDODA1 are homologs of betalain-specific ADHα, CYP76AD1-α, and DODAα clades, respectively . In this study, we found that HuCYP76AD1–1 and HuDODA1 were key structural genes, while HuADH1 encoded a limitation enzyme in pitaya betalain biosynthesis (Fig. 1). These betalain biosynthesis genes were co-localized on chromosome 3, while HuCYP76AD1–1 and HuDODA1 showed a closer distance than HuADH1 in H. undatus, probably implicating a more efficient betalain biosynthesis in pitaya (Fig. 7). B. vulgaris also has betalain-specific CYP76AD1-α and DODAα homologs on chromosome 2 that is close to each other . So, HuADH1, HuCYP76AD1–1, and HuDODA1 are three important structural genes contributing to betalain biosynthesis in pitaya.
Based on the WGCNA of betalain biosynthesis genes and potential TFs, HuMYB132 was preferentially expressed in red-pulp pitaya at the mature stage , which is further verified by RT-qPCR analyses (Fig. 2). MYB-related subfamily members contain a single or a partial MYB repeat and are classified into five subgroups: CIRCADIAN CLOCK ASSOCIATED 1 (CCA1)-like, CAPRICE (CPC)-like, TELOMERIC DNA-BINDING PROTEIN (TBP)-like, I-box-binding-like and R-R-type. A root tree of HuMYB132 and 64 Arabidopsis 1R-MYB proteins was constructed using AT1G09770.1 (AtCDC5, an ‘unusual’ MYB gene), suggesting that HuMYB132 was a putative R-R-type MYB-like TF (Fig. 2B). The sequence analysis indicated that HuMYB132 has two separated DBD domains and a histidine region (Fig. 2C). The SHAQKYF amino acid signature motif in DBD II is highly conserved in MYB-related genes and can be used to distinguish between R2R3-MYB genes . Moreover, HuMYB132 is a nucleus protein with transcriptional activator activity, consistent with the previous report that DBD I was crucial for transcriptional activation activity (Fig. 3) . MYB TFs can regulate downstream gene expression by interacting with MYB-binding sites in promoters [33,34,35]. In our study, HuMYB132 promoted the expressions of HuADH1, HuCYP76AD1–1, and HuDODA1 by specifically binding to the MER and MBS elements of their promoters on pitaya chromosome 3, enhancing the reaction from arogenate to products of L-DOPA, cyclo-DOPA, and betalamic acid (Fig. 7A). Subsequently, those products were transferred into kinds of betacyanins via GT genes, which were located in other pitaya chromosome, and/or generate betaxanthins through spontaneous reactions (Fig. 7B). Thus, HuMYB132 is a positive regulator and has a transcription effect on betacyanin biosynthesis in pitaya. Similarly, HpWRKY44 and HmoWRKY40 were also positive regulators of pitaya betacyanin biosynthesis by promoting the expression of HuCYP76AD1–1 [36, 37]. However, more studies are necessary to determine whether HuMYB132 can coordinate with the other TFs in regulation of pitaya betalain biosynthesis. Further studies are necessary to explore the relationships between HuMYB132 or other TFs and other key betalain biosynthesis genes.
A previous study reported that a rice R-R-type MYB TF, MID1 (Os05g37060), was responsive to drought stresses during reproductive development by improving the contents of proline, soluble sugar, and the activities of peroxidase and superoxide dismutase . During seed germination of Arabidopsis, AtDIV2 (At5g04760, a R-R-type MYB TF) integrated ABA signaling to negatively regulate salt stress . At5g58900 (A Arabidopsis R-R MYB TF) was involved in the NO signaling pathway through binding to the NAC027 promoter, resulting in a response to aluminum stress . Moreover, AtMYBL (At1g49010) was reported to promote leaf senescence via decreasing chlorophyll content and regulating ABA or salinity signal transduction . These results suggested that R-R-type MYB TFs are responsive to abiotic stresses. In general, betalain-producing plants can adapt to regions with high temperatures, drought, and salty environmental conditions due to betalain metabolites playing important roles in plant abiotic and biotic stress resistance . However, more research is needed to determine whether HuMYB132 is responsible for betalain biosynthesis in response to stress in the environment.
In summary, HuADH1, HuCYP76AD1–1, and HuDODA1 genes were critical for pitaya betalain biosynthesis. HuMYB132 is a member of the R-R-type sub-family and it is a nuclear protein with transcriptional activation activity. HuMYB132 positively regulates the expression of betalain biosynthesis genes by binding to the MRE and MBS elements in the promoters of HuADH1, HuCYP76AD1–1, and HuDODA1. Silencing HuMYB132 reduced the expression levels of HuADH1, HuCYP76AD1–1, and HuDODA1, leading to the reduction of betalain content in red pitaya pulp. The present studies provide new insights into understanding the molecular mechanism of betalain biosynthesis in pitayas.
‘Guanhuahong’ (red peel with red pulp, H. monacanthus) and Nicotiana benthamiana were used as plant materials. N. benthamiana was grown in a greenhouse with a condition of 16 h/8 h day/night at 25 °C and was used for transient expression, subcellular localization, and dual luciferase reporter assay. The South China Agricultural University provided all plant materials used in this study, and no specific permissions were required for the collection of those samples for research purposes following institutional, national, and international guidelines. Fruits of ‘Guanhuahong’ pitaya were cultivated in the orchard of Jinsuinong (Zhongluotan Village, Guangzhou City, Guangdong Province, China) and collected on the 14th, 17th, 19th, 23rd, 25th, 27th, and 32nd day after flowering (DAF). All samples were frozen in liquid nitrogen immediately, and stored at − 80 °C until future analyses.
Expression and sequence analyses
Total RNA of infiltration N. benthamiana leaves, pitaya pulps (14th, 17th, 19th, 23rd, 25th, 27th, and 32nd DAF), and injection pitaya pulps were isolated using the EASYspin Plus Complex Plant RNA Kit (RN53) (Aidlab Biotechnology, Beijing, China), and then cDNA was synthesized using the PrimeScript™ RT Reagent Kit with gDNA Eraser (TaKaRa, Shiga, Japan).
Reverse transcription PCR (RT-PCR) was conducted using I-5™ 2 × High-Fidelity Master Mix (MCLAB, USA) and pictured by Gel Doc™ XR+ Imaging System (Bio-Rad, USA). Real-time quantitative PCR (RT-qPCR) was performed on an CFX384-Real-Time System (C1000 Touch Thermal Cycler, USA) using the RealUniversal Color PreMix (SYBR Green) (TIANGEN, China) and specific primers (Additional file 2). The Actin genes of pitaya and N. benthamiana are respectively referred to Chen et al.  and Zhao et al. . All determinations were performed in three biological repetitions with three technical replicates.
The full-length coding sequence of HuADH1, HuCYP76AD1–1, HuDODA1 and HuMYB132 were cloned from ‘Guanhuahong’ pitaya pulps using I-5™ 2 × High-Fidelity Master Mix (MCLAB, USA) with specific primers (Additional file 2). Besides, the phylogenetic tree of HuMYB132 and Arabidopsis 1R-MYB proteins was constructed by the Maximum likelihood method (ML) in MEGA 7 with 1000 bootstrap and displayed by EVOLVIEW online tool (https://www.evolgenius.info/evolview/). The protein sequences of HuMYB132, XP_021748976.1 (Chenopodium quinoa), XP_010673224.1 (Beta vulgaris), XP_021854462.1 (Spinacia oleracea), AT1G49010.1 (Arabidopsis thaliana), and RVX14086.1 (Vitis vinifera) were aligned by Clustal X software and shown by GeneDoc software.
Measurement of Betalain contents
Betalains were extracted and measured following our previously described method . In brief, 0.5 g of freeze-dried samples were ground into fine powder and extracted with 5 mL 80% aqueous methanol (v/v) solution by sonication for 10 min and stirred for 20 min in the dark at room temperature. After centrifuging at 5000 rpm for 15 min, the residues were subjected to a similar second extraction. The supernatants were measured by spectrophotometry (Infinite M200, Tecan Co.) at 478 nm for betaxanthins and 538 nm for betacyanins. All determinations were performed in three biological repetitions with three technical replicates.
The transient expression of Betalain pathway genes in N. benthamiana
The coding sequences of betalain biosynthesis genes were inserted into the pEAQ vector to create pEAQ-HuADH1, pEAQ-HuCYP76AD1–1 and pEAQ-HuDODA1 fusion constructs (primers are listed in Additional file 2), which were separately transformed into Agrobacterium tumefaciens strain GV3101 and then infiltrated into N. benthamiana leaves . 5 days after infiltration, N. benthamiana leaves were frozen in liquid nitrogen immediately and stored at − 80 °C for betalain content measurement and gene expression analyses. The experiments were repeated three biological replicates.
The full-length coding sequence of HuMYB132 was inserted into the pGreen-35S-GFP vector (primers are listed in Additional file 2), transformed into A. tumefaciens strain GV3101 (pSoup-p19), and infiltrated into N. benthamiana leaves. 2 days after infiltration, leaf protoplasts were isolated according to  and the GFP fluorescence was observed by a fluorescence microscope (ZEISS LCM-800, Germany).
Transactivation activity in yeast cells and N. benthamiana
The coding sequence of HuMYB132 was fused with the GAL4 DNA-binding domain (DBD) in the pGBKT7 vector. Then, pGBKT7-HuMYB132, negative control (pGBKT7) and positive control (pGBKT7-P53 and pGADT7-T) were separately transformed into yeast cells and cultured on the synthetic dropout medium without tryptophan (SD/−Trp) or without tryptophan, histidine, and adenine (SD/−Trp -His -Ade). The transactivation activity of HuMYB132 in yeast cells was assessed based on their growth status in SD/−Trp -His -Ade medium after 3 d in 30 °C, followed by x-α-galactosidase (X-ɑ-Gal) incubation for 30 min.
The full-length of HuMYB132 was fused with GAL4BD in the BD-62-SK vector driven by the Cauliflower Mosaic Virus (CaMV) 35S promoter. Then, the effector of BD-62-SK-HuMYB132, internal control (BD-62-SK), positive control (BD-62-SK-VP16), and reporter plasmids were separately introduced into the A. tumefaciens strain GV3101 (pSoup), followed by the co-infiltration of effector and reporter into N. benthamiana leaves. The transactivation activity of HuMYB132 in N. benthamiana was assessed using the Dual-Luciferase Reporter Assay System (Promega, USA) after 3 d. The experiments were repeated three biological replicates.
Yeast one-hybrid assay
‘Guanhuahong’ pitaya DNA was isolated using Plant DNA Extraction Kit (DN14) (Aidlab, China) and Ribonuclease A (RNase A) (TaKaRa, Shiga, Japan). The DNA fragments of HuADH1, HuCYP76AD1–1, and HuDODA1 promoters were cloned using Seq Amp DNA Polymerase (TaKaRa, Shiga, Japan) and inserted into the pABAi vector with specific primers (Additional file 2). Then, the recombination plasmids of pABAi-HuADH1, pABAi-HuCYP76AD1–1, and pABAi-HuDODA1 were transformed into yeast cells and cultured on the SD medium lacking Ura with or without 200 ng/mL AbA. After 5 days of culture on a 30 °C incubator, the yeast strains harboring pABAi-HuADH1, pABAi-HuCYP76AD1–1, and pABAi-HuDODA1 plasmid sequences were verified via Seq Amp DNA Polymerase (TaKaRa, Shiga, Japan). Besides, the cDNA sequence of HuMYB132 was inserted into the pGADT7 vector with specific primers (Additional file 2). Then, the plasmid of pGADT7-HuMYB132 was transformed into the yeast strains containing pABAi-HuADH1, pABAi-HuCYP76AD1–1, and pABAi-HuDODA1 and cultured on the SD medium lacking Leu with or without 200 ng/mL AbA. After 5 days of culture on a 30 °C incubator, the abilities of HuMYB132 to bind HuADH1, HuCYP76AD1–1 and HuDODA1 promoters were evaluated by the growing status on SD/−Leu.
Dual luciferase reporter assay
The full-length of HuMYB132 was inserted into the effector vector while the DNA sequences of HuADH1, HuCYP76AD1–1, and HuDODA1 promoters were inserted into the reporter vector (Additional file 2). After individually transforming them into the A. tumefaciens strain GV3101 (pSoup), the effector and reporter were infiltrated together into N. benthamiana leaves. The ratio of LUC to Ren activity was measured using the Dual-Luciferase Reporter Assay System (Promega, USA) after 3 d. The experiments were repeated three biological replicates.
Protein purification and electrophoretic mobility shift assay
The coding region of HuMYB132 was fused with GST in the pGEX-4 T-1 vector and expressed in Escherichia coli strain BM Rosetta (DE3) (Additional file 2). The GST-HuMYB132 protein was further purified using Glutathione-Superflow Resin (Yeasen, China) (Additional file 1: Fig. S2). Oligonucleotide probes were synthesized and labelled with biotin at their 3′-ends (Sangon Biotech, China) (Additional file 1: Fig. S3). Electrophoretic mobility shift assay (EMSA) was conducted using the LightShift Chemiluminescent EMSA Kit (Thermo Scientific, USA). 50 ng GST-HuMYB132 protein was incubated with biotin-labeled probes with the addition of biotin-unlabeled probes, or mutant probes, or nothing. The shift bands of free and protein-DNA complexes were separated using a 6% native polyacrylamide gel. Then, the blots were cut after hybridisation with antibodies and transferred to a 5 × 9 cm nylon membrane. The blots were detected using the Chemiluminescent Biotin-labeled Nucleic Acid Detection Kit (Beyotime, China) and pictured by the ChemiDoc MP imaging system (Bio-Rad, USA) with the chemiluminescence method. Probe sequences used for EMSA assay are listed in Additional file 2.
Gene silencing assay of HuMYB132 in pulps of ‘Guanhuahong’ pitaya
The coding region of HuMYB132 was inserted into the pTRV2 vector and transformed into A. tumefaciens strain GV3101 (Additional file 2). The pTRV2-HuMYB132 and pTRV1 were co-infiltrated in the pulps of ‘Guanhuahong’ pitaya and the empty vector was used as a negative control according to the method of Xie et al. . Each pulp was injected with 1 mL of A. tumefaciens. 14 days after infiltration, pulps were frozen in liquid nitrogen immediately and stored at − 80 °C for betalain content measurement and gene expression analyses. The experiments were repeated three biological replicates.
cytochrome P450 enzyme
DNA binding domains
4,5-DOPA extradiol dioxygenase
electrophoretic mobility shift assay
weighted correlation network analysis
Chen J, Xie F, Cui Y, Chen C, Lu W, Hu X, et al. A chromosome-scale genome sequence of pitaya (Hylocereus undatus) provides novel insights into the genome evolution and regulation of betalain biosynthesis. Hortic Res. 2021;8(1):164. https://doi.org/10.1038/s41438-021-00612-0.
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. https://doi.org/10.1016/j.tplants.2010.06.005.
Xie F, Hua Q, Chen C, Zhang Z, Zhang R, Zhao J, et al. Genome-wide characterization of R2R3-MYB transcription factors in pitaya reveals a R2R3-MYB repressor HuMYB1 involved in fruit ripening through regulation of betalain biosynthesis by repressing betalain biosynthesis-related genes. Cells. 2021;10(8):1949. https://doi.org/10.3390/cells10081949.
Chen Y, Yang X, He K, Liu M, Li J, Gao Z, et al. The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2006;60(1):107–24. https://doi.org/10.1007/s11103-005-2910-y.
Stracke R, Werber M, Weisshaar B. The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2006;4(5):447–56. https://doi.org/10.1016/s1369-5266(00)00199-0.
Hatlestad G, Akhavan N, Sunnadeniya R, Elam L, Cargile S, Hembd A, et al. The beet Y locus encodes an anthocyanin MYB-like protein that activates the betalain red pigment pathway. Nat Genet. 2015;47(1):92–6. https://doi.org/10.1038/ng.3163.
Ibrahim S, Mohamed G, Khedr A, Zayed M, El-Kholy A. Genus Hylocereus: beneficial phytochemicals, nutritional importance, and biological relevance-a review. J Food Biochem. 2018;42(2):e12491. https://doi.org/10.1111/jfbc.12491.
Timoneda A, Feng T, Sheehan H, Walker-Hale N, Pucker B, Lopez-Nieves S, et al. The evolution of betalain biosynthesis in Caryophyllales. New Phytol. 2019;224(1):71–85. https://doi.org/10.1111/nph.15980.
Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54(4):733–49. https://doi.org/10.1111/j.1365-313X.2008.03447.x.
Lopez-Nieves S, Yang Y, Timoneda A, Wang M, Feng T, Smith S, et al. Relaxation of tyrosine pathway regulation underlies the evolution of betalain pigmentation in Caryophyllales. New Phytol. 2018;217(2):896–908. https://doi.org/10.1111/nph.14822.
Chen Q, Man C, Li D, Tan H, Xie Y, Huang J. Arogenate dehydratase isoforms differentially regulate anthocyanin biosynthesis in Arabidopsis thaliana. Mol Plant. 2016;9(12):1609–19. https://doi.org/10.1016/j.molp.2016.09.010.
Moghe G, Smith S. The push and pull of plant specialized metabolism underlies a long-standing, colorful mystery. New Phytol. 2018;217(2):471–3. https://doi.org/10.1111/nph.14914.
Gandía-Herrero F, García-Carmona F. Characterization of recombinant Beta vulgaris 4,5-DOPA-extradiol-dioxygenase active in the biosynthesis of betalains. Planta. 2012;236(1):91–100. https://doi.org/10.1007/s00425-012-1593-2.
Hatlestad G, Sunnadeniya R, Akhavan N, Gonzalez A, Goldman I, McGrath J, et al. The beet R locus encodes a new cytochrome P450 required for red betalain production. Nat Genet. 2012;44(7):816–20. https://doi.org/10.1038/ng.2297.
Heuer S, Vogt T, Boehm H, Strack D. Partial purification and characterization of UDP-glucose: Betanidin 5-O- and 6-O-glucosyltransferases from cell suspension cultures of Dorotheanthus bellidiformis (Burm. f.) N.E.Br. Planta. 1996;199:244–50. https://doi.org/10.1007/BF00196565.
Sasaki N, Abe Y, Wada K, Koda T, Goda Y, Adachi T, et al. Amaranthin in feather cockscombs is synthesized via glucuronylation at the cyclo-DOPA glucoside step in the betacyanin biosynthetic pathway. J Plant Res. 2005;118(6):439–42. https://doi.org/10.1007/s10265-005-0237-z.
Polturak G, Heinig U, Grossman N, Battat M, Leshkowitz D, Malitsky S, et al. Transcriptome and metabolic profiling provides insights into betalain biosynthesis and evolution in Mirabilis Jalapa. Mol Plant. 2018;11(1):189–204. https://doi.org/10.1016/j.molp.2017.12.002.
Imamura T, Takagi H, Miyazato A, Ohki S, Mizukoshi H, Mori M. Isolation and characterization of the betalain biosynthesis gene involved in hypocotyl pigmentation of the allotetraploid Chenopodium quinoa. Biochem Biophys Res Commun. 2018;496(2):280–6. https://doi.org/10.1016/j.bbrc.2018.01.041.
Hua Q, Chen C, Chen Z, Chen P, Ma Y, Wu J, et al. Transcriptomic analysis reveals key genes related to betalain biosynthesis in pulp coloration of Hylocereus polyrhizus. Front Plant Sci. 2016;6:1179. https://doi.org/10.3389/fpls.2015.01179.
Hua Q, Zhou Q, Gan S, Wu J, Chen C, Li J, et al. Proteomic analysis of Hylocereus polyrhizus reveals metabolic pathway changes. Int J Mol Sci. 2016;17(10):1606. https://doi.org/10.3390/ijms17101606.
Suh D, Lee S, Heo D, Kim Y, Cho S, Lee S, et al. Metabolite profiling of red and white pitayas (Hylocereus polyrhizus and Hylocereus undatus) for comparing betalain biosynthesis and antioxidant activity. J Agric Food Chem. 2014;62(34):8764–71. https://doi.org/10.1021/jf5020704.
Chen C, Xie F, Hua Q, Tel-Zur N, Zhang L, Zhang Z, et al. Integrated sRNAome and RNA-Seq analysis reveals miRNA effects on betalain biosynthesis in pitaya. BMC Plant Biol. 2020;20(1):437. https://doi.org/10.1186/s12870-020-02622-x.
Xie F, Hua Q, Chen C, Zhang L, Zhang Z, Chen J, et al. Transcriptomics-based identification and characterization of glucosyltransferases involved in betalain biosynthesis in Hylocereus megalanthus. Plant Physiol Biochem. 2020;152:112–24. https://doi.org/10.1016/j.plaphy.2020.04.023.
Zhang X, Ju H, Chung M, Huang P, Ahn S, Kim C. The R-R-type MYB-like transcription factor, AtMYBL, is involved in promoting leaf senescence and modulates an abiotic stress response in Arabidopsis. Plant Cell Physiol. 2011;52(1):138–48. https://doi.org/10.1093/pcp/pcq180.
Lopez-Nieves S, El-Azaz J, Men Y, Holland C, Feng T, Brockington S, et al. Two independently evolved natural mutations additively deregulate TyrA enzymes and boost tyrosine production in planta. Plant J. 2022;109(4):844–55. https://doi.org/10.1111/tpj.15597.
Brockington S, Yang Y, Gandia-Herrero F, Covshoff S, Hibberd J, Sage R, et al. Lineage-specific gene radiations underlie the evolution of novel betalain pigmentation in Caryophyllales. New Phytol. 2015;207(4):1170–80. https://doi.org/10.1111/nph.13441.
Bean A, Sunnadeniya R, Akhavan N, Campbell A, Brown M, Lloyd A. Gain-of-function mutations in beet DODA2 identify key residues for betalain pigment evolution. New Phytol. 2018;219(1):287–96. https://doi.org/10.1111/nph.15159.
Polturak G, Breitel D, Grossman N, Sarrion-Perdigones A, Weithorn E, Pliner M, et al. Elucidation of the first committed step in betalain biosynthesis enables the heterologous engineering of betalain pigments in plants. New Phytol. 2016;210(1):269–83. https://doi.org/10.1111/nph.13796.
Sheehan H, Feng T, Walker-Hale N, Lopez-Nieves S, Pucker B, Guo R, et al. Evolution of L-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales. New Phytol. 2020;227(3):914–29. https://doi.org/10.1111/nph.16089.
Hua Q, Chen C, Xie F, Zhang Z, Zhang R, Zhao J, et al. A genome-wide identification study reveals that HmoCYP76AD1, HmoDODAα1 and HmocDOPA5GT involved in betalain biosynthesis in Hylocereus. Genes (Basel). 2021;12(12):1858. https://doi.org/10.3390/genes12121858.
Chang Y, Chiu Y, Tsao N, Chou Y, Tan C, Chiang Y, et al. Elucidation of the core betalain biosynthesis pathway in Amaranthus tricolor. Sci Rep. 2021;11(1):6086. https://doi.org/10.1038/s41598-021-85486-x.
Rose A, Meier I, Wienand U. The tomato I-box binding factor LeMYBI is a member of a novel class of myb-like proteins. Plant J. 1999;20(6):641–52. https://doi.org/10.1046/j.1365-313x.1999.00638.x.
Meneses E, Cárdenas H, Zárate S, Brieba L, Orozco E, López-Camarillo C, et al. The R2R3 Myb protein family in Entamoeba histolytica. Gene. 2010;455(1–2):32–42. https://doi.org/10.1016/j.gene.2010.02.004.
Wang Y, Liu W, Jiang H, Mao Z, Wang N, Jiang S, et al. The R2R3-MYB transcription factor MdMYB24-like is involved in methyl jasmonate-induced anthocyanin biosynthesis in apple. Plant Physiol Biochem. 2019;139:273–82. https://doi.org/10.1016/j.plaphy.2019.03.031.
Huang D, Tang Z, Fu J, Yuan Y, Deng X, Xu Q. CsMYB3 and CsRuby1 form an ‘activator-and-repressor’ loop for the regulation of anthocyanin biosynthesis in citrus. Plant Cell Physiol. 2020;61(2):318–30. https://doi.org/10.1093/pcp/pcz198.
Cheng M, Huang Z, Hua Q, Shan W, Kuang J, Lu W, et al. The WRKY transcription factor HpWRKY44 regulates CytP450-like1 expression in red pitaya fruit (Hylocereus polyrhizus). Hortic Res. 2017;4:17039. https://doi.org/10.1038/hortres.2017.39.
Zhang L, Chen C, Xie F, Hua Q, Zhang Z, Zhang R, et al. A novel WRKY transcription factor HmoWRKY40 associated with betalain biosynthesis in pitaya (Hylocereus monacanthus) through regulating HmoCYP76AD1. Int J Mol Sci. 2021;22(4):2171. https://doi.org/10.3390/ijms22042171.
Guo C, Yao L, You C, Wang S, Cui J, Ge X, et al. MID1 plays an important role in response to drought stress during reproductive development. Plant J. 2018;88(2):280–93. https://doi.org/10.1111/tpj.13250.
Fang Q, Wang Q, Mao H, Xu J, Wang Y, Hu H, et al. AtDIV2, an R-R-type MYB transcription factor of Arabidopsis, negatively regulates salt stress by modulating ABA signaling. Plant Cell Rep. 2018;37(11):1499–511. https://doi.org/10.1007/s00299-018-2321-6.
Agrahari R, Enomoto T, Ito H, Nakano Y, Yanase E, Watanabe T, et al. Expression GWAS of PGIP1 identifies STOP1-dependent and STOP1-independent regulation of PGIP1 in aluminum stress signaling in Arabidopsis. Front Plant Sci. 2021;12:774687. https://doi.org/10.3389/fpls.2021.774687.
Li G, Meng X, Zhu M, Li Z. Research progress of betalain in response to adverse stresses and evolutionary relationship compared with anthocyanin. Molecules. 2019;24(17):3078. https://doi.org/10.3390/molecules24173078.
Chen C, Wu J, Hua Q, Tel-Zur N, Xie F, Zhang Z, et al. Identification of reliable reference genes for quantitative real-time PCR normalization in pitaya. Plant Methods. 2019;15:70. https://doi.org/10.1186/s13007-019-0455-3.
Zhao J, Chen L, Ma A, Wang D, Lu H, Chen L, et al. R3-MYB transcription factor LcMYBx from Litchi chinensis negatively regulates anthocyanin biosynthesis by ectopic expression in tobacco. Gene. 2022;812:146105. https://doi.org/10.1016/j.gene.2021.146105.
Lai B, Du L, Liu R, Hu B, Su W, Qin Y, et al. Two LcbHLH transcription factors interacting with LcMYB1 in regulating late structural genes of anthocyanin biosynthesis in Nicotiana and Litchi chinensis during anthocyanin accumulation. Front Plant Sci. 2016;7:166. https://doi.org/10.3389/fpls.2016.00166.
The authors would also like to thank fellow members of the orchard of Jinsuinong for management of experimental materials.
This work was supported by the National Natural Science Foundation of China (31972367 and 31960578), Provincial Rural Revitalization Strategy Special Project of Guangdong in 2022 (No. 61), Science and Technology Program of Guangzhou (202002020060) and Yangjiang (Yangketong50).
Availability of data and materials
The genome data of H. undatus were deposited in the NCBI’s SRA (Sequence Read Archive) data resource with the Bioproject ID PRJNA691451 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA691451), while the RNA-Seq data of ‘Guanhuahong’ pitaya pulps on the 17th, 23rd, 25th, and 32nd DAF were the Bioproject ID PRJNA704510 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA704510). The coding sequences of HuADH1, HuCYP76AD1–1, HuDODA1 and HuMYB132, as well as their promoter sequences, are listed in the Supplementary information 3.
FFX, JYC (Jianye Chen) and YHQ conceived the study, and designed the experiments. FFX conducted the experiment. FFX, CBC, JYC (Jiayi Chen), JXC, QZH, ZKZ, JTZ and GBH analyzed data. FFX, KS and YHQ wrote the manuscript. All Authors read and approved the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original images of RT-PCR and EMSA assays. Fig. S1. Full-length gel of Fig. 1B. Fig. S2. The SDS-PAGE gel stained with coomassie brilliant blue, presenting affinity purification of the GST (A) amd GST-HuMYB132 (B) proteins. Fig. S3. The unprocessed blots of the biotin probe of HuADH1, HuCYP76AD1–1 and HuDODA1 promoters. Fig. S4-S6. The original images of HuMYB132 binding to the HuADH1 (Fig. S4), HuCYP76AD1–1 (Fig. S5) and HuDODA1 (Fig. S6) promoters in EMSA assay.
The primers used in this study.
Coding sequences of HuADH1, HuCYP76AD1–1, HuDODA1 and HuMYB132, and promoter sequences of HuADH1, HuCYP76AD1–1 and HuDODA1.
About this article
Cite this article
Xie, F., Chen, C., Chen, J. et al. Betalain biosynthesis in red pulp pitaya is regulated via HuMYB132: a R-R type MYB transcription factor. BMC Plant Biol 23, 28 (2023). https://doi.org/10.1186/s12870-023-04049-6