Genome-Wide Identification, Expression Analysis, and Potential Roles under Abiotic Stress of the YUCCA Gene Family in Mungbean (Vigna radiata L.)

YUCCA, belonging to the class B flavin-dependent monooxygenases, catalyzes the rate-limiting step for endogenous auxin synthesis and is implicated in plant-growth regulation and stress response. Systematic analysis of the YUCCA gene family and its stress response benefits the dissection of regulation mechanisms and breeding applications. In this study, 12 YUCCA genes were identified from the mungbean (Vigna radiata L.) genome and were named based on their similarity to AtYUCCAs. Phylogenetic analysis revealed that the 12 VrYUCCAs could be divided into 4 subfamilies. The evidence from enzymatic assays in vitro and transgenetic Arabidopsis in vivo indicated that all the isolated VrYUCCAs had biological activity in response to IAA synthesis. Expression pattern analysis showed that functional redundancy and divergence existed in the VrYUCCA gene family. Four VrYUCCAs were expressed in most tissues, and five VrYUCCAs were specifically highly expressed in the floral organs. The response toward five stresses, namely, auxin (indole-3-acetic acid, IAA), salinity, drought, high temperatures, and cold, was also investigated here. Five VrYUCCAs responded to IAA in the root, while only VrYUCCA8a was induced in the leaf. VrYUCCA2a, VrYUCCA6a, VrYUCCA8a, VrYUCCA8b, and VrYUCCA10 seemed to dominate under abiotic stresses, due to their sensitivity to the other four treatments. However, the response modes of the VrYUCCAs varied, indicating that they may regulate different stresses in distinct ways to finely adjust IAA content. The comprehensive analysis of the VrYUCCAs in this study lays a solid foundation for further investigation of VrYUCCA genes’ mechanisms and applications in breeding.


Introduction
Auxin/indole-3-acetic acid (IAA) is one of the major endogenous phytohormones, playing a critical role in regulating plant growth and development as well as the stress response [1], and it results from the coordinated actions of auxin biosynthesis, metabolism, transport, and signaling [2]. During the plant life cycle, IAA adjusts almost all aspects of biological processes, including embryogenesis [3], shoot and root formation [4], leaf development [5], floral development and fertility [6,7], and fruit initiation [8]. Additionally, auxin homeostasis is essential for plants to adapt to biotic and abiotic stresses, such as pathogens, herbivory, temperature, salt, and drought [1,9,10].
Mungbean is an important warm-season grain legume with rich nutrition and short duration, facilitating environmental sustainability [36]. Our previous study showed that a YUC4-like protein regulated the chasmogamous flower trait in mungbean [37], indicating the crucial roles of YUCCA genes in developmental regulation. In this study, we systematically identified 12 YUCCA family genes (named as VrYUCCA) in the mungbean genome and investigated their gene structures, expression patterns, biological activity in vitro and in vivo, and potential roles under abiotic stresses. This research can provide valuable information about the functional explanation of the YUCCA gene family in mungbean, contributing to the application of these genes in mungbean genetic improvement.

