Genome-Wide Identification and Expression Profiling Analysis of WOX Family Protein-Encoded Genes in Triticeae Species

The WOX family is a group of plant-specific transcription factors which regulate plant growth and development, cell division and differentiation. From the available genome sequence databases of nine Triticeae species, 199 putative WOX genes were identified. Most of the identified WOX genes were distributed on the chromosomes of homeologous groups 1 to 5 and originated via the orthologous evolution approach. Parts of WOX genes in Triticum aestivum were confirmed by the specific PCR markers using a set of Triticum. durum-T. aestivum genome D substitution lines. All of these identified WOX proteins could be grouped into three clades, similar to those in rice and Arabidopsis. WOX family members were conserved among these Triticeae plants; all of them contained the HOX DNA-binding homeodomain, and WUS clade members contained the characteristic WUS-box motif, while only WUS and WOX9 contained the EAR motif. The RNA-seq and qPCR analysis revealed that the TaWOX genes had tissue-specific expression feature. From the expression patterns of TaWOX genes during immature embryo callus production, TaWOX9 is likely closely related with the regulation of regeneration process in T. aestivum. The findings in this study could provide a basis for evolution and functional investigation and practical application of the WOX family genes in Triticeae species.


Introduction
The Triticeae tribe belongs to the Poaceae family and is made up of more than 350 plant species in 35 genera [1]. In Triticeae plants, 111 annual species in 19 genera such as Triticum aestivum (bread wheat), Hordeum vulgare (barley), Secale cereale (rye), Avena sativa (Oat), Triticum urartu, Triticum dicoccoides, Triticum turgidum, and so on, have been cultivated as crops, which provide necessary nutrition for more than two billion people in the world [1,2]. For a long time, people are interested in understanding the origination, genetic basis and evolution of Triticeae plants for their improvement. The assembling of wheat genome is a milestone in interpreting the genetic information of Triticeae plants. However, due to the large genome size of greater than 16 Gb, the genomic study on wheat has lagged behind that of rice and maize [3]. The application of modern biotechnology tools such as transgene and gene editing in plant breeding can help us to increase yield, improve quality, and enhance biotic or abiotic resistance of major crops, but the realization of these aims depends on genetic transformation. The ability of regenerating new plantlets from in vitro tissues is an obstacle that restricts the application of genetic transformation and gene editing systems [4,5].
Regeneration ability is one of many important genetically physiological traits for most plants, which enables plants to recover from wound tissues and form new organs. For modifying plants using a genetic-engineering strategy, shoot or somatic embryo production from isolated tissues or cells is an indispensable step to achieve transgenic plants. However, it is still difficult to obtain regenerated plants in the process of genetic transformation from most genotypes (especially the extensively cultivated varieties) of wheat and other Triticeae species [5][6][7][8][9]. During plant regeneration, a series of genes are expressed in an orderly manner under the regulation of auxin and cytokinin. These regeneration-related genes include WUSCHEL-RELATED HOMEOBOX (WOX), AUXIN RESPONSE FACTOR (ARF), BABY BOOM (BBM), GROWTH-REGULATING FACTOR(GRF), SCARECROW (SCR), SHORT ROOT (SHR), PLETHORA (PLT), CUP-SHAPED COTYLEDON (CUC), and YUCCA (YUC) expressed during the progress of embryonic patterning, somatic embryogenesis, cell differentiation, wound reparation and epigenetic reprogramming [5,[10][11][12][13][14][15]. An in-depth understanding of regeneration-related genes at the molecular level will make it possible to break through the bottleneck in genetic transformation and build a more efficient transformation system with less genotype dependence. By overexpressing a fusion protein of TaGRF4 and its cofactor GRF-INTERACTING FACTOR 1 (TaGIF1), the wheat regeneration progress was greatly speeded up with an enhanced efficiency [15]. Up-regulation of TaTCP-1 in wheat also resulted in a higher callus regeneration rate [9]. The application of regeneration-associated genes including WUS2 and BBM in crop transformation has achieved a great success, by which various maize (Zea mays) inbred lines and tissues, and recalcitrant genotypes of Indica rice (Oryza sativa ssp indica), sugarcane (Saccharum officinarum) and sorghum (Sorghum bicolor) were efficiently transformed for stable transgenic plants [16,17].
The WOX family is a group of plant-specific transcription factors and belongs to the homeobox (HB) transcription factor family [18]. All the identified WOX genes contain a conserved sequence of amino acids (60-66 residues), which is called a homeodomain (HD) encoded by the HB DNA sequence [19,20]. The distinctive WUS-box motif forms as T-L-X-L-F-P-X-X(T-L-[DEQP]-L-F-P-[GITVL]-[GSKNTCV]), of which the consensus structure is TLELFPLH [18]. These homolog sequences fold into a DNA-binding domain. Published data suggest that WOX genes act as pivotal regulators during the progress of embryonic development and polarization, plant growth and development, stem cell differentiation, embryo patterning, and flower development [21][22][23][24][25]. There are 15 WOX genes in Arabidopsis thaliana, 13 in rice, and 21 in maize [18,26,27]. In Arabidopsis, as a stem cell regulator, AtWUS is expressed in the organizing-center (OC) cells in the shoot apical meristem and regulates plant growth and shoot stem cell maintenance [28,29]. Ectopic overexpression of WUS genes promotes cell dedifferentiation in shoot meristem, somatic embryo formation, adventitious shoot and lateral leaf origination [29][30][31]. WUSCHEL-related homeobox proteins also protect plants from virus infection by repressing the expression of plant S-adenosyl-Lmethionine-dependent methyltransferases, especially protecting plant stem cells against viral intrusion [32].
It is found that AtWOX1 possibly regulates the activity of S-adenosylmethionine decarboxylase polyamine homeostasis and/or the expression of CLAVATA3(CLV3) and has an important function in meristem development in Arabidopsis. Overexpression of AtWOX1 leads to abnormal meristem development and polyamine homeostasis [33]. Normally, AtWOX2 expresses in the zygote and early embryogenesis formation, and performs functions in correcting the apical domain development of the embryos [26]. AtWOX2 triggers the expression of PINFORMED1 (PIN1), which is an auxin transport and localizes auxin to the cotyledonary tips of early embryo and root pole [21]. AtWOX3 (PRESSED FLOWER1, PRS1) expresses in the peripheral layer of shoot meristem and regulates cells to form the lateral domain in vegetative and floral organs [34]. The expression of AtWOX2 and AtWOX3 are regulated by Leafy Cotyledon2 (LEC2), and AtWOX2 and AtWOX3 play essential roles in somatic embryogenesis [35]. AtWOX4 is expressed in a narrow domain in cambial cells, and AtWOX4, coordinating with PHLOEM INTERCALATED WITH XYLEM (PXY), acts as a key regulator for cambium activity in the main stem [36]. AtWOX5 is expressed in the QC of meristematic zone in root tips, regulates the columella stem cell (CSC) identity, and helps to maintain the root stem cell niche [37]. AtWOX6 (PRETTY FEW SEEDS2, PFS2) is expressed in developing ovules and primordial and differentiating organs, regulates ovule development, and affects differentiation and maturation of leaves, outer integuments and floral primordial [38]. AtWOX7 is expressed during all development stages of lateral root but is primarily involved in the initiation of lateral root [39]. AtWOX8 (STIMPY-LIKE) and AtWOX9 (STIMPY) are sister homologs [40,41] and responsible for maintaining the normal development of both basal and apical embryo lineages at the early development stage [21]. The expression of AtWOX8 is induced by AtWRKY2 in the basal cell lineage at the initiation stage of embryogenesis [42]. AtWOX11 plays a key role in the course of vascular cambium differentiation to new lateral root founder cells. AtWOX11 is strongly induced and expressed in de novo root organogenesis, which is the same as its homologous AtWOX12 [43,44]. The expression of AtWOX11/12 is regulated by auxin, which activates the transcription of AtWOX5/7 during the transition from root founder cells to root primordium [45]. AtWOX13 expresses mainly in meristematic tissues to promote replum development and orchestrate fruit patterning [46]. AtWOX14 is regulated by the CLAVATA3/ESRLIKE41/PHLOEM INTERCALATED WITH XYLEM (CLE41/PXY) pair, expressed in the procambium during stem maturation, and promotes xylem differentiation, vascular cell differentiation and lignification in inflorescence stems [47,48].
Based on the phylogenetic analysis in Arabidopsis, plant WOX proteins are naturally divided into three clades: WUS and WOX1 to WOX7 in the WUS clade; WOX8, 9, 11, and 12 in the intermediate clade; and WOX10, 13, and 14 in the ancient clade [18]. However, the WOX genes in Triticeae species have not been fully identified and characterized yet. Therefore, the objectives of this study are (1) identifying WOX genes in nine Triticeae species including T. aestivum, H. vulgare, S. cereale, A. sativa, Thinopyrum elongatum, T. turgidum, T. dicoccoides, Aegilops tauschii, and T. urartu, and aligning them onto chromosomes; (2) dividing all of the WOX proteins in these nine Triticeae species into groups by phylogenetic analysis using deduced protein sequences from all the WOX genes and the sequences of OsWOX genes from rice and AtWOX genes from Arabidopsis; and (3) analyzing the differential expression of TaWOX genes in five different tissues by RNA sequencing (RNAseq) and eight tissues by quantitative real-time PCR (qPCR). Our results provide insights for further understanding the functions and evolution clarification of WOX family genes in Triticeae plants, and facilitate their application in gene transformation for the improvement of Triticeae plants.

Identification of WUS Homoeologous Genes in Triticeae Species
In these nine Triticeae species, only one transcript of WUS gene was annotated as TaWUSa on chromosome 2A in wheat in the database (Table 1). We found the homoeologous fragments of TaWUSa on chromosomes 2B and 2D in T. aestivum (Table 1), 2H in H. vulgare (Table 2), 2A, 2C and 2D in A. sativa (Table 4), 2E in T. elongatum (Table 5), 2D in A. tauschii (Table S2), 2A and 2B in T. dicoccoides and T. turgidum (Tables S3 and S4), 2A in T. urartu (Table S5), and chromosome 2 in S. cereale (Table 3). According to the results of multiple sequence alignment, the full length of the open reading frame (ORF) of these homologous genes can be achieved, and their deduced amino acid sequences were highly consistent with TaWUS ( Figure 1A). To understand if these genes can normally transcribe and express, promoter analysis was performed. It was shown that the promoter region of the WUS genes in six Triticeae plant species all contained core promoter elements, including transcription start TATA-box and AT~TATA-box, indicating they possessed potential transcriptional activity ( Figure 1B). In the promoter region of TaWUSa, TdWUSa, TtWUSa, and TuWUS, a fragment of GGTCCAT existed, which is a cis-acting regulatory element involved in auxin responsiveness. Nevertheless, this element was not detected in the promoter of AtaWUS, TaWUSb, TaWUSd, TdWUSb, and TtWUSb ( Figure 1B).

Chromosomal Location of WOX Genes in Triticeae Species
In general, except ScWOX7, AsWOX10a, AsWOX10.2c, and AsWOX11a genes in S. cereale, no WOX gene was found on homologous groups 6 and 7 in the genomes of these nine Triticeae plant species (Table 1 and Tables S2-S5). In T. aestivum, all the TaWOX genes had three copies in its genomes A, B, and D. Three homologous alleles of TaWUS were located on chromosomes 2A, 2B, and 2D. The homologous genes of TaWOX2 or TaWOX12 were located on chromosomes 1A, 1B, and 1D. Three copies of TaWOX4 or TaWOX11 were located on chromosomes 2A, 2B, and 2D. The three homologous genes of TaWOX7 to TaWOX10, TaWOX13 and TaWOX14 were all located on chromosomes 3A, 3B, and 3D. The three alleles of TaWOX6 were located on chromosomes 4A, 4B, and 4D. The three alleles of TaWOX3 or TaWOX5 were located on chromosomes 5A, 5B, and 5D. Further investigation would be needed for the unknown chromosomal location of an incomplete transcript of TaWOX8. No WOX gene was found on homologous groups 6 and 7 in T. aestivum (Table 1, Figure 2A). The HvWOX genes in H. vulgare showed the similar chromosomal localization to the TaWOX genes in T. aestivum and AtaWOX genes in A. tauschii. HvWOX2 and HvWOX12 were located on chromosome 1H; HvWOX4 and HvWOX11 were located on chromosome 2H; HvWOX7 to HvWOX10, HvWOX13, and HvWOX14 were located on chromosome 3H; HvWOX6 was located on chromosome 4H, and HvWOX3 and HvWOX5 were located on chromosome 5H. (Table 2; Figure 2B). There are additional copies of HvWOX8 and HvWOX10 on chromosome 3H. The HvWOX10.1 and HvWOX10.2 showed complete sequence consistency, but HvWOX8.2 was shortened compared with HvWOX8.1 (Table 2; Figure 2B).
A similar situation was observed in A. tauschii. AtaWOX2 and AtaWOX12 were located on chromosome 1D. AtaWOX4 and AtaWOX11 were located on chromosome 2D. AtaWOX7 to AtaWOX10, AtaWOX13, and AtaWOX14 were all located on chromosome 3D. AtaWOX6 was located on chromosome 4D, AtaWOX3 and AtaWOX5 were located on chromosome 5D (Table S2, Figure S1A). Similar results were also obtained in T. elongatum, T. dicoccoides, T. turgidum, and T. urartu. As expected, all the TeWOX, TdWOX, TtWOX, and TuWOX genes were located on the corresponding chromosomes of their genomes A and B because the two species only have the two genomes (Table 4 and Tables S3-S5, Figure 2D and Figure S1B,C,E). Additional copies of TdWOX8a and TtWOX14a also existed on the corresponding chromosomes.
The chromosomal location of AsWOX genes in A. sativa was complicated. The three copies of AsWOX3, AsWOX4, AsWOX7, and AsWOX12 were distributed on the homoeologous chromosome groups 5, 2, 3, and 1 in order (Table 5, Figure S1D). However, the three copies of AsWOX6 were located on chromosomes 5A, 4C, and 5D; the three copies of AsWOX8 were located on chromosomes 4A, 3C, and 3D; the three copies of AsWOX9 were located on chromosomes 4A, 3C, and 4D; the three copies of AsWOX11 were located on chromosomes 6A, 4C, and 5D. Moreover, only a few fragments of WUS, WOX2, WOX5, WOX10 genes were found in the avaliable genome sequences of A. sativa (Table 5).
It is interestingly found that each WOX or WUS homoeologous gene was collinearly located on the corresponding chromosome among the nine Triticeae species. For example, WOX2 and 12 were located on chromosome group 1 in T. aestivum, H. vulgare, S. cereale, T. elongatum, and A. sativa; WOX4 and WUS were located on chromosome group 2 in the five species; WOX9 was located on chromosome group 3 in the five species; WOX3 was located on chromosome group 5 in the five species (Figure 2 and Figure S1, Tables 1-5 and  Table S1). The results indicated that the WOX or WUS homoeologous genes in Triticeae species were originated via orthologous evolution approach. To verify the chromosomal locations of those WOX genes in these nine Triticeae species, partial sequences of some TaWOX genes were amplified by their specific primers using a set of T. durum-T. aestivum genome D substitution lines (Figure 3). The TaWUSa and its two homologs (named as TaWUSb and TaWUSd) were detected in T. aestivum L. cv CS (ABD genome), T. durum cv Langdon (AB genome), and other substitution lines except 2D(2A), indicating that the two copies TaWUSa and TdWUSa were located on chromosome 2A. TaWUSb was amplified in CS, Langdon, and other substitution lines except 2D(2B), indicating that TaWUSb was located on chromosome 2B. TaWUSd only appeared in CS, 2D(2A) and 2D(2B), indicating that it was located on chromosome 2D (Figure 3). Similarly, WOX2a, WOX2b, WOX6a, and WOX6b were absent in 1D(1A), 1D(1B), 6D(6A), and 6D(6B), respectively. WOX2d and WOX6d were only detected in CS and the substitution lines which contain chromosome 1D or 4D (Figure 3).
The phylogenetic tree of the WOX family proteins from nine Triticeae species was also constructed via maximum likelihood method ( Figure 5). Based on the tree, it was clearly seen that the WOX proteins with the same names from these nine Triticeae species were clustered together ( Figure 5), indicating that the WOX proteins were conserved in these plant species.

Analysis of the Conserved Motifs of WOX Proteins in Triticeae Species
All the amino acid sequences of WOX proteins in nine Triticeae species were deduced from their transcripts mentioned above. Each member contained HOX homeodomain, which were the most noteworthy symbol and defining feature of this protein family ( Figure 6, Figure 7 and Figure S3). Sequences of HOX homeodomain of the three clades of WOX proteins were conserved in these nine Triticeae species ( Figure 7A and Figure S3).

The carboxy-terminal ERFassociated amphiphilic repression (EAR) domain of L-[ED]-L-[RST]-L only exists in WUS
and WOX9 (Figure 6A), and the EAR domain of WOX9 in these Triticeae species was highly conserved ( Figure 7C).

Expression Patterns of TaWOX Genes in Various Wheat Tissues
The WOX genes were mainly expressed in the meristematic region, and played a regulatory role in the process of plant growth and tissue differentiation. We retrieved the data from expVIP website (http://wheat-expression.com, accessed date: 5 September 2020) and sketched the contours of expression pattern of TaWOX genes. It is shown that TaWUS is expressed in the root during the seedling stage, in spike during the vegetative stage, and in spike and leave/shoot during the productive stage. Its expression level was higher in spike than other organs ( Figure S4A). All the three are homologous of TaWOX2 to 4, 7, 8, and TaWOX12 showed a higher expression level in developing spike than other organs, and even higher at the vegetative stage than the reproductive stage (Figure S4B-D,G,H,L). The expression level of TaWOX5 was higher in grain than that in other organs at the reproductive stage ( Figure S4E). TaWOX6, 9 to 11 showed a high transcriptional activity in root ( Figure S4F,I-K). The transcripts of TaWOX10 and TaWOX11 mainly accumulated in root at seedling stage while the expression level of TaWOX9 was high in root at vegetative stage ( Figure S4I-K). The transcript levels of TaWOX6b and TaWOX6d in root were increased at productive stage compared with vegetative stage ( Figure S4F).
Furthermore, we used wheat root, stem, leave, spike at the booting stage, and anther at the heading stage as well as immature embryo and callus derived from the immature embryos at proliferative and differential stages as materials to perform expression profiling analysis of TaWOX genes by qPCR assay. The results indicated that expression patterns of TaWOX genes changed greatly in different organs at different stages ( Figure 8). The expression levels of TaWUS and TaWOX6 to 8 were relative high in spike ( Figure 8A,B), and the expression levels of TaWOX9 and TaWOX11 were high in root ( Figure 8B,C). Additionally, TaWOX2 showed high activity in embryo, and TaWOX3 and TaWOX4 showed high expression levels in embryogenic callus and differential callus, respectively ( Figure 8A).

Expression Patterns of TaWOX Genes during Wheat Callus Proliferation
The expression profiles of 14 TaWOX genes in the immature embryo and the calluses cultured for one, two, and three weeks of wheat were detected by qPCR assay. All the tested TaWOX genes showed a constitutive expression pattern in the calluses. In the early stage of the callus proliferation, the expression of most TaWOX genes was activated and upregulated, and then repressed after two weeks ( Figure 9). The expression level of TaWUS, TaWOX10, TaWOX13, and TaWOX14 in the first week stage was higher than that in other two stages ( Figure 9); the expression level of TaWOX3, TaWOX6, TaWOX8, TaWOX9, and TaWOX12 was increased after the callus initiation and constitutively expressed in the calluses ( Figure 9); the expression level of TaWOX2, TaWOX7, and TaWOX11 was always low during the callus production period ( Figure 9); meanwhile, the expression level of TaWOX4 and TaWOX5 was down-regulated during the callus proliferation course ( Figure 9A). Particularly, TaWOX9 was highly expressed in the calluses of Fielder ( Figure 9B).

Discussion
In Triticeae species, wheat and barley are two important crops globally which account for a large proportion of food production in the world. With the completion of wheat genome assembly and annotation, a great progress on functional genomic study in Triticeae plants, especially in wheat, has been achieved [51][52][53][54]. It is well-known that the wheat genome was originated from the natural hybridization of its three ancestor species. Therefore, the wheat genome consisting of three genomes of A, B, and D has a large number of repeated gene sequences, and most wheat genes have three or more copies [55]. In the present study, we identified 43 putative WOX gene copies in the genome of T. aestivum, 42 of which were consistent with the result reported by Li et al. [56], and a new locus of TaWOX8 was added to the results of TaWOX family. Particularly, we firstly identified 17 putative WOX genes in H. vulgare, 23 putative WOX genes in S. cereale, 24 putative WOX genes in A. sativa, 14 putative WOX genes in T. elongatum, 13 in A. tauschii, 30 in T. turgidum, 25 in T. dicoccoides, and 16 in T. urartu. There were still several duplicated copies of the WOX gene such as TaWOX14a, TaWOX14d, HvWOX10, TdWOX14, ScWOX11, andScWOX13. A few WOX-like pseudo genes were found to be scattered over Triticeae genomes, which might be a duplication of WOX genes or the other genes losing transcriptional activity during their evolution progress.
WUS plays an indispensable role on stem cell niche maintenance in shoot apical meristem (SAM), lateral primordia differentiation and other diverse cellular processes [29]. The deficiency of the WUS gene will lead to the loss of function of SAM and terminated plant growth [28]. However, only the allele of TaWUS located on chromosome 2A was annotated as a transcript. TdWOX12a, TdWOX12b, TdWOX7b, and TdWOX13b, which have a high sequence identity with their homologous genes from wheat, were also not annotated as transcripts in the database. The DNA sequences and deduced protein sequences of four genes TdWOX12a, TdWOX12b, TdWOX7b, and TdWOX13b were added into the WOX members in the six Triticeae species (Table S2). In barley, the annotation of HORVU1Hr1G087940 and HORVU1Hr1G087950 and their deduced protein sequences A0A287GM87 and A0A287GM65 are actually originated from HvWOX12 (Table 2).
In previous studies, the classification and naming of WOX genes in wheat were confusing to some extent. This might be attributed to the different naming scheme of WOX genes in Arabidopsis and rice [18,26,27]. For example, the TaWOX5 reported by Zhao et al. [57] was regarded as TaWOX9 due to its high similarity to OsWOX9, even though it showed a close similarity to AtWOX5 in all the WOX members in Arabidopsis (Figure 4). Several reported TaWOX members such as TraesCS3A02G358100, TraesCS3B02G391100, TraesCS3D02G352500, TraesCS3A02G358200, TraesCS3A02G358400, TraesCS3B02G391200, TraesCS3D02G352600, and TraesCS3D02G352700 on chromosomes 3A, 3B, and 3D, respectively, were named TaWOX13 and TaWOX14 [51] according to new nomination regulations. However, TaWOX13 was not similar to AtWOX13 or OsWOX13, and TaWOX14 was also not similar to AtWOX14 in transcripts, while TaWOX13 and TaWOX14 were similar to the homologs of TaWOX5 according to phylogenic analysis ( Figure 4). The WOX13 and WOX14 in other Triticeae species showed the similar phylogenetic relationship with WOX5 members ( Figure 5).
All the TaWOX genes in wheat have three or more copies. Due to their sequence similarity, it is difficult to distinguish the expression level of each copy of TaWOX genes. A feasible approach was applied to estimate the amount of mRNA by calculating transcript amount of each copy. Zhao et al. indicated that the transcriptional level of the individual TaWOX5 allele was varied during the period of callus growth in wheat [57]. Based on the results in the present investigation, the expression profiles of other WOX alleles were also changed in different wheat organs, which need to be justified by further research.
In the plant regeneration system through somatic embryogenesis from somatic tissues, the activation of an embryonic developmental program is an essential step [58], which is widely adopted in plant genetic engineering. As a group of important transcription factors, WOX family proteins are involved in fate transition of stem cells, cell growth regulation, and somatic embryogenesis process. Overexpression of WOX genes or co-expression of WOX genes and other regeneration-related genes is a potential strategy to resolve the difficult transformation situation in many plant species. It was demonstrated that the overexpression of ZmBBM and ZmWUS2 produced high transformation frequencies in maize, sorghum, sugarcane, and indica rice [16]. A recent study revealed that the overexpression of ZmBBM and ZmWUS2 was driven by the specific promoters of ZmPLTP and ZmAxig1 to embryo development, respectively, initiated abundant somatic embryos, and further developed into plantlet directly without a callus phase [17]. Besides the WUS genes, other WOX gene members were also proved to play a part role in the regeneration process. The AtWOX5 regulated wound healing after laser ablation and root tip excision, in which the molecular program of callus induction and development is similar to those in root establishment, even if the callus is derived from the non-root tissue [59,60]. Following an auxin maximum after wounding in leaf explants, WOX11 and 12 which are involved in the specification of a hypophysis-like root founder cell activate WOX5 and 7 and promote the root primordia initiation [44,45].
Based on the findings achieved in this investigation, TaWUS, TaWOX8, TaWOX10, and TaWOX14 were greatly up-expressed in the differentiated callus ( Figure 8A-C), which indicated that these genes may be related to callus differentiation process. In the callus proliferation process of the immature embryos of four wheat cultivars used in this study, the expression of the 14 TaWOX genes showed a similar pattern (Figure 9). However, in the model spring wheat genotype Fielder with high regeneration and transformation efficiency [7,61], TaWOX9, TaWOX6, TaWOX10, and TaWOX12 all had a significant upregulation during callus production process (Figure 9). Considering that TaWOX9 was a homologous gene of AtWOX5/AtWOX7, and TaWOX6, TaWOX10, and TaWOX12 showed a high homology with AtWOX11/AtWOX12 (Figure 4), we speculate that the expression pattern of these TaWOX genes might be related to the regulation of regeneration process in T. aestivum.

Database Used for Searching WOX Family Genes in Triticeae Plants
Twenty-six predicted WOX family protein sequences of T. aestivum were obtained from Plant TFDB database (http://planttfdb.gao-lab.org/, accessed date: 1 September 2020) and retrieved Genbank (https://www.ncbi.nlm.nih.gov/genbank, accessed date: 1 September 2020) with AtWOX of Arabidopsis and OsWOX of rice (Table S1). Using all of the protein sequences above as queries to conduct a TBLASTN search on Gramene (http://ensembl.gramene.org/Tools/Blast, accessed date: 2 September 2020) and URGI (https://urgi.versailles.inra.fr accessed date: 2 September 2020) for the identification of WOX protein-encoded genes in the genome of T. aestivum. Then, a BLASTN search with sequences of TaWOX genes was performed in the genomes of H. vulgare, T. dicoccoides, T. turgidum, and A. tauschii. Using all the sequences of TaWOX genes as queries to conduct a TBLASTN search on GrainGenes (https://wheat.pw.usda.gov, accessed date: 29 July 2021) for the identification of WOX protein-encoded genes in the genomes of A. sativa and T. urartu, a TBLASTN search on WheatOmics (https://http://202.194.139.32/, accessed date: 29 July 2021) for the identification of WOX protein-encoded genes in the genomes of S. cereale and T. elongatum All the genetic analysis was carried out using these protein sequences of these nine Triticeae plants. Based on the BLAST results from Gramene and URGI, the WOX genes from the Triticeae plants were located on exact chromosomes. The location chart was created by MapGene2Chrom web v2.1 (http://mg2c.iask.in/mg2c_v2.1/ accessed date: 29 July 2021).

Phylogenetic Trees Construction
The full-length of the WOX proteins of Triticeae species were aligned by ClustalW algorithm. Phylogenetic analysis and phylogenetic tree construction were performed by the MEGA X program (version 10.0.5) [62] using the Maximum Likelihood method, JTT matrix-based model [63], and 1000 bootstrap replicates. The initial tree for the heuristic search was obtained automatically by employing Neighbor-Join and BioNJ algorithms to a matrix of pairwise distance estimated using a JTT model, and then selecting the topology with a superior log likelihood value. In total, 130 amino acid sequences were used for this analysis. All the positions with less than 95% site coverage were eliminated, and alignment gaps fewer than 5%, missing data, and ambiguous bases were allowed at any position (partial deletion option). Consequently, 93 positions were remained in the final dataset. Sequences of TaWOX proteins were aligned with AtWOX and OsWOX proteins, and a phylogenetic tree was constructed to confirm classification and phylogenetic relationship of the identified TaWOX members. Then, taking OsWOX proteins as model, the phylogenetic trees were constructed between OsWOX and HvWOX, OsWOX and TdWOX, OsWOX and TtWOX, and OsWOX and AtaWOX members to name and classify the WOX members in the six Triticeae species.

Conserved Protein Motif Analysis
The conserved domain HD was identified by SMART software (http://smart.emblheidelberg.de/, accessed date: were both defined in a strict sense. TEXshade software [64] was employed to perform the multiple sequence alignments for HD domains, WUS-box motifs, and EAR motifs. The logo diagrams were drawn by canonical conserved residues including HD domains, WUS-box motifs, and EAR motifs by SeqLOGO in TBTools (version 1.075) [65].

DNA Isolation and PCR Analysis
Wheat genomic DNA was isolated by NuClean Plant Genomic DNA kit (Cwbio, CW0531M, Beijing, China) from the leaf samples at the three-leaf stage. The PCR reaction system (20 µL) contained 10 µL of 2× Taq Master Mix (containing Mg 2+ and dNTP, Vazyme), 0.5 µL of each forward primer and reverse primer (10 mM), and 1 µL of gDNA (1 µg·µL −1 ), adding ddH 2 O up to 20 • C. Sequences of all the primers used for the detection are shown in Table S6. The thermal cycling conditions were 94 • C for 5 min, 35 cycles of 94 • C for 20 s, 60 • C for 20 s, 72 • C for 30 s, and then 72 • C for 10 min.

Callus Induction from Wheat Immature Embryos
The immature seeds were collected from wheat cultivars Fielder, CB037, CS, and Ningchun4 15 days post anthesis (DPA), and their immature embryos were isolated after surface sterilization as described previously [7,61]. The scutella with the flat side up after removing the embryonic axes were inoculated onto an MS medium containing 2,4-D 2.0 mg L −1 and cultured at 25 • C in darkness for 14 days. Then, the primary calluses were sliced vertically into halves and cultured on the same medium under the same conditions for another 14 days. Then, the embryonic calluses were cultured for differentiation on 1/2 MS medium in a photoperiod of 14 h light and 10 h darkness [7,61]. Callus samples were collected when they were cultured on callus induction medium for one to three weeks, and on calluses differentiation medium for one week.

RNA Isolation and qPCR Analysis
Wheat samples were collected from the seedlings at the three-leaf stage for roots, stems, and leaves, from the adult plants at the booting stage for young spikes, and at the heading stage for anthers. Wheat immature embryo samples were collected at 15 DPA, and callus samples were collected after culturing for one week, followed by two or three weeks on callus production medium, and one week on differentiation medium.

Expression Analysis of TaWOX Genes Using RNA-seq Data
RNA-seq data of 43 TaWOX genes was downloaded from expVIP (http://wheatexpression.com/ accessed date: 5 September 2020). Their expression levels in roots and leaves/shoots at the seedling stage, spikes at the vegetative stage, and grains at the reproductive stage were analyzed.

Statistical Analysis
The SPSS 19.0 software package was employed to statistically analyze the expression data of the target genes achieved by qPCR. Statistical comparisons of multiple sets of data was carried out by Duncan's multiple range test. The histogram was made using the Excel software.

Conclusions
To our knowledge, this is the first study on genome-wide and contrastive analysis on WOX family genes in annual Triticeae plant species. In total, 199 WOX genes were identified, including 43 in T. aestivum, 28 in T. turgidum, 23 in T. dicoccoides, 23 in S. cereale, 24 in A. sativa, 14 in T. elongatum, 15 in H. vulgare, 13 in A. tauschii, and 16 in T. urartu. The homoeologous genes of TaWUSb, TaWUSd, and WUS in the other five Triticeae species were annotated, which were predicted to express normally according to the promoter element analysis. Four novel homologous alleles of TaWOX genes including TdWOX12a, TdWOX12b, TdWOX7b, and TdWOX13b were also identified in T. dicoccoides. All of these WOX members showed similar chromosomal location arrangement and a collinearly orthologous evolution relationship in Triticeae species. Based on the RNA-seq data in the wheat-expression database and qPCR array results, TaWOX genes were found to have a tissue-specific expression feature, and the expression of part-TaWOX genes were closely associated with the regulation of the regeneration process in T. aestivum. The results obtained in this study would be helpful to further understand the molecular function and evolutionary relationship of WOX family genes in Triticeae plants, and potentially apply them in plant genetic transformation in the future.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/ijms22179325/s1, Figure S1. Chromosomal locations of WOX genes in A. tauschii, T. dicoccoides and T. turgidum. Figure S2. Phylogenetic relationships between WOX proteins in rice and barley, T. dicoccoides, T. turgidum, A. tauschii or T. urartu. Figure S3. Alignment of WOX homeodomains from intermediate and ancient clade of 6 Triticeae plant species. Figure S4. Expression profiling of TaWOX genes in different organs at different stages in wheat. Author Contributions: The experiment was conceived by X.Y. and H.L. L.S. and H.L. analyzed the data, K.W. and Y.S. assisted with bioinformatics analysis. L.S. and L.D. performed the PCR and qPCR experiments. The manuscript was drafted by L.S., X.Y. and H.L., and corrected and approved by all authors. All authors have read and agreed to the published version of the manuscript.
Funding: This research was financially supported by the National Natural Science Foundation of China (31971946), and the Chinese Academy of Agricultural Sciences in China (2060302-2-19). The funders had no role in the designing and conducting of this study and collection, analysis, and interpretation of data and in writing the manuscript.
Institutional Review Board Statement: Not applicable.