Over-expression of CcMYB24, encoding a R2R3-MYB transcription factor from a high-leaf-number mutant of Cymbidium, increases the number of leaves in Arabidopsis

Ornamental foliage plants have long been cultivated for their attractive leaves. Variation in leaf traits of ornamental foliage plants is one of the goals in breeding. MYB transcription factors regulate many aspects of leaf development, and thus influence morphological traits of leaves. However, little is known about the function of MYB transcription factors in leaf development of Cymbidium, one of the most economically important ornamental plants in the world. In the present study, a MYB transcription factor, CcMYB24, was identified and the corresponding gene cloned from a new orchid mutant, TRIR-2, which produces more leaves than control plants. The CcMYB24 showed a higher expression level in ‘TRIR-2’ than in control plants, and the protein was located in the nucleus. The sequence of CcMYB24 showed a high similarity with RAX2-like genes which belong to the R2R3-MYB gene family in other Cymbidium plants. Overexpression of CcMYB24 resulted in a phenotype with an increased number of leaves, elevated chlorophyll content, and decreased contents of carotenoids and flavonoids in Arabidopsis. These results provide functional evidence for the role of CcMYB24 in promoting the production of leaves in ‘TRIR-2’. Understanding the role of CcMYB24 in Cymbidium will be beneficial for the molecular breeding of ornamental foliage plants.

INTRODUCTION transcription factors may influence leaf traits in multiple pathways. With the development of gene function verification technology, it can be expected that more functions of MYB transcription factors in relation to the regulation of leaf traits will be revealed in different plant taxa in the future.
Cymbidium is one of the most economically important genera of the Orchidaceae and is highly popular around the world for its beautiful flowers and leaf variegation . A mutagenesis-based approach is the main way to produce various flower and leaf color mutants in different Cymbidium species . Leaf traits, such as color, variegation patterning, shape and number of leaves, are important in determining the quality of ornamental foliage plants. In China, the variation in leaf traits in ornamental foliage plants has been called 'leaf art', and novel variation in leaf traits is one of the goals in ornamental plant breeding. In our previous study, we obtained a new orchid mutant, temporarily named 'TRIR-2', through radiation mutation breeding on a hybrid between two Cymbidium species (Su et al., 2016). A comparative cytological study between the non-mutant Cymbidium inter-specific hybrid 'TRIR-1' and 'TRIR-2', which was not induced by irradiation, revealed significant differences. The study found that 'TRIR-2' had fewer chloroplasts, mitochondria, and other organelles, and some leaf cells of 'TRIR-2' contained only minimal cytoplasm. Consequently, the leaves of 'TRIR-2' appeared albino, with green color only visible at the edges and tips (Su et al., 2016). Furthermore, it was obvious that 'TRIR-2' generated more leaves than 'TRIR-1'. According to our unpublished data, the result of RNA sequencing had revealed that expression of a MYB gene CcMYB24 was significantly upregulated in 'TRIR-2', compared with 'TRIR-1' as the control, indicating that this gene may contribute to the morphological differences in 'TRIR-2'.
In the present study, we describe the cloning and characterization of CcMYB24 from the ornamental foliage mutant line 'TRIR-2'. The sequence analysis and protein structure prediction showed that CcMYB24 belongs to the R2R3-MYB subfamily. Overexpression of CcMYB24 in transgenic Arabidopsis seedlings significantly increased the number of leaves. These results indicate that this gene plays a positive role in increasing the number of leaves in 'TRIR-2', providing a basis for exploring the molecular mechanism of CcMYB24 in leaf development, and identifying a potential gene target for the molecular breeding of new cultivars in ornamental foliage plants.