Identification of the YUCCA Gene Family in Mungbean and Structure Analysis
Based on the output of the Gcorn plant, BLAST homology search, and conserved domain analysis, 12 VrYUCCA genes were isolated from the mungbean genome (Table S1), named according to their similarity to Arabidopsis YUCCA genes [16]. No homologs for AtYUC1, AtYUC5, AtYUC7, or AtYUC9 were found in mungbean. Meanwhile, there were three homologs for AtYUC8 and two each for AtYUC2, AtYUC4, and AtYUC6. The deduced amino acid lengths of the 12 VrYUCCA proteins ranged between 385 and 440 aa, with the pI between 7.96 and 9.54 and the molecular weight between 43.14 and 49.89 kD.
To investigate the diversity of the gene structures of the 12 VrYUCCAs, the exon and intron configurations were determined through visualized analysis using the GSDS tool. VrYUCCA3, VrYUCCA8a, VrYUCCA8b, and VrYUCCA8c, which were clustered in the same group, each contained three exons and two introns. The other VrYUCCAs all had four exons and three introns ( Figure 1A). The conserved protein motifs were analyzed by MEME, and a total of 12 motifs are exhibited here ( Figure 1A). All six of the recognized conserved motifs, including the FAD-binding motif, GC motif, ATG-containing motif1, FMO-identifying sequence, NADPH-binding motif, and ATG-containing motif2, are shown in a both schematic diagram and sequence alignment ( Figure 1B and S1). duration, facilitating environmental sustainability [36]. Our previous study showed that a YUC4-like protein regulated the chasmogamous flower trait in mungbean [37], indicating the crucial roles of YUCCA genes in developmental regulation. In this study, we systematically identified 12 YUCCA family genes (named as VrYUCCA) in the mungbean genome and investigated their gene structures, expression patterns, biological activity in vitro and in vivo, and potential roles under abiotic stresses. This research can provide valuable information about the functional explanation of the YUCCA gene family in mungbean, contributing to the application of these genes in mungbean genetic improvement.

Identification of the YUCCA Gene Family in Mungbean and Structure Analysis
Based on the output of the Gcorn plant, BLAST homology search, and conserved domain analysis, 12 VrYUCCA genes were isolated from the mungbean genome (Table  S1), named according to their similarity to Arabidopsis YUCCA genes [16]. No homologs for AtYUC1, AtYUC5, AtYUC7, or AtYUC9 were found in mungbean. Meanwhile, there were three homologs for AtYUC8 and two each for AtYUC2, AtYUC4, and AtYUC6. The deduced amino acid lengths of the 12 VrYUCCA proteins ranged between 385 and 440 aa, with the pI between 7.96 and 9.54 and the molecular weight between 43.14 and 49.89 kD.
To investigate the diversity of the gene structures of the 12 VrYUCCAs, the exon and intron configurations were determined through visualized analysis using the GSDS tool. VrYUCCA3, VrYUCCA8a, VrYUCCA8b, and VrYUCCA8c, which were clustered in the same group, each contained three exons and two introns. The other VrYUCCAs all had four exons and three introns ( Figure 1A). The conserved protein motifs were analyzed by MEME, and a total of 12 motifs are exhibited here ( Figure 1A). All six of the recognized conserved motifs, including the FAD-binding motif, GC motif, ATG-containing motif1, FMO-identifying sequence, NADPH-binding motif, and ATG-containing motif2, are shown in a both schematic diagram and sequence alignment ( Figures 1B and S1).

Phylogenetic Analysis of YUCCAs from Multiple Species
To further describe the evolutionary relationship, we performed an unrooted phylogenetic tree of YUCCA proteins from mungbean, Arabidopsis, soybean, rice, and Medicago ( Figure 2). A total of 75 YUCCA proteins from five plant species were classified into four subfamilies (I to IV), which was consistent with previous research on Arabidopsis [16] and soybean [19]. VrYUCCA2a, VrYUCCA2b, VrYUCCA6a, and VrYUCCA6b were grouped together in subfamily I; VrYUCCA4a and VrYUCCA4b were placed in subfamily II; VrYUCCA3, VrYUCCA8a, VrYUCCA8b, and VrYUCCA8c were located in subfamily III; and VrYUCCA10 and VrYUCCA11 were placed in subfamily IV. At the same time, four rice YUCCA proteins (OsYUCCA1, OsYUCCA4, OsYUCCA5, and OsYUCCA8) did not cluster well with the YUCCA proteins in subfamily II. Therefore, these four proteins could also be classified into another new subfamily, subfamily V. Among these five analyzed plant species, the evolutionary relationship of VrYUCCA proteins was closest to GmYUCCAs. The sequences information of all the YUCCAs is listed in Supplementary Data S1. and soybean [19]. VrYUCCA2a, VrYUCCA2b, VrYUCCA6a, and VrYUCCA6b w grouped together in subfamily I; VrYUCCA4a and VrYUCCA4b were placed in subfa ily II; VrYUCCA3, VrYUCCA8a, VrYUCCA8b, and VrYUCCA8c were located in s family III; and VrYUCCA10 and VrYUCCA11 were placed in subfamily IV. At the sa time, four rice YUCCA proteins (OsYUCCA1, OsYUCCA4, OsYUCCA5, and OsYU CA8) did not cluster well with the YUCCA proteins in subfamily II. Therefore, these f proteins could also be classified into another new subfamily, subfamily V. Among th five analyzed plant species, the evolutionary relationship of VrYUCCA proteins w closest to GmYUCCAs. The sequences information of all the YUCCAs is listed in S plementary Data S1.

