Next Article in Journal
Detection and Analysis of Genes Affecting the Number of Thoracic Vertebrae in Licha Black Pigs
Previous Article in Journal
Advances in Understanding Fusarium graminearum: Genes Involved in the Regulation of Sexual Development, Pathogenesis, and Deoxynivalenol Biosynthesis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 15
Issue 4
10.3390/genes15040476
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Identification and Expression Analysis of the WOX Transcription Factor Family in Foxtail Millet (Setaria italica L.)

by
Lizhang Nan
1,
Yajun Li
1,
Cui Ma
1,
Xiaowei Meng
1,
Yuanhuai Han
1,
Hongying Li
1,
Mingjing Huang
1,
Yingying Qin
2,* and
Xuemei Ren
2,*
1
College of Agriculture, Shanxi Agricultural University, Taigu, Jinzhong 030800, China
2
College of Life Sciences, Shanxi Agricultural University, Taigu, Jinzhong 030800, China
*
Authors to whom correspondence should be addressed.
Genes 2024, 15(4), 476; https://doi.org/10.3390/genes15040476
Submission received: 8 March 2024 / Revised: 29 March 2024 / Accepted: 7 April 2024 / Published: 10 April 2024
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
WUSCHEL-related homeobox (WOX) transcription factors are unique to plants and play pivotal roles in plant development and stress responses. In this investigation, we acquired protein sequences of foxtail millet WOX gene family members through homologous sequence alignment and a hidden Markov model (HMM) search. Utilizing conserved domain prediction, we identified 13 foxtail millet WOX genes, which were classified into ancient, intermediate, and modern clades. Multiple sequence alignment results revealed that all WOX proteins possess a homeodomain (HD). The SiWOX genes, clustered together in the phylogenetic tree, exhibited analogous protein spatial structures, gene structures, and conserved motifs. The foxtail millet WOX genes are distributed across 7 chromosomes, featuring 3 pairs of tandem repeats: SiWOX1 and SiWOX13, SiWOX4 and SiWOX5, and SiWOX11 and SiWOX12. Collinearity analysis demonstrated that WOX genes in foxtail millet exhibit the highest collinearity with green foxtail, followed by maize. The SiWOX genes primarily harbor two categories of cis-acting regulatory elements: Stress response and plant hormone response. Notably, prominent hormones triggering responses include methyl jasmonate, abscisic acid, gibberellin, auxin, and salicylic acid. Analysis of SiWOX expression patterns and hormone responses unveiled potential functional diversity among different SiWOX genes in foxtail millet. These findings lay a solid foundation for further elucidating the functions and evolution of SiWOX genes.

1. Introduction

Transcription factors encompass a wide range of classifications and functions. They possess the capacity to bind to the promoters of target genes, thereby either activating or repressing their transcriptional activity. These factors orchestrate the specific expression of target genes across diverse tissues, cells, and under various external conditions, governing vital life processes through cascades of transcriptional regulation. The WUSCHEL-related homeobox (WOX) transcription factor family, unique to plants, constitutes a distinct subfamily within the eukaryotic transcription factor homeobox family. Its members exhibit diverse biological functions, primarily contributing to the dynamic equilibrium of plant stem cell division and differentiation, embryonic and post-embryonic development, hormone signal transduction, and responses to crop stress [1,2,3]. Proteins within the WOX family feature a conserved homeodomain (HD) structure composed of 60 amino acid residues, with WUS additionally harboring a tyrosine (Y) residue. These amino acids are arranged into a “helix-loop-helix-turn-helix-loop-helix” structure, where the “helix-loop-helix-turn-helix” formed by the second and third helices can bind to specific DNA sequences, thereby regulating the transcription of downstream genes [4]. Research on the classification of WOX family members primarily focuses on their evolutionary origins. By reconstructing phylogenetic trees of plant WOX using homologous HD domain protein sequences, the phylogenetic tree of plant WOX can be categorized into three clades: Modern, intermediate, and ancient. Bioinformatics analysis of 350 WOX genes in 50 plants revealed that WOX members in lower plants originate solely from the ancient clade, while WOX members in higher plants are exclusively found in the modern clade, with the intermediate clade emerging in vascular plants. This study suggests that the intermediate clade and the evolutionary clade may have evolved from the duplication and differentiation of ancient clade members [5,6,7]. The conserved HD structure of the WOX transcription factor family remains highly conserved across different species. In addition to the HD domain, it also contains a WUS-box domain, exclusive to modern clade members, positioned at the carboxy-terminal end of the homologous domain. The amino acid sequence is in the form of T-L-X-L-F-P-X-X, with X representing any amino acid. Non-modern clade WOX family members exhibit variations at this position. The WUS-box domain plays a critical role in regulating the WUS gene’s control of stem cell characteristics in shoot apical meristems and floral meristem morphology [3,8,9]. Some family members feature carboxy-terminal ethylene-responsive element binding factors (ERF)-related domains, known to be involved in transcriptional repression [10].
In the model plant Arabidopsis thaliana, a total of 16 WOX family members have been identified, comprising AtWUS and AtWOX1 to AtWOX15. Previously documented WOX transcription factor family members are closely associated with regulating stem cell fate, as well as initiating and developing tissues and organs [7]. Among the 16 WOX members, the modern clade encompasses WUS and WOX1 to WOX7; the intermediate clade comprises WOX8, WOX9, WOX11, WOX12, and WOX15; while the ancient clade includes WOX10, WOX13, and WOX14 [7,8]. WUS was the pioneering gene discovered in the WOX family, acting as a pivotal regulator of stem cell regulation in shoot apical meristems. Through EMS mutagenesis, the wus mutant exhibited defective shoot apical meristems, reduced floral organ numbers, and missing pistils [11]. Subsequent investigations unveiled that in mutants with WUS loss of function, WUS plays a crucial role in maintaining stem cell characteristics and quantity, while inhibiting stem cell differentiation [12], Moreover, WUS is implicated in ovule and anther development in Arabidopsis [13,14], with similar roles observed in rice and maize [15]. In Arabidopsis plants overexpressing TaWUS, mutants exhibit early flowering and increased numbers of floral organs, including petals, sepals, and flowers, while the number of stamens and pistils remains unchanged [16]. Genetic and molecular studies have revealed that in the regulatory network of shoot apical meristems in Arabidopsis, WUSCHEL (WUS) interacts with CLAVATA (CLV). Assisted by related transport proteins, WUS relocates from the tissue center to the stem cell area via intercellular connections, stimulating the expression of the CLV3 gene to uphold stem cell characteristics. Once CLV3 binds to its receptor protein CLV1, it dampens WUS transcription levels, facilitating stem cell differentiation and thus forming a CLV-WUS feedback loop to regulate the dynamic balance of stem cells in Arabidopsis shoot apical meristems [17,18,19,20]. This feedback loop can also interact with plant hormones such as cytokinins and auxins to collectively regulate meristem development [21,22]. Furthermore, AtWOX5 exhibits specific expression in the root apical meristem, where it suppresses the differentiation of root apical stem cells, ensuring a constant stem cell count and fostering the growth of root apical meristems [23], CLE40, acting through receptors ACR4 and CLV1, inhibits the expression of AtWOX5, thereby promoting stem cell differentiation, suggesting the presence of a regulatory pathway akin to the shoot apical meristem (SAM) in Arabidopsis’ root apical meristem. During seed development in Arabidopsis [23,24]. AtWOX2 and AtWOX8 are co-expressed, with AtWOX2 playing a crucial role in initiating shoot apical stem cells in embryos [25,26]. AtWOX3 participates in regulating the lateral axis development of flowers, with mutants displaying loss or narrowing of petals [27], AtWOX11 and AtWOX12 collectively influence root meristem development [28], Similarly, AtWOX13 and AtWOX14 play roles in regulating both flower and root development, with AtWOX13 primarily expressed in the stigma and pistil [29]. Recent research has unveiled AtWOX13 as a pivotal regulatory factor in Arabidopsis callus tissue formation. Following explant injury, AtWOX13 is rapidly induced, partially reliant on the activity of the AP2/ERF transcription factor WOUND-INDUCED DEDIFFERENTIATION 1 (WIND1). Subsequently, AtWOX13 directly upregulates WIND2 and WIND3, thereby promoting cell proliferation and organ regeneration [30].
In the growth and development of crops such as rice, maize, and tomato, WOX transcription factors also play pivotal roles. For instance, OsWOX1/MOC3 serves as a homolog of the Arabidopsis WUS gene in rice, whereas MOC1 encodes a plant-specific GAI-RGA-SCR (GRAS) family protein. FON1, on the other hand, acts as a homolog of Arabidopsis CLV1 in rice. MOC1 activates FON1 expression by binding to its promoter, and further investigations have demonstrated that MOC1 can interact with MOC3, acting as a co-activator to enhance FON1 expression [31,32,33,34]. In maize, ZmWUS1 and ZmWUS2 function as homologs of Arabidopsis WUS. It has been observed that the maize mutant Bif3 contains an enhancer region with multiple binding sites for type B response regulators (B-type RR), which significantly upregulates the expression of ZmWUS1-B, resulting in excessive development of meristematic stem cells in the inflorescence [35]. In rice, OsNS, and in maize, ZmNS1/ZmNS2, serve as homologs of AtWOX3. OsNS exhibits predominant expression in rice leaf primordia, young leaves, and flowers, while ZmNS1/ZmNS2 ensure the proper phenotype of maize leaves and floral organs [36,37,38,39]. Overexpression of the tomato gene SlWOX13 significantly influences the expression of genes associated with leaf development, meristem development, and trichome development. Additionally, SlWOX13 is implicated in tomato plant gibberellin synthesis and response, thereby affecting normal plant growth and development [40]. Recent studies have unveiled that SlWOX13 directly activates the expression of multiple genes involved in ethylene synthesis and signaling, as well as carotenoid biosynthesis, thereby positively regulating tomato fruit ripening [41].
With the completion of foxtail millet whole-genome sequencing [42], numerous gene families of foxtail millet have been identified and subjected to bioinformatics analysis. As a model crop of C4 plants [43], foxtail millet has recently been found to possess YUC [44], CLE [45], and GRF [46] families associated with its growth and development. Despite extensive studies on the WOX gene family in various species such as Arabidopsis [8], rice [47], sorghum [47], cotton [48], wheat [16], maize [47,49], mung bean [50], sunflower [51], cucumber [52], melon [53], and Medicago sativa [54], there have been no reports on the WOX gene family in foxtail millet. This study identified 13 foxtail millet WOX genes and constructed a phylogenetic tree comprising foxtail millet (13), Arabidopsis (16), rice (13), maize (21), tomato (10), and green foxtail (13). Through RNA-seq data analysis, the study scrutinized the protein spatial structure, gene structure and conserved motifs, cis-acting elements, chromosome location, tissue-specific expression patterns, and response to different plant hormones of the WOX gene family in foxtail millet. These findings suggest that the WOX gene family plays a crucial role in the growth and development of foxtail millet and its response to various plant hormones, laying the groundwork for further elucidating the functions of WOX genes in foxtail millet.