Plant materials and growth conditions
The mutant plant material 'TRIR-2' used in this study was produced by 60 Co γ-irradiation on 'TRIR-1', the latter being the hybrid offspring of two Cymbidium species, Cymbidium tracyanum × Cymbidium iridioides. In brief, the two Cymbidium species, C. tracyanum and C. iridioides, were first cross-fertilized to generate 'TRIR-1', an orchid hybrid with a stable phenotype. Then, 'TRIR-1' was exposed to 60 Co γ-rays to induce mutations. Several individuals expressing variation in leaf traits were identified. After subculture and screening for 8 years, we isolated the stable mutant named 'TRIR-2'. 'TRIR-2' showed distinct phenotypic variation in leaf traits compared with the parental hybrid individuals. Pure-breeding seedlings of 'TRIR-1' and 'TRIR-2' were grown in a climate-controlled incubator (MGC-350HP-2; Yiheng Technical Co., Ltd, Shanghai, China) under long-day conditions (16 h light/8 h dark) at 22 C in the growth chamber. Several traits were measured and compared between 'TRIR-2' and non-irradiated 'TRIR-1' individuals grown under the same conditions, including plant height (the distance between the highest leaf tip and soil), number of leaves, length of leaves (the distance from tip to base of mature leaf) and breadth of leaves (the width of the widest part of mature leaf). Plant material cultivation was performed in the College of Horticulture and Landscape, Yunnan Agricultural University, Kunming, China.