Expression Pattern of VrYUCCAs in Mungbean
To explore the potential biological functions of VrYUCCA genes during mungbean development, we analyzed their expression patterns in different tissues, including seedling root (SR), seedling true leaf (STL), compound leaf (CL), stem (St), young inflorescence (YI), flower bud (FB), opening flower (OF), 7-day pod (P7), and seed (Se), using real-time PCR ( Figure 3). Several VrYUCCAs were constitutively expressed in all tissues examined, including VrYUCCA2a, VrYUCCA6b, VrYUCCA8a, and VrYUCCA10. Most YUCCAs had an extremely low expression in the stems except for VrYUCCA2a, VrYUCCA4a, and VrYUCCA8a; among them, VrYUCCA4a was specifically expressed. VrYUCCA3 was preferentially expressed in seed. There were five VrYUCCAs specifically enriched in flower organs, including VrYUCCA4b (IY), VrYUCCA6a (OF), VrYUCCA8b (OF), VrYUCCA8c (OF), and VrYUCCA11 (FB and OF). These results suggest that functional redundancy and divergence exist in the VrYUCCA gene family, and during different stages of growth or organs, distinct VrYUCCAs might play a major role. VrYUCCA4a, and VrYUCCA8a; among them, VrYUCCA4a was specifically expressed VrYUCCA3 was preferentially expressed in seed. There were five VrYUCCAs specifically enriched in flower organs, including VrYUCCA4b (IY), VrYUCCA6a (OF), VrYUCCA8 (OF), VrYUCCA8c (OF), and VrYUCCA11 (FB and OF). These results suggest that func tional redundancy and divergence exist in the VrYUCCA gene family, and during dif ferent stages of growth or organs, distinct VrYUCCAs might play a major role. . Expression pattern of VrYUCCAs in different tissues. SR, seeding root; STL, seedling tru leaf; CL, compound leaf; St, stem; YI, young inflorescence; FB, flower bud; OF, opening flower; P7 7-day pod; Se, seed. The VrACTIN3 gene (Vradi03g00210) was used as the internal control to nor malize the real-time PCR data. The 2 −△CT method was used to analyze relative expression leve Error bars indicated SDs (standard deviations) from three biological repetitions.

VrYUCCA Proteins Showed Biological Activity Both In Vitro and In Vivo
It was reported that YUCCA proteins might function in two possible routes to pro duce IAA, with IPA or TAM as the substrates, respectively [13,25]. To investigate th catalytic substrate of VrYUCCA proteins, enzyme activity detection was conducted in vitro. The results indicated that all 12 VrYUCCAs could catalyze IPA to produce IAA (Figures 4 and S2), but they showed no activity toward TAM ( Figure S3). Thus, in mungbean, VrYUCCAs participated in the IPA pathway to synthesize IAA. Further, t examine their biological activity in vivo, VrYUCCAs were overexpressed in Arabidopsi (Col-0, wild type; WT) and driven by the 35S strong promoter. The T1 generation of mos VrYUCCAs transgenic plants showed typical auxin-related phenotypes ( Figure 5), such as narrow rosette leaves, downward curling of the leaf-blade margin, twisted stem, and shorter or no siliques. What was noteworthy was that 35S::VrYUCCA10-transgenic plant showed a distinct phenotype, as their leaves were broader and slightly curled, with denser siliques on the top of the stem. These results indicated that all 12 of the isolated VrYUCCA genes had a veritable role in IAA synthesis.