2. Materials and Methods

2.1. Identification and Phylogenetic Analysis of WOXs in Foxtail Millet

The genomic, proteomic, coding sequence, and GFF annotation data of foxtail millet were downloaded from the foxtail millet multi-omics database MDSI (http://foxtail-millet.biocloud.net/home, accessed on 4 December 2023). Whole-genome data of Arabidopsis, tomato, rice, maize, and green foxtail were obtained from the Phytozome database (https://phytozome-next.jgi.doe.gov/, accessed on 4 December 2023). The HMMER model (PF00046) of the WOX gene family was downloaded from the Pfam database (http://pfam.xfam.org/, accessed on 4 December 2023). TBtools (v2.034) [55] software was used to perform HMM searches in the total protein sequences of foxtail millet, filtering for SiWOX genes. Protein sequences of WOX genes in Arabidopsis were downloaded from the TAIR database (https://www.arabidopsis.org/, accessed on 4 December 2023), and TBtools software was used to perform BLAST to compare Arabidopsis WOX gene family protein sequences with foxtail millet (threshold: 1 × 10−5), identifying SiWOX genes. SiWOX genes obtained from both methods were subjected to domain prediction using the Pfam website, selecting genes containing the HD domain as candidate genes, and finally determining the SiWOX gene family members. Published gene numbers were used to obtain WOX protein sequences of maize, rice, tomato, and Brachypodium distachyon. The MEGA (v11.0) [56] software was used to align WOX protein sequences of Arabidopsis (16), tomato (10), rice (13), maize (21), foxtail millet (13), and green foxtail (13) and construct a phylogenetic tree using the neighbor-joining method (bootstrap = 1000). The resulting tree was beautified using iTOL (https://itol.embl.de/, accessed on 26 December 2023).

2.2. Protein Properties and Sequence Analyses of SiWOX Genes

The physicochemical properties of SiWOX proteins, including the number of amino acids, molecular weight, isoelectric point, instability index, aliphatic index, and hydropathy index, were analyzed using the ExPASy website (https://web.expasy.org/protparam/, accessed on 5 December 2023). Subcellular localization prediction of SiWOX protein sequences was performed using the website (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/#, accessed on 5 December 2023). Protein spatial structure prediction of SiWOX was conducted using the protein structure database website (https://alphafold.ebi.ac.uk/, accessed on 5 December 2023).

2.3. Sequence Alignment, Gene Structure, and Conserved Motif Analysis of SiWOX Gene Family Members

SiWOX protein sequences were aligned using the CLUSTALW website (https://www.ebi.ac.uk/Tools/msa/clustalo/, accessed on 18 December 2023). The alignment results were downloaded and visualized using ESPript (https://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi, accessed on 18 December 2023). The exon and intron structures of SiWOX were analyzed using GSDS2.0 (http://gsds.gao-lab.org/, accessed on 19 December 2023). Conserved motifs of SiWOX genes were predicted using the MEME online website (https://meme-suite.org/meme/, accessed on 25 December 2023), with the number of motifs set to 12, the minimum amino acid residues set to 3, the maximum amino acid residues set to 70, and other parameters set to default.

2.4. Analysis of Promoter Cis-Acting Elements of SiWOX Gene Family Members

The upstream 2000 bp sequences of SiWOX genes were extracted using TBtools software and submitted to the PlantCARE online website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 5 December 2023) for prediction of promoter cis-acting elements.

2.5. Chromosomal Localization and Collinearity Analysis of SiWOX Gene Family Members

Chromosomal position data of SiWOX were obtained from the gene annotation file. Using TBtools’ One Step MCScanX with default parameters, genes undergoing tandem duplication were identified, and Advanced Circos was used for visualization. Similarly, One Step MCScanX was used to identify collinear genes between foxtail millet and green foxtail, rice, and maize. The results were integrated using File Merge For MCScanX, and Multiple Synteny Plot was used to draw inter-species collinearity diagrams.

2.6. Analysis of Expression Patterns of SiWOX Gene Family Members

Transcriptome data of different developmental stages and tissues of Yugu1 were obtained from the Setaria-db database (www.setariadb.com/millet, accessed on 8 January 2023) [57], The tissue naming method followed the study by Meng et al. Heatmaps were generated using TBtools to compare differences in SiWOX expression. In the experimental field of Shanxi Agricultural University’s Minor Crops molecular Breeding Team, Jingu21 (provided by the SXAU Minor Crops Functional Genomics Center) was planted. At the 11–12 leaf stage of foxtail millet, the first period of cone panicles (approximately 1.0–1.5 mm) was sampled, and at the 13–14 leaf stage, the second period of panicles (approximately 2.5–3.0 mm) was sampled. Each sample had 3 biological replicates. The Yugu1 (provided by SXAU Minor Crops Functional Genomics Center) were hydroponically cultivated in an artificial climate chamber at Shanxi Agricultural University (light intensity: 50,000 LX; 16 h of daylight at 28 °C, 8 h of night at 22 °C) with ABA (2 μM) and distilled water as controls. After 9 days, the roots of seedlings were sampled, with 3 biological replicates per sample. After freezing with liquid nitrogen, the panicles and roots were sent to Novogene Bioinformatics Technology Co., Ltd. (Beijing, China) for high-throughput sequencing on the Illumina Hiseq platform. The sequencing results were analyzed using the R language to generate bar charts.

2.7. Response of SiWOX to Plant Hormones

Jingu21 seedlings were hydroponically cultivated in a growth chamber at the College of Life Sciences, Shanxi Agricultural University, using Hoagland’s solution (Beijing Coolaber Technology Co., Ltd., Beijing, China). The light intensity was set to 50,000 LX, with a photoperiod of 16 h of light at 28 °C during the day and 8 h of darkness at 22 °C during the night. After 28 days of foxtail millet germination, seedlings were treated with six hormones: 6-benzylaminopurine (6-BA) (100 µM), abscisic acid (ABA) (100 µM), salicylic acid (SA) (100 µM), gibberellic acid (GA3) (100 µM), jasmonic acid (JA) (100 µM), and indole-3-acetic acid (IAA) (100 µM). Leaf samples were taken at 0 h, 0.5 h, 2 h, 6 h, and 12 h after treatment. Each sample had three biological replicates and was frozen in liquid nitrogen and stored at −80 °C. Primers were designed using the foxtail millet multi-omics database MDSi (http://foxtail-millet.biocloud.net/home, accessed on 4 December 2023) and synthesized by Sangon Biotech Co., Ltd. (Shanghai, China) (Table S1) with specificity checked. The gene Si9g37480 was used as an internal reference [43]. Total RNA was extracted from the leaves using RNA extraction reagent (Beijing Coolaber Technology Co., Ltd., Beijing, China) following the Trizol method. Real-time quantitative polymerase chain reaction (qPCR) was performed using the UnionScript First-strand cDNA Synthesis Mix for qPCR (with dsDNase) for reverse transcription. Before amplification, the cDNA template was diluted fivefold. The qPCR reaction system was prepared using ChamQTM Universal SYBR qPCR Master Mix (Vazyme Biotech Co., Ltd., Nanjing, China). All qPCR reactions were performed on a Bio-Rad CFX96 Touch qPCR instrument, with each reaction repeated three times. The relative expression levels were calculated using the 2−ΔΔCT method [58], and line graphs were generated using R language.

3. Results

3.1. Identification and Phylogenetic Analysis of SiWOX Gene Family

To identify WOX genes in foxtail millet, we conducted BLAST and HMM searches using the 13 WOX genes from Arabidopsis. A total of 13 putative WOX genes were identified in foxtail millet, the same number as in rice and green foxtail [47]. Structural domain prediction using Pfam revealed that all 13 foxtail millet WOX genes contain the HD domain. Based on their chromosomal distribution order, they were named SiWOX1 to SiWOX13. To analyze the evolutionary pattern of the WOX gene family in foxtail millet, we performed cluster analysis and constructed a phylogenetic tree using protein sequences of WOX genes from Arabidopsis (16), tomato (10), rice (13), maize (21), foxtail millet (13), and green foxtail (13) (Figure 1). The 86 WOX proteins from these six species were divided into three evolutionary clades, consistent with previous classification results [3,7,47]. Among them, the modern clade had the most members (46), followed by the intermediate clade (29), and the ancient clade (11). In foxtail millet, SiWOX2, SiWOX3, SiWOX7, SiWOX8, SiWOX10, SiWOX11, and SiWOX12 belonged to the modern clade; SiWOX1, SiWOX4, SiWOX5, SiWOX9, and SiWOX13 belonged to the intermediate clade; only SiWOX6 belonged to the ancient clade. This distribution indicates conserved functional characteristics in foxtail millet growth and development. Further analysis revealed that SiWOX2, along with WUS from maize, Arabidopsis, and tomato, formed a separate cluster, indicating the conservation of WUS protein function in these species. SiWOX2 is speculated to play a role in SAM differentiation during foxtail millet growth and development. Additionally, because foxtail millet, green foxtail, maize, and rice are monocotyledonous plants, while Arabidopsis and tomato are dicotyledonous plants, the genes of foxtail millet always clustered with those of green foxtail, maize, and rice, indicating a closer phylogenetic relationship with them than with Arabidopsis and tomato. Notably, the number of WOX genes in foxtail millet is the same as that in green foxtail, and their relationship in the evolutionary tree is the closest, suggesting that each pair of genes in foxtail millet and green foxtail may be orthologous genes, supporting the hypothesis that foxtail millet was domesticated from green foxtail [59].

3.2. Structure Characterization of the WOX Proteins of Foxtail Millet

Analysis of the characteristics of foxtail millet WOX proteins (Table 1) reveals that the 13 foxtail millet WOX proteins consist of 212 to 531 amino acids, with relative molecular weights ranging from 24.01 to 54.80 kDa. The theoretical isoelectric points range from 6 to 9.26, with most proteins being alkaline, while only 3 are acidic (SiWOX4, SiWOX6, and SiWOX13). The instability coefficients range from 50.26 to 78.12, with SiWOX6 being the only member belonging to the ancient clade and its protein being the most stable. The protein hydropathy coefficients are all negative, indicating that these family members are hydrophilic proteins.Since WOX proteins are transcription factors, subcellular localization results predict that all SiWOX proteins are located in the cell nucleus. Analysis of the predicted protein spatial structures of SiWOX proteins reveals that all SiWOX proteins have similar spatial structures (Figure 2). Among them, SiWOX1 and SiWOX13, SiWOX4 and SiWOX5, and SiWOX11 and SiWOX12 have the most similar protein structures, and they are also clustered into the same subfamily in the evolutionary tree. Detailed information about the amino acids of SiWOX proteins is provided in Text S1.

3.3. Sequence Alignment, Gene Structure, and Conserved Motif Analysis of SiWOX Gene Family Members

Visualization of the multiple sequence alignment results of identified SiWOX protein amino acids (Figure S1) reveals that all SiWOX proteins contain a conserved homologous domain composed of 60 amino acid residues, known as the HD. This number of amino acid residues in the homologous domain is consistent with that found in WOX proteins from rice, Arabidopsis, maize, and sorghum [47]. SiWOX2 has an additional tyrosine (Y) residue in its homologous domain compared to other SiWOX proteins. Previous studies have shown that the additional tyrosine residues in the homologous domains of WUS proteins from Arabidopsis, maize, rice, and sorghum are highly conserved [47], further indicating that SiWOX2 likely functions similarly to the WUS gene in foxtail millet, Arabidopsis, maize, rice, and sorghum, primarily maintaining the number and characteristics of stem cell populations during foxtail millet SAM differentiation [12]. This also suggests that this residue may play an important role in the functionality of WUS transcription factors. Within the homologous domains of SiWOX, there are 17 highly conserved amino acid residues, consistent with previous studies [47]. All modern clade members of SiWOX proteins contain a WUS-box domain, with the amino acid sequence present in the form of T-L-X-L-F-P-X-X (where X represents any amino acid). Previous studies have shown that this domain is unique to modern clade members [3]. In Arabidopsis WOX transcription factors, WUS, WOX5, and WOX7 contain a conserved ethylene-responsive element binding factor-associated amphiphilic repression (EAR)-like motif at their C-terminus, the main function of which is to suppress downstream gene expression [3], SiWOX2 and SiWOX8 also contain an EAR-like motif at their carboxyl terminus, with the amino acid sequence L-E-L-X-L-X (where X represents any amino acid). In the phylogenetic tree, SiWOX2 clusters with AtWUS, while SiWOX8 clusters with AtWOX5 and AtWOX7, suggesting that SiWOX2 and SiWOX8 in foxtail millet may have similar functions to AtWUS, AtWOX5, and AtWOX7, which are homologous to Arabidopsis.
In the phylogenetic evolution tree of foxtail millet WOX, SiWOX genes with similar gene structures and identical conserved motifs often cluster into one subfamily (Figure 3A). Gene structure analysis results (Figure 3B) show that SiWOX2 and SiWOX10 each contain only one intron, while the remaining SiWOX genes contain one intron in both the N-terminal and C-terminal regions. SiWOX3, SiWOX7, SiWOX8, SiWOX11, and SiWOX12 contain two exons, while SiWOX1, SiWOX2, SiWOX5, SiWOX6, SiWOX9, SiWOX10, and SiWOX13 contain three exons, and only SiWOX4 contains four exons. Conserved motif prediction results (Figure 3C) identified 12 conserved motifs among the 13 SiWOX genes, with SiWOX1 and SiWOX13 containing the most conserved motifs, with a total of 7. All SiWOX proteins contain motifs 1 and 2, indicating that motifs 1 and 2 are conserved in SiWOX proteins, with the ancient clade SiWOX6 being the most conserved protein, containing only motifs 1 and 2. Motif 1 and motif 2 correspond to the HD domain of SiWOX. Motif 4 includes motif 11, with motif 11 corresponding to the WUS-box domain. Motif 11 is only present in SiWOX2, SiWOX3, SiWOX7, SiWOX8, and SiWOX10, while motif 4 is only present in SiWOX11 and SiWOX12. The WUS-box is only present in the modern clade, and these proteins belong to the modern clade in the evolutionary tree. Intermediate clade SiWOX proteins all contain motif 3, with SiWOX4 and SiWOX5 also containing motif 10, and SiWOX1, SiWOX9, and SiWOX13 all containing motifs 5, 7, and 9. In addition, SiWOX11 and SiWOX12 uniquely contain motif 8.

3.4. Promoter Cis-Acting Element Analysis of SiWOX Gene Family Members

To further elucidate the potential functions of the SiWOX gene family in foxtail millet growth, development, and response to abiotic stresses, upstream 2000 bp sequences of WOX genes in foxtail millet were extracted for analysis using the PlantCare database (Table S3). The results indicated an uneven distribution of various functional cis-acting elements on each SiWOX gene, primarily including two major categories: Stress response and plant hormone response cis-acting regulatory elements. SiWOX9 contained the highest number of cis-acting elements, with 114, while SiWOX2 only contained 30 cis-acting elements. All SiWOX promoters contained light-responsive elements, except for SiWOX13, suggesting that the remaining SiWOX genes may be associated with environmental stresses, as they contained cis-acting elements related to low temperature, drought, defense, and stress. Cis-acting elements related to hormone responses in SiWOX included methyl jasmonate (50), abscisic acid (36), gibberellin (13), auxin (6), and salicylic acid (3). The ancient clade SiWOX6 had all types of cis-acting elements, while SiWOX2 had only one salicylic acid-responsive cis-acting element, with the other two distributed in SiWOX6 and SiWOX7. SiWOX1, SiWOX6, SiWOX7, SiWOX9, SiWOX10, and SiWOX11 each had one auxin-responsive cis-acting element, suggesting that these genes are involved in auxin regulatory feedback during foxtail millet growth and development. Additionally, SiWOX1, SiWOX2, SiWOX4, SiWOX5, SiWOX9, SiWOX10, SiWOX12, and SiWOX13 all contained cis-acting elements involved in seed-specific regulation, with SiWOX2 and SiWOX4 also containing cis-acting elements involved in endosperm-specific expression. These results suggest that SiWOX genes may be associated with foxtail millet stress resistance and hormone response pathways.

3.5. Chromosomal Localization and Collinearity Analysis of SiWOX Gene Family Members

To further analyze the evolutionary relationships of foxtail millet WOX family members, the identified foxtail millet WOX genes were mapped onto chromosomes based on gene annotation information, and the results were visualized using Advanced Circos (Figure 4). Except for chromosome1 and chromosome4, WOX genes were distributed on all other chromosomes, with chromosome5 having the most with four WOX genes, and chromosome2, chromosome6, and chromosome9 each having one WOX gene. Gene tandem duplication is an important internal driving force for crop domestication and breeding, playing a significant role in plant phenotypic variation [60,61,62]. Using MCScanX, the tandem repeat relationship of WOX genes in foxtail millet was analyzed, revealing that there are three pairs of tandem repeats in foxtail millet WOX genes, namely SiWOX1 and SiWOX13, SiWOX4 and SiWOX5, and SiWOX11 and SiWOX12. Tandem repeats of SiWOX11 and SiWOX12 occurred in the modern clade, while the other two pairs occurred in the middle clade, with the ancient clade being the most conservative and having no gene tandem repeats. Gene tandem repeats can lead to changes in gene structure [63], such as SiWOX4 having one more intron than SiWOX5, indicating that gene tandem duplication enriches the diversity of the SiWOX gene family.
Using the multiple synteny plot, a collinearity graph between green foxtail, foxtail millet, maize, and rice was generated (Figure 5). From the graph, it can be observed that the collinearity of WOX genes between green foxtail and foxtail millet is the highest, further suggesting that they may share a common ancestor. The collinearity of WOX genes between foxtail millet and maize is higher than that between foxtail millet and rice, possibly because foxtail millet and maize both belong to C4 crops, while rice is a C3 crop.

3.6. Analysis of the Expression Patterns of SiWOX Genes

Transcriptome data from various tissues at different growth and development stages of Yugu1 were obtained from the Setaria-db database, and a heatmap of gene expression was generated (Figure 6). The ancient clade member SiWOX6 exhibited the most conservative function, with expression observed across different tissues at various stages of foxtail millet development. Among the intermediate clade members, SiWOX1 was expressed in all other tissues at various stages except for anthers, mainly during the vegetative growth and seed maturation stages of foxtail millet, suggesting that its evolutionary rate is slower compared to other intermediate clade genes. Genes clustered in the same clade of the phylogenetic tree often share similar functions, but differences exist. For instance, SiWOX11 and SiWOX12 were primarily expressed in the later stages of stem apex meristem development, with SiWOX12 also showing expression during the spikelet primordia and flowering stages. Additionally, expression of SiWOX4 decreased initially and then increased with seed maturation, while SiWOX5 exhibited the opposite trend, indicating distinct functions during foxtail millet seed maturation. SiWOX1, SiWOX6, SiWOX8, and SiWOX13 were highly expressed in germinating seeds, suggesting their involvement in seed germination. SiWOX1 and SiWOX13 showed high expression levels in various parts of foxtail millet such as roots, stems, and leaves, indicating their association with cell proliferation, division, differentiation, and tissue elongation. SiWOX10 exhibited high expression levels in the stems and nodes of foxtail millet and contained cis-acting elements responsive to auxin, suggesting its involvement in cytokinin and auxin synthesis. SiWOX1 and SiWOX6 were highly expressed in primary branching during flowering, indicating their joint influence on inflorescence development in foxtail millet. As foxtail millet seeds matured, expression of SiWOX9 gradually decreased, indicating its involvement in early seed development. Only SiWOX2 exhibited tissue specificity, with high expression only during the development of stem apical meristems, suggesting its role in controlling panicle development in foxtail millet. These results suggest that SiWOXs may play important roles in the vegetative and reproductive growth processes of foxtail millet.
The development of the foxtail millet panicle involves a continuous differentiation process, starting from the SAM forming a “smooth round growth cone”, which elongates and develops small protrusions to form primary branches. As primary branches grow, they initiate the differentiation of secondary branches, and after tertiary branch differentiation, spikelets and florets begin to differentiate at the tip of tertiary branches, completing panicle development [64]. Analysis of transcriptome sequencing results of foxtail millet panicle development (Figure 7A) revealed higher expression levels of SiWOX2, SiWOX3, SiWOX4, SiWOX5, and SiWOX6 compared to SiWOX1, SiWOX10, SiWOX11, and SiWOX12. Among them, the ancient clade member SiWOX6 exhibited the highest expression level, followed by SiWOX3 and SiWOX5, while SiWOX7, SiWOX8, SiWOX9, and SiWOX13 showed almost negligible expression levels. During the second stage of panicle development, the expression levels of SiWOX3 and SiWOX10 significantly increased, while SiWOX5 decreased, and the expression levels of other genes remained almost unchanged. In foxtail millet seedling roots, fewer SiWOX genes were expressed compared to panicle development, but their expression levels were relatively high (Figure 7B). Among them, SiWOX6 exhibited the highest expression level, followed by SiWOX13, SiWOX8, SiWOX1, and SiWOX9. ABA plays a crucial role in regulating plant root growth and development [65]. Compared to the control, the expression levels of SiWOX6 and SiWOX13 were significantly upregulated in seedling roots treated with ABA, while SiWOX1 and SiWOX9 also showed upregulation, and SiWOX8 was downregulated. Overall, many members of the SiWOX gene family are specifically involved in the development of panicle and root in foxtail millet.

3.7. SiWOX Gene Response to Plant Hormones

To investigate the response of SiWOX genes to various hormones, 28-day-old foxtail millet seedlings were treated with six plant hormones: 6-BA (100 μM), ABA (100 μM), GA3 (100 μM), SA (100 μM), MeJA (100 μM), and IAA (100 μM). The expression patterns of SiWOX1, SiWOX2, SiWOX6, SiWOX7, SiWOX10, and SiWOX13 genes, distributed across different subfamilies in foxtail millet seedling leaves, were analyzed (Figure 8). Results showed that after treatment with MeJA, all six SiWOX genes exhibited an up-down-up-down expression pattern, with expression peaking at 0.5 h after treatment. By 12 h after treatment, except for SiWOX7, the expression levels of the other five SiWOX genes decreased to levels observed at 2 h after treatment. Except for SiWOX7 and SiWOX10, the expression levels of the remaining four SiWOX genes remained nearly unchanged in response to GA treatment. Compared to the other four hormones, the ancient clade member SiWOX6 showed minimal response to ABA and GA. SiWOX1 demonstrates analogous trends in response to JA and SA. The relative expression level at 0.5 h post-hormone treatment surpasses that at 12 h by more than double. Conversely, under IAA treatment, the relative expression level at 0.5 h is approximately half of that observed at 12 h. Additionally, SiWOX10 and SiWOX13 exhibited similar trends in response to the six hormones. These results suggest that the SiWOX gene family may have different potential functions during foxtail millet leaf development.

4. Discussion

Foxtail millet holds significance as an important cereal crop due to its resilience to drought, tolerance to poor soil, and salinity, characterized by well-developed roots [66,67,68]. WOX transcription factors, unique to plants, play a pivotal and conserved role in plant growth, development, response to various abiotic stresses, and plant hormone signaling [69,70,71]. While extensively studied in plants such as Arabidopsis, maize, rice, and tomato [8,47,72], the WOX gene family in foxtail millet has not been previously reported. This study employed two methods and conservative domain prediction to identify 13 WOX genes in foxtail millet (Table 1), consistent with the number in green foxtail and rice, as evidenced by phylogenetic analysis (Figure 1), underscoring the highly conserved distribution of unique clades (ancient, intermediate, and modern) within this family in foxtail millet. Previous studies on WOX family members have primarily focused on evolution. The foxtail millet WOX gene family exhibits only one member in the ancient clade, with the maximum number of members in the modern clade being seven, akin to findings in other species [47]. Unlike the AtWUS protein sequence, which features a non-conservative acidic amino acid sequence between the HD homologous domain and the WUS-box domain [47], the multiple sequence alignment results (Figure S1) indicate the absence of such an amino acid sequence in the SiWOX protein sequence, suggesting that the acidic region domain may predominantly influence Arabidopsis WOX transcription factors.
In the phylogenetic tree (Figure 1), SiWOX2 clusters with WUS in maize, Arabidopsis, and tomato, forming a distinct cluster. Through multiple sequence alignment of the foxtail millet WOX gene family (Figure S1), we identified additional tyrosine (Y) residues in the homologous domain of SiWOX2, which are highly conserved with the additional tyrosine (Y) residues in the WUS domain of Arabidopsis, rice, maize, and tomato [47]. The main function of AtWUS is to maintain the number and characteristics of SAM stem cells [12]. In the gene expression pattern (Figure 6), SiWOX2 shows specific expression in the SAM, suggesting a similar function in regulating the size of the SAM in foxtail millet. In maize BIF3 mutants, overexpression of ZmWUS1 leads to spherical spikes and a significant reduction in yield [73]. Controlling the expression of meristem genes can alter plant yield without affecting the growth of explants [74], highlighting the potential importance of SiWOX2 in improving foxtail millet yield traits. Furthermore, collinearity analysis (Figure 5) reveals a high degree of collinearity between the WOX family members of foxtail millet and green foxtail. The differences between the inflorescences of foxtail millet and green foxtail, such as changes in the number, density, and order of primary branches [75], may be related to differences in the WUS gene between the two species, warranting further investigation.
In Arabidopsis, rice, maize, and tomato, different members of the WOX gene family have distinct functions, and phylogenetic tree studies provide references for predicting the functions of WOX genes in foxtail millet. The WOX gene family plays a crucial role in regulating plant growth and development processes [76,77,78]. The CLV-WUS feedback pathway ensures the self-renewal of stem cells while forming new plant tissues and organs [79]. Recently, the foxtail millet CLE gene family has also been identified and reported [45]. In the evolutionary tree (Figure 1), SiWOX2, together with AtWUS, ZmWUS1, ZmWUS2, SvWUS, and SlWUS, forms a subcluster and is specifically expressed during foxtail millet inflorescence meristem development, indicating the possible existence of a CLV-WUS feedback pathway in the SAM of foxtail millet, with SiWOX2 playing a crucial role in the signal transduction of the SAM CLV-WUS feedback. AtWOX5 is specifically expressed in the root apical meristem, and similar regulatory pathways to those in the SAM also exist in Arabidopsis root apical meristem [23,24]. In the evolutionary tree, SiWOX8 clusters with AtWOX5 in the same subcluster; in the gene expression pattern (Figure 6), SiWOX8 shows higher expression levels in the roots of foxtail millet seedlings and at the one-leaf-two-heart stage, indicating that SiWOX8 may play an important role in the development of the root apical meristem in foxtail millet. OsWOX11 is a key regulator of crown root development in rice, and its expression is induced by auxin and cytokinin. Recently, a dynamic feedback regulation mechanism formed by WOX11-JMJ706-LBD16 was discovered in the growth process of rice crown roots, ensuring that LBD16 maintains an appropriate expression level to ensure the growth of rice crown roots [80]. In the evolutionary tree (Figure 1), SiWOX1 and OsWOX11 cluster together, and combined with the analysis of SiWOX gene expression patterns in different tissues (Figure 6) and the expression pattern of SiWOX genes in the roots of foxtail millet (Figure 7B), we found that SiWOX1 is highly expressed in the roots of foxtail millet at different stages, indicating its important role in the root development of foxtail millet.

5. Conclusions

This study provides a comprehensive analysis of the foxtail millet SiWOX gene family, with 13 SiWOX genes divided into three subclusters. All SiWOX proteins contain a conserved HD domain, with SiWOX2 featuring an additional tyrosine (Y) residue in this domain. SiWOX genes clustered in the same clade of the phylogenetic tree exhibit similar protein spatial structures, gene structures, and conserved motifs. Tandem repeats analysis reveals three pairs of repeats among the SiWOX genes, while collinearity analysis shows the highest collinearity between foxtail millet and green foxtail, suggesting a common ancestor. Analysis of SiWOX expression patterns and hormone responses reveals potential functional diversity among different SiWOX genes in foxtail millet. Overall, this study provides a foundation for further elucidating the functions of SiWOX genes in foxtail millet.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15040476/s1.

Author Contributions

H.L., Y.Q. and X.R. designed the study; L.N., Y.L. and C.M. wrote the manuscript; L.N., X.M. and M.H. performed the bioinformatics analysis; L.N. and Y.L. performed the experiment; Y.H., X.R. and Y.Q. revised the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) (32001515, to Ren Xuemei) and (32201721, to Qin Yingying), Basic Research Program of Shanxi Province (20210302124503, to Qin Yingying); Scientific and Technological Innovation Programs of Higher Education (2021L114, to Qin Yingying); Innovation Project of Shanxi Foxtail Millet Industry (202104010911013-1, to Han Yuanhuai); Science and Technology Innovation Fund Project of Shanxi agricultural university (2020BQ58, to Qin Yingying) and (2018YJ32, to Ren Xuemei).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhou, Y.; Liu, X.; Engstrom, E.M.; Nimchuk, Z.L.; Pruneda-Paz, J.L.; Tarr, P.T.; Sun, Y.; Zhang, Y.; Oh, E.M.; Jacobsen, S.E. Meyerowitz, Control of plant stem cell function by conserved interacting transcriptional regulators. Nature 2015, 517, 377–380. [Google Scholar] [CrossRef]
  2. Li, X.; Hamyat, M.; Liu, C.; Ahmad, S.; Gao, X.; Guo, C.; Li, R.; Wang, Y.; Abdullah, M.; Gao, X. Identification and Characterization of the WOX Family Genes in Five Solanaceae Species Reveal Their Conserved Roles in Peptide Signaling. Genes 2018, 9, 260. [Google Scholar] [CrossRef]
  3. Van der Graaff, E.; Laux, T.; Rensing, S.A. Rensing, The WUS homeobox-containing (WOX) protein family. Genome Biol. 2009, 10, 248. [Google Scholar] [CrossRef] [PubMed]
  4. Gehring, W.J.; Affolter, M.; Bürglin, T. Homeodomain proteins. Annu. Rev. Biochem. 1994, 63, 487–526. [Google Scholar] [CrossRef] [PubMed]
  5. Richardt, S.; Lang, D.; Reski, R.; Frank, W.; Rensing, S.A. Rensing, PlanTAPDB, a Phylogeny-Based Resource of Plant Transcription-Associated Proteins. Plant Physiol. 2007, 143, 1452–1466. [Google Scholar] [CrossRef]
  6. Palovaara, J.; Hakman, I. Conifer WOX-related homeodomain transcription factors, developmental consideration and expression dynamic of WOX2 during Picea abies somatic embryogenesis. Plant Mol. Biol. 2008, 66, 533–549. [Google Scholar] [CrossRef] [PubMed]
  7. Lian, G.; Ding, Z.; Wang, Q.; Zhang, D.; Xu, J. Origins and Evolution of WUSCHEL-Related Homeobox Protein Family in Plant Kingdom. Sci. World J. 2014, 2014, 534140. [Google Scholar] [CrossRef] [PubMed]
  8. Haecker, A.; Groß-Hardt, R.; Geiges, B.; Sarkar, A.; Breuninger, H.; Herrmann, M.; Laux, T. Expression dynamics of WOX genes mark cell fate decisions during early embryonic patterning in Arab. Thaliana Dev. 2004, 131, 657–668. [Google Scholar] [CrossRef]
  9. Ikeda, M.; Mitsuda, N.; Ohme-Takagi, M. Arabidopsis WUSCHEL Is a Bifunctional Transcription Factor That Acts as a Repressor in Stem Cell Regulation and as an Activator in Floral Patterning. Plant Cell 2009, 21, 3493–3505. [Google Scholar] [CrossRef]
  10. Paponov, I.A.; Teale, W.; Lang, D.; Paponov, M.; Reski, R.; Rensing, S.A.; Palme, K. The evolution of nuclear auxin signalling. BMC Evol. Biol. 2009, 9, 126. [Google Scholar] [CrossRef]
  11. Laux, T.; Mayer, K.F.X.; Berger, J.; Jürgens, G. The WUSCHEL gene is required for shoot and floral meristem integrity in Arabidopsis. Development 1996, 122, 87–96. [Google Scholar] [CrossRef]
  12. Mayer, K.F.X.; Schoof, H.; Haecker, A.; Lenhard, M.; Jürgens, G.; Laux, T. Role of WUSCHEL in Regulating Stem Cell Fate in the Arabidopsis Shoot Meristem. Cell 1998, 95, 805–815. [Google Scholar] [CrossRef] [PubMed]
  13. Gross-Hardt, R.; Lenhard, M.; Laux, T. WUSCHEL signaling functions in interregional communication during Arabidopsis ovule development. Genes Dev. 2002, 16, 1129–1138. [Google Scholar] [CrossRef] [PubMed]
  14. Deyhle, F.; Sarkar, A.K.; Tucker, E.J.; Laux, T. WUSCHEL regulates cell differentiation during anther development. Dev. Biol. 2007, 302, 154–159. [Google Scholar] [CrossRef]
  15. Nardmann, J.; Werr, W. The Shoot Stem Cell Niche in Angiosperms: Expression Patterns of WUS Orthologues in Rice and Maize Imply Major Modifications in the Course of Mono- and Dicot Evolution. Mol. Biol. Evol. 2006, 23, 2492–2504. [Google Scholar] [CrossRef] [PubMed]
  16. Li, Z.; Liu, D.; Xia, Y.; Li, Z.; Jing, D.; Du, J.; Niu, N.; Ma, S.; Wang, J.; Song, Y.; et al. Identification of the WUSCHEL-Related Homeobox (WOX) Gene Family, and Interaction and Functional Analysis of TaWOX9 and TaWUS in Wheat. Int. J. Mol. Sci. 2020, 21, 1581. [Google Scholar] [CrossRef] [PubMed]
  17. Clark, S.E. Cell signalling at the shoot meristem. Nat. Rev. Mol. Cell Biol. 2001, 2, 276–284. [Google Scholar] [CrossRef]
  18. Schoof, H.; Lenhard, M.; Haecker, A.; Mayer, K.F.X.; Jürgens, G.; Laux, T. The Stem Cell Population of Arabidopsis Shoot Meristems Is Maintained by a Regulatory Loop between the CLAVATA and WUSCHEL Genes. Cell 2000, 100, 635–644. [Google Scholar] [CrossRef] [PubMed]
  19. Su, Y.H.; Zhou, C.; Li, Y.J.; Yu, Y.; Tang, L.P.; Zhang, W.J.; Yao, W.J.; Huang, R.; Thomas, L.; Zhang, X.S. Integration of pluripotency pathways regulates stem cell maintenance in the Arabidopsis shoot meristem. Proc. Natl. Acad. Sci. USA 2020, 117, 22561–22571. [Google Scholar] [CrossRef]
  20. Daum, G.; Medzihradszky, A.; Suzaki, T.; Lohmann, J.U. A mechanistic framework for noncell autonomous stem cell induction in Arabidopsis. Proc. Natl. Acad. Sci. USA 2014, 111, 14619–14624. [Google Scholar] [CrossRef]
  21. Lopes, F.L.; Galvan-Ampudia, C.; Landrein, B. WUSCHEL in the shoot apical meristem: Old player, new tricks. J. Exp. Bot. 2021, 72, 1527–1535. [Google Scholar] [CrossRef] [PubMed]
  22. Leibfried, A.; To, J.P.; Busch, W.; Stehling, S.; Kehle, A.; Demar, M.; Kieber, J.J.; Lohmann, J.U. WUSCHEL controls meristem function by direct regulation of cytokinin-inducible response regulators. Nature 2005, 438, 1172–1175. [Google Scholar] [CrossRef] [PubMed]
  23. Sarkar, A.K.; Luijten, M.; Miyashima, S.; Lenhard, M.; Hashimoto, T.; Nakajima, K.; Scheres, B.; Heidstra, R.; Laux, T. Conserved factors regulate signalling in Arabidopsis thaliana shoot and root stem cell organizers. Nature 2007, 446, 811–814. [Google Scholar] [CrossRef] [PubMed]
  24. Berckmans, B.; Kirschner, G.; Gerlitz, N.; Stadler, R.; Simon, R. CLE40 Signaling Regulates Root Stem Cell Fate. Plant Physiol. 2019, 182, 1776–1792. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, Z.; Tucker, E.; Hermann, M.; Laux, T. A Molecular Framework for the Embryonic Initiation of Shoot Meristem Stem Cells. Dev. Cell 2017, 40, 264–277.e4. [Google Scholar] [CrossRef] [PubMed]
  26. Verma, S.; Attuluri, V.P.S.; Robert, H.S. Transcriptional control of Arabidopsis seed development. Planta 2022, 255, 90. [Google Scholar] [CrossRef] [PubMed]
  27. Matsumoto, N.; Okada, K. A homeobox gene, PRESSED FLOWER, regulates lateral axis-dependent development of Arabidopsis flowers. Genes Dev. 2001, 15, 3355–3364. [Google Scholar] [CrossRef]
  28. Liu, J.; Sheng, L.; Xu, Y.; Li, J.; Yang, Z.; Huang, H.; Xu, L. WOX11 and 12 are involved in the first-step cell fate transition during de novo root organogenesis in Arabidopsis. Plant Cell 2014, 26, 1081–1093. [Google Scholar] [CrossRef] [PubMed]
  29. Deveaux, Y.; Toffano-Nioche, C.; Claisse, G.; Thareau, V.; Morin, H.; Laufs, P.; Moreau, H.; Kreis, M.; Lecharny, A. Genes of the most conserved WOX clade in plants affect root and flower development in Arabidopsis. BMC Evol. Biol. 2008, 8, 291. [Google Scholar] [CrossRef]
  30. Ikeuchi, M.; Iwase, A.; Ito, T.; Tanaka, H.; Favero, D.S.; Kawamura, A.; Sakamoto, S.; Wakazaki, M.; Tameshige, T.; Fujii, H.; et al. Wound-inducible WUSCHEL-RELATED HOMEOBOX 13 is required for callus growth and organ reconnection. Plant Physiol. 2022, 188, 425–441. [Google Scholar] [CrossRef]
  31. Liao, Z.; Yu, H.; Duan, J.; Yuan, K.; Yu, C.; Meng, X.; Kou, L.; Chen, M.; Jing, Y.; Liu, G.; et al. SLR1 inhibits MOC1 degradation to coordinate tiller number and plant height in rice. Nat. Commun. 2019, 10, 2738. [Google Scholar] [CrossRef] [PubMed]
  32. Li, X.; Qian, Q.; Fu, Z.; Wang, Y.; Xiong, G.; Zeng, D.; Wang, X.; Liu, X.; Teng, S.; Fujimoto, H.; et al. Control of tillering in rice. Nature 2003, 422, 618–621. [Google Scholar] [CrossRef] [PubMed]
  33. Lu, Z.; Shao, G.; Xiong, J.; Jiao, Y.; Wang, J.; Liu, G.; Meng, X.; Liang, Y.; Xiong, G.; Wang, Y.; et al. MONOCULM 3, an Ortholog of WUSCHEL in Rice, Is Required for Tiller Bud Formation. J. Genet. Genom. 2015, 42, 71–78. [Google Scholar] [CrossRef]
  34. Takuya, S.; Makoto, S.; Motoyuki, A.; Masahiro, M.; Yasuo, N.; Hiro-Yuki, H. The gene FLORAL ORGAN NUMBER1 regulates floral meristem size in rice and encodes a leucine-rich repeat receptor kinase orthologous to Arabidopsis CLAVATA1. Development 2004, 131, 5649–5657. [Google Scholar]
  35. Chen, Z.; Cortes, L.; Gallavotti, A. Genetic dissection of cis-regulatory control of ZmWUSCHEL1 expression by type-B RESPONSE REGULATORS. Plant Physiol. 2023, 194, 2240–2248. [Google Scholar] [CrossRef] [PubMed]
  36. Dai, M.; Hu, Y.; Zhao, Y.; Liu, H.; Zhou, D.X. A WUSCHEL-LIKE HOMEOBOX gene represses a YABBY gene expression required for rice leaf development. Plant Physiol. 2007, 144, 380–390. [Google Scholar] [CrossRef] [PubMed]
  37. Li, J.; Yuan, Y.; Lu, Z.; Yang, L.; Gao, R.; Lu, J.; Li, J.; Xiong, G. Glabrous Rice 1, encoding a homeodomain protein, regulates trichome development in rice. Rice 2012, 5, 32. [Google Scholar] [CrossRef] [PubMed]
  38. Zhang, H.; Wu, K.; Wang, Y.; Peng, Y.; Hu, F.; Wen, L.; Han, B.; Qian, Q.; Teng, S. A WUSCHEL-like homeobox gene, OsWOX3B responses to NUDA/GL-1 locus in rice. Rice 2012, 5, 30. [Google Scholar] [CrossRef] [PubMed]
  39. Scanlon, M.J.; Freeling, M. The narrow sheath leaf domain deletion: A genetic tool used to reveal developmental homologies among modified maize organs. Plant J. 1998, 13, 547–561. [Google Scholar]
  40. Chen, Y. Study on Transcription Factor SlWOX13 Regulating Tomato Growth and Development through Gibberellin Pathway. Master’s Thesis, Chongqing University, Chongqing, China, 2022. [Google Scholar]
  41. Jiang, G.; Li, Z.; Ding, X.; Zhou, Y.; Lai, H.; Jiang, Y.; Duan, X. WUSCHEL-related homeobox transcription factor SlWOX13 regulates tomato fruit ripening. Plant Physiol. 2023, 194, 2322–2337. [Google Scholar] [CrossRef]
  42. Zhang, G.; Liu, X.; Quan, Z.; Cheng, S.; Xu, X.; Pan, S.; Xie, M.; Zeng, P.; Yue, Z.; Wang, W.; et al. Genome sequence of foxtail millet (Setaria italica) provides insights into grass evolution and biofuel potential. Nat. Biotechnol. 2012, 30, 549–554. [Google Scholar] [CrossRef]
  43. Yang, Z.; Zhang, H.; Li, X.; Shen, H.; Gao, J.; Hou, S.; Zhang, B.; Mayes, S.; Bennett, M.; Ma, J.; et al. A mini foxtail millet with an Arabidopsis-like life cycle as a C(4) model system. Nat. Plants 2020, 6, 1167–1178. [Google Scholar] [CrossRef]
  44. Meng, Q.; Zhang, R.; Wang, Y.; Zhi, H.; Tang, S.; Jia, G.; Diao, X. Genome-Wide Characterization and Haplotypic Variation Analysis of the YUC Gene Family in Foxtail Millet (Setaria italica). Int. J. Mol. Sci. 2023, 24, 15637. [Google Scholar] [CrossRef]
  45. Ren, X.; Chen, J.; Chen, S.; Zhang, H.; Li, L. Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica, L.). Genes 2023, 14, 2046. [Google Scholar] [CrossRef]
  46. Chen, H.; Ge, W. Identification, Molecular Characteristics, and Evolution of GRF Gene Family in Foxtail Millet (Setaria italica, L.). Front. Genet. 2022, 12, 727674. [Google Scholar] [CrossRef]
  47. Zhang, X.; Zong, J.; Liu, J.; Yin, J.; Zhang, D. Genome-wide analysis of WOX gene family in rice, sorghum, maize, Arabidopsis and poplar. J. Integr. Plant Biol. 2010, 52, 1016–1026. [Google Scholar] [CrossRef]
  48. Sun, R.; Zhang, X.; Ma, D.; Liu, C. Identification and Evolutionary Analysis of Cotton (Gossypium hirsutum) WOX Family Genes and Their Potential Function in Somatic Embryogenesis. Int. J. Mol. Sci. 2023, 24, 11077. [Google Scholar] [CrossRef]
  49. Chen, X.; Hou, Y.; Cao, Y.; Wei, B.; Gu, L. A Comprehensive Identification and Expression Analysis of the WUSCHEL Homeobox-Containing Protein Family Reveals Their Special Role in Development and Abiotic Stress Response in Zea mays, L. Int. J. Mol. Sci. 2023, 25, 441. [Google Scholar] [CrossRef]
  50. Guo, X.; Zhang, H.; Wang, Z.; Li, S. Molecular characterization and transcriptional analysis of VrWOX genes in mungbean [Vigna radiate (L.) Wilczek]. Sheng Wu Gong. Cheng Xue Bao 2023, 39, 566–585. [Google Scholar]
  51. Riccucci, E.; Vanni, C.; Vangelisti, A.; Fambrini, M.; Giordani, T.; Cavallini, A.; Mascagni, F.; Pugliesi, C. Genome-Wide Analysis of WOX Multigene Family in Sunflower (Helianthus annuus, L.). Int. J. Mol. Sci. 2023, 24, 3352. [Google Scholar] [CrossRef]
  52. Yin, S.; Zhao, L.; Liu, J.; Sun, Y.; Li, B.; Wang, L.; Ren, Z.; Chen, C. Pan-genome Analysis of WOX Gene Family and Function Exploration of CsWOX9 in Cucumber. Int. J. Mol. Sci. 2023, 24, 17568. [Google Scholar] [CrossRef] [PubMed]
  53. Tang, L.; He, Y.; Liu, B.; Xu, Y.; Zhao, G. Genome-Wide Identification and Characterization Analysis of WUSCHEL-Related Homeobox Family in Melon (Cucumis melo, L.). Int. J. Mol. Sci. 2023, 24, 12326. [Google Scholar] [CrossRef] [PubMed]
  54. Yang, T.; Gao, T.; Wang, C.; Wang, X.; Chen, C.; Tian, M.; Yang, W. In silico genome wide identification and expression analysis of the WUSCHEL-related homeobox gene family in Medicago sativa. Genom. Inf. 2022, 20, e19. [Google Scholar] [CrossRef] [PubMed]
  55. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  56. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  57. He, Q.; Wang, C.; He, Q.; Zhang, J.; Liang, H.; Lu, Z.; Xie, K.; Tang, S.; Zhou, Y.; Liu, B.; et al. A complete reference genome assembly for foxtail millet and Setaria-db, a comprehensive database for Setaria. Mol. Plant 2024, 17, 219–222. [Google Scholar] [CrossRef] [PubMed]
  58. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  59. Li, H.W.; Li, C.H.; Pao, W.K. Cytological and Genetical Studies of the Interspecific Cross of the Cultivated Foxtail Millet, Setaria Italica (L.) Beauv., and the Green Foxtail Millet, S. Viridis, L.1. Agron. J. 1945, 37, 32–54. [Google Scholar] [CrossRef]
  60. Soyk, S.; Lemmon, Z.H.; Sedlazeck, F.J.; Jiménez-Gómez, J.M.; Alonge, M.; Hutton, S.F.; Van Eck, J.; Schatz, M.C.; Lippman, Z.B. Duplication of a domestication locus neutralized a cryptic variant that caused a breeding barrier in tomato. Nat. Plants 2019, 5, 471–479. [Google Scholar] [CrossRef]
  61. Alonge, M.; Wang, X.; Benoit, M.; Soyk, S.; Pereira, L.; Zhang, L.; Suresh, H.; Ramakrishnan, S.; Maumus, F.; Ciren, D.; et al. Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato. Cell 2020, 182, 145–161.e23. [Google Scholar] [CrossRef]
  62. Lye, Z.N.; Purugganan, M.D. Copy Number Variation in Domestication. Trends Plant Sci. 2019, 24, 352–365. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, Y.; Tan, X.; Paterson, A.H. Paterson, Different patterns of gene structure divergence following gene duplication in Arabidopsis. BMC Genom. 2013, 14, 652. [Google Scholar] [CrossRef] [PubMed]
  64. Doust, A.N.; Devos, K.M.; Gadberry, M.D.; Gale, M.D.; Kellogg, E.A. The genetic basis for inflorescence variation between foxtail and green millet (poaceae). Genetics 2005, 169, 1659–1672. [Google Scholar] [CrossRef] [PubMed]
  65. De Smet, I.; Signora, L.; Beeckman, T.; Inzé, D.; Foyer, C.H.; Zhang, H. An abscisic acid-sensitive checkpoint in lateral root development of Arabidopsis. Plant J. 2003, 33, 543–555. [Google Scholar] [CrossRef] [PubMed]
  66. He, L.; Zhang, B.; Wang, X.; Li, H.; Han, Y. Foxtail millet: Nutritional and eating quality, and prospects for genetic improvement. Front. Agric. Sci. Eng. 2015, 2, 124–133. [Google Scholar] [CrossRef]
  67. Hou, D.; Chen, J.; Ren, X.; Wang, C.; Diao, X.; Hu, X.; Zhang, Y.; Shen, Q. A whole foxtail millet diet reduces blood pressure in subjects with mild hypertension. J. Cereal Sci. 2018, 84, 13–19. [Google Scholar] [CrossRef]
  68. Ren, X.; Yin, R.; Hou, D.; Xue, Y.; Zhang, M.; Diao, X.; Zhang, Y.; Wu, J.; Hu, J.; Hu, X.; et al. The Glucose-Lowering Effect of Foxtail Millet in Subjects with Impaired Glucose Tolerance: A Self-Controlled Clinical Trial. Nutrients 2018, 10, 1509. [Google Scholar] [CrossRef] [PubMed]
  69. Brandt, R.; Cabedo, M.; Xie, Y.; Wenkel, S. Wenkel, Homeodomain leucine-zipper proteins and their role in synchronizing growth and development with the environment. J. Integr. Plant Biol. 2014, 56, 518–526. [Google Scholar] [CrossRef]
  70. Costanzo, E.; Trehin, C.; Vandenbussche, M. The role of WOX genes in flower development. Ann. Bot. 2014, 114, 1545–1553. [Google Scholar] [CrossRef]
  71. Wang, M.M.; Liu, M.M.; Ran, F.; Guo, P.C.; Ke, Y.Z.; Wu, Y.W.; Wen, J.; Li, P.F.; Li, J.N.; Du, H. Global Analysis of WOX Transcription Factor Gene Family in Brassica napus Reveals Their Stress- and Hormone-Responsive Patterns. Int. J. Mol. Sci. 2018, 19, 3470. [Google Scholar] [CrossRef] [PubMed]
  72. Li, X.X.; Liu, C.; Li, W.; Zhang, Z.L.; Gao, X.M.; Zhou, H.; Guo, Y.F. Genome-wide identification, phylogenetic analysis and expression profiling of the WOX family genes in Solanum lycopersicum. Hereditas 2016, 38, 444–460. [Google Scholar] [PubMed]
  73. Chen, Z.; Li, W.; Gaines, C.; Buck, A.; Galli, M.; Gallavotti, A. Structural variation at the maize WUSCHEL1 locus alters stem cell organization in inflorescences. Nat. Commun. 2021, 12, 2378. [Google Scholar] [CrossRef] [PubMed]
  74. Lindsay, P.; Swentowsky, K.W.; Jackson, D. Cultivating potential: Harnessing plant stem cells for agricultural crop improvement. Mol. Plant 2024, 17, 50–74. [Google Scholar] [CrossRef] [PubMed]
  75. Doust, A.N.; Kellogg, E.A. Inflorescence diversification in the panicoid “bristle grass” clade (Paniceae, Poaceae): Evidence from molecular phylogenies and developmental morphology. Am. J. Bot. 2002, 89, 1203–1222. [Google Scholar] [CrossRef] [PubMed]
  76. Cheng, S.; Zhou, D.X.; Zhao, Y. WUSCHEL-related homeobox gene WOX11 increases rice drought resistance by controlling root hair formation and root system development. Plant Signal. Behav. 2016, 11, e1130198. [Google Scholar] [CrossRef] [PubMed]
  77. Nakata, M.T.; Tameshige, T.; Takahara, M.; Mitsuda, N.; Okada, K. The functional balance between the WUSCHEL-RELATED HOMEOBOX1 gene and the phytohormone auxin is a key factor for cell proliferation in Arabidopsis seedlings. Plant Biotechnol. 2018, 35, 141–154. [Google Scholar] [CrossRef] [PubMed]
  78. Minh-Thu, P.T.; Kim, J.S.; Chae, S.; Jun, K.M.; Lee, G.S.; Kim, D.E.; Cheong, J.J.; Song, S.I.; Nahm, B.H.; Kim, Y.K. A WUSCHEL Homeobox Transcription Factor, OsWOX13, Enhances Drought Tolerance and Triggers Early Flowering in Rice. Mol. Cells 2018, 41, 781–798. [Google Scholar] [PubMed]
  79. Rodriguez-Leal, D.; Xu, C.; Kwon, C.T.; Soyars, C.; Demesa-Arevalo, E.; Man, J.; Liu, L.; Lemmon, Z.H.; Jones, D.S.; Van Eck, J.; et al. Lippman, Evolution of buffering in a genetic circuit controlling plant stem cell proliferation. Nat. Genet. 2019, 51, 786–792. [Google Scholar] [CrossRef]
  80. Geng, L.; Tan, M.; Deng, Q.; Wang, Y.; Zhang, T.; Hu, X.; Ye, M.; Lian, X.; Zhou, D.X.; Zhao, Y. Transcription factors WOX11 and LBD16 function with histone demethylase JMJ706 to control crown root development in rice. Plant Cell 2024, koad318. [Google Scholar] [CrossRef]
Genes 15 00476 g001
Figure 1. Phylogenetic analysis of WOX proteins of foxtail millet, green foxtail, tomato, maize, rice, and Arabidopsis. A phylogenetic evolutionary tree was constructed using MEGA-11, where clades of the same color represent that they belong to the same subfamily. The WOX gene IDs of the six species are supplemented in Table S2. Node numbers: bootstrap values, a put back sampling statistical method used to test the credibility of evolutionary tree branches.
Figure 1. Phylogenetic analysis of WOX proteins of foxtail millet, green foxtail, tomato, maize, rice, and Arabidopsis. A phylogenetic evolutionary tree was constructed using MEGA-11, where clades of the same color represent that they belong to the same subfamily. The WOX gene IDs of the six species are supplemented in Table S2. Node numbers: bootstrap values, a put back sampling statistical method used to test the credibility of evolutionary tree branches.
Genes 15 00476 g001
Genes 15 00476 g002
Figure 2. Prediction of the spatial structure of SiWOX proteins. Modern clade: SiWOX2, SiWOX3, SiWOX7, SiWOX8, SiWOX10, SiWOX11, and SiWOX12; Intermediate clade: SiWOX1, SiWOX4, SiWOX5, SiWOX9, and SiWOX13; Ancient clade: SiWOX6. The spatial structure predictions of SiWOX proteins links: https://alphafold.ebi.ac.uk/ (accessed on 5 December 2023).
Figure 2. Prediction of the spatial structure of SiWOX proteins. Modern clade: SiWOX2, SiWOX3, SiWOX7, SiWOX8, SiWOX10, SiWOX11, and SiWOX12; Intermediate clade: SiWOX1, SiWOX4, SiWOX5, SiWOX9, and SiWOX13; Ancient clade: SiWOX6. The spatial structure predictions of SiWOX proteins links: https://alphafold.ebi.ac.uk/ (accessed on 5 December 2023).
Genes 15 00476 g002
Genes 15 00476 g003
Figure 3. Analysis of SiWOX gene structure and conserved motifs. (A) Phylogenetic tree. (B) Gene structure. Coding sequences (CDS) and untranslated regions (UTRs) are represented by different colored boxes, while introns are indicated by lines. (C) Conserved motifs. Conserved motifs within the SiWOX genes are represented by boxes of different colors. The weblogo diagram of the twelve conserved motifs is provided in Figure S2.
Figure 3. Analysis of SiWOX gene structure and conserved motifs. (A) Phylogenetic tree. (B) Gene structure. Coding sequences (CDS) and untranslated regions (UTRs) are represented by different colored boxes, while introns are indicated by lines. (C) Conserved motifs. Conserved motifs within the SiWOX genes are represented by boxes of different colors. The weblogo diagram of the twelve conserved motifs is provided in Figure S2.
Genes 15 00476 g003
Genes 15 00476 g004
Figure 4. Chromosome distribution and gene tandem repeats of SiWOX genes. Gray lines represent tandem repeats between foxtail millet genomes, while red lines represent tandem repeats between SiWOX genes.
Figure 4. Chromosome distribution and gene tandem repeats of SiWOX genes. Gray lines represent tandem repeats between foxtail millet genomes, while red lines represent tandem repeats between SiWOX genes.
Genes 15 00476 g004
Genes 15 00476 g005
Figure 5. Collinearity analysis of WOX genes in Setaria viridis, Setaria italica, Zea mays, and Oryza sativa. Gray lines indicate collinearity between genes from different species, while purple lines indicate collinearity between WOX genes from different species. Sv: Setaria viridis, Si: Setaria italica, Zm: Zea mays, Os: Oryza sativa.
Figure 5. Collinearity analysis of WOX genes in Setaria viridis, Setaria italica, Zea mays, and Oryza sativa. Gray lines indicate collinearity between genes from different species, while purple lines indicate collinearity between WOX genes from different species. Sv: Setaria viridis, Si: Setaria italica, Zm: Zea mays, Os: Oryza sativa.
Genes 15 00476 g005
Genes 15 00476 g006
Figure 6. Expression pattern analysis of SiWOX genes in different tissues. The heatmap of the expression profiles of the SiWOX gene in different developmental stages is represented by normalized values using RNA-seq data, with the color from red to blue indicating the expression levels from high to low. Tissue information for SiWOX expression analysis is provided in Table S4.
Figure 6. Expression pattern analysis of SiWOX genes in different tissues. The heatmap of the expression profiles of the SiWOX gene in different developmental stages is represented by normalized values using RNA-seq data, with the color from red to blue indicating the expression levels from high to low. Tissue information for SiWOX expression analysis is provided in Table S4.
Genes 15 00476 g006
Genes 15 00476 g007
Figure 7. Expression patterns of SiWOX genes in foxtail millet panicles and roots. (A) Expression patterns of SiWOX genes in panicles at different developmental stages: Stage 1, panicle length approximately 1.0–1.5 mm; stage 2, panicle length approximately 2.5–3.0 mm. (B) Expression patterns of SiWOX genes in 9-day-old seedling roots treated with ABA (2 μM) and without ABA (CK). Bar graphs represent differences in transcriptome sequencing duplicates between spike samples (A) and root samples (B). Statistical significance between panicle (A) and root (B) stages was determined by t-test (* p < 0.05). TPM: transcripts per million.
Figure 7. Expression patterns of SiWOX genes in foxtail millet panicles and roots. (A) Expression patterns of SiWOX genes in panicles at different developmental stages: Stage 1, panicle length approximately 1.0–1.5 mm; stage 2, panicle length approximately 2.5–3.0 mm. (B) Expression patterns of SiWOX genes in 9-day-old seedling roots treated with ABA (2 μM) and without ABA (CK). Bar graphs represent differences in transcriptome sequencing duplicates between spike samples (A) and root samples (B). Statistical significance between panicle (A) and root (B) stages was determined by t-test (* p < 0.05). TPM: transcripts per million.
Genes 15 00476 g007
Genes 15 00476 g008
Figure 8. Expression pattern of SiWOX genes in response to plant hormones. Six plant hormones, 6-BA (100 μM), ABA (100 μM), GA3 (100 μM), SA (100 μM), MeJA (100 μM), and IAA (100 μM), were used to treat 28-day-old seedlings, and leaf samples were taken at 0 h, 0.5 h, 2 h, 6 h, and 12 h post-treatment. The gene expression levels at different time intervals were detected using the qPCR method. Line graphs represent changes in gene expression levels.
Figure 8. Expression pattern of SiWOX genes in response to plant hormones. Six plant hormones, 6-BA (100 μM), ABA (100 μM), GA3 (100 μM), SA (100 μM), MeJA (100 μM), and IAA (100 μM), were used to treat 28-day-old seedlings, and leaf samples were taken at 0 h, 0.5 h, 2 h, 6 h, and 12 h post-treatment. The gene expression levels at different time intervals were detected using the qPCR method. Line graphs represent changes in gene expression levels.
Genes 15 00476 g008
Table 1. The properties of SiWOX proteins.
Table 1. The properties of SiWOX proteins.
Gene NameGene IDNumber of Amino AcidsMolecular WeightIsoelectric PointInstability IndexAliphatic IndexGRAVYSubcellular Localization
SiWOX1Si2g4319026327.507.8462.0663.57−0.264Nucleus
SiWOX2Si3g0366032533.498.6160.862.06−0.306Nucleus
SiWOX3Si3g0520027328.597.6666.9269.3−0.137Nucleus
SiWOX4Si3g1456048350.916.9267.9282.11−0.19Nucleus
SiWOX5Si5g2713053154.807.1657.2670.96−0.2Nucleus
SiWOX6Si5g3634027830.91650.2659.42−0.733Nucleus
SiWOX7Si5g3774033135.198.9378.1262.51−0.419Nucleus
SiWOX8Si5g3863021224.078.3652.1268.07−0.703Nucleus
SiWOX9Si6g0943025727.897.2264.865.41−0.352Nucleus
SiWOX10Si7g2698022825.399.2666.3274.87−0.534Nucleus
SiWOX11Si7g3351022024.018.8368.0452.95−0.673Nucleus
SiWOX12Si8g0112022524.558.9871.0751.33−0.732Nucleus
SiWOX13Si9g4113031232.006.746963.27−0.249Nucleus
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nan, L.; Li, Y.; Ma, C.; Meng, X.; Han, Y.; Li, H.; Huang, M.; Qin, Y.; Ren, X. Identification and Expression Analysis of the WOX Transcription Factor Family in Foxtail Millet (Setaria italica L.). Genes 2024, 15, 476. https://doi.org/10.3390/genes15040476

AMA Style

Nan L, Li Y, Ma C, Meng X, Han Y, Li H, Huang M, Qin Y, Ren X. Identification and Expression Analysis of the WOX Transcription Factor Family in Foxtail Millet (Setaria italica L.). Genes. 2024; 15(4):476. https://doi.org/10.3390/genes15040476

Chicago/Turabian Style

Nan, Lizhang, Yajun Li, Cui Ma, Xiaowei Meng, Yuanhuai Han, Hongying Li, Mingjing Huang, Yingying Qin, and Xuemei Ren. 2024. "Identification and Expression Analysis of the WOX Transcription Factor Family in Foxtail Millet (Setaria italica L.)" Genes 15, no. 4: 476. https://doi.org/10.3390/genes15040476

APA Style

Nan, L., Li, Y., Ma, C., Meng, X., Han, Y., Li, H., Huang, M., Qin, Y., & Ren, X. (2024). Identification and Expression Analysis of the WOX Transcription Factor Family in Foxtail Millet (Setaria italica L.). Genes, 15(4), 476. https://doi.org/10.3390/genes15040476

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

Article Metrics

Back to TopTop