Cloning and sequence analysis of CcMYB24
The expressed sequence tag of CcMYB24 (File S1) was aligned to NCBI using BLAST. As a result, the RAX2-like gene (query coverage: 100%; percentage identity: 96.81%) of Cymbidium ensifolium showed the greatest sequence similarity and was used to design degenerate primers to obtain the full-length sequence of CcMYB24 in 'TRIR-1' and 'TRIR-2'. The sequences of primers were listed in Table S1. The sequence of 'TRIR-1' and 'TRIR-2' was aligned to identify their sequence differences by using Multiple Sequence Alignment function in DNAMAN 8.0. Then, we used ORFfinder (https://www.ncbi.nlm.nih.gov/ orffinder) to obtain the coding sequence (CDS) and amino acid sequence of the corresponding protein, and aligned the latter to NCBI using BlastP to screen for homologs. The homology was compared using DNAMAN 8.0 software. Additionally, phylogenetic trees were constructed based on full-length amino acid sequences by the Maximum Likelihood method in MEGA 7.0 software, with bootstraps set to 1,000.

Prediction of physicochemical properties and structure of CcMYB24 protein
The conservation of CcMYB24 domains was examined using NCBI Conserved Domains (www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) (Yang et al., 2020). The online tool ExPASy-ProtParam (https://web.expasy.org/protparam/) and ExPASy Protscale analysis program (https://web.expasy.org/protscale/) was used to estimate the physicochemical parameters (including molecular weight, theoretical pI (isoelectric point), instability index, aliphatic index and grand average of hydropathy (GRAVY)), hydrophobicity and hydrophilicity values of the CcMYB24 protein, by providing sequence information and analyzing in the website following the user guide (Gasteiger et al., 2003). The TMHMM-2.0 (https://services.healthtech.dtu.dk/service.php?TMHMM-2.0) tool was adopted to predict the transmembrane helical region, which gives some statistics and a list of the location of the predicted transmembrane helices and the predicted location of the intervening loop regions. To predict the presence of a signal peptide, SignalP-5.0 software (https://services. healthtech.dtu.dk/service.php?SignalP-5.0) was used by uploading protein sequence and analyzing by online tools. The SOPMA software (http://npsa-pbil.ibcp.fr/cgi-bin/npsa_ automat.pl?page=npsa_sopma.html) and SWISS-MODEL (http://www.expasy.ch/ swissmod/SWISS-MODEL.html) software were used to predict the secondary structure and tertiary crystal structure following the online guideline of programs, respectively (Geourjon & Deléage, 1995;Arnold et al., 2006).

Subcellular localization
The subcellular localization of CcMYB24 was initially predicted using the pSORT web server (http://www.genscript.com/psort/wolf_psort.html). Then, the coding sequence of CcMYB24 was amplified by using the primers listed in Table S1 and cloned into the vector pCAMBIA1300-35S-N-GFP, which expresses the green fluorescent protein (GFP), using EcoR. and SalI restriction enzymes. The recombinant vector was introduced into competent Escherichia coli DH5-Alpha cells and then transformed into Agrobacterium tumefaciens EHA105. The Agrobacterium was injected into the leaves of four-week-old seedlings of Nicotiana benthamiana. Following the subsequent transient gene expression, subcellular localization of the CcMYB24-GFP fusion protein was observed using confocal laser scanning microscopy (LSM710; Carl Zeiss, Oberkochen, Germany) for 72 h after the injection (Olympus, Tokyo, Japan).

Vector construction and transformation of CcMYB24 in Arabidopsis
DNA fragments containing the full-length coding sequence of CcMYB24 were amplified from the DNA of 'TRIR-2' using gene-specific primers (Table S1) and then cloned into the pCAMBIA1300-35S vector between SacI and BamH. restriction sites. Then, the recombinant vector was introduced into A. tumefaciens strain EHA105 and transformed into Arabidopsis wild-type (WT) Columbia ecotype (Col-0) via the floral dip method (Clough & Bent, 1998). The seeds of transformed Arabidopsis were surface sterilized with ClO 2 for 8 min, rinsed thoroughly and selected for growth on selection medium containing 40 mg/L hygromycin. The Arabidopsis plants, including WT and genetically transformed individuals, were grown in incubator under long-day conditions (16 h light/8 h dark) at 23 C during the light time and 21 C during dark time, respectively. The presence of the gene insert was confirmed by PCR using 2×T5 DIRECT PCR Kit (Tsingke Biotech Co., Ltd., Beijing, China) with the primers listed in Table S1.

Quantification of CcMYB24 gene expression by qRT-PCR analysis
The expression of CcMYB24 gene was quantified by quantitative reverse transcription-PCR (qRT-PCR) in leaf of forty-day-old 'TRIR-2' and non-irradiated 'TRIR-1' Cymbidium seedlings, as well as in stem, leaf, and flower of genetically transformed and WT Arabidopsis. The RNAs were extracted by using the RNAprep Pure Plant Plus Kit (Tiangen, Beijing, China), following the manufacturer's instructions. For each sample, 1.0 µg of total RNA was reverse transcribed into cDNA using a PrimeScriptTM RT reagent kit (TaKaRa, Dalian, China). The qPCR was performed using SYBRÒ Premix Ex Taq Kit (TaKaRa, Dalian, China) and ABI7900 Fast Real-Time PCR System (ABI, Foster City, CA, USA). The 18S rRNA gene and UBQ5 were selected as the internal reference gene for TRIR-1/2 and Arabidopsis, respectively. The expression value was calculated by using the 2 −ΔΔCt method based on three biological replications for each sample. The relative expression was further calculated in Arabidopsis using log2 fold change of genetically transformed vs. WT plants (Li et al., 2017). Primers are listed in Table S1, and the Cq value for each sample was listed in Table S2.

Trait measurement in Arabidopsis
Several traits were measured and compared in genetically transformed and WT Arabidopsis. Morphological traits: number of leaves and branching were measured after 4 and 5 weeks of growth, respectively; physiological traits were measured after 4 weeks of growth: the concentration of photosynthetic pigments including chlorophyll a, chlorophyll b, total chlorophyll, and total carotenoids, were determined by the spectrophotometric method that homogenizing plant tissue in a solvent, centrifuging the homogenate, and analyzing the supernatant's absorbance at specific wavelengths to calculate the concentration of photosynthetic pigments (Jodłowska & Latała, 2011); the soluble sugar concentration was measured by the anthrone method, which involves adding anthrone reagent to the sample and heating it in a water bath. The colored complex formed between the anthrone and the sugar is then measured using a spectrophotometer, and the concentration of soluble sugars is calculated based on a standard curve (Clegg, 1956); the soluble protein concentration was determined by the Coomassie Brilliant Blue G-250 (Snyder & Desborough, 1978), which involves adding the Coomassie Brilliant Blue G-250 reagent to a protein sample, allowing it to bind to the protein, and measuring the absorbance of the resulting complex at a specific wavelength using a spectrophotometer. The protein concentration is then determined by comparing the absorbance of the sample to a standard curve of known protein concentrations. The total flavonoid concentration was measured by the Al(NO 3 ) 3 -NaNO 2 -NaOH colorimetric assay (Zheng et al., 2011), which involves the formation of a colored compound specific to flavonoids, and its absorbance is measured at a specific wavelength using a spectrophotometer.

Statistical analysis
All the data are presented as the mean ± standard error of at least six biological replicates for morphological traits and three biological replicates for physiological traits to confirm reproducibility. Statistical analyses were conducted using SPSS 22.0 (IBM, Armonk, NY, USA). Differences between two groups were calculated using Student's t-test. Statistical significance was defined as p < 0.05 ( Ã ) or p < 0.01 ( ÃÃ ). The raw data are listed in Tables  S3-S5.
The complete sequence of CcMYB24 was amplified and sequenced from 'TRIR-1' and 'TRIR-2'. The full-length sequence of CcMYB24 was 1389 bp, and the length of the open reading frame (ORF) was 906 bp, encoding 301 amino acids (File S1). The amplified sequence of 'TRIR-1' was determined to be 861 base pairs in length. Upon sequence alignment with the corresponding sequence of 'TRIR-2', it was observed that 'TRIR-1' exhibited a segmental variation of 95 base pairs, including three gaps with a deletion of 45 base pairs, and 50 single nucleotide polymorphisms (SNPs) when compared to 'TRIR-2'. These observed variations led to a loss of 15 amino acids and a substitution of 24 amino acids in the peptide chain of ccMYB24 in 'TRIR-1'. The alignment results are provided in File S2. Analysis of conserved domains showed that the protein encoded by CcMYB24 had a conserved Myb DNA-binding domain ( Fig. 2A). According to the results of BLAST, the ten sequences with the highest amino acid similarity to CcMYB24 were obtained, including CeRAX2-like (from Cymbidium ensifolium, 95.99% similarity), DcRAX2 (from Dendrobium catenatum, 86.15% similarity), PeRAX2-like (from Phalaenopsis equestris, 76.49% similarity), EgRAX2-like (from Elaeis guineensis, 60.80% similarity), PdMYB36 (from Phoenix dactylifera, 64.17% similarity), MeRAX2 (from Manihot esculenta, 52.51% similarity), JrRAX2-like (from Juglans regia, 51.39% similarity), QlRAX2-like (from Quercus lobata, 45.12% similarity), QsRAX2-like (from Quercus suber 50.53% similarity), NtRAX (from Nicotiana tabacum, 45.85% similarity) and AtRAX2 (from A. thaliana, 39.60% similarity). All of these genes belong to the R2R3-MYB family which includes most MYB TFs (Kelemen et al., 2015). The result of sequence alignments showed that all these sequences had conserved domains of the Myb superfamily proteins, such as the PLN03091 superfamily, REB1 superfamily, and two MYB DNA-binding domains (Fig. 2B). Phylogenetic analysis showed that the CcMYB24 protein is closely related to the CeRAX2like protein and the DcRAX2 protein (Fig. 2C).
The physicochemical properties and structure of CcMYB24 protein The ProtParam ExPaSy program was used to analyze the amino acid sequence and chemical properties of the CcMYB24 protein. The result showed that CcMYB24 possesses the molecular formula C 1457 H 2263 N 407 S 15 , and the molecular weight was calculated to be 33,242.35 Da. The theoretical isoelectric point (pI) was 8.59, and the instability index was 63.02, indicating that the protein is unstable. The grand average value of the hydropathy index (GRAVY) was calculated as the sum of the hydropathy values of all amino acids, and the estimated value of GRAVY was −0.483 (Fig. 3A). This finding suggests that the CcMYB24 protein is likely to be largely hydrophilic. The result of the transmembrane helix region and signal peptide prediction analysis showed that neither transmembrane structure nor signal peptide could be found in the CcMYB24 protein, indicating that it is not a transmembrane protein, nor a secretory protein. The secondary structure of the protein was predicted using the SOPMA server and shown in Fig. 3B. The result showed that CcMYB24 protein was mainly composed of random coil (52.49%), alpha-helix (24.58%), extended strand (16.28%) and beta turn (6.64%). The prediction of the tertiary protein structure of CcMYB24 (Fig. 3C) by the SWISS-MODEL analysis platform showed a similar structure with a R2R3-MYB transcription factor from A. thaliana which has six alpha-helices belonging to R2 and R3 repeats (Wang et al., 2020).

Subcellular localization of CcMYB24
Most R2R3-MYB proteins are exclusively localized in the nucleus, where they function as TFs (Stracke, Werber & Weisshaar, 2001). For CcMYB24, the result of localization prediction using the pSORT web server predicted that it was located in the nucleus. We further determined the subcellular localization of CcMYB24 via transient expression of a translation fusion with GFP in tobacco leaves. The result showed that the fluorescence signal of the GFP control was localized to the cytoplasm and the nucleus in N. benthamiana leaves, whereas the CcMYB24-GFP was located exclusively in the nucleus (Fig. 4), a finding similar to that of other R2R3-MYB TFs, suggesting that CcMYB24 also functions in the nucleus as a TF.

Expression pattern of CcMYB24 in 'TRIR-1', 'TRIR-2' and transgenic Arabidopsis
The expression pattern of CcMYB24 was examined in wild-type Cymbidium 'TRIR-1', mutant Cymbidium 'TRIR-2' and transgenic Arabidopsis by using quantitative real-time PCR (qPCR). As shown in Fig. 5A, the relative expression of CcMYB24 transcripts was about 1.4 time higher in 'TRIR-2' than in 'TRIR-1', indicating that more transcripts were generated in leaves of 'TRIR-2' than in those of the wild-type parent. In the case of Arabidopsis, the expression level was evaluated in terms of log2 fold change between transgenic and WT individuals. The result showed that CcMYB24 expression was much higher in the stem, leaf and flower of transgenic Arabidopsis than in the corresponding organs of the WT (Fig. 5B).

CcMYB24 induced production of more leaves when overexpressed in Arabidopsis
To investigate the function of CcMYB24 in 'TRIR-2', the coding region was expressed in A. thaliana Col-0. We observed some differences between the transgenic and WT individuals (Figs. 6A-6D). Despite a p-value of 0.06, statistical analysis revealed no significant difference in the number of branches between WT and transgenic individuals (Fig. 6E). The mean number of leaves in transgenic individuals was significantly higher than that in WT plants (p < 0.01, Fig. 6F).

CcMYB24 markedly affected the accumulation of photosynthetic pigments and total flavonoids
The influence of CcMYB24 on a number of physiological traits, namely photosynthetic pigments, total flavonoids, and soluble sugar and protein, were examined and compared between transgenic and WT Arabidopsis. For photosynthetic pigments, the results showed that chlorophyll a (p < 0.01), chlorophyll b (p = 0.02), and total chlorophyll concentrations (p < 0.01) were significantly higher in transgenic Arabidopsis (Fig. 7A), whereas total carotenoid concentration was significantly lower in transgenic Arabidopsis plants than in WT individuals (p < 0.01, Fig. 7A). The total flavonoid concentration was also significantly lower in transgenic Arabidopsis than in WT individuals (p < 0.01; Fig. 7B), but soluble sugar and protein concentrations did not differ significantly between transgenic and WT plants (Figs. 7C and 7D).

DISCUSSION
Cymbidium plants enjoy a prime position in the floriculture industry. The ornamental value of Cymbidiums is conferred not only by the large and colorful flowers and extended flower longevity, but also by its foliage and floral fragrance . Some Cymbidium varieties have long been cultivated as foliage plants which are grown for attractive leaves, and the combined character of leaf variegation and colorful flowers have an added advantage in the ornamental plant market. Variations in leaf traits, including color, shape and number of leaves, are vital for the quality of foliage plant. In this study, a new foliage variant of Cymbidium was generated by combining inter-specific hybridization with radiation mutation breeding. The new variant, 'TRIR-2', has significantly broader but shorter leaves compared with the parent, 'TRIR-1', and showed an increased number of leaves and a part-albino leaf phenotype, with a large area of white in the middle of the leaves with green only at the edges and tip of the leaf (Fig. 1A), making it an attractive ornamental plant. Several researches have explored mutation in leaf development including shape development, photosynthetic mechanisms, chloroplast development and other mechanisms of genetic control on color variation of leaves. These studies mainly focused on model plants or crops such as A. thaliana (Moon et al., 2008), Zea mays (Lonosky et al., 2004, Nicotiana tabacum (Wu et al., 2020), and Oryza sativa (Zheng et al., 2019). The knowledge of genetic control of number of leaves is still limited, especially in ornamental foliage plants which commercial value strongly influenced by the number of leaves. Recent studies have shown that transcription factors control many aspects of leaf developmental processes (Liu et al., 2018;Zhao et al., 2020;Heisler et al., 2022). The MYB transcription factors are versatile in regulation of development, stress responses as well as in leaf development (Jiang & Rao, 2020;Wang, Niu & Zheng, 2021;Wu et al., 2022), but the knowledge about the role it plays in control of number of leaves is still limited.
In the present study, a R2R3-MYB gene, CcMYB24, was reported to participate in producing more leaves in 'TRIR-2', a high-leaf-number mutant of Cymbidium.
The expression analysis showed that CcMYB24 exhibited a higher expression level in 'TRIR-2' than in 'TRIR-1'. The mechanism of expression level increasing is complex and remain elusive. It may be due the mutation of upstream and/or regulatory gene, or the difference in epigenetic control between 'TRIR-1' and 'TRIR-2'. More works are needed in the future to discover the reason why there were more transcripts of CcMYB24 in 'TRIR-2'. We plan to conduct comparative transcriptome and gene co-expression network analysis between 'TRIR-1' and 'TRIR-2' and further integrated the results with the analysis of epigenetic control of CcMYB24.
According to the results of sequence similarity analysis, CcMYB24 appeared to be a RAX2-like gene which belongs to the R2R3-MYB transcription factor gene family. Many studies have proved that the RAX homolog genes play important roles in determination of axillary meristems by generating a tissue environment conducive to the establishment of meristems to control the spatial pattern of axillary meristems development (reviewed by Zhang et al., 2022), thus participate in branching and leaf development. In Arabidopsis, transcripts of RAX genes accumulated in the axils of young leaf primordia and controlled a very early step of axillary meristem initiation, with rax mutants showing new phenotypes characterized by defects in lateral bud formation in overlapping zones along the shoot axis (Müller, Schmitz & Theres, 2006). Much evidence has also been reported that RAX genes control axillary meristem formation and shoot branching in Arabidopsis by collaborating in gene modules with other transcriptional factors, such as bHLH and WRKY (Yang et al., 2012;Guo & Qin, 2016). In Chrysanthemum morifolium, CmRAX2 was expressed in the lateral branches and roots, with overexpression of CmRAX2 leading to increases in lateral branch number and plant height . RAX genes were also be found to participate in shoot branching and leaf dissection in tomato (Busch et al., 2011;Zhang et al., 2022). Based on these evidences, it could be inferred that the CcMYB24 may be responsible for the variation in foliage traits shown by 'TRIR-2'.
The current studies showed the regulatory role played by RAX2 genes in shoot branching and leaf development. The results of transgenic experiment in this study showed that overexpression of CcMYB24 in Arabidopsis stimulated production of more leaves (Fig. 6E) and more branches (Fig. 6F), although the latter effect was not significant (but close with p-value of 0.06), indicating that CcMYB24 may contribute to the phenotype of more leaves in 'TRIR-2' compared with 'TRIR-1'. Orchid plants grow by only one or two growth patterns, namely sympodial or monopodial. Both types grow by leaves rather than obvious branches. Therefore, it can be inferred that the phenotype of more leaves in 'TRIR-2' may be related to the branching regulation nature of RAX2 genes. However, the factors contributing to the higher expression level of MYB24 in 'TRIR-2' than that in 'TRIR-1' are still unknown. More research focusing on the expression pattern and regulatory role at different leaf developmental stages is needed to explore the relationship between CcMYB24 and the formation of novel leaf traits in 'TRIR-2'.
The content of pigments strongly influences the color of leaves, and in this study, carotenoid concentration was found to be significantly lower in transgenic Arabidopsis than in WT plants (Fig. 7A). It is known that carotenoid mutants in Arabidopsis can cause developmental defects in plastids, resulting in albino patches on leaves (Bartley & Scolnik, 1995). However, the higher concentrations of chlorophyll a, chlorophyll b, and total chlorophyll in transgenic Arabidopsis indicated that the expression of CcMYB24 was not the direct cause of the albino phenotype in 'TRIR-2' leaves, as chlorophylls a and b are the primary sources of green color.
Several studies have shown that MYB TFs work as regulators of primary and secondary metabolism in plants, and are essential for regulating pathways involved in plant growth, development, and responses to biotic and abiotic stress (Wang, Niu & Zheng, 2021). Plant metabolites include many types of compounds, such as sugars, proteins, flavonoids etc. Evidence has shown that the MYB TFs participate in flavonoid biosynthesis, and that the expression of some MYB transcription factor genes could promote flavonoid accumulation in plants, which could result in increased tolerance to environmental stresses (Jin et al., 2022). In the present study, accumulation of total flavonoids in transgenic Arabidopsis was significantly lower than in WT individuals (Fig. 7B), indicating that CcMYB24 may negatively influence the accumulation of total flavonoids. Some MYB genes could also influence sugar transport and promote accumulation of soluble sugars and proteins. For instance, two MYB proteins, CSA2 and CSA, were proved to regulate sugar transport in rice (Wang et al., 2021). Overexpression of MYB1 in Arabidopsis resulted in accumulation of greater concentrations of soluble sugars in transgenic individuals than in the WT plants (Cheng et al., 2013). However, the concentration of soluble sugars and proteins was not affected by overexpression of CcMYB24 in Arabidopsis (Figs. 7C and 7D), indicating that CcMYB24 may not be involved in the determination of the concentrations of soluble sugars and proteins. In addition, there were differences in traits exhibited by other plant organs between 'TRIR-1' and 'TRIR-2', such as root number and plant height ( Fig. 1A). As MYB TFs play roles in many aspects of plant development (Ambawat et al., 2013;Li et al., 2019;Jiang & Rao, 2020), more research is needed to explore the function of CcMYB24 in 'TRIR-2'.
Overexpression of the CcMYB24 gene significantly increased the number of leaves in transgenic Arabidopsis plants, indicating it may play a positive role in controlling the number of leaves in 'TRIR-2'. The number of leaves is an important trait that directly affects the quality of ornamental foliage plants, with the ornamental value of foliage plants usually being increased by increasing the number of leaves. Therefore, CcMYB24 can be considered as a potential target gene to regulate the number of leaves in molecular breeding of ornamental foliage plants.
Longjie Cheng performed the experiments, analyzed the data, prepared figures and/or tables, and approved the final draft. Zhilin Li conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft. Yiran Zhao performed the experiments, analyzed the data, prepared figures and/or tables, and approved the final draft. Yuying Wang conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft.

DNA Deposition
The following information was supplied regarding the deposition of DNA sequences: The CcMYB24 sequence is available in the Supplemental Files and at GenBank: OQ341202.