VrYUCCA Proteins Showed Biological Activity Both In Vitro and In Vivo
It was reported that YUCCA proteins might function in two possible routes to produce IAA, with IPA or TAM as the substrates, respectively [13,25]. To investigate the catalytic substrate of VrYUCCA proteins, enzyme activity detection was conducted in vitro. The results indicated that all 12 VrYUCCAs could catalyze IPA to produce IAA (Figure 4 and S2), but they showed no activity toward TAM ( Figure S3). Thus, in mungbean, VrYUCCAs participated in the IPA pathway to synthesize IAA. Further, to examine their biological activity in vivo, VrYUCCAs were overexpressed in Arabidopsis (Col-0, wild type; WT) and driven by the 35S strong promoter. The T 1 generation of most VrYUCCAs transgenic plants showed typical auxin-related phenotypes ( Figure 5), such as narrow rosette leaves, downward curling of the leaf-blade margin, twisted stem, and shorter or no siliques. What was noteworthy was that 35S::VrYUCCA10-transgenic plants showed a distinct phenotype, as their leaves were broader and slightly curled, with denser siliques on the top of the stem. These results indicated that all 12 of the isolated VrYUCCA genes had a veritable role in IAA synthesis.

VrYUCCAs Expression in Response to Abiotic Stress Treatments
To investigate the response of VrYUCCAs toward surroundings stresses, 12-day-old mungbean seedlings were exposed to IAA (1 µM), high salinity (200 mM NaCl), drought (15% PEG6000), heat (35 • C), and chilling (4 • C) for some time. The expression of VrYUCCA11 could not be detected well in the root and leaf, even under these stresses, due to its specific expression in flower organ. As shown in Figure 6, VrYUCCA2b, VrYUCCA6b, VrYUCCA8a, and VrYUCCA10 were inhibited markedly in the root by IAA, while VrYUCCA8b was induced. In the leaf, only the expression of VrYUCCA8a showed a highly significant increase. long-term stress, such as VrYUCCA4a and VrYUCCA8b in the leaf. Additionally, VrYUCCA2a, VrYUCCA6b, VrYUCCA8c, and VrYUCCA10 mainly responded to salinity in the root. Under drought stress (Figures 7 and S5), VrYUCCA6a was inhibited, while VrYUCCA8a and VrYUCCA8b were induced in the leaf. In the root, however, VrYUC-CA2a, VrYUCCA8a, VrYUCCA8c, and VrYUCCA10 might play a crucial part. When the mungbean seedlings were exposed to high temperatures, the expression of VrYUCCA2a, VrYUCCA6a, and VrYUCCA6b markedly increased, while that of VrYUCCA8b and VrYUCCA8c was inhibited. Cold prominently induced the expression of VrYUCCA8a and VrYUCCA8b. Meanwhile the expression of VrYUCCA2a and VrYUCCA10 declined (Figures 8 and S6). Overall, VrYUCCA2a, VrYUCCA6a, VrYUCCA8a, VrYUCCA8b, and VrYUCCA10 proved to be sensitive to all four treatments and might play a major role under abiotic stresses.  As for the other abiotic stresses, all the analyzed VrYUCCA genes showed a visible increase and decrease under high-salinity treatment (Figure 7 and S4). Both at 6 h and 24 h after NaCl treatment, VrYUCCA4b and VrYUCCA8a had high expression in the leaf; meanwhile, the expression of VrYUCCA2b and VrYUCCA6a was inhibited significantly, indicating their main roles under salt stress. Several genes seemed to only function under long-term stress, such as VrYUCCA4a and VrYUCCA8b in the leaf. Additionally, VrYUCCA2a, VrYUCCA6b, VrYUCCA8c, and VrYUCCA10 mainly responded to salinity in the root. Under drought stress (Figure 7 and S5), VrYUCCA6a was inhibited, while VrYUCCA8a and VrYUCCA8b were induced in the leaf. In the root, however, VrYUCCA2a, VrYUCCA8a, VrYUCCA8c, and VrYUCCA10 might play a crucial part. When the mungbean seedlings were exposed to high temperatures, the expression of VrYUCCA2a, VrYUCCA6a, and VrYUCCA6b markedly increased, while that of VrYUCCA8b and VrYUCCA8c was inhibited. Cold prominently induced the expression of VrYUCCA8a and VrYUCCA8b. Meanwhile the expression of VrYUCCA2a and VrYUCCA10 declined (Figure 8 and S6). Overall, VrYUCCA2a, VrYUCCA6a, VrYUCCA8a, VrYUCCA8b, and VrYUCCA10 proved to be sensitive to all four treatments and might play a major role under abiotic stresses.

