
==== Front
Int J Genomics
Int J Genomics
IJG
International Journal of Genomics
2314-436X
2314-4378
Wiley

10.1155/2024/2708223
Research Article
Genome-Wide Characterization and Analysis of the SPL Gene Family in Eucalyptus grandis
https://orcid.org/0009-0004-4057-1825
An Lijun 1
Ma Jiasi 1
Fan Chunjie 2
Li Huiling 1
https://orcid.org/0000-0001-8322-8833
Wu Aimin wuaimin@scau.edu.cn
1
1 Guangdong Key Laboratory for Innovative Development and Utilization of Forest Plant Germplasm College of Forestry and Landscape Architectures South China Agricultural University, Guangzhou 510642, China
2 State Key Laboratory of Tree Genetics and Breeding Key Laboratory of State Forestry and Grassland Administration on Tropical Forestry Research Institute of Tropical Forestry Chinese Academy of Forestry, Guangzhou 510520, China
Academic Editor: Antonio Ferrante

2024
11 9 2024
2024 270822331 3 2024
13 6 2024
13 8 2024
Copyright © 2024 Lijun An et al.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
SQUAMOSA promoter-binding protein-like (SPL) gene family, a group of plant-specific transcription factors, played crucial roles in regulating plant growth, development, signal transduction, and stress response. This study focuses on the SPL gene family in the fast-growing Eucalyptus grandis, employing bioinformatics approaches to identify and analyze the gene physiochemical characteristics, conserved domains, structural composition, chromosomal distribution, phylogenetic relationships, cis-acting elements, and their expression patterns in various tissues and stress treatments. Twenty-three SPL genes were identified in E. grandis, which uneven distributed across seven chromosomes and classified into five groups. Prediction of cis-acting elements revealed that these genes might be related to light, hormone, and stress responses. Furthermore, EgSPL9 and EgSPL23, mainly expressed in the stem apex and lateral branches, seem to be involved in hormone stress resistance. Our study provides insights into the potential functions of the EgSPL genes in plant growth, stress response, and hormone transduction, offering valuable perspectives for subsequent research into their biological roles.

Keywords

bioinformatics analysis
E. grandis
gene expression
SPL gene family
Natural Science Foundation of Guangdong Province2023A1515010379 China Postdoctoral Science Foundation2022M721198 Guangzhou Science and Technology Program key projects202201011282
==== Body
pmc1. Introduction

The SQUAMOSA promoter-binding protein-like (SPL) contained a highly conserved squamous promoter-binding protein (SBP) domain consisting of 76 amino acid residues and two specific zinc finger motifs (composed of 8 cysteine (Cys) or histidine (His) residues) [1]. As the plant-specific transcription factor family, it was first discovered in Arabidopsis thaliana [2] and was able to specifically recognize the SQUAMOSA (SQUA) promoter. Consequently, the SPL transcription factors play important roles in multiple stages of plant growth and development including initiation of flowering [3, 4], stem and leaf development [5, 6], flower and fruit formation [7], male sterility, and response to abiotic stress [8–10].

In plant, SPL genes showed diverse and critical roles for various developmental and physiological processes. For instance, AtSPL2, AtSPL10, AtSPL11, and AtSPL13 are involved in regulating transition from vegetative growth to reproduction growth [11]. AtSPL3 and its homologous genes AtSPL4 and AtSPL5 are involved in flowering [12, 13]. AtSPL9 and AtSPL15 mainly regulated leaf growth and the initiation of flowering [14] while they also were explored involving in the transition from juvenile to adult developmental stages [15–17]. On the contrary, AtSPL14 seems to have a function in delaying phase transition [18]. Moreover, comprehensive studies have revealed that most SPL family members harbor the response target element of miR156, a type of noncoding RNA predominantly involved in regulating gene expression at the posttranscriptional level [19]. Specifically, in A. thaliana, the mutation of the miR156 binding site impedes its negative regulation of AtSPL13, significantly delaying the appearance of the leaf primordium and therefore delaying the growth of true leaves during germination [20]. AtSPL8 cooperated with other SPL genes with miR156 target sequence, which further regulate early anther development and pistil differentiation, as well as sepal epidermal trichome formation, root growth, microsporogenesis, and megasporogenesis [21–24].

Additionally, SPLs also played an important role in stress resistance. AtSPL7 regulates copper (Cu) homeostasis and cadmium (Cd) stress [25, 26]. The AtSPL10 binds to and activates key defense genes, thereby enhancing resistance in A. thaliana leaves [27]. Under conditions of NaCl stress, the expression of the BpSPL9 gene in birch roots and leaves is upregulated, facilitating the scavenging of reactive oxygen species [9]. Tea CsSPL1 and CsSPL12 were inhibited by abscisic acid (ABA), GA, and salt treatment and induced by ethylene [28].

Since the functional importance of SPL genes, researchers have developed genome-wide identification of SPL genes in numerous plant species. There was 17 SPL genes identified in A. thaliana [29], 28 SPL genes in Populus trichocarpa [30], 21 SPL genes in hybrid poplar [31], and 18 SPL genes in Betula luminifera [32]. Although SPL genes in the same evolutionary branch have certain similarities, many SPL genes from the same branch might have different functions in different plant species as well.

