Next Article in Journal
Enzymatic Spermine Metabolites Induce Apoptosis Associated with Increase of p53, caspase-3 and miR-34a in Both Neuroblastoma Cells, SJNKP and the N-Myc-Amplified Form IMR5
Next Article in Special Issue
Role of Glycoproteins during Fruit Ripening and Seed Development
Previous Article in Journal
2′O-Ribose Methylation of Ribosomal RNAs: Natural Diversity in Living Organisms, Biological Processes, and Diseases
Previous Article in Special Issue
Characterization of RIN Isoforms and Their Expression in Tomato Fruit Ripening
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 10
Issue 8
10.3390/cells10081949
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Fangfang Xie
,
Qingzhu Hua
,
Canbin Chen
,
Zhike Zhang
,
Rong Zhang
,
Jietang Zhao
,
Guibing Hu
,
Jianye Chen
* and
Yonghua Qin
*
State Key Laboratory for Conservation and Utilization of Subtropical Agrobioresources/Guangdong Provincial Key Laboratory of Postharvest Science of Fruits and Vegetables/Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture and Rural Affairs, College of Horticulture, South China Agricultural University, Guangzhou 510642, China
*
Authors to whom correspondence should be addressed.
Cells 2021, 10(8), 1949; https://0-doi-org.brum.beds.ac.uk/10.3390/cells10081949
Submission received: 14 June 2021 / Revised: 26 July 2021 / Accepted: 28 July 2021 / Published: 31 July 2021
(This article belongs to the Special Issue Molecular Regulation of Fruit-Set, Development, Maturation, Ripening and Shelf-Life)
Browse Figures
Versions Notes

Abstract

:
The MYB (myeloblastosis) superfamily constitutes one of the most abundant transcription factors (TFs) regulating various biological processes in plants. However, the molecular characteristics and functions of MYB TFs in pitaya remain unclear. To date, no genome-wide characterization analysis of this gene family has been conducted in the Cactaceae species. In this study, 105 R2R3-MYB members were identified from the genome data of Hylocereus undatus and their conserved motifs, physiological and biochemical characteristics, chromosome locations, synteny relationship, gene structure and phylogeny were further analyzed. Expression analyses suggested that three up-regulated HuMYBs and twenty-two down-regulated HuMYBs were probably involved in fruit ripening of pitaya. Phylogenetic analyses of R2R3-MYB repressors showed that seven HuMYBs (HuMYB1, HuMYB21, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101) were in clades containing R2R3-MYB repressors. HuMYB1 and HuMYB21 were significantly down-regulated with the betalain accumulation during fruit ripening of ‘Guanhuahong’ pitaya (H. monacanthus). However, only HuMYB1 had R2 and R3 repeats with C1, C2, C3 and C4 motifs. HuMYB1 was localized exclusively to the nucleus and exhibited transcriptional inhibition capacities. Dual luciferase reporter assay demonstrated that HuMYB1 inhibited the expression of betalain-related genes: HuADH1, HuCYP76AD1-1 and HuDODA1. These results suggested that HuMYB1 is a potential repressor of betalain biosynthesis during pitaya fruit ripening. Our results provide the first genome-wide analyses of the R2R3-MYB subfamily involved in pitaya betalain biosynthesis and will facilitate functional analysis of this gene family in the future.

1. Introduction

Transcription factors (TFs) play important roles in plant growth and development by controlling gene expression, suppressing translation or activating transcription of target genes, or interacting with other proteins. Generally, TFs are composed of at least four discrete domains: DNA-binding domain, nuclear localization signal, transcription-activation domain and oligomerization site [1]. The MYB family is one of the largest TF families known in eukaryotes and is defined by the presence of one to four highly conserved MYB repeats. The MYB repeat contains approximately 52 amino acids with three regularly spaced tryptophan, encoding three α-helices. The third helix of each repeat plays a recognition role in direct binding to DNA sequence and intercalating in the major groove [2]. Based on the MYB repeat number and the identity of the MYB repeats, MYB superfamily members are generally classified into four major subfamilies, namely, MYB-related (1R-MYB), R2R3-MYB (2R-MYB), R1R2R3-MYB (3R-MYB) and 4R-like MYB (4R-MYB) proteins. All four subfamilies are found in plants, representing the taxon with the highest diversity of MYB proteins. The smallest subfamily is the 4R-MYB group, whose members contain four R1/R2-like repeats in several plant genomes [3]. The 3R-MYB subfamily contains R1, R2 and R3 repeats and typically encoded by five genes in higher plant genomes. The 1R-MYB subfamily comprises proteins with a single or a partial MYB repeat and falls into several subclasses, including R3-type, R1/R2-type and the GARP family [4].
Most plant MYB genes encode proteins of the 2R-MYB class, which contain an R2 and an R3 domain at N terminal and an activation or repression motif at the C terminus. The R2R3-MYB subfamily probably evolved from an R1-MYB genes through duplication of an R1 repeat or from an R1R2R3-MYB gene through loss of the R1 repeat [4,5]. Subsequently, genome-wide duplication and tandem duplication events play important roles in the rapid expansion of R2R3-MYB members and its subfamily/clade-asymmetric and lineage-specific in land plants [6,7]. The first MYB gene, named v-MYB, was identified from avian myeloblastosis virus (AMV) [8]. Meanwhile, the first plant MYB gene, named C1, was isolated from Zea mays and encoded a c-MYB-like TF involved in anthocyanin pigmentation [9]. With the growing number of fully sequenced plant genomes, great progress has been made in identification of R2R3-MYB genes in recent years. Based on their well conserved MYB domains, R2R3-MYB families have been annotated genome-wide in many plant such as Arabidopsis (Arabidopsis thaliana, 126 members) [3], grape (Vitis vinifera, 117 members) [10], apple (Malus × domestica, 222 members) [11], kiwifruit (Actinidia chinensis, 93 members) [12], pineapple (Ananas comosus, 87 members) [13], watermelon (Citrullus lanatus, 162 members) [14], banana (Musa acuminata, 285 members) [15] and beet (Beta vulgaris, 70 members) [16]. These R2R3-MYB subfamilies are involved in the regulation of plant development, cell differentiation, primary and secondary metabolism and stress responses. R2R3-MYBs associated with flavonoid production or suppression can be identified by phylogenetic analyses. In Arabidopsis, the members of subgroup 4–7 could regulate the flavonoid biosynthesis [3]. BvMYB1 was found to activate the expression levels of CYP76AD1 and DODA1 resulting in betalain accumulation in beet [17]. However, BvMYB1 is clustered in subgroup 6 responsible for anthocyanin biosynthesis [7]. Moreover, there is no study about the R2R3-MYB subfamily in pitaya and no information is available about repressor involved in betalain biosynthesis.
Pitaya, belonging to Hylocereus in the Cactaceae family of Caryophyllales, is mainly classified into three groups including H. undatus (red peel with white pulp), H. monacanthus or H. polyrhizus (red peel with red pulp) and H. megalanthus or Selenicereus megalanthus (yellow peel with white pulp) [18]. Pitaya is the only edible fruit with abundant betalains that has been commercially produced at a large scale. Betalains are vacuole-localized, water-soluble and nitrogen-containing pigment derived from the L-tyrosine and divided into red-violet betacyanins and orange-yellow betaxanthins [19]. Four key enzymes are involved in betalain biosynthesis: (i) ADH (arogenate dehydrogenase), catalyzing the tyrosine (Tyr) biosynthesis [20]; (ii) CYP76AD1 (cytochrome P450), catalyzing the hydroxylation of tyrosine to DOPA and the subsequent oxidation of DOPA to cyclo-DOPA [21]; (iii) DODA (4,5-DOPA dioxygenase extradiol), catalyzing the L-DOPA to seco-DOPA and then quickly converting to betalamic acid spontaneously [22]; and (iv) GT (glucosyltransferase), catalyzing cyclo-DOPA to betanin [23]. As important natural colorants, betalains provide nutrients and health benefits for humans by increasing the antioxidant properties [24,25]. Betalain accumulation is an important index and parameter for measuring harvest time and fruit quality in pitaya. The color of pitaya peels and pulps is determined by betalain biosynthesis and accumulation. In this study, the genome-wide identification and characterization of R2R3-MYB genes were performed to screen candidate genes related to betalain biosynthesis using available genomic information. The aim of the present study is to identify and validate the R2R3-MYB involved in pitaya betalain biosynthesis and candidate R2R3-MYB genes can be used for genetic improvement of pitaya.

2. Materials and Methods

2.1. Identification and Sequence Analyses of Pitaya R2R3-MYB Subfamily