Discussion
YUCCA, a key rate-limiting enzyme for IAA synthesis, has been widely characterized in quite a few higher plants [11], existing in forms of families in different organisms, with numbers ranging from a few (e.g., 6 in tomato [38]) to dozens (e.g., 57 in Brassica napus [24] and 63 in wheat [14]), except for petunia, which only has one YUCCA gene in its genome [39]. Reports about YUCCA genes in mungbean are limited. Previously, a YUC4-like protein (Vradi06g12650) was reported to be responsible for the mungbean chasmogamous mutant (CM) [37], indicating the vital role of VrYUCCAs during plant growth and development. In this study, we identified 12 VrYUCCAs in the mungbean genome, named based on their similarity to AtYUCCAs. Homologous genes of AtYUC1, AtYUC5, AtYUC7, and AtYUC9 were absent, but two each for AtYUC2, AtYUC4, and AtYUC6 and three for AtYUC8 were isolated in mungbean (Figure 1). This phenomenon also exists in apple and cucumber [2,22], which seems normal for gene family studies, and chromosomal fragment deletion and genetic duplication might be the causative factors of this phenomenon. Further, multiple copies of the VrYUCCAs displayed distinct expression patterns (Figure 3), indicating their functional differentiation during evolution.
The 12 VrYUCCA proteins could be divided into four groups (Figure 2), possessing the conserved motifs of the YUCCA gene family. Subfamily I (YUC2-YUC6 group) and subfamily II (YUC1-YUC4 group) were reported to function in the inflorescence apex and flower in Arabidopsis, and AtYUC4 was specifically expressed in the floral meristems [16]. We previously obtained one mungbean mutant related to VrYUCCA4b that had developmental defects in its flower organs [37]. VrYUCCA4b (Vradi06g12650) seems to have a similar role in floral development, which is also highly expressed in young inflorescence (Figure 3). In Arabidopsis, the Atyuc4 single mutant displayed the normal phenotype, but Atyuc1yuc4 showed severe defects in its floral organs and fertility [16]; in mungbean, the homologous genes of AtYUC1 were vacant, and other VrYUCCAs including VrYUCCA4a seemed not to complement the function of VrYUCCA4b, resulting in the petal-loss trait of CM (the Vryucca4b mutant) [37]. Thus, YUCCA paralogous genes have functional redundancy to ensure normal plant development. Additionally, in Isatis indigotica Fort, IiYUCCA1, and IiYUCCA4 were also mainly detected in reproductive organs [40]. In Brassica napus, BnYUC4 was also highly expressed in early flower buds [24]. Thus, the expression pattern and function of YUCCAs seem to be conserved. It was reported that subfamily III (YUC3-7-8-5-9 group) showed distinct expression patterns during root development in Arabidopsis [41,42]. Interestingly, in the root defects of these five YUCs (YUC3-7-8-5-9), a multiple mutant (yucQ) could be rescued by expressing a YUCCA gene in subfamily III rather than in other subfamily members [41]. At the same time, in mungbean, a high transcript of VrYUCCA3 was detected in its seeds (Figure 3), implying a likely vital role during seed formation. VrYUCCA8b and VrYUCCA8c seemed like gene clusters located adjacently on chromosome 8, and both were expressed extremely highly during flower opening. However, VrYUCCA8a had a relatively high expression in the root. VrYUCCA10 and VrYUCCA11 were found in subfamily IV (YUC10-11 group), which was assumed to be the most conserved group of YUCCAs in the plant kingdom [42]. AtYUC10 and AtYUC11 mainly expressed in siliques and apices and were suggested to have overlapping functions with AtYUC1 and AtYUC4 during embryogenesis [26]. In mungbean, the relative transcript level of VrYUCCA10 was the highest in all tissues (Figure 3), hinting at its essential role during the whole life cycle. At the same time, the transcript of VrYUCCA11 was examined only in the floral organs and pods (Figure 3). In plants, YUCCA11 seemed to be specifically involved in the floral organs, embryogenesis, and fruit development. In cucumber, CsYUC11 was specifically enriched in the opened male flower (mainly detected in the tapetum cell layer and the microspores in the anthers), and transgenic Arabidopsis lines of 35S::CsYCU11 displayed defective pollens and reduced fertility [2]. In peach, the transcript of PpYUC11 increased during the late ripening stage in melting flesh peaches, implying its vital function during fruit ripening [43]. The number of YUCCA11 homologous genes increased in several Rosaceae species, such as apple (six YUCCA11) and white pear (seven YUCCA11). MdYUCCA11a, MdYUCCA11b, and MdYUCCA11d showed high expression in young fruit, the receptacle, and the pistil, respectively [22]. Thus, like YUCCA in other species, distinct VrYUCCA participated in specific biological processes with redundancy and differential functions.
In this study, we verified 12 VrYUCCA proteins' enzyme activity both in vivo and in vitro. The substrate of YUCCA was controversial. Initially, it was demonstrated that the Arabidopsis YUCCA1 protein could catalyze the N-hydroxylation of TAM to form N-hydroxy-TAM in vitro [25]. Meanwhile, later in vivo research showed that YUCCA1 also catalyzed the oxidative decarboxylation of IPA to form IAA [44,45]. Subsequently, a growing number of studies have suggested that the real substrate of YUCCA was IPA rather than TAM [13,46,47]. Here, VrYUCCA exhibited explicit activity toward IPA, but no novel compound was produced toward TAM in vitro (Figure 4 and S3). The 35S::VrYUCCAs transgenetic Arabidopsis lines all showed typical IAA over-accumulation defects ( Figure 5), indicating the authentic function of the VrYUCCAs for IAA synthesis in vivo. No seeds were harvested from the 35S::VrYUCCA11 line, which displayed the most severe abnormity during T 1 generation. Notable was, the 35S::VrYUCCA10 line displayed less of an effect in the vegetative growth stage (even stronger than WT), but the siliques showed apex density. Whether this line has special potential advantages during developmental processes or under stress need to be studied further.
Auxin, a central regulator involved in almost all aspects of the plant life cycle, helps plants to adaptive external stresses [1,48]. YUCCA protein catalyzes IPA to produce IAA, and IAA can also feedback modulate YUCCA expression. Here, five VrYUCCA genes responded to IAA stimulation, which mainly functions in the root. Notably, only VrYUCCA8a was induced in the leaf, but it was significantly inhibited in the root (Figure 6), indicating that it might control gene expression by regulating hormone signaling through complex crosstalk, as reported in soybean [19]. In apple and soybean, the homologous gene of VrYUCCA8a also displayed sensitivity to hormone treatments [19,22], indicating its functional conservation in the plant kingdom. Additionally, YUCCA was also proven to be involved in IAA-mediated stress response [2,9,31]. In this study, we investigated the potential roles of VrYUCCAs in stress response. More VrYUCCAs were sensitive under salt stress than drought treatment (Figure 7, S5 and S6). VrYUCCAs showed different expression patterns to respond to the same stress; for example, leaf VrYUCCA2b and VrYUCCA6a were inhibited markedly under high salinity, while the expression of VrYUCCA4b and VrYUCCA8a increased significantly. These results hinted that the regulation mechanism of the stress response by YUCCA-mediated IAA synthesis is complex and can even antagonize to create a buffering system for endogenous auxin synthesis, as reported in cucumber [2]. Regarding heat and chilling, VrYUCCA2a, VrYUCCA6a, VrYUCCA6b, VrYUCCA8a, and VrYUCCA8b might act as the main regulators. In Arabidopsis and apple, PIF4 could specifically bind to the G-box motif (CACGTG) in the promoter of YUCCA8 to activate its expression, which was deemed a conserved molecular mechanism [22,30]. In cucumber, CsYUC8 and CsYUC9 were selectively upregulated after being exposed to 38 • C for 3 h [2]. Remarkably, VrYUCCA8a and VrYUCCA8b were markedly induced by low temperatures rather than heat, similar to that in Populus [9], suggesting that high-temperature regulatory mechanisms might be diverse in mungbean. It is possible that several VrYUCCAs have a dual role, like AtYUC6 [29], and transmit stress signals through another functional domain mediated by ROS or hormone crosstalk, such as ABA and SA, which needs further exploration.
Overall, the YUCCA gene family is relatively conservative among plants, and the homologous genes in different species tend to act similarly in their roles during plant growth and stress responses. The crucial role of YUCCAs in the auxin pathway is beyond doubt, so further in-depth research should be conducted to investigate its functional mechanism and application value.

