Next Article in Journal
Generation of Mutants from the 57B Region of Drosophila melanogaster
Previous Article in Journal
Ovine KRTAP36-2: A New Keratin-Associated Protein Gene Related to Variation in Wool Yield
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 11
10.3390/genes14112046
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

Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica L.)

by
Xuemei Ren
*,†,
Jinjie Chen
†,
Shuwan Chen
,
Hui Zhang
and
Li Li
College of Life Science, Shanxi Agricultural University, Jinzhong 030801, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2023, 14(11), 2046; https://doi.org/10.3390/genes14112046
Submission received: 18 September 2023 / Revised: 31 October 2023 / Accepted: 3 November 2023 / Published: 6 November 2023
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
The CLAVATA3/EMBRYO-SURROUNDING REGION (CLE) genes encode signaling peptides that play important roles in various developmental and physiological processes. However, the systematic identification and characterization of CLE genes in foxtail millet (Setaria italica L.) remain limited. In this study, we identified and characterized 41 SiCLE genes in the foxtail millet genome. These genes were distributed across nine chromosomes and classified into four groups, with five pairs resulting from gene duplication events. SiCLE genes within the same phylogenetic group shared similar gene structure and motif patterns, while 34 genes were found to be single-exon genes. All SiCLE peptides harbored the conserved C-terminal CLE domain, with highly conserved positions in the CLE core sequences shared among foxtail millet, Arabidopsis, rice, and maize. The SiCLE genes contained various cis-elements, including five plant hormone-responsive elements. Notably, 34 SiCLE genes possessed more than three types of phytohormone-responsive elements on their promoters. Comparative analysis revealed higher collinearity between CLE genes in maize and foxtail millet, which may be because they are both C4 plants. Tissue-specific expression patterns were observed, with genes within the same group exhibiting similar and specific expression profiles. SiCLE32 and SiCLE41, classified in Group D, displayed relatively high expression levels in all tissues except panicles. Most SiCLE genes exhibited low expression levels in young panicles, while SiCLE6, SiCLE24, SiCLE25, and SiCLE34 showed higher expression in young panicles, with SiCLE24 down-regulated during later panicle development. Greater numbers of SiCLE genes exhibited higher expression in roots, with SiCLE7, SiCLE22, and SiCLE36 showing the highest levels and SiCLE36 significantly down-regulated after abscisic acid (ABA) treatment. Following treatments with ABA, 6-benzylaminopurine (6-BA), and gibberellic acid 3 (GA3), most SiCLE genes displayed down-regulation followed by subsequent recovery, while jasmonic acid (JA) and indole-3-acetic acid (IAA) treatments led to upregulation at 30 min in leaves. Moreover, identical hormone treatments elicited different expression patterns of the same genes in leaves and stems. This comprehensive study enhances our understanding of the SiCLE gene family and provides a foundation for further investigations into the functions and evolution of SiCLE genes in foxtail millet.

1. Introduction

Multicellular organisms often rely on intercellular communication to coordinate growth and development, regulate physiological activities, and ensure efficient and orderly operation. In plants, intercellular communication is primarily achieved through hormones, peptide ligands, small RNAs, and even gases [1]. Peptide ligands, specifically small peptides with less than 100 amino acids in length, play diverse roles in plant growth, development, and reproduction, as well as symbiotic interactions and stress responses [2]. The CLE (CLAVATA3/Endosperm surrounding region-related) gene family encodes the largest known peptides in plants. Research has demonstrated that each CLE gene encodes a precursor peptide containing a highly conserved CLE domain. The CLE conserved motif undergoes proteolytic cleavage from the precursor peptide and then experiences further post-translational modifications, such as hydroxylation and glycosylation, to produce mature CLE peptides [3]. The mature peptide can be secreted into the extracellular space and subsequently perceived by receptors located on the plasma membrane of adjacent cells, thereby mediating intercellular communication [4].
The CLE peptide family exhibits widespread distribution in various plant tissues and organs, displaying tissue-specific expression patterns [5]. Specific CLE peptides have been found to exert influence on plant development in specific tissues [6]. The maintenance of shoot apical meristem (SAM) mainly relies on the receptor kinases mediated CLV-WUS (CLAVATA-WUSCHEL) negative feedback regulatory loop [7,8]. Within this regulatory loop, the transcription factor WUS functions as a central regulatory element in maintaining the activity of stem cells within the SAM [9]. The CLV-WUS negative feedback pathway that maintains the homeostasis of the SAM was first discovered in Arabidopsis thaliana. However, this pathway is highly conserved in higher plants. In rice (Oryza sativa), the FON2 (FLORAL ORGAN NUMBER 2) and FCP1 (FON2-like CLE PROTEIN 1) genes encode homologous proteins of CLV3, both containing conserved CLE domains. However, unlike in Arabidopsis, where CLV3 simultaneously regulates the SAM both in vegetative and reproductive phases, FON2 and FCP1, respectively, maintain the homeostasis of the SAM during the reproductive and vegetative growth phases [10,11]. The rice gene FON1, the maize gene TD1 (THICK TASSEL DWARF 1), and the tomato (Solanum lycopersicum L.) gene FAB (FASCIATED AND BRANCHED) encode homologous proteins of CLV1 [12,13,14]. The maize gene FEA2 (FASCIATED EAR 2) and the tomato gene SlCLV2 encode homologous proteins of CLV2 [14,15]. Mutations in these genes lead to enlarged SAM and the formation of additional floral organs.
The regulation of plant development by CLE peptides often involves the production and participation of phytohormones, and the interaction between phytohormone signaling and CLE peptides has been widely studied. The cle42 mutant displayed earlier senescence phenotypes, while overexpression of CLE42 or application of the synthesized 12-amino-acid peptide (CLE42p) delayed leaf senescence under natural and dark conditions [16]. Further research found that CLE42 can delay leaf senescence by suppressing ethylene biosynthesis in Arabidopsis. The classical phytohormone cytokinin and the small CLE peptides are potent regulators of cell division and cell differentiation. CLE and cytokinin signaling are highly intertwined developmental regulators with antagonistic functions in shoots and synergistic functions in roots [17]. The application of gibberellic acid 4 (GA4) to the shoot apex of the GA-deficient mutant ga3ox1/ga3ox2 induces substantial accumulation of CLE6 in the root and rescued the short-root phenotype, which suggests that CLE6 plays a systemic role in shoot growth under the influence of GA in Arabidopsis [18]. AtCLE9 is expressed in stomata and acts as an essential regulator in the induction of stomatal closure, while in abscisic acid (ABA)-deficient mutants, CLE9-induced stomatal closure was impaired, indicating that ABA is required for CLE9-mediated guard cell signaling [19]. Arabidopsis CLE5 and CLE6 are expressed specifically at the base of developing leaves and floral organs, and their transcript levels are regulated by auxin to modulate the final leaf morphology [20]. The mutual cooperation or antagonistic regulation between CLE peptides and plant hormones in various plant tissues is gradually being discovered.
The CLE gene family is gradually being discovered in various species. So far, the CLE gene family has been identified and analyzed in several species, such as Arabidopsis [1], tomato [21,22,23], cucumber [24], soybean [25], Populus trichocarpa [26], grape [27], wheat [28], Brassica napus [29], cotton [30], and potato [31]. Foxtail millet is an important C4 cereal crop that has outstanding drought resistance and barren resistance. In addition, foxtail millet is a diploid (2n = 2X = 18) self-pollinated crop, which has the advantages of a small genome (about 400 Mb), few repeat sequences, and easy mutagenic and mutant screening, and it is developing into a new model plant for C4 grasses [32]. However, no studies related to the CLE gene family in foxtail millet have been reported. In this study, we identified 41 SiCLE genes by whole-genome analysis in foxtail millet and conducted a phylogenetic tree of the CLE genes from four species, including foxtail millet, Arabidopsis, maize, and rice. Focusing on SiCLE genes in foxtail millet, we analyze the chromosome location, gene structure and conserved motifs, cis-acting elements, and tissue-specific expression pattern based on RNA-seq data. Further, we examine the response of SiCLE genes to different plant hormones in foxtail millet. Our study provides a comprehensive understanding of the SiCLE gene family that opens the gate to further detailed studies on the gene function and evolution of SiCLE genes in foxtail millet.

2. Materials and Methods

2.1. Genome-Wide Identification of SiCLE Genes in Foxtail Millet