Eucalyptus grandis, distinguished by its rapid growth rate and widespread geographical distribution, is highly valued for its superior wood and fiber qualities. Completeness of genome sequence of E. grandis provided an opportunity to elucidate the molecular mechanisms governing its growth, development, and response to environmental stressors. Currently, the biological role and precise molecular mechanism of the SPL genes were elusive. Comprehensive annotation of EgSPLs is the first step towards elucidation of the underlying molecular mechanisms. This study conducted a genome-wide analysis of the E. grandis SPLs, including the systematic identification of conserved structure, phylogeny, chromosomal location, collinearity analysis, promoter homeopathic element analysis, and its expression in various tissues and different stress conditions. Our study lays the foundation for exploring functions of SPLs in transcriptional regulation of E. grandis.

2. Materials and Methods

2.1. Identification and Physicochemical Properties

We acquired genome-wide data of E. grandis from E. grandis Genome Database (https://phytozome-next.jgi.doe.gov/info/Egrandis_v2_0), and downloaded 17 A. thaliana AtSPL gene sequences from the TAIR database (https://www.Arabidopsisthaliana.org/browse/genefamily/sbp_box_genefamily.jsp). Through blast comparison, initial screening is conducted to obtain EgSPL family–related genes. Utilizing the NCBI CDD online tool to predict SPL protein domains, excluding atypical domains and incomplete redundant sequences, the EgSPL family members are identified. Analysis of the physicochemical properties of EgSPL proteins is performed using the online software ExPASy ProtParam (https://web.expasy.org/protparam/), the subcellular localization of EgSPL genes was predicted using the Plant-mPLoc platform (http://www.csbio.sjtu.edu.cn/cgi-bin/PlantmPLoc.cgi), and potential miR156 target sites within the EgSPL genes were identified through predictions made on the psRNA target online platform (http://plantgm.noble.org/psrnatarget/).

2.2. Phylogenetic Tree Analysis

To elucidate the evolutionary relationships among SPL genes across diverse plant groups, including algae, ferns, dicots, and woody plants, the protein sequences of Physcomitrium patens, Selaginella moellendorffii, A. thaliana, and E. grandis were selected for comparison analysis. The resulting data were analyzed using the maximum likelihood (ML) method implemented in the TBtools software [33], with bootstrap analysis set to 2000 replicates, to construct a comprehensive phylogenetic tree. Lastly, the phylogenetic tree file was input into the Evolview website (https://www.evolgenius.info/evolview/) for annotation and optimization.

2.3. Conservative Domain and Genetic Structure Analysis

The conserved motifs of the EgSPL proteins were analyzed using the MEME online software (http://memesuite.org/tools/meme), with the number of identified motifs set to 10. The results of this analysis were saved in an XML file format, which was subsequently visualized using TBtools [33]. Additionally, the gene structure of the EgSPL gene family was elucidated through analysis in TBtools software, utilizing the gff3 file of the E. grandis genome.

2.4. Chromosome Positioning and Collinearity Analysis

The gff3 file in the E. grandis genome was used in TBtools software to obtain the staining and positioning information of the EgSPL family members and draw the distribution map of the genes on the chromosome. Genome sequences of E. grandis, A. thaliana, and P. alba were downloaded from the Phytozome (https://phytozome-next.jgi.doe.gov/). Genome collinearity analysis was conducted using the One Step MCScanX function in TBtools software, focusing on the species E. grandis and its comparative analyses with A. thaliana and P. alba. The parameter “Num of BlastHits” was set to 4 for the interspecific comparisons, and the E-value was set to 1e−10.

2.5. Cis-Element Analysis of EgSPLs Promoters

The upstream 2000 bp sequence of SPL genes as promoter was extracted by using TBtools, and PlantCARE (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/) was applied to predict the cis-acting elements. The predicted cis-elements were processed and visualized by using Excel and TBtools, respectively.

2.6. Plant Material and Gene Expression Pattern Analysis and Growth Conditions

This tissue from E. grandis trees was from 6-month-old saplings providing roots, leaves, xylem, phloem, and stem internodes, while mature trees (3 and 6 years old) supplied organs such as flowers, xylem, phloem, and cambium. Hormonal treatments with salicylic acid (SA) or methyl jasmonate (MeJA) and salinity stress were applied to 2-month-old shoots. Samples were collected at 0, 1, 6, 24, and 168 h posttreatment. All experiments were performed triple replicates at least, ensuring a robust dataset for subsequent analysis. The detail transcriptomic data for E. grandis, detailing gene expression across tissues and following treatments with SA, MeJA, and salt [34], were retrieved from a dedicated eucalyptus database (https://bar.utoronto.ca/%7Easher/eplant_eucalyptus/). Subsequent to data extraction, heat maps were generated using TBtools to visually delineate and interpret the expression profiles, offering a comprehensive overview of the gene activity under various physiological and stress conditions.

2.7. Prediction of Protein Structure

The three-dimensional structure of EgSPL protein was predicted through the SWISS-MODEL website (https://swissmodel.expasy.org/). Upload the amino acid sequence to the website, select the model with the highest matching degree, among which the similarity of the selected models is above 40%, and output the file for analysis.

3. Results

3.1. Identification and Physicochemical Properties

A total of 23 members of the SPLs were identified by removing proteins without conserved domains from the genome of E. grandis. According to the result from two-way blast comparison, these genes were named EgSPL1-23 and the details are listed in Table 1. And their physiochemical properties are also shown in Table 1. The lengths of EgSPL proteins are 147–1078 aa (no. of amino acids), with a relative molecular weight (MW) ranging from 19.95 to 119.29 kDa (molecular mass). The isoelectric points (pI) of the EgSPL proteins range from 5.79 to 9.50, with EgSPL1, EgSPL4, EgSPL9, EgSPL16, and EgSPL22 having values less than 7, classifying them as acidic proteins. The instability index of EgSPL protein ranged from 50.07 to 88.37, of which EgSPL2 was 88.37, was most instable protein. In addition, the grand average of hydropathicity of all EgSPL protein was negative, which ranged from −1.44 to −0.41. Among the 23 EgSPL genes, 16 genes contain target sites for miR156. However, there are still seven genes including EgSPL3, EgSPL9, EgSPL10, EgSPL12, EgSPL20, EgSPL22, and EgSPL23, which were not identified targeting sites of miR156. Most of the EgSPL genes were located in the nucleus, and only a few genes, including EgSPL1, EgSPL4, EgSPL14, EgSPL16, and EgSPL22, appeared in both the nucleus and the cytoplasm (Table 1).

3.2. Phylogenetic Analysis

To elucidate the evolutionary relationships of SPL genes in plants, we constructed a phylogenetic tree using SPL protein sequences from E. grandis, A. thaliana, S. moellendorffii, and P. patens, which were divided into five groups (Figure 1). Notably, no EgSPL genes existed in Group I, raising questions about potential gene loss in the evolutionary trajectory of the E. grandis genome. In Group III, EgSPL genes form a close cluster with SPL genes from S. moellendorffii and P. patens, suggesting a shared common ancestor during evolutionary processes. In Group IV, a pronounced clustering of SPL genes from both A. thaliana and E. grandis is observed, suggesting a more recent common evolutionary lineage or a closer phylogenetic relationship between these genes from the two species. Specifically, EgSPL3 and EgSPL5 align with AtSPL9 and AtSPL15 from A. thaliana, respectively. The convergence of EgSPL22, AtSPL7, PpSPL9, and SmSPL9 on the same branch is particularly intriguing and warrants further investigation to elucidate any potential functional similarities or shared evolutionary trajectories. Further analysis revealed that the number of SPL family genes has increased and gene duplications have been significantly enhanced from lower to higher plants, reflecting the evolutionary success of flowering plants in adapting to diverse and fluctuating environmental conditions [35].

3.3. Gene Structure and Conserved Motif Analysis

Evolutionary analysis of EgSPL genes suggests that genes from the same evolutionary branch often have similar conserved domains, regardless of the closeness of their relationship (Figure 2(a)). Conserved domain analysis of the EgSPL proteins indicates that all members include Motifs 1, 2, and 3, suggesting that these are fundamental to SPLs function (Figure 2(b)). Specific structural domains were also identified, such as Motifs 4 and 8 existed in EgSPL10, EgSPL12, and EgSPL20 suggesting the possibility of their functional attributes and performing similar biological roles. Gene structure analysis illustrates variation in exon–intron configurations among the EgSPL family members, with the number of exons ranging from 2 to 11 (Figure 2(c)). Furthermore, EgSPL9 has an 11-exon structure, the most among all EgSPL genes. Moreover, genes that are phylogenetically closed often show similar intron–exon structures.

3.4. Chromosomal Locations and Synteny Analysis

Based on the E. grandis genome data, chromosome position information was obtained. Twenty-three EgSPL genes are unevenly distributed on seven chromosomes (Figure 3). Three sets of genes were found in clusters on Chromosome 11 (EgSPL5/6/11, EgSPL1/16, and EgSPL7/17/19), and a set of genes were found in clusters on Chromosome 05 (EgSPL10/12/20) (Figure 3). This suggests that clustering may enhance their functional capabilities.

The synteny analysis conducted on A. thaliana, E. grandis, and P. alba revealed that the EgSPL genes exhibited a conserved collinear relationship with homolog genes within A. thaliana and P. alba. Significantly, the assortments of collinear gene pairs vary among A. thaliana, E. grandis, and P. alba, underscoring the intricacy of their phylogenetic interrelations (Figure 4). There are 10 collinear gene pairs identified between A. thaliana and E. grandis and 18 pairs between E. grandis and P. alba, speculating that the disparity and presence of gene clusters on the chromosomes probably caused by gene duplication and rearrangement events throughout evolutionary history. Furthermore, the closer genetic affinity between the SPLs of E. grandis and P. alba indicates a higher degree of conservation in woody plants.

3.5. Cis-Element Analysis of EgSPLs Promoters

A total of 27 cis-elements were obtained by PlantCARE from upstream 2 kb sequence of EgSPLs, and a large number of cis elements were related to light (9/27); stress (6/27) including low temperature, hypoxia, drought, and defense; growth and development (4/27); and hormone response (7/27) (Figure 5). Hormone-responsive components including ABA, MeJA, auxin, GA, and SA were also identified in promoter of EgSPL genes, which was in accordance with important roles in flowering [3], defense response [8], and stress homeostasis of SPL genes in previous studies [36]. These cis-elements exist in EgSPLs, indicating that the EgSPL genes in E. grandis may play the similar roles.

3.6. Analysis of Tissue Expression Patterns of SPLs in E. grandis at Different Stages

The gene expression pattern often implied their function. Based on the transcriptome data of E. grandis, differential expression patterns of EgSPLs were observed across various tissues. Meanwhile, 23 EgSPL genes displayed low expression in adult leaves and young roots, while showing high expression in the stem apex and mature xylem and phloem (Figure 6(a)). For instance, EgSPL5, EgSPL6, and EgSPL11 showed high expression levels in young leaves, whereas EgSPL9 and EgSPL22 were predominantly expressed in xylem. Additionally, EgSPL23 demonstrated significant expression in both xylem and flowers, alongside EgSPL9 and EgSPL21, suggesting a potential role in the flowering process (Figure 6(b)). In 6-year-old E. grandis, EgSPL9, EgSPL22, and EgSPL23 were highly expressed in both xylem and phloem. Interestingly, with increasing age, there was an upsurge in the expression of EgSPL3, EgSPL10, EgSPL12, EgSPL15, and EgSPL20 genes in these tissues (Figure 6(c)). This pattern suggests that EgSPL9, EgSPL22, and EgSPL23 may play a substantial role in the overall growth and stem development of E. grandis.

3.7. Analysis of Induced Expression Patterns of SPLs in E. grandis Under Different Stress

To further explore the expression patterns of different EgSPL family genes under MeJA, SA, and salt stress, transcriptome data were obtained. Under various stress conditions, the EgSPL genes exhibit a spectrum of expression profiles, underscoring the multifaceted role in the plant's stress response mechanism. Specifically, under SA treatment, elevated expression levels were recorded for EgSPL4, EgSPL9, and EgSPL23, whereas the rest of EgSPLs did not show significant changes. Concurrently, EgSPL13 showed a short-term increase, whereas EgSPL7, EgSPL17, EgSPL18, and EgSPL19 demonstrated a short-term decrease (Figure 7(a)). Under MeJA treatment, EgSPL13 exhibited a decrease expression, contrasting with the general increase observed in other EgSPL genes. This differential response suggests intricate regulatory mechanisms in MeJA treatment, potentially reflecting the diverse roles of these genes in MeJA-mediated pathway in E. grandis (Figure 7(b)). Under salt stress conditions, a majority of EgSPL genes exhibited a decrease expression, either in the short or long term. Conversely, EgSPL8, EgSPL14, EgSPL22, and EgSPL23 showed a short-term increase expression, which suggests that these genes could play crucial roles in plant adaptation salt stress, highlighting their potential importance in stress response mechanisms (Figure 7(c)). Differently, EgSPL9 consistently maintained higher expression levels.

3.8. Prediction of Protein Structure

In order to study the function of the protein, we predicted the structure of EgSPLs from the SWISS-MODEL website. The results showed that these proteins had certain similarities overall, suggesting that they may have similar functions or originate from the same protein family (Figure 8). In addition, the tertiary structures of EgSPL4, EgSPL13, and EgSPL23 contain multiple functional domains and show high conservation. In comparison, EgSPL1-3 and EgSPL5-12 are mainly random coils, containing a small amount of α-helix and β-sheet structures, and the overall structure is relatively simple. EgSPL9 and EgSPL18 exhibit moderately complex tertiary structures. Overall, EgSPL proteins have high conservation with small structural diversities. These results provide important basis for further functional research and evolutionary analysis.

4. Discussion

The completion of genome sequencing of E. grandis [37] enabled researchers to study the classification, evolutionary characteristics, and functional prediction of gene families at the genome level. Twenty-three EgSPL genes were identified, and then, their gene and protein character, phylogenetic relationships, collinearity, and expression patterns were performed. This provides a basis for in-depth research into the functional roles of EgSPL genes in E. grandis.

Phylogenetic and evolutionary analysis showed that the SPLs of E. grandis are closely related to A. thaliana, and there are some orthologous genes. Notably, EgSPL22 and AtSPL7 are orthologous genes, and EgSPL22 is highly expressed in xylem and phloem. Previous studies have shown that AtSPL7 played a major part in copper homeostasis and cadmium stress [25, 26], suggesting that EgSPL22 may also be involved in these processes and contribute to the development of vascular tissue. Additionally, EgSPL3 and EgSPL15 are on the same branch with AtSPL9 and AtSPL15, with high expression in stem apex and lateral branches, which might also mirror their roles in leaf growth regulation and flowering, as well as defense and immunity [14]. Furthermore, EgSPL14, positioned alongside AtSPL10 and AtSPL11 known for their involvement in nutrient transduction to reproductive growth [11], may serve a comparable function. In A. thaliana, SOC1 regulates AtSPL3, AtSPL4, and AtSPL5 by directly binding to the promoter, thereby integrating photoperiodic signals to promote flowering [15]. Interestingly, these genes are all part of the phylogenetic Group IV. Correspondingly, several EgSPL genes in E. grandis, particularly EgSPL7, EgSPL17, EgSPL19, EgSPL21, EgSPL1, EgSPL16, and EgSPL18, are also classified in this same subgroup and have higher expression in flower and stem apex. Given the close phylogenetic relationship to AtSPL3/4/5 and the abundance of light-responsive elements in their promoters, these EgSPL genes in E. grandis are likely to have similar functions in photoperiodic control, potentially influencing flowering and other light-dependent processes.

MiR156, which is known for its highly conserved role across the plant kingdom, and its target genes SPL are considered to be the main regulators of plant growth and development [38, 39] including mediating crucial pathways related to flowering and aging [29]. This mediation occurs either through the modulation of transcription in downstream target genes or by engaging in diverse biological processes via protein interactions [40–43]. In A. thaliana, 11 out of 17 SPL genes are regulated by miR156. In this study, of the 23 EgSPL genes identified in E. grandis, 16 EgSPL genes were found to contain miR156 binding sites. Notably, these genes exhibited elevated expression levels in specific tissues such as the stem apex, flower, and xylem, indicating a potential tissue-specific regulatory role. The concurrent presence of miR156 binding sites in these highly expressed genes suggests a possible intricate interaction with miR156, potentially modulating spatial and developmental gene expression within these key tissues. Notably, EgSPL9 and EgSPL23 showed high expression levels in various tissues and under hormone-induced conditions and are not regulated by miR156. This suggests their significant involvement in distinct developmental processes of E. grandis.

The analysis of gene expression patterns facilitates the identification of key EgSPLs, which are instrumental in mediating plant responses to both biotic and abiotic stresses. Given the observed upregulation of EgSPL9, EgSPL13, and EgSPL23 under SA stress, coupled with the presence of SA-responsive cis-regulatory elements in their promoters, it is plausible to hypothesize that these genes may play a crucial role in the SA-mediated pathways. Similarly, the increased expression of EgSPL9, EgSPL23, EgSPL4, EgSPL14, EgSPL21, and EgSPL22 under MeJA stress, in conjunction with the presence of MeJA-responsive cis-regulatory elements in their promoters, suggests a potential involvement in MeJA-mediated pathways. Interestingly, under SA treatment, EgSPL13 exhibited a transient increase in expression, whereas its expression decreased under MeJA treatment, accentuating its potential role in the SA-mediated pathway. Further experimental validation is necessary to confirm the functional implications of these observations.

5. Conclusion

This study identified 23 SPL genes of E. grandis through bioinformatics analysis and analyzed their physical and chemical properties, evolutionary tree relationships, conserved domains, promoter elements, tissue-specific expression, and hormone stress expression. EgSPL9, EgSPL22, and EgSPL23 genes were found to have strong specific expression, and their functions were predicted by a comprehensive phylogenetic tree and promoter functional elements, which are of guiding significance for future generations to study the light response, growth and development, and stress resistance of E. grandis.

Acknowledgments

This work was supported by the Science and Technology Program of Guangzhou, China (202201011282); China Postdoctoral Science Foundation (2022M721198); and Natural Science Foundation of Guangdong Province (2023A1515010379).

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflicts of interest.

Author Contributions

The authors confirm contribution to the paper as follows: study conception and design: A.W., L.A., and J.M.; data collection: C.F., L.A., J.M., and H.L.; analysis and interpretation of results: L.A. and J.M.; draft manuscript preparation: L.A. All authors reviewed the results and approved the final version of the manuscript.

Funding

This work was supported by the Science and Technology Program of Guangzhou, China (202201011282); China Postdoctoral Science Foundation (2022M721198); and Natural Science Foundation of Guangdong Province (2023A1515010379).

Figure 1 Phylogenetic analysis of SPLs from E. grandis, A. thaliana, S. moellendorffii, and P. patens. The ML tree was constructed using TBtools software with 2000 bootstrap replicates. Signs of different shapes represent SPL genes from E. grandis (green round, Eg), A. thaliana (red star, At), S. moellendorffii (yellow triangle, Sm), and P. patens (blue triangle, Pp).

Figure 2 Analysis of conserved motifs and gene structure of SPL gene family members in E. grandis. (a) TBtools are used to construct the evolutionary tree of 23 EgSPL. (b) The conserved motifs of the EgSPL proteins predicted on MEME. (c) The structure of the EgSPL genes is visualized based on gff file. Exons, introns, and UTRs are represented by yellow square rectangles, black lines, and green square rectangles, respectively.

Figure 3 Chromosome location of EgSPL genes in E. grandis chromosomes. The scale on the left denotes chromosome length in megabases (Mb). Chromosome numbers are indicated by the black numerals at the top of each chromosomal representation. Locations of EgSPL genes are highlighted in red text along the chromosomes. The gradient of color from blue to red represents the relative gene density, with warmer colors corresponding to regions of higher EgSPL gene density.

Figure 4 Collinearity analysis between E. grandis with A. thaliana and P. alba. The horizontal bars represent the chromosomes of each species, labeled at the top of A. thaliana (Chr1 to Chr5), in the middle for E. grandis (Chr1 to Chr11), and at the bottom for P. alba (Chr01 to Chr 19). Syntenic relationships are indicated by blue lines connecting the chromosomes, suggesting evolutionary conservation of gene order between the species.

Figure 5 Analysis of homeopathic elements in promoter sequences of EgSPL family 266 genes by PlantCARE. Distribution of homeopathic elements 2000 bp upstream of 267 EgSPL on the left and the type of homeopathic elements on the right.

Figure 6 Tissue expression pattern of EgSPL family genes. (a) Expression patterns of EgSPL genes in various tissues of semiannual E. grandis, including young leaves, adult leaves, xylem, phloem, roots, and stem apex. (b) Expression patterns of EgSPL genes in different tissues of 3-year-old E. grandis, encompassing young leaves, adult leaves, xylem, flowers, lateral branches, and phloem. (c) Expression patterns of EgSPL genes in xylem and phloem of 6-year-old E. grandis, highlighting age-related changes.

Figure 7 Expression patterns of EgSPL genes under (a) MeJA, (b) SA, and (c) NaCl treatment in E. grandis.

Figure 8 Three-dimensional structure of EgSPL predicted by the SWISS model. The α-helix and β-sheet regions are shown in blue and red, respectively.

Table 1 Physiochemical properties of the SPLs in E. grandis.

Name	Sequence ID	Length (aa) no. of amino acids	MW (kD) molecular mass	pI isoelectric points	Instability index	Grand average of hydropathicity	miR156 target site	Subcellular localization	
EgSPL1	Eucgr.K02545.2.v2.0	551	59.58332	6.83	50.07	−0.48	+	Cytoplasm nucleus	
EgSPL2	Eucgr.E03260.1.v2.0	147	16.97751	6.00	88.37	−1.44	+	Nucleus	
EgSPL3	Eucgr.K01828.1.v2.0	376	39.87973	9.12	59.39	−0.75	−	Nucleus	
EgSPL4	Eucgr.H04114.2.v2.0	869	97.26905	5.92	50.80	−0.44	+	Cytoplasm nucleus	
EgSPL5	Eucgr.K01046.2.v2.0	189	21.07562	9.21	67.96	−1.12	+	Nucleus	
EgSPL6	Eucgr.K01046.1.v2.0	236	26.05625	8.66	64.51	−0.78	+	Nucleus	
EgSPL7	Eucgr.K02708.2.v2.0	363	39.80499	8.96	55.21	−0.76	+	Nucleus	
EgSPL8	Eucgr.D02505.1.v2.0	348	38.35489	8.45	64.98	−0.64	+	Nucleus	
EgSPL9	Eucgr.A01019.1.v2.0	983	109.05482	5.79	56.97	−0.41	−	Nucleus	
EgSPL10	Eucgr.E01600.1.v2.0	367	38.75971	9.07	72.79	−0.64	−	Nucleus	
EgSPL11	Eucgr.K01046.3.v2.0	178	19.94519	9.16	71.57	−1.23	+	Nucleus	
EgSPL12	Eucgr.E01600.2.v2.0	321	34.08069	8.96	65.75	−0.62	−	Nucleus	
EgSPL13	Eucgr.F03303.1.v2.0	186	21.27769	8.86	69.09	−1.28	+	Nucleus	
EgSPL14	Eucgr.B03500.1.v2.0	488	53.53203	8.67	53.61	−0.66	+	Cytoplasm nucleus	
EgSPL15	Eucgr.A02441.1.v2.0	390	41.27372	9.12	54.92	−0.61	+	Nucleus	
EgSPL16	Eucgr.K02545.1.v2.0	551	59.58332	6.83	50.07	−0.48	+	Cytoplasm nucleus	
EgSPL17	Eucgr.K02708.3.v2.0	315	34.66033	9.41	62.03	−0.78	+	Nucleus	
EgSPL18	Eucgr.B00631.1.v2.0	524	57.17946	8.74	55.34	−0.64	+	Nucleus	
EgSPL19	Eucgr.K02708.1.v2.0	363	39.80499	8.96	55.21	−0.76	+	Nucleus	
EgSPL20	Eucgr.E01600.3.v2.0	257	26.71659	7.37	59.71	−0.43	+	Nucleus	
EgSPL21	Eucgr.B03518.1.v2.0	318	34.53342	9.50	53.22	−0.68	+	Nucleus	
EgSPL22	Eucgr.B01228.1.v2.0	821	91.22122	5.77	55.35	−0.42	−	Cytoplasm nucleus	
EgSPL23	Eucgr.F01828.1.v2.0	1078	119.28528	7.11	58.38	−0.53	−	Nucleus	
Note: In the penultimate titled “miR156 Target Site,” the presence of a miR156 target site within the SPL protein sequence is indicated by a plus sign (+), while its absence is denoted by a minus sign (−).
==== Refs
1 Yamasaki K. Kigawa T. Inoue M. An Arabidopsis SBP-domain fragment with a disrupted C-terminal zinc-binding site retains its tertiary structure FEBS Letters 2006 580 8 2109 2116 10.1016/j.febslet.2006.03.014 2-s2.0-33645231123 16554053
2 Salinas M. Xing S. Höhmann S. Berndtgen R. Huijser P. Genomic organization, phylogenetic comparison and differential expression of the SBP-box family of transcription factors in tomato Planta 2012 235 6 1171 1184 10.1007/s00425-011-1565-y 2-s2.0-84861644291 22160465
3 Wang J.-W. Czech B. Weigel D. miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana Cell 2009 138 4 738 749 10.1016/j.cell.2009.06.014 2-s2.0-68749095252 19703399
4 Gandikota M. Birkenbihl R. P. Höhmann S. Cardon G. H. Saedler H. Huijser P. The miRNA156/157 recognition element in the 3′ UTR of the Arabidopsis SBP box gene SPL3 prevents early flowering by translational inhibition in seedlings Plant Journal 2007 49 4 683 693 10.1111/j.1365-313X.2006.02983.x 2-s2.0-33846904347 17217458
5 Wang Z. Wang Y. Kohalmi S. E. Amyot L. Hannoufa A. SQUAMOSA promoter binding protein-like 2 controls floral organ development and plant fertility by activating asymmetric leaves 2 in Arabidopsis thaliana Plant Molecular Biology 2016 92 6 661 674 10.1007/s11103-016-0536-x 2-s2.0-84986327451 27605094
6 Shikata M. Koyama T. Mitsuda N. Ohme-Takagi M. Arabidopsis SBP-box genes SPL10, SPL11 and SPL2 control morphological change in association with shoot maturation in the reproductive phase Plant and Cell Physiology 2009 50 12 2133 2145 10.1093/pcp/pcp148 2-s2.0-72649096308 19880401
7 Cui L. Zheng F. Wang J. miR156a-targeted SBP-box transcription factor SlSPL13 regulates inflorescence morphogenesis by directly activating SFT in tomato Plant Biotechnology Journal 2020 18 8 1670 1682 10.1111/pbi.13331 31916387
8 Ren Y. Ma R. Fan Y. Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in quinoa (Chenopodium quinoa) BMC Genomics 2022 23 1 p. 773 10.1186/s12864-022-08977-9 36434504
9 Ning K. Chen S. Huang H. Jiang J. Yuan H. Li H. Molecular characterization and expression analysis of the SPL gene family with BpSPL9 transgenic lines found to confer tolerance to abiotic stress in Betula platyphylla Suk Plant Cell Tissue and Organ Culture 2017 130 3 469 481 10.1007/s11240-017-1226-3 2-s2.0-85023772189
10 Wang Y. Ruan Q. Zhu X. Wang B. Wei B. Wei X. Identification of alfalfa SPL gene family and expression analysis under biotic and abiotic stresses Scientific Reports 2023 13 1 10.1038/s41598-022-26911-7
11 Xu M. Hu T. Zhao J. Developmental functions of miR156-regulated SQUAMOSA promoter binding protein-like (SPL) genes in Arabidopsis thaliana PLoS Genetics 2016 12 8, article e1006263 10.1371/journal.pgen.1006263 2-s2.0-84984868368 27541584
12 Zimmermann P. Hirsch-Hoffmann M. Hennig L. Gruissem W. GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox Plant Physiology 2004 136 1 2621 2632 10.1104/pp.104.046367 2-s2.0-13444264729 15375207
13 Cardon G. H. Höhmann S. Nettesheim K. Saedler H. Huijser P. Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3: a novel gene involved in the floral transition The Plant Journal 1997 12 2 367 377 10.1046/j.1365-313X.1997.12020367.x 2-s2.0-0031214215 9301089
14 Wang J.-W. Schwab R. Czech B. Mica E. Weigel D. Dual effects of miR156-targeted SPL genes and CYP78A5/KLUH on plastochron length and organ size in Arabidopsis thaliana Plant Cell 2008 20 5 1231 1243 10.1105/tpc.108.058180 2-s2.0-54949093292 18492871
15 Jung J.-H. Ju Y. Seo P. J. Lee J.-H. Park C.-M. The SOC1-SPL module integrates photoperiod and gibberellic acid signals to control flowering time in Arabidopsis Plant Journal 2012 69 4 577 588 10.1111/j.1365-313X.2011.04813.x 2-s2.0-84856744161 21988498
16 Zheng J. Ma Y. Zhang M. Lyu M. Yuan Y. Wu B. Expression pattern of FT/TFL1 and miR156-targeted SPL genes associated with developmental stages in Dendrobium catenatum International Journal of Molecular Sciences 2019 20 11 p. 2725 10.3390/ijms20112725 2-s2.0-85067447534 31163611
17 Padmanabhan M. S. Ma S. Burch-Smith T. M. Czymmek K. Huijser P. Dinesh-Kumar S. P. Novel positive regulatory role for the SPL6 transcription factor in the N TIR-NB-LRR receptor-mediated plant innate immunity PLoS Pathogens 2013 9 3, article e1003235 10.1371/journal.ppat.1003235 2-s2.0-84875993687 23516366
18 Stone J. M. Liang X. Nekl E. R. Stiers J. J. Arabidopsis AtSPL14, a plant-specific SBP-domain transcription factor, participates in plant development and sensitivity to fumonisin B1 Plant Journal 2005 41 5 744 754 10.1111/j.1365-313X.2005.02334.x 2-s2.0-14844320081 15703061
19 Bartel D. P. MicroRNAs Cell 2004 116 2 281 297 10.1016/S0092-8674(04)00045-5 2-s2.0-0347444723 14744438
20 Preston J. C. Hileman L. C. Functional evolution in the plant SQUAMOSA-promoter binding protein-like (SPL) gene family Frontiers in Plant Science 2013 4 p. 80 10.3389/fpls.2013.00080 2-s2.0-84893492962 23577017
21 Unte U. S. Sorensen A. M. Pesaresi P. SPL8, an SBP-box gene that affects pollen sac development in Arabidopsis Plant Cell 2003 15 4 1009 1019 10.1105/tpc.010678 2-s2.0-0037392721 12671094
22 Xing S. Salinas M. Höhmann S. Berndtgen R. Huijser P. miR156-targeted and nontargeted SBP-box transcription factors act in concert to secure male fertility in Arabidopsis The Plant Cell 2010 22 12 3935 3950 10.1105/tpc.110.079343 2-s2.0-79551618858 21177480
23 Zhang Y. Schwarz S. Saedler H. Huijser P. SPL8, a local regulator in a subset of gibberellin-mediated developmental processes in Arabidopsis Plant Molecular Biology 2007 63 3 429 439 17093870
24 Xing S. Salinas M. Garcia-Molina A. Höhmann S. Berndtgen R. Huijser P. SPL8 and miR156-targeted SPL genes redundantly regulate Arabidopsis gynoecium differential patterning Plant Journal 2013 75 4 566 577 10.1111/tpj.12221 2-s2.0-84882449072 23621152
25 Araki R. Mermod M. Yamasaki H. Kamiya T. Fujiwara T. Shikanai T. SPL7 locally regulates copper-homeostasis-related genes in Arabidopsis Journal of Plant Physiology 2018 224-225 137 143 10.1016/j.jplph.2018.03.014 2-s2.0-85053501122 29635212
26 Gielen H. Remans T. Vangronsveld J. Cuypers A. Toxicity responses of Cu and Cd: the involvement of miRNAs and the transcription factor SPL7 BMC Plant Biology 2016 16 1 p. 145 10.1186/s12870-016-0830-4 2-s2.0-84975841254 27352843
27 Hu L. Molecular Mechanisms of the Age-Related Resistance in Arabidopsis thaliana 2023 University of Georgia ProQuest Dissertations & Theses
28 Teng R. M. Wang Y. X. Li H. Liu H. Wang Y. Zhuang J. Identification and expression analysis of the SBP-box gene family related to abiotic stress in tea plant (Camellia sinensis (L.) O. Kuntze) Plant Molecular Biology Reporter 2022 40 1 148 162 10.1007/s11105-021-01306-6
29 Cardon G. Höhmann S. Klein J. Nettesheim K. Saedler H. Huijser P. Molecular characterisation of the Arabidopsis SBP-box genes Gene 1999 237 1 91 104 10.1016/S0378-1119(99)00308-X 2-s2.0-0344418715 10524240
30 Li C. Lu S. Molecular characterization of the SPL gene family in Populus trichocarpa BMC Plant Biology 2014 14 1 p. 131 10.1186/1471-2229-14-131 2-s2.0-84901621246 24884654
31 Zhou Q. Zhang S. Chen F. Genome-wide identification and characterization of the SBP-box gene family in Petunia BMC Genomics 2018 19 1 p. 193 10.1186/s12864-018-4537-9 2-s2.0-85043461248 29703141
32 Li X. Y. Lin E. P. Huang H. H. Niu M. Y. Tong Z. K. Zhang J. H. Molecular characterization of SQUAMOSA promoter binding protein-like (SPL) gene family in Betula luminifera Frontiers in Plant Science 2018 9 p. 608 10.3389/fpls.2018.00608 2-s2.0-85046886697
33 Chen C. Chen H. Zhang Y. TBtools: an integrative toolkit developed for interactive analyses of big biological data Molecular Plant 2020 13 8 1194 1202 10.1016/j.molp.2020.06.009 32585190
34 Fan C. Lyu M. Zeng B. Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis Plant, Cell & Environment 2024 47 4 1363 1378 10.1111/pce.14814 38221855
35 Corrêa L. G. G. Riaño-Pachón D. M. Schrago C. G. Vicentini dos Santos R. Mueller-Roeber B. Vincentz M. The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes PLoS One 2008 3 8, article e2944 10.1371/journal.pone.0002944 2-s2.0-51449092810 18698409
36 Wang Z.-M. Xue W. Dong C.-J. A comparative miRNAome analysis reveals seven fiber initiation-related and 36 novel miRNAs in developing Cotton ovules Molecular Plant 2012 5 4 889 900 10.1093/mp/ssr094 2-s2.0-84866405189 22138860
37 Myburg A. A. Grattapaglia D. Tuskan G. A. The genome of Eucalyptus grandis Nature 2014 510 7505 356 362 10.1038/nature13308 2-s2.0-84903203459 24919147
38 Yu S. Lian H. Wang J.-W. Plant developmental transitions: the role of microRNAs and sugars Current Opinion in Plant Biology 2015 27 1 7 10.1016/j.pbi.2015.05.009 2-s2.0-84930208941 26042537
39 Fouracre J. P. Poethig R. S. The role of small RNAs in vegetative shoot development Current Opinion in Plant Biology 2016 29 64 72 10.1016/j.pbi.2015.11.006 2-s2.0-84952673589 26745378
40 Shikata M. Yamaguchi H. Sasaki K. Ohtsubo N. Overexpression of Arabidopsis miR157b induces bushy architecture and delayed phase transition in Torenia fournieri Planta 2012 236 4 1027 1035 22552637
41 Wang J.-W. Park M. Y. Wang L.-J. MiRNA control of vegetative phase change in trees PLoS Genetics 2011 7 2 p. e1002012 10.1371/journal.pgen.1002012 2-s2.0-79952261472 21383862
42 Wu G. Park M. Y. Conway S. R. Wang J.-W. Weigel D. Poethig R. S. The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis Cell 2009 138 4 750 759 10.1016/j.cell.2009.06.031 2-s2.0-68749100821 19703400
43 Fu C. Sunkar R. Zhou C. Overexpression of miR156 in switchgrass (Panicum virgatum L.) results in various morphological alterations and leads to improved biomass production Plant Biotechnology Journal 2012 10 4 443 452 10.1111/j.1467-7652.2011.00677.x 2-s2.0-84862800354 22239253