Plant Materials and Treatment
The seedlings of mungbean cultivar 'Sulyu 1' were grown in the illumination incubator at 25 • C with 16 h light/8 h dark photoperiod. Adult mungbean was planted in experimental plots at Jiangsu Academy of Agricultural Sciences, Nanjing, China (32 • 21 N, 118 • 53 E).
Seedling root, seedling true leaf, stem, compound leaf, young inflorescence, flower bud, opening flower, pod (7 d), and dry seed were sampled from mungbean seedlings (12 d) or adult plants. All the tissues were immediately frozen in liquid nitrogen, and each sample was collected independently three times for further VrYUCCA expression pattern analysis using real-time PCR.
For IAA treatment, 12-day-old mungbean seedlings were soaked in liquid MS (Murashige and Skoog) medium with 1 µM IAA and liquid MS medium only (as the control) for 4 h. The treated roots and leaves were sampled separately after flushing with ddH2O and immediate freezing in liquid nitrogen. Treatment was repeated three times. The expression of VrYUCCA genes was then examined by real-time PCR.
For high-salinity and drought treatments, 12-day-old mungbean seedlings were soaked in liquid MS medium with 200 mM NaCl and 15% (w/v) polyethylene glycol (PEG) 6000 for 6 h and 24 h, respectively. The seedlings soaked in liquid MS medium served as the control, and the treated roots and leaves were sampled separately after flushing with ddH2O and immediate freezing in liquid nitrogen. Treatment was repeated three times. Then, the expression of VrYUCCA genes was analyzed using RNA-seq (BGI, Shenzhen, China).
For high-and low-temperature treatments, 12-day-old mungbean seedlings were placed in an incubator at 35 • C and 4 • C, respectively, for 24 h. The seedlings grown under 25 • C were taken as a control. Additionally, the treated leaves were sampled and immediately frozen in liquid nitrogen. The expression of VrYUCCA genes was detected using RNA-seq (BGI, Shenzhen, China).

Expression Analysis Using Real-Time PCR
Total RNA was extracted using Polysaccharide polyphenol Plant total RNA Extraction Kit (PD Biotech, Shanghai, China), and cDNAs were synthesized using Goldenstar ® RT6 cDNA Synthesis Kit Ver.2 (TSINGKE, Beijing, China). Real-time PCR was performed using ChamQ SYBR qPCR Master Mix (Vazyme, Nanjing, China) on an ABI prism 7500 real-time PCR System. Mungbean VrACTIN3 (Vradi03g00210) and Arabidopsis ACTIN2 (AT3G18780) were selected as reference genes to normalize the expression data. The 2 −∆CT or 2 −∆∆CT methods were used to analyze relative expression levels. The results were calculated from three biological replicates and three technical replicates. Statistical analysis and plot drawing were conducted using Excel and GraphPad Prism 5 software. Heat maps were drawing by TBtools after the FPKM value taken logarithm (LOG 2). The primers used in this study are listed in Table S2.

Prokaryotic Expression, Purification, and Enzymatic Assays of VrYUCCA
The full-length CDS sequence of each VrYUCCA gene was ligated into a pMAL-c2x expression vector (NEB, MBP-tag). Then, the recombinant VrYUCCA proteins were expressed in E. coli (BL21, DE3) and purified by dextrin Sepharose affinity chromatography (Amylose resin, NEB), in accordance with the instructions described previously [50]. The protein concentration was determined using the Bradford method with BSA as the standard.
The enzymatic assays were conducted by referring to a previous study [47]. A total 40 µL reaction system containing 1~2 µg VrYUCCA protein, 50 mM Tris-HCl (pH 7.4), 10 mM NADPH/NADP, and 2 mM of the substrates IPA or TAM was incubated at 30 • C for 4 h. A mixture containing one boiled VrYUCCA protein was served as the control. The products were analyzed on silica-gel-coated TLC plates with chloroform/methanol/acetic acid (75:20:5) as the mobile phase. The TLC plates were stained using Ehrlich's reagent (Sigma, 02560-500ML, St. Louis, MO, USA). Meanwhile, the products with IPA as the substrate were also detected by ESI-HPLC-MS/MS conducted by Nanjing WEBiolotech Biotechnology Co., Ltd. (Nanjing, China).

Arabidopsis Transformation
The CDS sequence of each VrYUCCA gene was cloned into a pCAMBIA1305.1-GFP vector through XbaI and BamHI sites, to create a fusion construct under the control of the CaMV 35S promoter. The fusion vectors 35S::VrYUCCAs-GFP were transformed into Agrobacterium tumefaciens strain GV3101 and then were introduced into the Arabidopsis thaliana (Col-0) wild type (WT) using the floral dipping method described previously [51]. The seeds of transgenic Arabidopsis were screened on 1/2 MS culture media containing 15 mg/L hygromycin antibiotics. Phenotypic analysis was performed in the progeny-positive T 1 lines.

Conclusions
In this study, 12 VrYUCCAs were systematically identified in mungbean, containing the six conserved motifs of the YUCCA gene family. The 12 VrYUCCA genes could be divided into four subfamilies according to the phylogenetic analysis. The VrYUCCAs showed different tissue expression patterns, and five had a highly specific transcript in the floral organs, including VrYUCCA4b, VrYUCCA6a, VrYUCCA8b, VrYUCCA8c, and VrYUCCA11. VrYUCCA2a, VrYUCCA6b, VrYUCCA8a, and VrYUCCA10 seemed to regulate the IAA content throughout the whole growing period, due to the constitutive expression pattern in all the tissues. The VrYUCCA proteins participated in the IPA pathway for IAA synthesis in vitro and exhibited veritable biological activity to the IAA product in vivo, which was certified by 35S::VrYUCCAs transgenetic Arabidopsis. Five YUCCAs in the root responded to IAA, but only VrYUCCA8a responded in the leaf. Regarding the other four stresses, VrYUCCA2a, VrYUCCA6a, VrYUCCA8a, VrYUCCA8b, and VrYUCCA10 might be mainly responsible for auxin regulation, due to their sensitivity to all four treatments used in this study.