The foxtail millet genome, proteome, coding sequences, and GFF annotation file were downloaded from MDSi, Multi-omics Database for S. italica (http://foxtail-millet.biocloud.net/home, accessed on 26 November 2021) [33]. CLE peptide sequences of Arabidopsis (https://www.arabidopsis.org/, accessed on 26 November 2021), maize (https://maizegdb.org/, accessed on 26 November 2021), and rice (https://rapdb.dna.affrc.go.jp/, accessed on 26 November 2021) were downloaded and used as queries for the Basic Local Alignment Sequence Tool for Protein (BLASTP) [34] to search the proteomes of foxtail millet with E-value of 1 × 10−5. Then, the Hidden Markov Models (HMMs) of the CLE domain were used to perform HMM searches on the protein sequences of foxtail millet using the HMM-based search program Hmmer3 (http://hmmer.org, accessed on 26 November 2021) [35,36]. The retrieved candidate gene sequences were further filtered through BLASTP screening. Finally, based on important features of the CLE gene family, a more thorough screening was conducted to determine the final SiCLE genes in foxtail millet. The SiCLE peptide sequences were subjected to protein physical and chemical property analysis, including molecular weight, isoelectric point, and hydrophilicity, using the ExPASy Proteomics (http://www.expasy.org, accessed on 20 March 2022) [37].

2.2. Phylogenetic Analysis and Classification of SiCLE Genes in Foxtail Millet

MAFFT (https://mafft.cbrc.jp/alignment/server/, accessed on 21 March 2022) was used to perform multiple sequence alignment of the CLE protein sequences of Arabidopsis, rice, maize, and foxtail millet [38]. The alignment results were imported into MEGA7.0 software, and a neighbor-joining algorithm with 1000 bootstrap replicates was used to construct the phylogenetic tree [39]. The iTOL (https://itol.embl.de, accessed on 21 March 2022) online tools were used to beautify the phylogenetic tree [40].

2.3. Gene Structure and Conserved Motif Analysis of SiCLE Genes in Foxtail Millet

GSDS2.0 (http://gsds.cbi.pku.edu.cn, accessed on 28 November 2021 ) was used to visualize the gene structures of the candidate genes [41]. The online tool Multiple Expectation Maximization for Motif Elucidation suite (MEME, https://meme-suite.org/, accessed on 28 November 2021) was used to search for conserved motifs with a maximum parameter value of 5 and the optimum width was 6 to 200 residue width with an E-value < 1 × 10−10 was retained [42].

2.4. Cis-Regulatory Elements Analysis of SiCLE Genes in Foxtail Millet

Promoter sequences (2000 bp of sequence upstream of the transcription start site) for SiCLE genes were extracted by TBtools and scanned for cis-regulatory elements on the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 30 March 2022) [43]. The analysis results were plotted using R scripts.

2.5. Cene Duplication and Collinearity Analysis of the SiCLEs in Foxtail Millet

The chromosomal localization information for each SiCLE gene was obtained from the gene annotation files, and the results were visualized using the Graphics function in TBtools [44]. The SiCLE genes were analyzed for duplication events using MCScanX with default parameters. The Dual Synteny Plotter tool in TBtools was used to analyze the conserved synteny between foxtail millet and three other species (Arabidopsis, rice, and maize).

2.6. Expression Pattern Analysis of SiCLE Genes in Foxtail Millet

Transcriptome data of different tissues (roots, stems, leaves, grains, and panicles) from Jingu21 at various developmental stages were obtained from the MDSi database. Germinated seeds were sampled at 3 days. Plants were sampled at one-tip-two-leaf stage. Leaf 2 was sampled at heading stage. Neck–panicle internodes, flag–leaf, flag–leaf–sheath, stem, leaf 1, leaf sheath 1, and root were sampled during the filling stage. Immature seeds were sampled at early grain filling stage (S1), middle grain filling stage (S2), late grain filling stage (S3), final grain filling stage (S4), and grain maturation stage (S5) [45]. Immature spikelets were sampled at the S2 and S4 stages. Panicle 1 are primary branches, and panicle 2 are tertiary branches at panicle–branch differentiation stage. TPM values were calculated to assess gene expression levels, and a heatmap was generated using R scripts [46]. Young panicle samples were collected from Jingu21 at two distinct formation stages: stage 1 (approximately 1.0–1.5 mm) when branch meristems were specified and stage 2 (approximately 2.5–3.0 mm) when branch meristems were clearly formed. Root samples were collected from 9-day-old Yugu 1 seedlings treated with ABA (2 μM) and untreated control (CK) in a greenhouse under the following conditions: light intensity of 50,000 LX, 16 h daylight at 28 °C, and 8 h darkness at 22 °C. Subsequently, young panicle and root samples were rapidly frozen in liquid nitrogen and stored at −80 °C. Novogene Bioinformatics Technology Co., Ltd. (Beijing, China) was commissioned to construct the RNA libraries. High-throughput sequencing was performed on the Illumina Hiseq platform with three biological replicates.

2.7. Response of SiCLE Genes to Phytohormones in Foxtail Millet

The seedlings of Jingu 21 were cultured with Hoagland hydroponic nutrient solution (Beijing Coolaber Technology Co., Ltd., Beijing, China) and grown in greenhouse with light intensity of 50,000 LX for 16 h in the daylight at 28 °C and 8 h in the dark at 22 °C. Seedlings (28 days old) were treated with plant hormones (Beijing Coolaber Technology Co., Ltd., Beijing, China) including 6-benzylaminopurine (6-BA) (100 µM), abscisic acid (ABA) (100 µM), GA3 (100 µM), jasmonic acid (JA) (100 µM), and indole-3-acetic acid (IAA) (100 µM). Samples of leaves and stems were collected at 0 h, 0.5 h, 2 h, 6 h, and 12 h after the treatments. The samples were frozen in liquid nitrogen and stored at −80 °C. Three replicates were performed for each treatment. RNA from different tissues was extracted by Trizol method using Total RNA Extract Reagent and RNA Extraction solution (Beijing Coolaber Technology Co., Ltd., Beijing, China). Reverse transcription was performed using All-in-One First-Strand Synthesis MasterMix (with dsDNase) (BestEnzymes Biotech Co., Ltd., Lianyungang, China) for Real-time Quantitative Polymerase Chain Reaction (qPCR). F488 SYBR qPCR Mix (Universal) (BestEnzymes Biotech Co., Ltd. Lianyungang, China) was used as a fluorescent dye. Primers were designed using Primer 5.0 (Table S1). We used the gene Si9g37480 as an internal control, which was stably expressed at each growth stage in almost all tissues [32]. Each reaction was performed three times, and the 2−ΔΔCT method [47] was used to calculate the relative gene expression levels.

3. Results

3.1. Identification of the SiCLE Genes in Foxtail Millet

The whole genome was comprehensively analyzed using HMM search and multiple rounds of BLASTP to identify genes with a significance threshold of E < 10−5. Candidate genes with a sequence length exceeding 300 amino acids were excluded based on CLE gene characteristics [48]. Furthermore, candidate peptide sequences that contained the conserved CLE motif at the C-terminus were retained after motif prediction using MEME [42]. A total of 41 candidate CLE genes were successfully identified in the foxtail millet genome, and they were named SiCLE1-SiCLE41 based on their respective chromosomal locations (Table 1). The CLE genes are found in diploid Arabidopsis (32 genes) [1], tomato (52 genes) [22], and cucumber (26 genes) [24], and the number of SiCLE genes in foxtail millet is between Arabidopsis and tomato. Nevertheless, this number was higher in polyploid species soybean (84 genes) [25], wheat (104 genes) [28], and B. napus (116 genes) [29], which may be attributed to complex gene duplication events throughout evolution.
Physicochemical analysis revealed that the 41 SiCLE peptides ranged from 77 to 127 amino acids, with predicted molecular weights ranging from 7.93 to 13.42 kDa and isoelectric points ranging from 5.45 to 12.16 (Table 1). The majority of SiCLE peptides displayed alkaline properties, with only six being acidic (SiCLE1, SiCLE9, SiCLE14, SiCLE20, SiCLE27, and SiCLE36). Signal peptide prediction identified the presence of signal peptides in 36 out of the 41 CLE peptides. Based on the average hydrophobicity index (grand average of hydropathicity, GRAVY) principle, which classifies proteins as biphasic with indices ranging from −0.5 to 0.5, hydrophobic with positive values, and hydrophilic with negative values, all SiCLE proteins were categorized as biphasic proteins. Despite variations in the number of CLE genes among different species, the physicochemical properties of the peptides have not significantly changed [48].

3.2. Phylogenetic Analysis and Classification of SiCLE Peptides in Foxtail Millet

To evaluate the evolutionary relationship of the SiCLE peptides in foxtail millet, Arabidopsis, rice, and maize, a phylogenetic tree using the neighbor-joining (NJ) method was constructed based on the full-length sequences of 41 SiCLE peptides with other CLE peptides identified in the three species (Table S2). All the CLE members clustered into four groups, named Groups A-D, in the phylogenetic analysis (Figure 1). Group C and Group D contain a relatively larger number of CLE genes, with eleven and twenty-one SiCLE genes, respectively. Group C CLE peptides include AtCLV3, OsFCP1, OsFCP2, OsFON2/OsFON4, ZmFCP1, and ZmCLE7 [1,10,11,49,50], which all participate in the stem cell differentiation and maintenance in SAM. The eleven SiCLE genes in Group C may play a conserved function in the stem cell differentiation and maintenance in the SAM of foxtail millet. Almost all the AtCLE peptides in Group D are called root-active CLE peptides, and they all led to short-root, suppressed protophloem differentiation and consumption of root apical meristem (RAM) phenotypes by applications of individual peptides in vitro [51,52]. The twenty-one SiCLE genes may be involved in the stem cell homeostasis in the root of foxtail millet. Group B has the least number of CLE genes and contains only two SiCLE genes, namely SiCLE6 and SiCLE24. AtCLE41 and AtCLE44 belong to Group B. It has been identified that AtCLE41p and AtCLE44p have the ability to inhibit the differentiation of mesophyll cells into xylem cells in an in vitro xylogenic culture system [53]. SiCLE6 and SiCLE24 may have similar functions in vascular bundle development in foxtail millet. The remaining seven SiCLE genes belong to Group A, and the mechanism of CLE genes in this group has not yet been resolved clearly. The results indicate a high degree of conservation and similar functions within each group, suggesting that SiCLE peptides in Groups C and D may play similarly important roles in maintaining homeostasis of the SAM and RAM, respectively. Additionally, Group B may be involved in vascular bundle development.

3.3. Gene Structure and Conserved Motif Analysis of SiCLE Peptides in Foxtail Millet

Based on genome annotation, we conducted a comprehensive analysis of the gene structure of SiCLE genes (Figure 2). The majority of SiCLE genes exhibited either intron-less structures or had only a few introns. Out of the total 41 SiCLE genes, 34 were single-exon genes, 6 were two-exon genes, and 1 was a three-exon gene. Moreover, more than half (24/41) of the SiCLE genes lacked untranslated regions (UTRs), while the remaining 17 SiCLE genes possessed either 3′UTRs or both 5′ and 3′UTRs. The lack of introns and UTRs in SiCLE genes may be attributed to genomic simplification, which minimizes non-coding regions to enhance genome efficiency and replication speed [54]. Notably, SiCLE genes within the same phylogenetic cluster demonstrated similar features despite displaying diverse gene structures.
To identify conserved motifs in SiCLE peptides, we utilized the MEME website. A total of five conserved motifs were identified, ranging in length from 7 to 36 amino acids (Figure S1). The CLE precursor peptide is characterized by three main domains: the N-terminal signal peptide region, the central variable domain, and the C-terminal conserved CLE domain, which is highly conserved and approximately 12–13 amino acids in length. The predicted motifs obtained from the MEME online software aligned well with the reported features of CLE precursors [22]. Motif 1, located near the C-terminus, encompassed the typical CLE domain and was present in all 41 members. Motif 2, positioned close to the N-terminus, was found in the majority of members. Motifs 3, 4, and 5 were scattered throughout the middle region of different members, with only a few members containing them.
To investigate the conservation of the CLE core sequences comprising 12–13 amino acids within Motif 1, we performed a comparative analysis among foxtail millet, Arabidopsis, rice, and maize (Figure S2). The results revealed a high degree of similarity between foxtail millet and the other three species. Specifically, positions 2, 4, 5, 7, 8, 10, and 12 exhibited notable conservation, indicating their significance in CLE peptide maturation and functionality [26,55].

3.4. Cis-Regulatory Elements in SiCLE Gene Promoters

Cis-elements play a crucial role as molecular switches in the transcriptional regulation of genes during plant growth and responses to abiotic stress [56,57,58]. To investigate the potential regulatory mechanisms of SiCLE genes in the growth and development of foxtail millet, the 2-kb promoter sequences upstream from the translation start codon of 41 SiCLEs were retrieved from foxtail millet genome sequences and analyzed using the PlantCARE website. We identified eight types of cis-acting elements grouped into three major categories (Figure 3, Table S3). These elements included plant hormone-responsive elements, including auxin, gibberellin, SA, ABA, and Methyl-jasmonic acid (Me-JA); stress-responsive elements contain low temperature and defense and stress; and elements involved in meristem expression regulation. Notably, multiple types of cis-acting elements were unevenly distributed across each SiCLE gene. In accordance with previous findings [24], it was observed that each SiCLE gene in foxtail millet contained at least one phytohormone-related cis-acting element. Most (34/41) SiCLE genes possessed more than three types of phytohormone-responsive elements on their promoters, with the exception of SiCLE5, SiCLE10, SiCLE20, SiCLE24, SiCLE36, SiCLE38, and SiCLE39. Notably, six SiCLE genes (SiCLE4, SiCLE8, SiCLE9, SiCLE16, SiCLE18, and SiCLE26) exhibited all five types of plant hormone-responsive elements, while SiCLE26 had all eight types of elements. Among the various types of cis-acting elements, ABA responsiveness elements were found to be the most abundant in the promoters of SiCLE genes in foxtail millet, with an average of 3.66 elements per gene. Notably, SiCLE16 exhibited the highest number of ABA responsiveness elements, with a total of 11 ABA elements present in its promoter region. Only SiCLE28 and SiCLE38 contained no elements involved in ABA responsiveness. On the other hand, the number of cis-acting elements involved in defense and stress responsiveness was the lowest, and only nine SiCLE genes (SiCLE9, SiCLE17, SiCLE20, SiCLE26, SiCLE29, SiCLE30, SiCLE34, SiCLE37, and SiCLE41) had no more than two elements involved in defense and stress responsiveness. These findings suggest that SiCLE genes are not only involved in tissue development processes but also exhibit responsiveness to various plant hormones.

3.5. Chromosomal Distribution, Gene Duplication, and Synteny Analysis of the SiCLE Genes

The chromosomal localization analysis of the 41 SiCLE genes revealed an uneven distribution across the nine chromosomes of foxtail millet (Figure 4). Among them, Chr3, Chr5, and Chr9 exhibited an equal number of genes, each hosting eight SiCLE genes, accounting for 19.51% of the total SiCLE genes. Following closely were Chr1 with seven genes (17.07%), Chr4 with four genes (9.76%), and Chr2 and Chr7 with two genes each (4.88%). Chr6 and Chr8 had the least distribution, each containing only one SiCLE gene. We then defined tandem duplication events as the presence of two or more genes within a 200 kb chromosomal region [59]. In the SiCLE gene family, we identified five tandem duplication events. These included three events on chromosomes Chr3, Chr4, and Chr9, namely SiCLE13 and SiCLE14, SiCLE19 and SiCLE20, and SiCLE36 and SiCLE37, respectively. Additionally, two events were found on Chr5, namely SiCLE26 and SiCLE27, and SiCLE27 and SiCLE28. Among the five pairs of genes with duplication events, four pairs of SiCLE genes belonged to Group C, and one pair belonged to Group D. These duplication events are the primary driving force behind the expansion of SiCLE genes, and it is possible that Group C and Group D, which possess a relatively large number of SiCLE genes, underwent expansion during the whole-genome duplication process. Gene duplication events also led to changes in gene structure. For instance, SiCLE19 has both 5′UTR and 3′UTR, while SiCLE20 only possesses a 3′UTR. SiCLE37 has an additional 3′UTR compared to SiCLE36.
To further study the evolutionary mechanism of the SiCLE genes in foxtail millet, a collinear map of foxtail millet and three representative species was constructed, including two monocotyledons (rice and maize) and one dicotyledon (Arabidopsis) (Figure 5). As depicted in Figure 5, no collinearity was observed between CLE genes in Arabidopsis and foxtail millet, which could be attributed to the distant phylogenetic relationship between dicotyledonous Arabidopsis and monocotyledonous foxtail millet. In monocotyledonous plants, 22 pairs of collinear CLE genes were identified between foxtail millet and rice, while the highest number of collinear CLE genes was found between foxtail millet and maize, with 30 pairs of collinear CLE genes. Notably, the CLE genes in maize and foxtail millet exhibited higher collinearity, which may be due to the fact that they are all C4 plants.

3.6. Expression Pattern Analysis of SiCLE Genes in Different Tissues

To comprehensively understand the functions of SiCLE genes in foxtail millet, we conducted an analysis of the transcript levels of the 41 SiCLE genes in various tissues (root, stem, leaf, panicle, and spikelet) at different growth stages (three days after germination, one-tip-two-leaf stage, heading stage, filling stage, and panicle differentiation stage) in the Jingu 21 variety. We utilized publicly available transcriptome data from MDSi [33] for our analysis. The TPM values of the SiCLE genes are listed in Table S4, and a heatmap was generated to display the expression patterns of the SiCLE genes (Figure 6). First of all, the expression levels of the majority of SiCLE genes remained low throughout the entire period, while a subset of SiCLE genes exhibited distinct spatiotemporal specificity. For instance, SiCLE19 showed higher expression exclusively in roots during the filling stage, whereas SiCLE14 displayed higher expression in immature seed S5 at the same stage. Additionally, genes clustered within the same group demonstrated similar expression patterns but also exhibited differences. Both SiCLE32 and SiCLE41 were classified in Group D, and their expression levels were relatively high in all the tissues except for the panicles. SiCLE3 and SiCLE17, also belonging to Group D, exhibited relatively high expression levels in immature spikelets and immature seeds, while their expression levels in other tissues were low. SiCLE25, also classified in Group D, displayed a similar expression pattern to SiCLE3 but exhibited higher expression in the panicles, germinated seeds, and one-tip-two-leaf plants. SiCLE23, another member of Group D, demonstrated high expression in the neck–panicle internodes and stems during the filling stage. Both SiCLE6 and SiCLE24 were clustered in Group B, with higher expression levels observed in the panicles, germinated seeds, and one-tip-two-leaf plants, while their expression levels were relatively lower in other tissues. These observations highlight the tissue-specific expression of genes within this family, as well as the similar and specific expression patterns within the same group.
CLE peptides play a critical role in regulating stem cell maintenance in plants [60]. Stem cells located within the SAM, RAM, and procambium are established during embryogenesis and continue to generate new tissues and organs during post-embryonic growth and development. The SAM is responsible for the formation of inflorescence meristems, which give rise to flowers and panicles, while the RAM is involved in root growth and development. Transcriptome sequencing was conducted to investigate the expression levels of SiCLE genes in young panicles at the formation stage and roots of seedlings (Figure 7). Our analysis revealed that the expression levels of most SiCLE genes were relatively low in the young panicle (Figure 7A). Notably, SiCLE6, SiCLE24, SiCLE25, and SiCLE34 exhibited higher expression levels compared to SiCLE16, SiCLE33, SiCLE35, and SiCLE38. Additionally, the remaining SiCLE genes showed lower expression levels compared to these eight genes. Based on their expression patterns, these eight genes are included in four groups of the phylogenetic tree: Group A (SiCLE34, SiCLE35, and SiCLE38), Group B (SiCLE6 and SiCLE24), Group C (SiCLE33), and Group D (SiCLE16 and SiCLE25). In particular, the expression of SiCLE24 and SiCLE32 was significantly reduced in the young panicle at stage 2 compared to stage 1. Conversely, SiCLE9, SiCLE23, SiCLE28, and SiCLE38 exhibited higher expression in the later stage. Compared to the young panicle, a greater number of SiCLE genes exhibited higher expression levels in roots (Figure 7B). Notably, SiCLE7, SiCLE22, and SiCLE36 showed the highest expression in roots, followed by SiCLE41, SiCLE1, SiCLE17, SiCLE3, SiCLE34, and SiCLE35, while the remaining genes displayed lower expression levels. These nine SiCLE genes are classified into four groups within the phylogenetic tree: Group A (SiCLE34 and SiCLE35), Group C (SiCLE1), and Group D (SiCLE3, SiCLE7, SiCLE16, SiCLE17, SiCLE22, SiCLE36, and SiCLE41). Moreover, ABA plays a crucial role in regulating root development and growth. It influences root meristem activity, elongation, and branching, thereby modulating overall root architecture [61]. In order to detect the response of SiCLE genes after ABA treatment, transcriptome sequencing analysis was performed on the roots of foxtail millet seedlings treated with ABA. The expression of eleven SiCLE genes (SiCLE4, SiCLE10, SiCLE12, SiCLE20, SiCLE21, SiCLE25, SiCLE26, SiCLE27, SiCLE28, SiCLE30, and SiCLE38) was significantly up-regulated, while six SiCLE genes (SiCLE6, SiCLE7, SiCLE14, SiCLE16, SiCLE34, and SiCLE36) were significantly down-regulated. The results presented above indicate that the expression of SiCLE genes exhibits temporal and spatial specificity, and many members of this gene family are specifically involved in the development of young panicles and roots.

3.7. Response of SiCLEs to Plant Hormones in Foxtail Millet

The interplay between plant hormones and CLE peptides has been extensively investigated in Arabidopsis, but the understanding of foxtail millet remains limited. To explore the response of SiCLE genes to different hormones in foxtail millet, we subjected 28-day-old seedlings to five plant hormone treatments and analyzed the expression patterns of 11 SiCLE genes distributed across different groups in leaves and stems (Figure 8). Previous studies have suggested antagonistic effects between cytokinins and CLE peptides in stems and synergistic effects in roots [17]. A similar trend was observed when treating the plants with the cytokinin 6-BA. Except for SiCLE6, the expression levels of the other genes were significantly suppressed with prolonged treatment time. Interestingly, within a short duration of 6-BA treatment, these genes exhibited significant changes, with different sensitivities observed among genes in different subgroups. Group C and a subset of Group D genes (SiCLE3 and SiCLE17) showed the lowest expression levels at 2 h of treatment, followed by a gradual recovery. Group A and another subset of Group D genes (SiCLE7 and SiCLE41) maintained a lower expression level even after 2 h of treatment. In contrast, SiCLE6 from Group B demonstrated clear tissue specificity, with significantly increased expression in stems and a pattern of initial upregulation followed by downregulation in leaves. In ABA treatment, no notable tissue or functional differences were observed among genes in different subgroups. Most genes exhibited their lowest expression levels at 2 or 6 h after treatment, except for SiCLE6, SiCLE20, SiCLE31, and SiCLE33, which showed significant upregulation in leaves, potentially associated with stomatal closure [19,62]. Upon short-term application of GA3, with the exception of SiCLE24, the other SiCLE genes exhibited similar expression patterns across different tissues. Further analysis revealed an antagonistic interaction between SiCLE genes and gibberellin during the seedling stage, as their expression was significantly suppressed by GA3, except for SiCLE6, which showed increased expression. However, in leaves, the expression of SiCLE6 remained suppressed with prolonged treatment. JA plays a critical role in plant defense [63]. Treatment of foxtail millet with JA resulted in increased expression of several SiCLE genes in leaves within 30 min, including SiCLE1, SiCLE3, SiCLE6, SiCLE7, SiCLE20, and SiCLE24. Notably, SiCLE6 and SiCLE31 maintained high expression levels with prolonged treatment. Conversely, in stems, the majority of SiCLE genes exhibited downregulation, such as SiCLE1, SiCLE3, SiCLE17, SiCLE20, SiCLE24, SiCLE31, SiCLE33, and SiCLE38. Hence, there is significant tissue specificity in the response of stems and leaves to JA. Similarly to JA, treatment with IAA caused a significant decrease in the expression of SiCLE genes in stems (except for SiCLE7), with the strongest inhibitory effect observed on SiCLE33 and SiCLE36. In leaves, some SiCLE genes showed upregulation, such as SiCLE7, SiCLE20, and SiCLE41, while others exhibited downregulation, including SiCLE1, SiCLE3, SiCLE24, SiCLE33, and SiCLE38.
We conducted further analysis to examine the varying response of SiCLE genes to different hormones (Figure S3), and notable disparities in SiCLE gene expression were observed under different hormone treatments. The majority of SiCLE genes exhibited a down-regulation followed by a subsequent recovery in expression levels after treatments with 6-BA, ABA, and GA3. However, certain SiCLE genes exhibited exceptions to this trend. For instance, SiCLE31 in leaves demonstrated sustained downregulation following 6-BA treatment, while SiCLE6 in stems displayed sustained upregulation after GA3 treatment. In comparison to other hormones, JA and IAA treatments resulted in a higher number of genes exhibiting upregulation at 30 min in leaves. Moreover, the expression patterns of the same genes differed between leaves and stems under identical hormone treatments. For instance, SiCLE6 showed an upregulation in leaf expression after 30 min of 6-BA treatment, followed by a continuous downregulation as the treatment duration extended. In contrast, in stems, 6-BA treatment initially led to downregulation at 30 min and 2 h, followed by a continuous upregulation with prolonged treatment time. Hence, it can be inferred that different SiCLE genes exhibit distinct hormone responses in stems and leaves, suggesting potential functional differences within this gene family across different tissues.

4. Discussion

Intercellular communication is crucial for coordinating growth and development in eukaryotes. In plants, various signaling molecules such as plant hormones, polypeptides, and small RNAs transmit information between cells [1,64,65]. Polypeptides, in particular, play a significant role in cell signaling within plants [66]. Many small peptides act as secretory signaling molecules that are sensed by specific receptor kinases located on the cell membrane. This interaction mediates complex intercellular communication and induces responses from neighboring cells, including the regulation of meristem homeostasis [67,68]. The CLE gene family represents one of the largest known polypeptide families in plants. These genes are present throughout the evolutionary lineage of terrestrial plants [36,69,70] and generally play a conserved role in regulating plant development and physiology [71,72]. Consistent with the findings of this study, the expression patterns of SiCLE genes exhibited distinct spatiotemporal specificity (Figure 6), indicating their diverse functions in plant development. For instance, during the filling stage, SiCLE19 exhibited higher expression exclusively in roots, while SiCLE14 displayed higher expression in immature seed S5. Within the same gene groups, SiCLE genes showed similar but distinct expression patterns. In Group B, SiCLE6 and SiCLE24 demonstrated a similar expression pattern, with higher expression levels observed in panicles, germinated seeds, and one-tip-two-leaf plants while exhibiting relatively lower expression levels in other tissues. Similarly, in Group D, SiCLE32 and SiCLE41 showed relatively high expression levels in all tissues except for the panicles. On the other hand, members such as SiCLE3 and SiCLE17 in Group D displayed relatively high expression levels in immature spikelets and immature seeds. Furthermore, transcriptome sequencing was conducted to investigate the expression levels of SiCLE genes in young panicles during the formation stage and in seedling roots (Figure 7). Most SiCLE genes exhibited relatively low expression levels in young panicles. SiCLE6, SiCLE24, SiCLE25, and SiCLE34 showed higher expression levels in young panicles, with SiCLE24 showing a significant reduction in expression during the later stage 2 of panicle development. In comparison, a greater number of SiCLE genes exhibited higher expression levels in roots. SiCLE7, SiCLE22, and SiCLE36 displayed the highest expression levels in roots, and the expression of SiCLE7 and SiCLE36 was significantly down-regulated after ABA treatment. Our study indicates that the expression of SiCLE genes demonstrates temporal and spatial specificity, with many members of this gene family playing specific roles in the development of young panicles and roots. The predicted functions of SiCLE genes in foxtail millet are shown in Figure 9.
Numerous CLE loci have been identified across different species [69], and a recent study uncovered 104 CLE loci in bread wheat [28]. In the case of foxtail millet, a total of 41 SiCLE genes were identified in our study (Table 1). All the SiCLEs were distributed on nine chromosomes and divided into four groups, Group A-D (Figure 1). Five pairs of genes with duplication events, four pairs of SiCLE genes belonged to Group C, and one pair belonged to Group D (Figure 4). SiCLE genes within the same phylogenetic group exhibited similar features and displayed a diverse gene structure (Figure 2). Although the overall sequence conservation of SiCLE peptides is limited, they have a molecular weight of less than 15 kDa, and most of the polypeptides are highly basic and possess signal peptides at the N-terminal. This molecular characteristic is highly similar to CLE polypeptides found in other species [48,72]. Although these SiCLE genes exhibit diverse structural features, they all possess a conserved CLE domain located at the C-terminal, which is crucial for the functional specificity of the CLE family [73]. In species like rice, wheat, and Medicago truncatula, genes producing multiple conserved CLE domains have been found due to series repetition [36]. This structural arrangement promotes the generation of more effective CLE peptides [74]. However, such a phenomenon was not observed in foxtail millet, suggesting that the SiCLE genes in foxtail millet are highly conserved and have undergone minimal changes throughout evolution. Each CLE gene encodes a pre-propeptide consisting of a signal sequence that guides the protein through the secretory pathway, a central variable region, and a highly conserved CLE domain located at the C-terminal [71]. The full-length propeptide undergoes proteolysis [75], resulting in the production of a mature polypeptide comprising 12–14 amino acids from the CLE domain. This mature polypeptide can undergo additional post-translational modifications [14,76,77]. In our study, five conserved motifs were identified (Figure 2). Motif 1 contained the typical CLE domain and was located at the C-terminus of all 41 SiCLE peptides. Motif 2 was found in the majority of SiCLE peptides and located near the N-terminus. Motifs 3, 4, and 5 were in the central variable region, and only a few SiCLE peptides contained them. The comparative analysis of conserved CLE sequences revealed a high degree of similarity between foxtail millet and the other three species, including Arabidopsis, rice, and maize (Figure S2). The collinearity of foxtail millet and the three species indicates that CLE genes in maize and foxtail millet exhibited higher collinearity with 30 collinear gene pairs, which may be due to the fact that they are all C4 plants.
The CLE gene family in Arabidopsis is classified into type A and type B based on their functional characteristics, specifically whether they participate in the differentiation of root or stem meristem cells [78]. Type A CLE peptides include CLV3, CLE1-CLE27, CLE40, and CLE45. Based on clustering analysis, it is evident that most A-type CLE peptides are grouped in clusters C and D, indicating their important roles in the differentiation of plant root or stem meristem cells. On the other hand, B-type CLE peptides are clustered in Groups A and B, suggesting their potential involvement in plant vascular bundle development [79]. The homeostasis of plant stem cells is primarily maintained through a negative feedback pathway involving WUS and CLV3. Genes such as AtCLV3, OsFON2/FON4, OsFCP1, OsFCP2, ZmCLE7, ZmFCP1, and others have been extensively studied and demonstrated to be essential members of this pathway in their respective species [10,11,49,50]. These genes have been found to play crucial roles in extracellular signaling for meristem maintenance. In addition, the CLV3 homolog in other species also significantly influences the regulation of floral meristem size, the number of floral organs, and even the fruit size, such as Setaria viridis [80], cucumber [81], tomato [82], and B. napus [83,84]. The pre-propeptides encoded by these genes are all clustered in Group C (Figure 1), suggesting that the SiCLE genes clustered in Group C may also have similar functions in regulating floral meristem size and the number of floral organs in foxtail millet. Notably, we observed that the gene SiCLE33, located in Group C, exhibited higher expression levels during the transcriptome sequencing of the young panicle development stage. This indicates that SiCLE33 may play a vital role in CLV-WUS feedback signal transduction in foxtail millet. Furthermore, we found that genes clustered within the same group shared not only similar gene structures and relatively conserved motif compositions (Figure 2) but also displayed similar expression patterns of SiCLE genes throughout the entire life cycle of foxtail millet, as revealed by transcriptome data (Figure 6 and Figure 7). These findings are significant in comprehending the functions of SiCLE peptides in foxtail millet. Mature CLE peptides function by binding to membrane-bound receptors from various families, thereby transmitting CLE signals into cells and regulating downstream transcription factors and plant hormone signaling pathways [1]. The interaction between phytohormone signaling and CLE peptides has been extensively studied in Arabidopsis. For example, AtCLE42 can delay leaf senescence by suppressing ethylene biosynthesis [16], while AtCLE6 plays a systemic role in shoot growth under the influence of gibberellin in Arabidopsis [18]. AtCLE5 and AtCLE6 are specifically expressed at the base of developing leaves and floral organs, with their transcript levels regulated by auxin to modulate final leaf morphology [20]. In our study, we identified eight types of cis-acting elements in the upstream 2-kb region of the 41 SiCLE genes (Figure 3, Table S3). Five of these elements were related to plant hormone response, including auxin, GA3, SA, ABA, and Me-JA responsive elements. These findings suggest possible interactions between SiCLE genes and this phytohormone signaling in foxtail millet. The frequent occurrence of ABA and Me-JA cis-acting elements indicates their significant roles in regulating the expression of SiCLE genes. To investigate further, we treated foxtail millet seedlings with five different hormones during the early growth stage and observed that the expression of the majority of SiCLE genes was downregulated under different hormone treatments (Figure 8 and Figure S3). Interestingly, even under the same hormone treatment, different SiCLE genes displayed distinct expression patterns in leaves and stems, implying potential functional differences within this gene family across different tissues. The roles of AtCLE peptides have been well established in a range of developmental and physiological processes, including shoot stem cell homeostasis, root xylem development, root protophloem cell differentiation, vascular cambium activity, and stomatal formation and closure [1]. As the functions of AtCLE peptide signaling pathways have been extensively elucidated in A. thaliana, our study will contribute to exploring the diverse functions of SiCLE peptides in various developmental and physiological processes in foxtail millet.

5. Conclusions

The present study provides the first comprehensive identification and analysis of the SiCLE gene family in foxtail millet. A total of 41 SiCLE genes were distributed on nine chromosomes, with five pairs of gene duplication events, which were divided into four groups. SiCLE genes within the same phylogenetic group exhibited comparable gene structure and motif patterns. All the SiCLE proteins had the C-terminal conserved CLE domain and highly conserved positions of CLE core sequences between foxtail millet, Arabidopsis, rice, and maize. The SiCLE genes had eight types of cis-elements, and five of them were plant hormone-responsive elements, including auxin, gibberellin, salicylic acid, ABA, and Me-JA. Thirty-four SiCLE genes possessed more than three types of phytohormone-responsive elements on their promoters. CLE genes in maize and foxtail millet exhibited higher collinearity with 30 pairs of collinear CLE genes, which may be due to the fact that they are all C4 plants. The tissue-specific expression patterns and plant hormone responsiveness of SiCLE genes suggest their crucial involvement in various aspects of plant development and physiology. These findings provide a foundation for further in-depth investigations into the functional characterization and evolutionary aspects of SiCLE genes in foxtail millet.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes14112046/s1. Figure S1: Weblogo plots of the 5 conserved motifs of SiCLE peptides; Figure S2: A comparative analysis of CLE core sequences among foxtail millet, Arabidopsis, rice, and maize; Figure S3: Response of SiCLE genes to different plant hormones; Table S1: Primers used to amplify the SiCLE genes and reference gene for qPCR; Table S2: CLE genes used for phylogenetic tree in Arabidopsis, rice, and maize; Table S3: Major cis-regulatory elements in SiCLE gene promoters. Table S4: The Transcripts Per Million (TPM) values of the SiCLE genes in different tissues.

Author Contributions

S.C., H.Z. and L.L. collected the data and utilized software; X.R. and J.C. conducted the experiments and analyzed the data; X.R. and J.C. wrote the manuscript and reviewed the manuscript. 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), the Youth Fund Project on Application of Basic Research Project of Shanxi Province (201901D211363), the Scientific and Technological Innovation Programs of Shanxi Agricultural University (2018YJ32 and 2021BQ112), and the Excellent doctors come to Shanxi to reward scientific research projects (SXBYKY2022059 and SXYBKY2018035).

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 conflict of interest.

References

  1. Fletcher, J.C. Recent Advances in Arabidopsis CLE Peptide Signaling. Trends Plant Sci. 2020, 25, 1005–1016. [Google Scholar] [PubMed]
  2. Tavormina, P.; De Coninck, B.; Nikonorova, N.; De Smet, I.; Cammue, B.P. The Plant Peptidome: An Expanding Repertoire of Structural Features and Biological Functions. Plant Cell 2015, 27, 2095–2118. [Google Scholar] [PubMed]
  3. Matsubayashi, Y. Small post-translationally modified Peptide signals in Arabidopsis. Arab. Book 2011, 9, e0150. [Google Scholar]
  4. Betsuyaku, S.; Sawa, S.; Yamada, M. The function of the CLE peptides in plant development and plant-microbe interactions. Arab. Book 2011, 9, e0149. [Google Scholar]
  5. Rodriguez-Garcia, A.; Smertenko, A. CLE peptides: Emerging players in plant development and stress responses. J. Exp. Bot. 2022, 73, 3509–3524. [Google Scholar]
  6. Fiume, E.; Fletcher, J.C.; Bergmann, D.C. The CLE signaling peptides regulate diverse aspects of plant development. Curr. Opin. Plant Biol. 2021, 63, 102049. [Google Scholar]
  7. Brand, U.; Fletcher, J.C.; Hobe, M.; Meyerowitz, E.M.; Simon, R. Dependence of stem cell fate in Arabidopsis on a feedback loop regulated by CLV3 activity. Science 2000, 289, 617–619. [Google Scholar]
  8. Fletcher, J.C.; Brand, U.; Running, M.P.; Simon, R.; Meyerowitz, E.M. Signaling of cell fate decisions by CLAVATA3 in Arabidopsis shoot meristems. Science 1999, 283, 1911–1914. [Google Scholar]
  9. Schoof, H.; Lenhard, M.; Haecker, A.; Mayer, K.F.; 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]
  10. Chu, H.; Qian, Q.; Liang, W.; Yin, C.; Tan, H.; Yao, X.; Yuan, Z.; Yang, J.; Huang, H.; Luo, D.; et al. The FLORAL ORGAN NUMBER4 gene encoding a putative ortholog of Arabidopsis CLAVATA3 regulates apical meristem size in rice. Plant Physiol. 2006, 142, 1039–1052. [Google Scholar]
  11. Ohmori, Y.; Tanaka, W.; Kojima, M.; Sakakibara, H.; Hirano, H.Y. WUSCHEL-RELATED HOMEOBOX4 is involved in meristem maintenance and is negatively regulated by the CLE gene FCP1 in rice. Plant Cell 2013, 25, 229–241. [Google Scholar] [PubMed]
  12. Suzaki, T.; Sato, M.; Ashikari, M.; Miyoshi, M.; Nagato, Y.; Hirano, H.Y. 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] [PubMed]
  13. Bommert, P.; Lunde, C.; Nardmann, J.; Vollbrecht, E.; Running, M.; Jackson, D.; Hake, S.; Werr, W. Thick tassel dwarf1 encodes a putative maize ortholog of the Arabidopsis CLAVATA1 leucine-rich repeat receptor-like kinase. Development 2005, 132, 1235–1245. [Google Scholar]
  14. Xu, C.; Liberatore, K.L.; MacAlister, C.A.; Huang, Z.; Chu, Y.H.; Jiang, K.; Brooks, C.; Ogawa-Ohnishi, M.; Xiong, G.; Pauly, M.; et al. A cascade of arabinosyltransferases controls shoot meristem size in tomato. Nat. Genet. 2015, 47, 784–792. [Google Scholar] [CrossRef] [PubMed]
  15. Taguchi-Shiobara, F.; Yuan, Z.; Hake, S.; Jackson, D. The fasciated ear2 gene encodes a leucine-rich repeat receptor-like protein that regulates shoot meristem proliferation in maize. Gene Dev. 2001, 15, 2755–2766. [Google Scholar] [PubMed]
  16. Zhang, Y.; Tan, S.; Gao, Y.; Kan, C.; Wang, H.L.; Yang, Q.; Xia, X.; Ishida, T.; Sawa, S.; Guo, H.; et al. CLE42 delays leaf senescence by antagonizing ethylene pathway in Arabidopsis. New Phytol. 2022, 235, 550–562. [Google Scholar] [CrossRef] [PubMed]
  17. Cammarata, J.; Roeder, A.H.; Scanlon, M.J. Cytokinin and CLE signaling are highly intertwined developmental regulators across tissues and species. Curr. Opin. Plant Biol. 2019, 51, 96–104. [Google Scholar]
  18. Bidadi, H.; Matsuoka, K.; Sage-Ono, K.; Fukushima, J.; Pitaksaringkarn, W.; Asahina, M.; Yamaguchi, S.; Sawa, S.; Fukuda, H.; Matsubayashi, Y.; et al. CLE6 expression recovers gibberellin deficiency to promote shoot growth in Arabidopsis. Plant J. 2014, 78, 241–252. [Google Scholar]
  19. Zhang, L.; Shi, X.; Zhang, Y.; Wang, J.; Yang, J.; Ishida, T.; Jiang, W.; Han, X.; Kang, J.; Wang, X.; et al. CLE9 peptide-induced stomatal closure is mediated by abscisic acid, hydrogen peroxide, and nitric oxide in Arabidopsis thaliana. Plant Cell Environ. 2019, 42, 1033–1044. [Google Scholar]
  20. DiGennaro, P.; Grienenberger, E.; Dao, T.Q.; Jun, J.H.; Fletcher, J.C. Peptide signaling molecules CLE5 and CLE6 affect Arabidopsis leaf shape downstream of leaf patterning transcription factors and auxin. Plant Direct 2018, 2, e00103. [Google Scholar]
  21. Zhang, Y.; Yang, S.; Song, Y.; Wang, J. Genome-wide characterization, expression and functional analysis of CLV3/ESR gene family in tomato. BMC Genom. 2014, 15, 827. [Google Scholar] [CrossRef]
  22. Carbonnel, S.; Falquet, L.; Hazak, O. Deeper genomic insights into tomato CLE genes repertoire identify new active peptides. BMC Genom. 2022, 23, 756. [Google Scholar] [CrossRef] [PubMed]
  23. Wulf, K.; Sun, J.; Wang, C.; Ho-Plagaro, T.; Kwon, C.T.; Velandia, K.; Correa-Lozano, A.; Tamayo-Navarrete, M.I.; Reid, J.B.; García Garrido, J.M.; et al. The role of CLE peptides in suppression of mycorrhizal colonisation of tomato. Plant Cell Physiol. 2023, pcad124. [Google Scholar] [CrossRef] [PubMed]
  24. Qin, N.; Gao, Y.; Cheng, X.; Yang, Y.; Wu, J.; Wang, J.; Li, S.; Xing, G. Genome-wide identification of CLE gene family and their potential roles in bolting and fruit bearing in cucumber (Cucumis sativus L.). BMC Plant Biol. 2021, 21, 143. [Google Scholar] [CrossRef]
  25. Hastwell, A.H.; Gresshoff, P.M.; Ferguson, B.J. Genome-wide annotation and characterization of CLAVATA/ESR (CLE) peptide hormones of soybean (Glycine max) and common bean (Phaseolus vulgaris), and their orthologues of Arabidopsis thaliana. J. Exp. Bot. 2015, 66, 5271–5287. [Google Scholar] [CrossRef] [PubMed]
  26. Han, H.; Zhang, G.; Wu, M.; Wang, G. Identification and characterization of the Populus trichocarpa CLE family. BMC Genom. 2016, 17, 174. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, P.; Wang, Y.; Ren, F. Genome-wide identification of the CLAVATA3/EMBRYO SURROUNDING REGION (CLE) family in grape (Vitis vinifera L.). BMC Genom. 2019, 20, 553. [Google Scholar] [CrossRef] [PubMed]
  28. Li, Z.; Liu, D.; Xia, Y.; Li, Z.; Niu, N.; Ma, S.; Wang, J.; Song, Y.; Zhang, G. Identification and Functional Analysis of the CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Wheat. Int. J. Mol. Sci. 2019, 20, 4319. [Google Scholar] [CrossRef]
  29. Han, S.; Khan, M.H.U.; Yang, Y.; Zhu, K.; Li, H.; Zhu, M.; Amoo, O.; Khan, S.U.; Fan, C.; Zhou, Y. Identification and comprehensive analysis of the CLV3/ESR-related (CLE) gene family in Brassica napus L. Plant Biol. 2020, 22, 709–721. [Google Scholar] [CrossRef]
  30. Wan, K.; Lu, K.; Gao, M.; Zhao, T.; He, Y.; Yang, D.L.; Tao, X.; Xiong, G.; Guan, X. Functional analysis of the cotton CLE polypeptide signaling gene family in plant growth and development. Sci. Rep. 2021, 11, 5060. [Google Scholar] [CrossRef]
  31. Gancheva, M.S.; Lutova, L.A. Nitrogen-Activated CLV3/ESR-Related 4 (CLE4) Regulates Shoot, Root, and Stolon Growth in Potato. Plants 2023, 12, 3468. [Google Scholar] [CrossRef] [PubMed]
  32. 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 C4 model system. Nat. Plants 2020, 6, 1167–1178. [Google Scholar] [CrossRef] [PubMed]
  33. Li, X.; Hou, S.; Feng, M.; Xia, R.; Li, J.; Tang, S.; Han, Y.; Gao, J.; Wang, X. MDSi: Multi-omics Database for Setaria italica. BMC Plant Biol. 2023, 23, 223. [Google Scholar] [CrossRef]
  34. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.H.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [PubMed]
  35. Finn, R.D.; Clements, J.; Eddy, S.R. HMMER web server: Interactive sequence similarity searching. Nucleic Acids Res. 2011, 39, W29–W37. [Google Scholar] [CrossRef]
  36. Oelkers, K.; Goffard, N.; Weiller, G.F.; Gresshoff, P.M.; Mathesius, U.; Frickey, T. Bioinformatic analysis of the CLE signaling peptide family. BMC Plant Biol. 2008, 8, 1. [Google Scholar] [CrossRef]
  37. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef]
  38. Rozewicki, J.; Li, S.L.; Amada, K.M.; Standley, D.M.; Katoh, K. MAFFT-DASH: Integrated protein sequence and structural alignment. Nucleic Acids Res. 2019, 47, W5–W10. [Google Scholar] [CrossRef]
  39. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  40. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  41. Guo, A.Y.; Zhu, Q.H.; Chen, X.; Luo, J.C. GSDS: A gene structure display server. Hereditas 2007, 29, 1023–1026. [Google Scholar] [CrossRef] [PubMed]
  42. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.Y.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  43. Higo, K.; Ugawa, Y.; Iwamoto, M.; Korenaga, T. Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27, 297–300. [Google Scholar] [CrossRef]
  44. 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]
  45. He, L.; Cheng, L.; Wang, J.; Liu, J.; Cheng, J.; Yang, Z.; Cao, R.; Han, Y.; Li, H.; Zhang, B. Carotenoid Cleavage Dioxygenase 1 Catalyzes Lutein Degradation To Influence Carotenoid Accumulation and Color Development in Foxtail Millet Grains. J. Agric. Food Chem. 2022, 70, 9283–9294. [Google Scholar] [CrossRef] [PubMed]
  46. Shalabi, L.A.; Shaaban, Z.; Kasasbeh, B. Data mining: A preprocessing engine. J. Comput. Sci 2006, 2, 735–739. [Google Scholar] [CrossRef]
  47. 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]
  48. Khan, S.U.; Khan, M.H.U.; Ahmar, S.; Fan, C. Comprehensive study and multipurpose role of the CLV3/ESR-related (CLE) genes family in plant growth and development. J. Cell. Physiol. 2021, 236, 2298–2317. [Google Scholar] [CrossRef]
  49. Chu, H.; Liang, W.; Li, J.; Hong, F.; Wu, Y.; Wang, L.; Wang, J.; Wu, P.; Liu, C.; Zhang, Q.; et al. A CLE-WOX signalling module regulates root meristem maintenance and vascular tissue development in rice. J. Exp. Bot. 2013, 64, 5359–5369. [Google Scholar] [CrossRef]
  50. Je, B.I.; Xu, F.; Wu, Q.; Liu, L.; Meeley, R.; Gallagher, J.P.; Corcilius, L.; Payne, R.J.; Bartlett, M.E.; Jackson, D. The CLAVATA receptor FASCIATED EAR2 responds to distinct CLE peptides by signaling through two downstream effectors. eLife 2018, 7, e35673. [Google Scholar] [CrossRef]
  51. Jun, J.; Fiume, E.; Roeder, A.H.; Meng, L.; Sharma, V.K.; Osmont, K.S.; Baker, C.; Ha, C.M.; Meyerowitz, E.M.; Feldman, L.J.; et al. Comprehensive analysis of CLE polypeptide signaling gene expression and overexpression activity in Arabidopsis. Plant Physiol. 2010, 154, 1721–1736. [Google Scholar] [CrossRef]
  52. Kinoshita, A.; Nakamura, Y.; Sasaki, E.; Kyozuka, J.; Fukuda, H.; Sawa, S. Gain-of-function phenotypes of chemically synthetic CLAVATA3/ESR-related (CLE) peptides in Arabidopsis thaliana and Oryza sativa. Plant Cell Physiol 2007, 48, 1821–1825. [Google Scholar] [CrossRef] [PubMed]
  53. Qiang, Y.; Wu, J.; Han, H.; Wang, G. CLE peptides in vascular development. J. Integr. Plant Biol. 2013, 55, 389–394. [Google Scholar] [CrossRef] [PubMed]
  54. Lynch, M.; Richardson, A.O. The evolution of spliceosomal introns. Curr. Opin. Genet. Dev. 2002, 12, 701–710. [Google Scholar] [CrossRef] [PubMed]
  55. Imin, N.; Patel, N.; Corcilius, L.; Payne, R.J.; Djordjevic, M.A. CLE peptide tri-arabinosylation and peptide domain sequence composition are essential for SUNN-dependent autoregulation of nodulation in Medicago truncatula. New Phytol. 2018, 218, 73–80. [Google Scholar] [CrossRef] [PubMed]
  56. Nakashima, K.; Ito, Y.; Yamaguchi-Shinozaki, K. Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. Plant Physiol. 2009, 140, 88–95. [Google Scholar] [CrossRef] [PubMed]
  57. Khong, G.; Richaud, F.; Coudert, Y.; Pati, P.K.; Santi, C.; Périn, C.; Breitler, J.C.; Meynard, D.; Vinh, D.N.; Guiderdoni, E.; et al. Modulating rice stress tolerance by transcription factors. Biotechnol. Genet. Eng. Rev. 2008, 25, 381–404. [Google Scholar] [CrossRef] [PubMed]
  58. Baker, S.S.; Wilhelm, K.S.; Thomashow, M.F. The 5′-region of Arabidopsis thaliana cor15a has cis-acting elements that confer cold-, drought-and ABA-regulated gene expression. Plant Mol. Biol. 1994, 24, 701–713. [Google Scholar] [CrossRef]
  59. Casneuf, T.; De Bodt, S.; Raes, J.; Maere, S.; Van de Peer, Y. Nonrandom divergence of gene expression following gene and genome duplications in the flowering plant Arabidopsis thaliana. Genome Biol. 2006, 7, R13. [Google Scholar] [CrossRef]
  60. Song, X.F.; Hou, X.L.; Liu, C.M. CLE peptides: Critical regulators for stem cell maintenance in plants. Planta 2021, 255, 5. [Google Scholar] [CrossRef]
  61. 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]
  62. Yoshida, T.; Fernie, A.R.; Shinozaki, K.; Takahashi, F. Long-distance stress and developmental signals associated with abscisic acid signaling in environmental responses. Plant J. 2021, 105, 477–488. [Google Scholar] [CrossRef]
  63. Li, C.; Xu, M.; Cai, X.; Han, Z.; Si, J.; Chen, D. Jasmonate Signaling Pathway Modulates Plant Defense, Growth, and Their Trade-Offs. Int. J. Mol. Sci. 2022, 23, 3945. [Google Scholar] [CrossRef] [PubMed]
  64. Wang, Y.; Jiao, Y. Cell signaling in the shoot apical meristem. Plant Physiol. 2023, 193, 70–82. [Google Scholar] [CrossRef] [PubMed]
  65. Awale, P.; McSteen, P. Hormonal regulation of inflorescence and intercalary meristems in grasses. Curr. Opin. Plant Biol. 2023, 102451. [Google Scholar] [CrossRef] [PubMed]
  66. Selby, R.; Jones, D.S. Complex peptide hormone signaling in plant stem cells. Curr. Opin. Plant Biol. 2023, 75, 102442. [Google Scholar] [CrossRef] [PubMed]
  67. Shiu, S.H.; Bleecker, A.B. Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. Proc. Natl. Acad. Sci. USA 2001, 98, 10763–10768. [Google Scholar] [CrossRef] [PubMed]
  68. John, A.; Smith, E.S.; Jones, D.S.; Soyars, C.L.; Nimchuk, Z.L. A network of CLAVATA receptors buffers auxin-dependent meristem maintenance. Nat. Plants 2023, 9, 1306–1317. [Google Scholar] [CrossRef]
  69. Goad, D.M.; Zhu, C.; Kellogg, E.A. Comprehensive identification and clustering of CLV3/ESR-related (CLE) genes in plants finds groups with potentially shared function. New Phytol. 2016, 216, 605–616. [Google Scholar] [CrossRef]
  70. Hirakawa, Y.; Uchida, N.; Yamaguchi, Y.L.; Tabata, R.; Ishida, S.; Ishizaki, K.; Nishihama, R.; Kohchi, T.; Sawa, S.; Bowman, J.L. Control of proliferation in the haploid meristem by CLE peptide signaling in Marchantia polymorpha. PLoS Genet. 2019, 15, e1007997. [Google Scholar] [CrossRef]
  71. Somssich, M.; Je, B.I.; Simon, R.; Jackson, D. CLAVATA–WUSCHEL signalling in the shoot meristem. Development 2016, 143, 3238–3248. [Google Scholar] [CrossRef] [PubMed]
  72. Cock, J.M.; McCormick, S. A large family of genes that share homology with CLAVATA3. Plant Physiol. 2001, 126, 939–942. [Google Scholar] [CrossRef] [PubMed]
  73. Meng, L.; Ruth, K.C.; Fletcher, J.C.; Feldman, L. The roles of different CLE domains in Arabidopsis CLE polypeptide activity and functional specificity. Mol. Plant 2010, 3, 760–772. [Google Scholar] [CrossRef]
  74. Miyawaki, K.; Tabata, R.; Sawa, S. Evolutionarily conserved CLE peptide signaling in plant development, symbiosis, and parasitism. Curr. Opin. Plant Biol. 2013, 16, 598–606. [Google Scholar] [CrossRef]
  75. Ni, J.; Clark, S.E. Evidence for functional conservation, sufficiency, and proteolytic processing of the CLAVATA3 CLE domain. Plant Physiol. 2006, 140, 726–733. [Google Scholar] [CrossRef] [PubMed]
  76. Kondo, T.; Sawa, S.; Kinoshita, A.; Mizuno, S.; Kakimoto, T.; Fukuda, H.; Sakagami, Y. A plant peptide encoded by CLV3 identified by in situ MALDI-TOF MS analysis. Science 2006, 313, 845–848. [Google Scholar] [CrossRef]
  77. Ohyama, K.; Shinohara, H.; Ogawa-Ohnishi, M.; Matsubayashi, Y. A glycopeptide regulating stem cell fate in Arabidopsis thaliana. Nat. Chem. Biol. 2009, 5, 578–580. [Google Scholar] [CrossRef]
  78. Whitford, R.; Fernandez, A.; De Groodt, R.; Ortega, E.; Hilson, P. Plant CLE peptides from two distinct functional classes synergistically induce division of vascular cells. Proc. Natl. Acad. Sci. USA 2008, 105, 18625–18630. [Google Scholar] [CrossRef]
  79. Yamaguchi, Y.L.; Ishida, T.; Sawa, S. CLE peptides and their signaling pathways in plant development. J. Exp. Bot. 2016, 67, 4813–4826. [Google Scholar] [CrossRef]
  80. Li, S.; Pan, Y.; Wen, C.; Li, Y.; Liu, X.; Zhang, X.; Behera, T.K.; Xing, G.; Weng, Y. Integrated analysis in bi-parental and natural populations reveals CsCLAVATA3 (CsCLV3) underlying carpel number variations in cucumber. Theor. Appl. Genet. 2016, 129, 1007–1022. [Google Scholar] [CrossRef]
  81. Zhu, C.; Liu, L.; Crowell, O.; Zhao, H.; Brutnell, T.P.; Jackson, D.; Kellogg, E.A. The CLV3 Homolog in Setaria viridis Selectively Controls Inflorescence Meristem Size. Front. Plant Sci. 2021, 12, 636749. [Google Scholar] [CrossRef] [PubMed]
  82. Yuste-Lisbona, F.J.; Fernández-Lozano, A.; Pineda, B.; Bretones, S.; Ortíz-Atienza, A.; García-Sogo, B.; Müller, N.A.; Angosto, T.; Capel, J.; Moreno, V.; et al. ENO regulates tomato fruit size through the floral meristem development network. Proc. Natl. Acad. Sci. USA 2020, 117, 8187–8195. [Google Scholar] [CrossRef] [PubMed]
  83. Fan, C.; Wu, Y.; Yang, Q.; Yang, Y.; Meng, Q.; Zhang, K.; Li, J.; Wang, J.; Zhou, Y. A novel single-nucleotide mutation in a CLAVATA3 gene homolog controls a multilocular silique trait in Brassica rapa L. Mol. Plant 2014, 7, 1788–1792. [Google Scholar] [CrossRef] [PubMed]
  84. Yang, Y.; Zhu, K.; Li, H.; Han, S.; Meng, Q.; Khan, S.U.; Fan, C.; Xie, K.; Zhou, Y. Precise editing of CLAVATA genes in Brassica napus L. regulates multilocular silique development. Plant Biotechnol. J. 2018, 16, 1322–1335. [Google Scholar] [CrossRef]
Genes 14 02046 g001
Figure 1. Phylogenetic tree of CLE peptides from foxtail millet, Arabidopsis, rice, and maize. The phylogenetic tree was constructed with the neighbor-joining (NJ) method in MEGA 7.0 software and was divided into four subgroups.
Figure 1. Phylogenetic tree of CLE peptides from foxtail millet, Arabidopsis, rice, and maize. The phylogenetic tree was constructed with the neighbor-joining (NJ) method in MEGA 7.0 software and was divided into four subgroups.
Genes 14 02046 g001
Genes 14 02046 g002
Figure 2. Gene structure and conserved motif analysis of SiCLE peptides. (A) Phylogeny tree. (B) Gene structure. Coding sequence (CDS) and untranslated region (UTR) are represented by different colored boxes, and introns are represented by lines. (C) Motif patterns. Conserved motifs in the SiCLE peptides are represented by different colored boxes. Weblogo plots of the five conserved motifs are shown in Figure S1.
Figure 2. Gene structure and conserved motif analysis of SiCLE peptides. (A) Phylogeny tree. (B) Gene structure. Coding sequence (CDS) and untranslated region (UTR) are represented by different colored boxes, and introns are represented by lines. (C) Motif patterns. Conserved motifs in the SiCLE peptides are represented by different colored boxes. Weblogo plots of the five conserved motifs are shown in Figure S1.
Genes 14 02046 g002
Genes 14 02046 g003
Figure 3. Cis-regulatory elements in SiCLE gene promoters. The elements are displayed in differently colored boxes.
Figure 3. Cis-regulatory elements in SiCLE gene promoters. The elements are displayed in differently colored boxes.
Genes 14 02046 g003
Genes 14 02046 g004
Figure 4. Chromosomal distribution and gene duplication of the SiCLE genes. Two SiCLE genes of the same segmental duplicated gene pair are labeled with blue lines.
Figure 4. Chromosomal distribution and gene duplication of the SiCLE genes. Two SiCLE genes of the same segmental duplicated gene pair are labeled with blue lines.
Genes 14 02046 g004
Genes 14 02046 g005
Figure 5. Synteny analyses of the SiCLE genes between S. italica and three representative plant species. Gray lines on the background indicate the collinear blocks within S. italica and other plant genomes; red lines highlight the syntenic S. italica SiCLE gene pairs.
Figure 5. Synteny analyses of the SiCLE genes between S. italica and three representative plant species. Gray lines on the background indicate the collinear blocks within S. italica and other plant genomes; red lines highlight the syntenic S. italica SiCLE gene pairs.
Genes 14 02046 g005
Genes 14 02046 g006
Figure 6. Expression pattern analysis of SiCLE genes in different tissues. The heatmap of the expression profiles of the SiCLE gene in different developmental stages is represented by normalized values using RNA-seq data, with the color from blue to red indicating the expression levels from high to low.
Figure 6. Expression pattern analysis of SiCLE genes in different tissues. The heatmap of the expression profiles of the SiCLE gene in different developmental stages is represented by normalized values using RNA-seq data, with the color from blue to red indicating the expression levels from high to low.
Genes 14 02046 g006
Genes 14 02046 g007
Figure 7. Expression patterns of SiCLE genes in young panicles and roots of foxtail millet. (A) Expression patterns of SiCLE genes in young panicles at different formation stages: branch meristems specified (stage 1, approximately 1.0–1.5 mm) and later stage with clearly formed branch meristems (stage 2, approximately 2.5–3.0 mm). (B) Expression patterns of SiCLE genes in the roots of 9-day-old seedlings with ABA treatment (2 μM) and without ABA as the control (CK). The bars represent the variation between replicates in the transcriptome sequencing of young panicle samples (A) and root samples (B). Statistically significant differences between the two stages of the panicle (A) and the root with and without ABA treated (B) were determined according to t-test (* p < 0.05, ** p < 0.01, *** p < 0.001).
Figure 7. Expression patterns of SiCLE genes in young panicles and roots of foxtail millet. (A) Expression patterns of SiCLE genes in young panicles at different formation stages: branch meristems specified (stage 1, approximately 1.0–1.5 mm) and later stage with clearly formed branch meristems (stage 2, approximately 2.5–3.0 mm). (B) Expression patterns of SiCLE genes in the roots of 9-day-old seedlings with ABA treatment (2 μM) and without ABA as the control (CK). The bars represent the variation between replicates in the transcriptome sequencing of young panicle samples (A) and root samples (B). Statistically significant differences between the two stages of the panicle (A) and the root with and without ABA treated (B) were determined according to t-test (* p < 0.05, ** p < 0.01, *** p < 0.001).
Genes 14 02046 g007
Genes 14 02046 g008
Figure 8. Expression patterns of SiCLE genes response to plant hormones in foxtail millet. Seedlings that were 28 days old were treated with plant hormones, including 6-BA (100 µM) (A), ABA (100 µM) (B), GA3 (100 µM) (C), MeJA (100 µM) (D), and IAA (100 µM) (E). Samples of leaves and stems were collected at 0 h, 0.5 h, 2 h, 6 h, and 12 h after the treatments. The gene expression levels at different time intervals were detected using the qPCR method. The lowercase letters in the heatmap represent significant differences between different time intervals after treatment. The color from red to blue indicates the expression levels from high to low in the heatmap.
Figure 8. Expression patterns of SiCLE genes response to plant hormones in foxtail millet. Seedlings that were 28 days old were treated with plant hormones, including 6-BA (100 µM) (A), ABA (100 µM) (B), GA3 (100 µM) (C), MeJA (100 µM) (D), and IAA (100 µM) (E). Samples of leaves and stems were collected at 0 h, 0.5 h, 2 h, 6 h, and 12 h after the treatments. The gene expression levels at different time intervals were detected using the qPCR method. The lowercase letters in the heatmap represent significant differences between different time intervals after treatment. The color from red to blue indicates the expression levels from high to low in the heatmap.
Genes 14 02046 g008
Genes 14 02046 g009
Figure 9. The schematic diagram illustrates the predicted functions of SiCLE genes in foxtail millet. The diagram is based on the phylogenetic relationship and expression patterns of SiCLE genes, as well as the knowledge from reported CLE genes in other species. The font color in the diagram corresponds to the background color of different groups in the phylogenetic tree.
Figure 9. The schematic diagram illustrates the predicted functions of SiCLE genes in foxtail millet. The diagram is based on the phylogenetic relationship and expression patterns of SiCLE genes, as well as the knowledge from reported CLE genes in other species. The font color in the diagram corresponds to the background color of different groups in the phylogenetic tree.
Genes 14 02046 g009
Table 1. The properties of SiCLE genes.
Table 1. The properties of SiCLE genes.
Gene IDGene NameChromosome LocationSignal Peptide PredictionNumber of Amino AcidsMolecular Weight
(kD)		Isoelectric PoinGRAVY
Si1g02740SiCLE1Chr1:1701179–1702028+838.675.450.071
Si1g09260SiCLE2Chr1:8075455–8076779+11010.9810.87−0.06
Si1g15040SiCLE3Chr1:16870079–16870441+12013.4211.20−0.497
Si1g21450SiCLE4Chr1:28806486–28807552+919.6610.39−0.204
Si1g25600SiCLE5Chr1:32729889–32731715+858.689.490.232
Si1g36220SiCLE6Chr1:41284255–41285566−9810.4611.83−0.427
Si1g37440SiCLE7Chr1:42036658–42038549+859.709.99−0.359
Si2g07840SiCLE8Chr2:6541423–6542683−10410.9211.49−0.34
Si2g37520SiCLE9Chr2:45004975–45005241+889.126.03−0.016
Si3g03240SiCLE10Chr3:1462866–1464336+929.5011.01−0.049
Si3g10460SiCLE11Chr3:6732083–6733514+869.2611.83−0.2
Si3g14760SiCLE12Chr3:10562885–10563172+959.648.250.082
Si3g17810SiCLE13Chr3:13430959–13431192+778.029.180.112
Si3g17830SiCLE14Chr3:13472475–13472753+929.275.680
Si3g18030SiCLE15Chr3:13627896–13628249+11712.6611.25−0.368
Si3g18460SiCLE16Chr3:13936956–13938953−9910.2411.34−0.269
Si3g22990SiCLE17Chr3:19091857–19092216+11913.1711.50−0.403
Si4g07460SiCLE18Chr4:5221231–5221955+838.547.780.094
Si4g16140SiCLE19Chr4:25203756–25204438+11912.229.98−0.317
Si4g16170SiCLE20Chr4:25243962–25244464+828.565.980.126
Si4g23620SiCLE21Chr4:34831659–34833580+12712.8611.03−0.366
Si5g09450SiCLE22Chr5:7782941–7783884+939.7211.860.103
Si5g12400SiCLE23Chr5:10508454–10508813+11912.5011.89−0.3
Si5g16690SiCLE24Chr5:15497255–15497578+10710.9511.910.068
Si5g18600SiCLE25Chr5:21570492–21570803+10311.0012.16−0.194
Si5g32640SiCLE26Chr5:37962048–37962314+888.9311.65−0.086
Si5g32660SiCLE27Chr5:37979626–37979868+807.936.570.173
Si5g32670SiCLE28Chr5:37981960–37982229+899.149.10−0.128
Si5g33310SiCLE29Chr5:38510416–38510766+11612.879.75−0.422
Si6g16350SiCLE30Chr6:27510936–27511193+859.6111.42−0.351
Si7g13900SiCLE31Chr7:21729033–21729898+11311.4011.66−0.113
Si7g19100SiCLE32Chr7:25632473–25632760+9510.0710.74−0.065
Si8g17400SiCLE33Chr8:31420418–31420996+10010.338.97−0.068
Si9g05550SiCLE34Chr9:3007939–3009361−9910.2811.34−0.362
Si9g13150SiCLE35Chr9:8512203–8514840+10310.4610.24−0.087
Si9g19660SiCLE36Chr9:14748515–14748763+829.106.56−0.402
Si9g19670SiCLE37Chr9:14751454–14751981+10111.0910.220.02
Si9g38030SiCLE38Chr9:45088642–45088950+10210.4411.110.047
Si9g38660SiCLE39Chr9:45679290–45679866+9410.389.81−0.068
Si9g44430SiCLE40Chr9:50178960–50179208−828.5311.610.263
Si9g55410SiCLE41Chr9:58385979–58386266+9510.8410.00−0.287
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

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. https://doi.org/10.3390/genes14112046

AMA Style

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(11):2046. https://doi.org/10.3390/genes14112046

Chicago/Turabian Style

Ren, Xuemei, Jinjie Chen, Shuwan Chen, Hui Zhang, and Li Li. 2023. "Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica L.)" Genes 14, no. 11: 2046. https://doi.org/10.3390/genes14112046

APA Style

Ren, X., Chen, J., Chen, S., Zhang, H., & Li, L. (2023). Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica L.). Genes, 14(11), 2046. https://doi.org/10.3390/genes14112046

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