Plot hidden Markov model (HMM) profile of MYB DNA-binding domain (PF00249) was downloaded from Pfam (https://pfam.xfam.org/, accessed on 8 September 2019) [26] to identify MYB genes from pitaya genome [27]. To ensure the presence of the core MYB domains, the putative MYB sequences were further examined using Pfam and SMART (https://smart.embl-heidelberg.de/) (accessed on 27 March 2020). The sequences with two MYB domains were considered to be the final R2R3-MYB gene family members.
The DNA-binding domains of R2R3-MYB proteins were aligned by Cluster X software and shown by WEBLOGO online tool (https://weblogo.berkeley.edu/logo.cgi) (accessed on 27 March 2020). The conserved motifs of complete amino acid sequences were analyzed by Multiple Em for Motif Elicitation (https://meme-suite.org/) (accessed on 26 March 2020). The exon/intron organizations of 105 pitaya R2R3-MYB genes were drawn by TBtools software [28]. The physiological and biochemical characteristics of pitaya R2R3-MYB proteins were analyzed using the ExPASy (https://www.expasy.org/) (accessed on 27 March 2020).

2.2. Chromosomal Locations and Synteny Analyses of Pitaya R2R3-MYB Members

The MapChart software was used to draw the location sites of HuMYBs according to their physical positions. The duplication pattern for each HuMYBs was analyzed by Multiple Collonearity Scan toolkit (MCscanX) software following the operation manual [29]. The synteny relationship of orthologous R2R3-MYB genes between pitaya, beet and Arabidopsis was shown by TBtools software [28].

2.3. Phylogenetic Analyses of Pitaya R2R3-MYB Members

The evolutionary trees of pitaya, beet and Arabidopsis R2R3-MYB proteins were constructed using the Maximum likelihood method (ML) in MEGA 7 with 1000 bootstrap replicates. The phylogenetic trees of 105 pitaya R2R3-MYB proteins and R2R3-MYB repressors were performed by Neighbor-joining method (NJ) method. All phylogenetic trees were exhibited by EVOLVIEW online tool (https://www.evolgenius.info/evolview/) (accessed on 2 April 2020).

2.4. Plant Materials

‘Guanhuahong’ pitaya (red peel and green scales with red pulp, H. monacanthus) from the orchard of Jinsuinong (Zhongluotan Village, Guangzhou City, Guangdong Province, China) were used as materials. Fruits of ‘Guanhuahong’ pitaya were 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 ℃ until future analyses.

2.5. Measurement of Betalain Contents

Betalains were extracted and measured following our previously described method [30]. In brief, 0.5 g freeze-dried pitaya pulp powder were extracted with 5 mL 80% aqueous methanol (v/v) solution by sonication for 10 min and stirred for 20 min in dark at room temperature. After centrifuging at 5000× g for 15 min, the residues were subjected to a similar second extraction. The supernatants were measured through spectrophotometry (Infinite M200, Tecan Co., Ltd, Shanghai, China) at 478 nm for betaxanthins and 538 nm for betacyanins. All determinations were performed in three biological repetition with three technical replicates.

2.6. Gene Expression Analyses and Cloning

Transcriptome data of four fruit developmental stages (17th, 23rd, 25th and 32nd) (PRJNA704510) were used to draw the transcript abundance of 2R-MYB repressors by TBtools software. The primers were designed in NCBI (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/tools/primer-blast/, accessed on 19 May 2020). Total RNA was isolated using the EASYspin Plus Complex Plant RNA Kit (RN53) (Aidlab Biotechnology, Beijing, China) according to the manufacturer’s protocol. Single-stranded cDNA was synthesized using the PrimeScript™ RT Reagent Kit with gDNA Eraser (TaKaRa, Shiga, Japan). qRT-PCR was performed with an CFX384-Real-Time System (C1000 Touch Thermal Cycler, Bio-Rad, Irvine, CA, USA) using the RealUniversal Color PreMix (SYBR Green) (TIANGEN, Beijing, China) with specific primers (Table S1). The Actin (1) was used as internal control for gene expression analyses [31]. All experiments were repeated in technical replicates. The relative expression levels of each transcript were calculated using the comparative Ct method [32].
The full-length coding sequences of HuMYBs were cloned using I-5TM 2× High-Fidelity Master Mix (MCLAB, San Francisco, CA, USA) with specific primers (Table S1). Alignment of the deduced protein sequences was performed by Clustal X and GeneDoc software.

2.7. Subcellular Localization Analyses

Full-length coding sequences of HuMYB1 were inserted into the pGreen-35S-GFP vector (primers are listed in Table S1). Subsequently, pGreen-35S-HuMYB1-GFP and pGreen-35S-GFP were separately transformed into Agrobacterium tumefaciens strain GV3101 (pSoup-p19) (Coolaber, Beijing, China) and infiltrated into Nicotiana benthamiana leaves. Two days after infiltration, leaf protoplasts were isolated according to Lai et al. [33] and the transient expression of GFP was observed using a fluorescence microscope (ZEISS LCM-800, Carl Zeiss, Oberkochen, Germany).

2.8. Transcriptional Activation Analyses in Yeast Cells

Full-length of HuMYB1 were inserted into the pGBKT7 vector (primers are listed in Table S1). Then, the pGBKT7-HuMYB1, positive (pGBKT7-53 + pGADT7-T-antigen) and negative (pGBKT7) controls were separately transformed into yeast strain GoldY2H (Coolaber, Beijing, China) and grown on minimal medium plates without tryptophan (SD/-Trp) or without tryptophan, histidine and adenine (SD/-Trp-His-Ade). The transactivation activity of HuMYB1 protein was evaluated according to their growth status after 3 days in 30 ℃ and confirmed by incubating with X-α-galactosidase (X-α-Gal) for 30 min.

2.9. Dual-Luciferase Transient Expression Assay

DNA was isolated according to the protocol of Plant DNA Extraction Kit (DN14) (Aidlab, Beijing, China) and RNA was removed with Ribonuclease A (RNase A) (Takara Biomedical Technology Co., Ltd, Beijing, China). Promoter sequences were cloned using Seq Amp DNA Polymerase (Takara Biomedical Technology Co., Ltd, Beijing, China) with specific primers (Table S1) and analyzed using Plant-Care (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) (accessed on 23 April 2020).
The dual-luciferase transient expression system was conducted according to the previous study [34]. To assess the transactivation activity of HuMYB1, full-length coding sequences of HuMYB1 was inserted into the pGreenⅡ BD-62-SK as effectors. The reporter vector was modified from the pGreenII 0800-LUC vector [33]. Then, effectors of pGreenⅡ BD-62-SK-HuMYB1, pGreenⅡ BD-62-SK, pGreenⅡ BD-62-SK-VP16 (positive control) and reporter were separately transformed into Agrobacterium tumefaciens strain GV3101. The effector and reporter were infiltrated into N. benthamiana leaves with a mixture of 9:1, and LUC and REN activities were measured by dual-luciferase assay kit (Promega, Madison, WI, USA) after 72 h.
To assess the transcriptional effects of HuMYB1 on the promoters of ADH1, CYP76AD1-1 and DODA1, these promoters were inserted into the pGreenII 0800-LUC vector as reporters, while HuMYB1 were inserted into the pGreenII 62-SK vector as effectors. The reporter and effector were transformed into Agrobacterium tumefaciens strain GV3101 and then infiltrated into N. benthamiana leaves. LUC and REN activities were analyzed as above.

3. Results

3.1. Identification of R2R3-MYB Genes from H. undatus

The pitaya MYB encoding genes were identified by the alignment of the pitaya genome and MYB domain HMM profile (PF00249) using BLASTP software. A total of 231 candidate deduced amino acid sequences containing MYB or MYB-like repeats were obtained. The MYB domains were subsequently analyzed by Pfam and SMART. In total, 75 1R-MYB proteins, 105 R2R3-MYB proteins, four R1R2R3-MYB proteins and one 4R-MYB protein were obtained from the H. undatus genome (Table S2). Besides this, the theoretical pI and molecular weight ranged from 4.26 kDa (HuMYB137) to 10.96 kDa (HuMYB171) and 8.91 kDa (HuMYB143) to 213.32 kDa (HuMYB154), respectively (Table S2).
To investigate the homologous domain sequence features, conservation and divergence of 105, 70 and 126 R2 and R3 repeats from H. undatus (betalain-producing plant), B. vulgaris (betalain-producing plant) and A. thaliana (anthocyanin-producing plant) were assessed using multiple alignment analyses (Figure 1). The results showed that the basic regions of pitaya, beet and Arabidopsis MYB domains (about 104 amino acid residues) were rarely deleted or inserted, as previously reported in the other plants [35,36], while they were widely divergent in the length and amino acid outside of the MYB domains. However, the length of Chinese cabbage and rice were larger than that of Arabidopsis based on the space between neighboring tryptophan (Trp, W) residues [37]. In pitaya, six highly conserved triplet W residues were located in R2 repeat (position 6, 26 and 46), R3 repeat (position 78 and 97) (Figure 1A,B) and the first Trp residue in R3 repeat (position 59) was replaced by a hydrophobic amino acid Phe (F). These conserved W (F) residues are indispensable in maintaining the helix-turn-helix structure of MYB domains [38]. Similar localization existed in both beet and Arabidopsis, suggesting the evolutionary conservation of MYBs among plants (Figure 1C,F). Besides this, the 3’ region of R2 domain in pitaya R2R3-MYBs contained a highly conserved LRPD motif. Compared to Arabidopsis, the positions at 12, 19, 25, 33 and 36 in R2 repeats and 69, 80 and 83 in R3 repeats showed difference in pitaya (Figure 1). Difference was detected between pitaya and beet in the position 69 and 80 in R3 repeats, indicating the divergence of the MYB domain between pitaya and beet was less than that of pitaya and Arabidopsis.

3.2. Genomic Location of HuMYBs

Genome chromosomal location analyses showed that all the 105 R2R3-MYB genes were distributed throughout all 11 chromosomes (Chrs) in the pitaya genome (Figure 2) and renamed as HuMYB1 to HuMYB105 according to their location in 11 Chrs. The largest number of MYB genes (fifteen) mapped on Chr 6, followed by thirteen on Chr 5. However, there were only three HuMYBs on Chr 9. In general, MYB genes are absent in the central sections of chromosomes including centromere and pericentromeric regions. Relatively high densities of HuMYB genes were observed at the end of Chr 3, 6 and 7 as well as the bottom of Chr 1, 4, 5, 6, 8 and 11 (Figure 2). The same uneven distributions in all chromosomes were also reported in the other plants, such as A. comosus [13], Brassica rapa [37] and Eucalyptus grandis [39].

3.3. Synteny Analyses of Pitaya R2R3-MYB Subfamily

Gene duplication is recognized to occur throughout plant evolution and plays an important role in expanding the large gene families in plants [40]. The duplication patterns and synteny relationship of the HuMYBs were analyzed using BLASTP and MCScanX methods. In total, 40 segmental duplication events with 61 HuMYBs happened in the pitaya genome (Figure 3A; Table S3). HuMYBs were located within synteny blocks on almost all chromosomes except for Chr 9, indicating that HuMYBs underwent through expansion during genome evolution. If two or more MYB genes resided within 20 kb in the same chromosome, a tandem duplication event was defined [41]. In this study, three HuMYB tandem duplication pairs happened in Chr 5 (HuMYB48 and HuMYB49), Chr 6 (HuMYB67 and HuMYB68) and Chr 10 (HuMYB93 and HuMYB94). These results suggested that large-scale segmental duplication events were the major contributors to the expansion of the R2R3-MYB gene family in H. undatus.
To further investigate the potential evolutionary mechanism of R2R3-MYB subfamily, we constructed the comparative syntenic maps of H. undatus, B. vulgaris and A. thaliana in pairs (Figure 3B). Finally, 46 collinear R2R3-MYB gene pairs between pitaya and beet (Figure 3B(a); Table S4-1), 28 between pitaya and Arabidopsis (Figure 3B(b); Table S4-2) and 10 between beet and Arabidopsis (Figure 3B(c); Table S4-3) were identified. The number of orthologous events of HuMYB-BvMYB was much more than that of HuMYB-AtMYB and BvMYB-AtMYB, indicating that the evolutionary distance between pitaya and beet was closer than that of pitaya and Arabidopsis. In addition, the area of all synteny blocks between pitaya and beet was far greater than that of pitaya and Arabidopsis. These results suggested that the evolutionary distance between pitaya and beet is closer than that of pitaya and Arabidopsis.

3.4. Phylogenetic, Gene Structure and Motif Composition Analyses of R2R3-MYBs

To explore the putative function of the 105 HuMYBs, functional groups were performed using MYB genes from B. vulgaris and A. thaliana. A. thaliana was chosen because MYB genes have been extensively studied in this model plant while B. vulgaris was a vital plant material for betalain studies. Most MYB proteins shared similar functions cluster in the same phylogenetic clades, suggesting that most closely-related MYBs could recognize similar target genes and possess redundant, overlapping and/or cooperative functions [42]. By the ML method, an evolutionary tree of 105 HuMYBs, 70 BvMYBs and 126 AtMYBs was constructed (Figure 4). In the classification of the R2R3-MYB genes, the subgroup (S) categories from A. thaliana were also labeled in the phylogenetic tree [7]. According to the phylogenetic tree topology, 105 HuMYB proteins were classified into nine groups (designated as Group 1 to Group 9). HuMYBs were distributed with AtMYBs and BvMYBs in most subgroups. However, species-specific AtMYBs were also detected in S10, S12 and S15 (clades without HuMYBs and BvMYBs).
HpWRKY44 and BvMYB1 are involved in betalain biosynthesis but they share high homology with anthocyanin regulated TFs including AtWRKY44, AtMYB75, AtMYB90, AtMYB113 and AtMYB114 [17,43]. Thus, there is probably a link of transcription regulation between betalains and anthocyanins. Ten R2R3-MYB genes in Group 1 attracted our attention due to their gathering in S4 to S7 from Arabidopsis and Bv_ihfg (Figure 4). Particularly, HuMYB12 was clustered with Bv_iogq and AtMYB111 in S7, which controlled flavonol biosynthesis [3,16]. HuMYB83 was close to Bv_ralf, Bv_jkkr and AtMYB90 in the S6 group, which is involved in betalain and anthocyanin biosynthesis [3,16,17]. HuMYB26 was near Bv_crae and AtMYB123 in S5, responsible for proanthocyanidins (PAs) biosynthesis in the seed coat of Arabidopsis [3]. In addition, five genes, i.e., HuMYB1, HuMYB21, HuMYB72, HuMYB78 and HuMYB101, were gathered in S4 which encoded transcriptional repressors [3]. HuMYB48 and HuMYB49 were grouped with Bv_ihfg containing negative flavonoid regulator FaMYB1 [16]. These results suggested that HuMYB12, HuMYB26 and HuMYB83 were candidate R2R3-MYB activators while HuMYB1, HuMYB21, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101 were candidate transcriptional repressors.
Phylogenetic reconstruction of all 105 HuMYBs proteins were performed by the NJ method with 1000 bootstraps in MEGA7.0 (Figure 5A). The 105 HuMYBs were divided into nine groups, which was consistent with the result of the phylogenetic tree constructed by ML method in Figure 4. To understand the distribution of conserved motifs, the 105 R2R3-MYB protein sequences were analyzed by MEME software with five conserved motifs (named as motif 1 to 5) (Figure 5B and Figure S1). A total of 105 HuMYBs contained motifs 1, 2 and 3 which constituted R2 and R3 domains in the N-terminal, except for HuMYB42 in Group 6 and HuMYB9, HuMYB14, HuMYB38, HuMYB43, HuMYB75 and HuMYB89 in Group 9, which were found in the central region. Besides this, motif 4 and motif 5 were respectively conserved in Group 4 and Group 5, which is consistent with the phylogeny analyses of 105 HuMYBs proteins.
The exon–intron structures of 105 R2R3-MYB coding sequences were analyzed (Figure 5C). Most R2R3-MYBs were clustered in the same group with the similar exon-intron structures, particularly with the same number of introns, such as Group 2 and Group 7. Except for six genes (HuMYB2, HuMYB4, HuMYB20, HuMYB46, HuMYB52 and HuMYB55) in Group 9 without any intron, most of HuMYBs were disrupted by introns. In addition, more than 73% HuMYBs had three exons and two introns, which was also found in other plants [13,14].

3.5. Betalain Measurement and Expression Analyses of R2R3-MYBs during Pitaya Fruit Ripening

Betacyanin and betaxanthin contents were measured in seven ‘Guanhuahong’ pitaya fruit from 14th d to 32nd DAF (Figure 6A). The contents of betacyanins and betaxanthins in pulps increased throughout pitaya fruit ripening (Figure 6B). Betacyanin contents were increased gradually during fruit ripening and reached their maximum levels at the fully mature stage. ADHα, CYP76AD1α, DODAα and cDOPA5GT were essential structure genes involved in betalain biosynthesis (Figure 6C). Based on H. undatus genome data, homologous genes were obtained and renamed as HuADH1 (HU03G02979.1), HuCYP76AD1-1 (HU03G00480.1), HuDODA1 (HU03G01342.1), cDOPA5GT1 (HU07G00239.1) and cDOPA5GT2 (HU03G00240.1). HuADH1, HuCYP76AD1-1 and HuDODA1 showed increasing trends and predominantly expressed in pulps of the later period (23rd to 32nd DAF) of ‘Guanhuahong’ pitaya (Figure 6D). cDOPA5GT1 and cDOPA5GT2 were predominantly expressed on the 23rd DAF when pulps started to accumulate betacyanins. The expression of the HuADH1, HuCYP76AD1-1 and HuDODA1 were closely related to betalain accumulation in pulps during fruit ripening of ‘Guanhuahong’ pitaya.
In order to analyze the expression levels of R2R3-MYB genes, qRT-PCR were performed to screen candidate R2R3-MYB genes related to betalain biosynthesis during fruit ripening of ‘Guanhuahong’ pitaya. As shown in Figure 6D, three HuMYBs (HuMYB3, HuMYB25 and HuMYB30) showed upward trends while twenty-two HuMYBs (HuMYB1, HuMYB4, HuMYB9, HuMYB13, HuMYB16, HuMYB18, HuMYB21, HuMYB23, HuMYB31, HuMYB37, HuMYB45, HuMYB60, HuMYB80, HuMYB81, HuMYB82, HuMYB84, HuMYB85, HuMYB87, HuMYB89, HuMYB99, HuMYB100 and HuMYB103) showed significant downward trends in pulps during fruit ripening of ‘Guanhuahong’ pitaya (Figure 6D). The qRT-PCR results were consistent with RNA-Seq data (PRJNA704510) that most HuMYBs were down-regulated in pulps during fruit ripening (17th, 23rd, 25th and 32nd) of ‘Guanhuahong’ pitaya (Figure S2). Additionally, fourteen genes were predominantly expressed on the 15th DAF when the seeds were white (Figure S3). Five genes were highly expressed on the 17th DAF when the seeds became red (Figure S3). Fourteen genes were predominantly expressed on the 19th DAF when seeds were black (Figure S3). Eleven genes were highly expressed on the 23rd DAF when the betacyanin started to accumulate in pulps (Figure S3). Twenty-two genes showed irregular trends during pitaya fruit ripening (Figure S4). These results suggested that many HuMYBs may function as repressors during pitaya fruit ripening.

3.6. Phylogenetic and Sequence Analyses of Pitaya 2R-MYB Repressors

To further explore the R2R3-MYB repressors in pitaya, the phylogenetic relationship of seven HuMYBs, seventeen FaMYB1-like repressors and twenty-four AtMYB4-like repressors were constructed using the ML method with 1000 replications in MEGA7.0 (Figure 7A, Table S5). HuMYB1 was first clustered with AtrMYB4, SoMYB308-like, BvMYB308 and CqMYB308-like from betalain-producing plants, with ZmMYB31 and ZmMYB42, which regulated phenylpropanoid genes in maize [44]. HuMYB21 and HuMYB78 were first grouped with CqMYB6-like and BvMYB6, and then with NtMYB2 which repressed the regulation of anthocyanin biosynthesis [45]. HuMYB72 and HuMYB101 were clustered with PpMYB20 which had repressive effects on the flavonoid pathway [46]. FaMYB1-like repressor group was divided into anthocyanin and proanthocyanin repressor subgroups [47]. HuMYB48 and HuMYB49 were grouped in proanthocyanin repressors. Both gathered first with CqMYB3-like and BvMYB6-like and then with VvMYBC2-L3 which negatively regulated proanthocyanin biosynthesis [48]. These results suggested that seven HuMYBs, including five AtMYB4-like repressors (HuMYB1, HuMYB21, HuMYB72, HuMYB78 and HuMYB101) and two FaMYB1-like repressors (HuMYB48 and HuMYB49), were R2R3-MYB repressors.
Furthermore, the sequence alignment of the seven R2R3-MYB repressor proteins indicated that all of them had the conserved R2 and R3 domains in the N-terminal. Except for HuMYB21 and HuMYB78, the R3 domain of the other genes contained the bHLH binding motif ([D/E]Lx2[R/K]x3Lx6Lx3R) which was responsible for interaction with some bHLH TFs (Figure 7B, Figures S5 and S6). In addition to R2 and R3 domains, only HuMYB1 had four conserved motifs located in C-terminal in terms of C1 motif (KLIsrGIDPxT/SHRxI/L), C2 motif (pdLNLD/ELxiG/S or LxLxL), C3 motif (Cx1-2Cx7-12Cx2C) and C4 motif (FLGLx4-7V/LLD/GF/YR/Sx1LEMK) (Figure 7B). The C1 and C2 motifs were involved in bHLH interactions and considered as promoter repression domains [49,50]. AtMYB4-like repressors are characterized by C1 and C2 motifs in the C-terminal [51]. C4 motif was identified as negative regulation factor of floral volatile benzoid/phenylpropanoid compounds [52] and absent in FaMYB1-like repressors [53]. Besides this, in AtMYB4-like repressor group, HuMYB72, HuMYB78 and HuMYB101 only contained the C1 motif compared to HuMYB21 without conserved motif in C-terminal (Figure S5). In FaMYB1-like repressor group, HuMYB48 and HuMYB49 had the C1 and C2 motifs in their C-terminal (Figure S6). These results suggested that HuMYB1, HuMYB21, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101 were R2R3-MYB TFs due to conserving in R2 and R3 domains in N-terminal. Five genes (HuMYB1, HuMYB48, HuMYB49, HuMYB72 and HuMYB101) had a bHLH binding motif in R3 domain while six genes (HuMYB1, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101) had one to two conserved repression motifs in their C-terminal.

3.7. HuMYB1 Is a Nucleus-Localized Transcription Repressor

HuMYB1 probably plays an essential role in suppressing betalain biosynthesis in pitaya according to the above evolutionary, expression and sequence alignment analyses. To investigate the subcellular localizations of HuMYB1, the full-length coding sequence of HuMYB1 was fused with the GFP to construct 35S-HuMYB1-GFP vectors which were transiently transformed into N. benthamiana leaves. The nuclear signals were captured in protoplasts after overexpressing HuMYB1-GFP, whereas GFP alone showed cytoplasmatic and nuclear signals (Figure 8A).
Full-length coding regions of HuMYB1 were fused with the GAL4BD in the pGBKT7 vector to study the transcriptional activation abilities of HuMYB1 in yeast cells. As shown in Figure 8B, the transformed yeast cells of positive control (pGBKT753 + pGADT7-T) grew well in SD/-Trp-His-Ade and showed α-Gal activity, while yeast cells containing pGBKT7 (negative control) or pGBKT7-HuMYB1 did not, suggesting that HuMYB1 has no transactivation activities in yeast cells and probably functions as a transcriptional repressor in gene regulation.
The transcriptional repressor activities of HuMYB1 were confirmed in N. benthamiana leaves using the dual-luciferase reporter system (Figure 8C). Compared to the expression of negative control (BD-62SK), the expression of the positive control (BD-62SK-VP16) resulted in a higher value of LUC/REN ratio while the BD-62SK-HuMYB1 showed a significantly lower value (Figure 8D). It is consistent with the result from yeast cells. These results demonstrated that HuMYB1 was a nucleus-localized transcriptional repressor.

3.8. HuMYB1 Inhibited the Transcription of Three Betalain Biosynthesis-Related Genes

HuADH1, HuCYP76AD1-1 and HuDODA1 promoters were cloned from ‘Guanhuahong’ pitaya and the conserved cis-element motifs were predicted in Plant-CARE database. MYB binding sites were present in the promoters of HuADH1, HuCYP76AD1-1 and HuDODA1 (Table S6). Transient dual-luciferase assays in N. benthamiana leaves were performed to study whether HuMYB1 acted as the transcriptional repressor of HuADH1, HuCYP76AD1-1 and HuDODA1 (Figure 8E). As shown in Figure 8F, compared to the empty control (pGreenII 62-SK), LUC/REN ratio was significantly reduced when HuMYB1 was co-transformed with HuADH1 or HuCYP76AD1-1 or HuDODA1 pro-LUC reporter. These results indicated that HuMYB1 participated in betalain biosynthesis by repressing the transcription of betalain biosynthesis-related HuADH1, HuCYP76AD1-1 and HuDODA1 genes (Figure 8F).

4. Discussion

Although the genome-wide of MYB superfamily has been extensively studied in various plants, the identification and characterization of R2R3-MYB TFs based on the whole genome sequence of H. undatus were not reported yet. In our study, a total of 185 MYB genes (75 1R-MYB, 105 R2R3-MYB, four R1R2R3-MYB and one 4R-MYB) were identified by genome-wide search and distributed in all H. undatus Chrs (Figure 2). The R2R3-MYB proteins are the largest subfamily which is consistent with the findings in Arabidopsis [3], pineapple [13] and Salvia miltiorrhiza [36]. Compared with beet, the number of pitaya R2R3-MYB members (105) is fewer than Arabidopsis (126) [3] but more than beet (70) [16], indicating R2R3-MYB subfamilies in pitaya and Arabidopsis were expanded. Low tandem and high segmental duplications resulted in the expansion of the R2R3-MYB gene family in H. undatus, which also occurred in cotton [54] and pineapple [13]. Based on the comparison of DNA-binding domains (Figure 1), synteny (Figure 3) and phylogenetic relationship (Figure 4) between pitaya, beet and Arabidopsis, closer genetic distance was detected between betalain-producing plants of pitaya and beet compared to pitaya and the anthocyanin-producing plant Arabidopsis [27].
Among 105 R2R3-MYB genes, expression patterns of 91 HuMYBs were analyzed in pulps of the seven fruit ripening stages of ‘Guanhuahong’ pitaya. Some genes showed preferential expressions on the 15th, 17th, 19th or 23rd DAF, suggesting that they play multiple regulatory roles during pitaya fruit ripening (Figure S3). HuMYB3, HuMYB25 and HuMYB30 shared a similar expression pattern with betalain biosynthesis-related structure genes of HuADH1, HuCYP76AD1-1 and HuDODA1. However, 20.9% R2R3-MYB genes kept downward trends and were contrary to the expressions of HuADH1, HuCYP76AD1-1 and HuDODA1 during pitaya fruit ripening (Figure 6). Additionally, HuMYB93 and HuMYB101 shared similar expression patterns with HucDOPA5GT2 involved in the formation of betanin (Figure S3). ADH and ADT (arogenate dehydratase) can catalyze arogenate to form tyrosine and phenylalanine, respectively, and then form betalains and anthocyanins [55]. In the present study, HuMYB83 was clustered with BvMYB1 and AtMYB90. BvMYB1 is responsible for betalain biosynthesis [17], while AtMYB90 is involved in anthocyanin biosynthesis [3]. The transcription regulation of betalains and anthocyanins may exist a connection in their evolutionary relationship. Thus, HuMYB12, HuMYB26 and HuMYB83 may be involved in betalain biosynthesis since they were clustered with Arabidopsis subgroups 7, 5 and 6 involved in the regulation of flavonoid biosynthesis. Besides this, five AtMYB4-like repressors (HuMYB1, HuMYB21, HuMYB72, HuMYB78 and HuMYB101) and two FaMYB1-like repressors (HuMYB48 and HuMYB49) were identified according to the phylogenetic tree. In total, eight positive regulators and twenty-seven negative regulators were obtained by expression profiles and phylogenetic analyses (Figure 9).
According to expression and sequence analyses, HuMYB1 showed a downward trend during fruit ripening (Figure 7) and was closer to AtMYB4, which negatively regulated phenylpropanoid pathway [3]. Besides this, HuMYB1 acted as a nucleus-localized transcriptional repressor involved in betalain biosynthesis pathway by repressing the transcription of HuADH1, HuCYP76AD1-1 and HuDODA1 (Figure 8). BvMYB1 is involved in betalain biosynthesis by targeting the CYP76AD1 and DODA promoter while HpWRKY44 has a functional role in the regulation of CYP76AD1 [17,43]. Thus, HuMYB TFs are probably involved in betalain biosynthesis by binding to the HuADH1, HuCYP76AD1-1 and HuDODA1 promoters in pitaya. Therefore, a hypothetical model of betalain biosynthesis in pitaya was proposed according to these candidate HuMYBs possibly involved in betacyanin biosynthesis by regulating the expression levels of HuADH1, HuCYP76AD1-1, HuDODA1 and/or GTs (Figure 9). Besides this, ABRE, E-box, G-box and W-box cis-elements were detected in the sequences of HuADH1, HuCYP76AD1-1 and HuDODA1 promoters, which were the binding sites of TFs such as AREB/ABF, bZIP, bHLH and WRKY [43,56,57,58]. Therefore, these HuMYB TFs are possibly involved in betalain biosynthesis in pitaya. Moreover, the transcription activity of MYB proteins was influenced by the interaction with other TFs, such as bZIP, bHLH, NAC and WRKY TFs [59,60,61,62]. The activity of MYBs and bHLHs appears to be regulated by cytosolic WD40 repeat proteins by forming highly dynamic MBW complexes. WRKY can activate the transcription of anthocyanin-related genes by engaging with the MBW complex to form MBWW complex. Thus, further studies need to be carried out to determine whether the other TFs such as AREB/ABF, bZIP, bHLH and WRKY are involved in betalain biosynthesis, or HuMYB1 can coordinate with the other TFs to regulate betalain biosynthesis.

5. Conclusions

In conclusion, this study is the first report on identification and characterization of R2R3-MYB gene family based on the genome of H. undatus. A total of 105 R2R3-MYB genes were obtained from the H. undatus genome. These R2R3-MYB genes were distributed in all 11 chromosomes with conserved R2 and R3 repeats and expanded its subfamily members through segmental duplication with a closer distance of beet than Arabidopsis. They were functionally classified into nine groups, consisting of conserved motif and exon–intron clustering. Phylogeny, sequence alignment and expression patterns of seven candidate R2R3-MYB repressors were further analyzed. HuMYB1 was grouped in the AtMYB1-like repressor group with bHLH binding motif and C1-C4 motifs. Besides this, HuMYB1 showed a downward expression trend opposite to betalain accumulation in pulps during fruit ripening of ‘Guanhuahong’ pitaya. HuMYB1 was a nucleus protein and could reduce the transcription activity of HuADH1, HuCYP76AD1-1 and HuDODA1 with transcription repression activities. The present study provides valuable information for a better understanding of MYB TFs involved in fruit ripening and betalain biosynthesis in pitaya.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/cells10081949/s1, Figure S1: Consensus sequences of the group specific motifs. Figure S2: The heatmap of R2R3-MYB genes in pulps of ‘Guanhuahong’ pitaya. Figure S3: The expression analyses of R2R3-MYB genes in pulps during fruit ripening of ‘Guanhuahong’ pitaya. Figure S4: The expression analyses of R2R3-MYB genes in pulps during fruit ripening of ‘Guanhuahong’ pitaya. Figure S5: The sequence analyses of HuMYB21, HuMYB72, HuMYB78 and HuMYB101. Figure S6: The sequence analyses of HuMYB48 and HuMYB49. Table S1: Summary of primers used in this study. Table S2: Summary of MYB members information in pitaya genome. Table S3: Segmentally and tandemly duplicated HuMYB gene pairs. Table S4: One-to-one orthologous relationships between pitaya, beet, and Arabidopsis. Table S5: Proteins used in Figure 7. Table S6: The regulatory motifs found within the HuADH1, HuCYP76AD1-1, and HuDODA1 promoter.

Author Contributions

Conceptualization, F.X., J.C. and Y.Q.; methodology and validation, F.X., Q.H. and C.C.; formal analysis, investigation, resources, data curation and visualization, F.X., Q.H., C.C., Z.Z., R.Z., J.Z. and G.H.; writing—original draft preparation, F.X.; writing—review and editing, F.X. and Y.Q.; supervision, J.C. and Y.Q.; project administration, G.H., J.C. and Y.Q.; funding acquisition, J.C. and Y.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by grants from the National Natural Science Foundation of China (grant nos. 31972367 and 31960578), Key Science and Technology Planning Project of Guangzhou (grant no. 201904020015), Science and Technology Program of Zhanjiang (grant no. 2019A01003), Science and Technology Program of Guangzhou (grant nos. 202002020060, 201704020003 and 2014Y2-00164) and Key Realm R&D Program of Guangdong Province (grant no. 2018B020202011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in the NCBI PRJNA704510.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Liu, L.S.; White, M.J.; Macrae, T.H. Transcription factors and their genes in higher plants: Functional domains, evolution and regulation. Eur. J. Biochem. 1999, 262, 247–257. [Google Scholar] [CrossRef]
  2. Jia, L.; Clegg, M.T.; Jiang, T. Evolutionary dynamics of the DNA-binding domains in putative R2R3-MYB genes identified from rice subspecies indica and japonica genomes. Plant Physiol. 2004, 134, 575–585. [Google Scholar] [CrossRef] [Green Version]
  3. Dubos, C.; Stracke, R.; Grotewold, E.; Weisshaar, B.; Martin, C.; Lepiniec, L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef]
  4. Rosinski, J.A.; Atchley, W.R. Molecular evolution of the Myb family of transcription factors: Evidence for polyphyletic origin. J. Mol. Evol. 1998, 46, 74–83. [Google Scholar] [CrossRef]
  5. Jiang, C.; Gu, J.; Chopra, S.; Gu, X.; Peterson, T. Ordered origin of the typical two- and three-repeat Myb genes. Gene 2004, 326, 13–22. [Google Scholar] [CrossRef] [PubMed]
  6. Du, H.; Liang, Z.; Zhao, S.; Nan, M.G.; Tran, L.S.; Lu, K.; Huang, Y.B.; Li, J.N. The evolutionary history of R2R3-MYB proteins across 50 eukaryotes: New insights into subfamily classification and expansion. Sci. Rep. 2015, 5, 11037. [Google Scholar] [CrossRef] [Green Version]
  7. Jiang, C.K.; Rao, G.Y. Insights into the diversification and evolution of R2R3-MYB transcription factors in plants. Plant Physiol. 2020, 183, 637–655. [Google Scholar] [CrossRef] [PubMed]
  8. Klempnauer, K.H.; Gonda, T.J.; Bishop, J.M. Nucleotide sequence of the retroviral leukemia gene v-myb and its cellular progenitor c-myb: The architecture of a transduced oncogene. Cell 1982, 31, 453–463. [Google Scholar] [CrossRef]
  9. Paz-Ares, J.; Ghosal, D.; Wienand, U.; Peterson, P.A.; Saedler, H. The regulatory c1 locus of Zea mays encodes a protein with homology to myb proto-oncogene products and with structural similarities to transcriptional activators. EMBO J. 1987, 6, 3553–3558. [Google Scholar] [CrossRef]
  10. Matus, J.T.; Aquea, F.; Arce-Johnson, P. Analysis of the grape MYB R2R3 subfamily reveals expanded wine quality-related clades and conserved gene structure organization across Vitis and Arabidopsis genomes. BMC Plant Biol. 2008, 8, 83. [Google Scholar] [CrossRef] [Green Version]
  11. Cao, Z.H.; Zhang, S.Z.; Wang, R.K.; Zhang, R.F.; Hao, Y.J. Genome wide analysis of the apple MYB transcription factor family allows the identification of MdoMYB121 gene confering abiotic stress tolerance in plants. PLoS ONE 2013, 8, e69955. [Google Scholar] [CrossRef] [PubMed]
  12. Li, W.; Ding, Z.; Ruan, M.; Yu, X.; Peng, M.; Liu, Y. Kiwifruit R2R3-MYB transcription factors and contribution of the novel AcMYB75 to red kiwifruit anthocyanin biosynthesis. Sci. Rep. 2017, 7, 16861. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, C.; Xie, T.; Chen, C.; Luan, A.; Long, J.; Li, C.; Ding, Y.; He, Y. Genome-wide organization and expression profiling of the R2R3-MYB transcription factor family in pineapple (Ananas comosus). BMC Genom. 2017, 18, 503. [Google Scholar] [CrossRef] [Green Version]
  14. Xu, Q.; He, J.; Dong, J.; Hou, X.; Zhang, X. Genomic survey and expression profiling of the MYB gene family in watermelon. Hortic. Plant J. 2018, 4, 1–15. [Google Scholar] [CrossRef]
  15. Pucker, B.; Pandey, A.; Weisshaar, B.; Stracke, R. The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns. PLoS ONE 2020, 15, e0239275. [Google Scholar] [CrossRef]
  16. Stracke, R.; Holtgräwe, D.; Schneider, J.; Pucker, B.; Sörensen, T.R.; Weisshaar, B. Genome-wide identification and characterisation of R2R3-MYB genes in sugar beet (Beta vulgaris). BMC Plant Biol. 2014, 14, 249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Hatlestad, G.J.; A Akhavan, N.; Sunnadeniya, R.M.; Elam, L.; Cargile, S.; Hembd, A.; González, A.; McGrath, J.M.; Lloyd, A.M. The beet Y locus encodes an anthocyanin MYB-like protein that activates the betalain red pigment pathway. Nat. Genet. 2014, 47, 92–96. [Google Scholar] [CrossRef] [PubMed]
  18. Ibrahim, S.R.M.; Mohamed, G.A.; Khedr, A.; Zayed, M.; El-Kholy, A.A.E.S. Genus Hylocereus: Beneficial phytochemicals, nutritional importance, and biological relevance—A review. J. Food Biochem. 2018, 42, e12491. [Google Scholar] [CrossRef]
  19. Hatlestad, G.J.; Sunnadeniya, R.M.; A Akhavan, N.; González, A.; Goldman, I.L.; McGrath, J.M.; Lloyd, A.M. The beet R locus encodes a new cytochrome P450 required for red betalain production. Nat. Genet. 2012, 44, 816–820. [Google Scholar] [CrossRef] [PubMed]
  20. Lopez-Nieves, S.; Yang, Y.; Timoneda, A.; Wang, M.; Feng, T.; Smith, S.A.; Brockington, S.F.; Maeda, H.A. Relaxation of tyrosine pathway regulation underlies the evolution of betalain pigmentation in Caryophyllales. New Phytol. 2018, 217, 896–908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Brockington, S.F.; Yang, Y.; Gandia-Herrero, F.; Covshoff, S.; Hibberd, J.M.; Sage, R.F.; Wong, G.K.S.; Moore, M.; Smith, S. Lineage-specific gene radiations underlie the evolution of novel betalain pigmentation in Caryophyllales. New Phytol. 2015, 207, 1170–1180. [Google Scholar] [CrossRef] [PubMed]
  22. 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, 287–296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Timoneda, A.; Feng, T.; Sheehan, H.; Walker-Hale, N.; Pucker, B.; Lopez-Nieves, S.; Guo, R.; Brockington, S. The evolution of betalain biosynthesis in Caryophyllales. New Phytol. 2019, 224, 71–85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gengatharan, A.; Dykes, G.A.; Choo, W.S. Betalains: Natural plant pigments with potential application in functional foods. LWT-Food Sci.Technol. 2015, 64, 645–649. [Google Scholar] [CrossRef]
  25. Polturak, G.; Aharoni, A. Advances and future directions in betalain metabolic engineering. New Phytol. 2019, 224, 1472–1478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bateman, A.; Coin, L.; Durbin, R.; Finn, R.D.; Hollich, V.; Griffits-Jones, S.; Khanna, A.; Marshall, M.; Moxon, S.; Sonnhammer, E.L.L.; et al. The Pfam protein families database. Nucleic Acids Res. 2004, 32, D138–D141. [Google Scholar] [CrossRef]
  27. Chen, J.; Xie, F.; Cui, Y.; Chen, C.; Lu, W.; Hu, X.; Hua, Q.; Zhao, J.; Wu, Z.; Gao, D.; 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–15. [Google Scholar] [CrossRef]
  28. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  29. Wang, Y.; Tang, H.; DeBarry, J.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  30. Hua, Q.; Chen, C.; Tel-Zur, N.; Wang, H.; Wu, J.; Chen, J.; Zhang, Z.; Zhao, J.; Hu, G.; Qin, Y. Metabolomic characterization of pitaya fruit from three red-skinned cultivars with different pulp colors. Plant Physiol. Biochem. 2018, 126, 117–125. [Google Scholar] [CrossRef]
  31. Chen, C.; Wu, J.; Hua, Q.; Tel-Zur, N.; Xie, F.; Zhang, Z.; Chen, J.; Zhang, R.; Hu, G.; Zhao, J.; et al. Identification of reliable reference genes for quantitative real-time PCR normalization in pitaya. Plant Methods 2019, 15, 1–12. [Google Scholar] [CrossRef] [Green Version]
  32. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  33. Lai, B.; Du, L.; Liu, R.; Hu, B.; Su, W.; Qin, Y.; Zhao, J.; Wang, H.; Hu, G. 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. [Google Scholar] [CrossRef] [Green Version]
  34. Yang, Y.; Shan, W.; Kuang, J.; Chen, J.; Lu, W. Four HD-ZIPs are involved in banana fruit ripening by activating the transcription of ethylene biosynthetic and cell wall-modifying genes. Plant Cell Rep. 2019, 39, 351–362. [Google Scholar] [CrossRef]
  35. Rotman, N.; Durbarry, A.; Wardle, A.; Yang, W.C.; Chaboud, A.; Faure, J.E.; Berger, F.; Twell, D. A novel class of MYB factors controls sperm-cell formation in plants. Curr. Biol. 2005, 15, 244–248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Li, C.; Lu, S. Genome-wide characterization and comparative analysis of R2R3-MYB transcription factors shows the complexity of MYB-associated regulatory networks in Salvia miltiorrhiza. BMC Genom. 2014, 15, 277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wang, Z.; Tang, J.; Hu, R.; Wu, P.; Hou, X.; Song, X.; Xiong, A. Genome-wide analysis of the R2R3-MYB transcription factor genes in Chinese cabbage (Brassica rapa ssp. pekinensis) reveals their stress and hormone responsive patterns. BMC Genom. 2015, 16, 17–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Wilkins, O.; Nahal, H.; Foong, J.; Provart, N.J.; Campbell, M.M. Expansion and diversification of the Populus R2R3-MYB family of transcription factors. Plant Physiol. 2009, 149, 981–993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Soler, M.; Camargo, E.L.O.; Carocha, V.; Cassan-Wang, H.; San Clemente, H.; Savelli, B.; Hefer, C.A.; Paiva, J.A.P.; Myburg, A.A.; Grima-Pettenati, J. The eucalyptus grandis R2R3-MYB transcription factor family: Evidence for woody growth-related evolution and function. New Phytol. 2015, 206, 1364–1377. [Google Scholar] [CrossRef]
  40. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004, 4, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Yu, J.; Wang, J.; Lin, W.; Li, S.G.; Li, H.; Zhou, J.; Ni, P.; Dong, W.; Hu, S.; Zeng, C.; et al. The genomes of Oryza sativa: A history of duplications. PLoS Biol. 2005, 3, 266–281. [Google Scholar] [CrossRef] [Green Version]
  42. Romano, J.M.; Dubos, C.; Prouse, M.B.; Wilkins, O.; Hong, H.; Poole, M.; Kang, K.Y.; Li, E.; Douglas, C.J.; Western, T.L.; et al. AtMYB61, an R2R3-MYB transcription factor, functions as a pleiotropic regulator via a small gene network. New Phytol. 2012, 195, 774–786. [Google Scholar] [CrossRef]
  43. Cheng, M.; Huang, Z.; Hua, Q.; Shan, W.; Kuang, J.; Lu, W.; Qin, Y.; Chen, J. The WRKY transcription factor HpWRKY44 regulates CytP450-like1 expression in red pitaya fruit (Hylocereus polyrhizus). Hortic. Res. 2017, 4, 17039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Agarwal, T.; Grotewold, E.; Doseff, A.I.; Gray, J. MYB31/MYB42 syntelogs exhibit divergent regulation of phenylpropanoid genes in maize, sorghum and rice. Sci. Rep. 2016, 6, 28502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Anwar, M.; Wang, G.; Wu, J.; Waheed, S.; Allan, A.C.; Zeng, L. Ectopic overexpression of a novel R2R3-MYB, NtMYB2 from Chinese narcissus represses anthocyanin biosynthesis in tobacco. Molecules 2018, 23, 781. [Google Scholar] [CrossRef] [Green Version]
  46. Zhou, H.; Peng, Q.; Zhao, J.; Owiti, A.; Ren, F.; Liao, L.; Wang, L.; Deng, X.; Jiang, Q.; Han, Y. Multiple R2R3-MYB transcription factors involved in the regulation of anthocyanin accumulation in peach flower. Front. Plant Sci. 2016, 7, 1557. [Google Scholar] [CrossRef]
  47. Ma, D.; Reichelt, M.; Yoshida, K.; Gershenzon, J.; Constabel, C.P. Two R2R3-MYB proteins are broad repressors of flavonoid and phenylpropanoid metabolism in poplar. Plant J. 2018, 96, 949–965. [Google Scholar] [CrossRef] [Green Version]
  48. Cavallini, E.; Matus, J.T.; Finezzo, L.; Zenoni, S.; Loyola, R.; Guzzo, F.; Schlechter, R.; Ageorges, A.; Arce-Johnson, P.; Tornielli, G.B. The phenylpropanoid pathway is controlled at different branches by a set of R2R3-MYB C2 repressors in grapevine. Plant Physiol. 2015, 167, 1448–1470. [Google Scholar] [CrossRef]
  49. Jin, H.; Cominelli, E.; Bailey, P.; Parr, A.; Mehrtens, F.; Jones, J.; Tonelli, C.; Weisshaar, B.; Martin, C. Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis. EMBO J. 2000, 19, 6150–6161. [Google Scholar] [CrossRef] [Green Version]
  50. Kagale, S.; Links, M.; Rozwadowski, K. Genome-wide analysis of ethylene-responsive element binding factor-associated amphiphilic repression motif-containing transcriptional regulators in Arabidopsis. Plant Physiol. 2010, 152, 1109–1134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Kranz, H.D.; Denekamp, M.; Greco, R.; Jin, H.; Leyva, A.; Meissner, R.C.; Petroni, K.; Urzainqui, A.; Bevan, M.; Martin, C.; et al. Towards functional characterisation of the members of the R2R3-MYB gene family from Arabidopsis thaliana. Plant J. 1998, 16, 263–276. [Google Scholar] [CrossRef]
  52. Colquhoun, T.A.; Verdonk, J.; Schimmel, B.; Tieman, D.M.; Underwood, B.A.; Clark, D.G. Petunia floral volatile benzenoid/phenylpropanoid genes are regulated in a similar manner. Phytochemistry 2010, 71, 158–167. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, L.; Hu, B.; Qin, Y.; Hu, G.; Zhao, J. Advance of the negative regulation of anthocyanin biosynthesis by MYB transcription factors. Plant Physiol. Biochem. 2019, 136, 178–187. [Google Scholar] [CrossRef]
  54. He, Q.; Jones, D.C.; Li, W.; Xie, F.; Ma, J.; Sun, R.; Wang, Q.; Zhu, S.; Zhang, B. Genome-wide identification of R2R3-MYB genes and expression analyses during abiotic stress in Gossypium raimondii. Sci. Rep. 2016, 6, 22980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Tanaka, Y.; Sasaki, N.; Ohmiya, A. Biosynthesis of plant pigments: Anthocyanins, betalains and carotenoids. Plant J. 2008, 54, 733–749. [Google Scholar] [CrossRef] [PubMed]
  56. Fujita, Y.; Fujita, M.; Satoh, R.; Maruyama, K.; Parvez, M.M.; Seki, M.; Hiratsu, K.; Ohme-Takagi, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. AREB1 is a transcription activator of novel ABRE-dependent ABA signaling that enhances drought stress tolerance in Arabidopsis. Plant Cell 2005, 17, 3470–3488. [Google Scholar] [CrossRef] [Green Version]
  57. Ezer, D.; Shepherd, S.J.; Brestovitsky, A.; Dickinson, P.; Cortijo, S.; Charoensawan, V.; Box, M.S.; Biswas, S.; Jaeger, K.E.; Wigge, P.A. The G-Box transcriptional regulatory code in Arabidopsis. Plant Physiol. 2017, 175, 628–640. [Google Scholar] [CrossRef] [Green Version]
  58. Chakraborty, M.; Gangappa, S.N.; Maurya, J.P.; Sethi, V.; Srivastava, A.K.; Singh, A.; Dutta, S.; Ojha, M.; Gupta, N.; Sengupta, M.; et al. Functional interrelation of MYC2 and HY5 plays an important role in Arabidopsis seedling development. Plant J. 2019, 99, 1080–1097. [Google Scholar] [CrossRef]
  59. An, J.; Yao, J.; Xu, R.; You, C.; Wang, X.; Hao, Y. Apple bZIP transcription factor MdbZIP44 regulates abscisic acid-promoted anthocyanin accumulation. Plant Cell Environ. 2018, 41, 2678–2692. [Google Scholar] [CrossRef]
  60. Sun, Q.; Jiang, S.; Zhang, T.; Xu, H.; Fang, H.; Zhang, J.; Su, M.; Wang, Y.; Zhang, Z.; Wang, N.; et al. Apple NAC transcription factor MdNAC52 regulates biosynthesis of anthocyanin and proanthocyanidin through MdMYB9 and MdMYB11. Plant Sci. 2019, 289, 110286. [Google Scholar] [CrossRef]
  61. Wang, L.; Tang, W.; Hu, Y.; Zhang, Y.; Sun, J.; Guo, X.; Lu, H.; Yang, Y.; Fang, C.; Niu, X.; et al. A MYB/bHLH complex regulates tissue-specific anthocyanin biosynthesis in the inner pericarp of red-centered kiwifruit Actinidia chinensis cv. Hongyang. Plant J. 2019, 99, 359–378. [Google Scholar] [CrossRef] [PubMed]
  62. Peng, Y.; Thrimawithana, A.H.; Cooney, J.M.; Jensen, D.J.; Espley, R.V.; Allan, A.C. The proanthocyanin-related transcription factors MYBC1 and WRKY44 regulate branch points in the kiwifruit anthocyanin pathway. Sci. Rep. 2020, 10, 14161. [Google Scholar] [CrossRef] [PubMed]
Cells 10 01949 g001 550
Figure 1. Comparison of DNA-binding domains of R2R3-MYB transcription factor in H. undatus, B. vulgaris and A. thaliana. Sequence logos of the R2 (A,C,E) and R3 (B,D,F) repeats are based on conserved alignments from H. undatus (A,B), B. vulgaris (C,D) and A. thaliana (E,F). The overall height of each stack indicates the conservation of the sequence at the position, whereas the height of letters within each stack represents the relative frequency of the corresponding amino acid. Highly conserved tryptophan (W) and phenylalanine (F) residues are indicated by asterisks. The positions with different patterns between H. undatus, B. vulgaris and A. thaliana are indicated by arrows.
Figure 1. Comparison of DNA-binding domains of R2R3-MYB transcription factor in H. undatus, B. vulgaris and A. thaliana. Sequence logos of the R2 (A,C,E) and R3 (B,D,F) repeats are based on conserved alignments from H. undatus (A,B), B. vulgaris (C,D) and A. thaliana (E,F). The overall height of each stack indicates the conservation of the sequence at the position, whereas the height of letters within each stack represents the relative frequency of the corresponding amino acid. Highly conserved tryptophan (W) and phenylalanine (F) residues are indicated by asterisks. The positions with different patterns between H. undatus, B. vulgaris and A. thaliana are indicated by arrows.
Cells 10 01949 g001
Cells 10 01949 g002 550
Figure 2. Distribution of R2R3-MYB members on the eleven chromosomes of pitaya. R2R3-HuMYB genes are shown on the right of each chromosome. Gene positions and the size of each chromosome can be estimated using the scale on the right of the figure, the scale indicates 20 megabases (Mb).
Figure 2. Distribution of R2R3-MYB members on the eleven chromosomes of pitaya. R2R3-HuMYB genes are shown on the right of each chromosome. Gene positions and the size of each chromosome can be estimated using the scale on the right of the figure, the scale indicates 20 megabases (Mb).
Cells 10 01949 g002
Cells 10 01949 g003 550
Figure 3. Synteny analyses of the pitaya, beet and Arabidopsis R2R3-MYB genes. (A) Schematic representation of interchromosomal relationships of the pitaya R2R3-MYB genes. (B) Gene duplication and synteny relationship of R2R3-MYB genes between pitaya and beet (a), pitaya and Arabidopsis (b), and beet and Arabidopsis (c). Gray lines in the background indicate the collinear blocks between genomes, while orange lines highlight the syntenic MYB gene pairs.
Figure 3. Synteny analyses of the pitaya, beet and Arabidopsis R2R3-MYB genes. (A) Schematic representation of interchromosomal relationships of the pitaya R2R3-MYB genes. (B) Gene duplication and synteny relationship of R2R3-MYB genes between pitaya and beet (a), pitaya and Arabidopsis (b), and beet and Arabidopsis (c). Gray lines in the background indicate the collinear blocks between genomes, while orange lines highlight the syntenic MYB gene pairs.
Cells 10 01949 g003
Cells 10 01949 g004 550
Figure 4. Phylogenetic Maximum likelihood (ML) tree (1000 bootstraps) with R2R3-MYB proteins from H. undatus (Hu), B. vulgaris (Bv) and A. thaliana (At) by MEGA7.0. Green circles indicate the bootstrap value range from 81 to 100 in the tree, orange is from 60 to 80, gray is from 0 to 59. Groups and Subgroups (according to the A. thaliana) are labeled with different alternating tones of gray and black background, respectively.
Figure 4. Phylogenetic Maximum likelihood (ML) tree (1000 bootstraps) with R2R3-MYB proteins from H. undatus (Hu), B. vulgaris (Bv) and A. thaliana (At) by MEGA7.0. Green circles indicate the bootstrap value range from 81 to 100 in the tree, orange is from 60 to 80, gray is from 0 to 59. Groups and Subgroups (according to the A. thaliana) are labeled with different alternating tones of gray and black background, respectively.
Cells 10 01949 g004
Cells 10 01949 g005 550
Figure 5. Phylogenetic relationship, gene structure and architecture of conserved protein motifs in R2R3-MYB genes from pitaya. (A) The neighbor-joining (NJ) tree on the left includes 105 R2R3-MYB proteins from pitaya with 1000 bootstrap. The MYB proteins were clustered into 9 groups. (B) Architecture of conserved protein motifs in 9 groups. Each motif is represented by a number on the colored box and their sequence were shown in Figure S1. (C) Exon–intron structures of R2R3-MYB proteins from pitaya. Exon(s) and intron(s) are represented by purple boxes and black lines, respectively. Untranslated region(s) are indicated by blue boxes.
Figure 5. Phylogenetic relationship, gene structure and architecture of conserved protein motifs in R2R3-MYB genes from pitaya. (A) The neighbor-joining (NJ) tree on the left includes 105 R2R3-MYB proteins from pitaya with 1000 bootstrap. The MYB proteins were clustered into 9 groups. (B) Architecture of conserved protein motifs in 9 groups. Each motif is represented by a number on the colored box and their sequence were shown in Figure S1. (C) Exon–intron structures of R2R3-MYB proteins from pitaya. Exon(s) and intron(s) are represented by purple boxes and black lines, respectively. Untranslated region(s) are indicated by blue boxes.
Cells 10 01949 g005
Cells 10 01949 g006 550
Figure 6. Expression analyses of R2R3-MYB genes in pulps during fruit ripening of ‘Guanhuahong’ pitaya. (A) Pulp coloration of fruit ripening of ‘Guanhuahong’ pitaya. Bars = 2 cm. (B) Contents of betacyanin and betaxanthin in pulps during fruit ripening of ‘Guanhuahong’ pitaya. DAF, day after flowering. (C) A brief pathway of betalain biosynthesis in pitaya. (D) Expression analyses of four structural genes of betalain biosynthesis pathway (labeled in pink), and candidate up-regulated (labeled in green) and down-regulated (labeled in blue) R2R3-MYB genes in pulps during fruit ripening of ‘Guanhuahong’ pitaya.
Figure 6. Expression analyses of R2R3-MYB genes in pulps during fruit ripening of ‘Guanhuahong’ pitaya. (A) Pulp coloration of fruit ripening of ‘Guanhuahong’ pitaya. Bars = 2 cm. (B) Contents of betacyanin and betaxanthin in pulps during fruit ripening of ‘Guanhuahong’ pitaya. DAF, day after flowering. (C) A brief pathway of betalain biosynthesis in pitaya. (D) Expression analyses of four structural genes of betalain biosynthesis pathway (labeled in pink), and candidate up-regulated (labeled in green) and down-regulated (labeled in blue) R2R3-MYB genes in pulps during fruit ripening of ‘Guanhuahong’ pitaya.
Cells 10 01949 g006
Cells 10 01949 g007 550
Figure 7. The sequence analyses of R2R3-MYB repressors. (A) The phylogenetic analyses of 2R-MYB repressor TFs based on the ML method with 1000 replications by MEGA7.0. Red background indicates the AtMYB4-like repressors while green and yellow background are anthocyanin and proanthocyanin repressors of FaMYB1-like repressors group, respectively. HuMYBs are labeled by the red star. (B) The alignment analyses of HuMYB1. R2 and R3 domains are labeled in green and orange shading, respectively. C1, C2, C3 and C4 motifs are labeled in purple shading. Blue asterisks indicate the motif ([D/E]Lx2[R/K]x3Lx6Lx3R) which can interact with a bHLH partner.
Figure 7. The sequence analyses of R2R3-MYB repressors. (A) The phylogenetic analyses of 2R-MYB repressor TFs based on the ML method with 1000 replications by MEGA7.0. Red background indicates the AtMYB4-like repressors while green and yellow background are anthocyanin and proanthocyanin repressors of FaMYB1-like repressors group, respectively. HuMYBs are labeled by the red star. (B) The alignment analyses of HuMYB1. R2 and R3 domains are labeled in green and orange shading, respectively. C1, C2, C3 and C4 motifs are labeled in purple shading. Blue asterisks indicate the motif ([D/E]Lx2[R/K]x3Lx6Lx3R) which can interact with a bHLH partner.
Cells 10 01949 g007
Cells 10 01949 g008 550
Figure 8. Subcellular localization and transcriptional activation analyses of HuMYB1. (A) Subcellular localization of HuMYB1 in protoplasts of Nicotiana benthamiana. Bars = 20 μm. (B) Transcriptional activation of HuMYB1 in yeast cells. pGBKT7 and pGBKT7-53 + pGADT7-T were used as negative and positive control, respectively. (C) Diagrams of the reporter and effector vectors. (D) Transcriptional activation of HuMYB1 in N. benthamiana leaves. The LUC/REN ratio of the empty BD-62SK vector was used as a calibrator (set as 1). BD-62SK-VP16 was used as a positive control. Asterisks represents highly significant differences at p value < 0.01 using two-tailed t-test, compared to BD-62SK. (E) Diagrams of the reporter and effector vectors. (F) HuMYB1 inhibited the transcription of betalain-related genes including HuADH1, HuCYP76AD1-1 and HuDODA1 by dual-luciferase transient expression assay in N. benthamiana leaves. The ratio of LUC/REN of the empty vector (62SK) plus promoter was used as a calibrator (set as 1). Asterisks represents highly significant differences at p value < 0.05 (one asterisks) and p value < 0.01 (two asterisks) using two-tailed t-test, compared to 62SK.
Figure 8. Subcellular localization and transcriptional activation analyses of HuMYB1. (A) Subcellular localization of HuMYB1 in protoplasts of Nicotiana benthamiana. Bars = 20 μm. (B) Transcriptional activation of HuMYB1 in yeast cells. pGBKT7 and pGBKT7-53 + pGADT7-T were used as negative and positive control, respectively. (C) Diagrams of the reporter and effector vectors. (D) Transcriptional activation of HuMYB1 in N. benthamiana leaves. The LUC/REN ratio of the empty BD-62SK vector was used as a calibrator (set as 1). BD-62SK-VP16 was used as a positive control. Asterisks represents highly significant differences at p value < 0.01 using two-tailed t-test, compared to BD-62SK. (E) Diagrams of the reporter and effector vectors. (F) HuMYB1 inhibited the transcription of betalain-related genes including HuADH1, HuCYP76AD1-1 and HuDODA1 by dual-luciferase transient expression assay in N. benthamiana leaves. The ratio of LUC/REN of the empty vector (62SK) plus promoter was used as a calibrator (set as 1). Asterisks represents highly significant differences at p value < 0.05 (one asterisks) and p value < 0.01 (two asterisks) using two-tailed t-test, compared to 62SK.
Cells 10 01949 g008
Cells 10 01949 g009 550
Figure 9. A hypothetical model of candidate R2R3-MYB genes involved in betalain biosynthesis of pitaya. Spon. indicates spontaneous.
Figure 9. A hypothetical model of candidate R2R3-MYB genes involved in betalain biosynthesis of pitaya. Spon. indicates spontaneous.
Cells 10 01949 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xie, F.; Hua, Q.; Chen, C.; Zhang, Z.; Zhang, R.; Zhao, J.; Hu, G.; Chen, J.; Qin, Y. 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, 1949. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10081949

AMA Style

Xie F, Hua Q, Chen C, Zhang Z, Zhang R, Zhao J, Hu G, Chen J, Qin Y. 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://0-doi-org.brum.beds.ac.uk/10.3390/cells10081949

Chicago/Turabian Style

Xie, Fangfang, Qingzhu Hua, Canbin Chen, Zhike Zhang, Rong Zhang, Jietang Zhao, Guibing Hu, Jianye Chen, and Yonghua Qin. 2021. "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 10, no. 8: 1949. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10081949

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop