Genome-wide characterization of the SPL gene family involved in the age development of Jatropha curcas

Transcriptomic and metabolomic profiling reveals molecular regulatory network involved in flower development and phenotypic changes in two Lonicera macranthoides varieties

Due to the medicinal importance of the flowers of Xianglei type (XL) Lonicera macranthoides, it is important to understand the molecular mechanisms that underlie their development. In this study, we elucidated the transcriptomic and metabolomic mechanisms that underlie the flower development mechanism of two L. macranthoides varieties. In this study, 3435 common differentially expressed unigenes (DEGs) and 1138 metabolites were identified. These common DEGs were mainly enriched in plant hormone signal transduction pathways. Metabolomic analysis showed that amino acids were the main metabolites of differential accumulation in wild-type (WT) L. macranthoides, whereas in XL, they were flavonoids and phenylalanine metabolites. Genes and transcription factors (TFs), such as MYB340, histone deacetylase 1 (HDT1), small auxin-up RNA 32 (SAUR32), auxin response factor 6 (ARF6), PIN-LIKES 7 (PILS7), and WRKY6, likely drive metabolite accumulation. Plant hormone signals, especially auxin signals, and various TFs induce downstream flower organ recognition genes, resulting in a differentiation of the two L. macranthoides varieties in terms of their developmental trajectories. In addition, photoperiodic, autonomous, and plant hormone pathways jointly regulated the L. macranthoides corolla opening. SAUR32, Arabidopsis response regulator 9 (ARR9), Gibberellin receptor (GID1B), and Constans-like 10 (COL10) were closely related to the unfolding of the L. macranthoides corolla. These findings offer valuable understanding of the flower growth process of L. macranthoides and the excellent XL phenotypes at the molecular level.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04019-1.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

Genome-wide analysis of the SPL transcription factor family and its response to water stress in sunflower (Helianthus annuus)

SPL (SQUAMOSA promoter binding proteins-like) are plant-specific transcription factors that play essential roles in a variety of developmental processes as well as the ability to withstand biotic and abiotic stresses. To date, numerous species have been investigated for the SPL gene family, but so far, no SPL family genes have been thoroughly identified and characterized in the sunflower (Helianthus annuus). In this study, 25 SPL genes were identified in the sunflower genome and were unevenly distributed on 11 chromosomes. According to phylogeny analysis, 59 SPL genes from H. annuus, O. sativa, and A. thaliana were clustered into seven groups. Furthermore, the SPL genes in groups-I and II were demonstrated to be potential targets of miR156. Synteny analysis showed that 7 paralogous gene pairs exist in HaSPL genes and 26 orthologous gene pairs exist between sunflower and rice, whereas 21 orthologous gene pairs were found between sunflower and Arabidopsis. Segmental duplication appears to have played a vital role in the expansion processes of sunflower SPL genes, and because of selection pressure, all duplicated genes have undergone purifying selection. Tissue-specific gene expression analysis of the HaSBP genes proved their diverse spatiotemporal expression patterns, which were predominantly expressed in floral organs and differentially expressed in stem, axil, and root tissues. The expression pattern of HaSPL genes under water stress showed broad involvement of HaSPLs in the response to flood and drought stresses. This genome-wide identification investigation provides detailed information on the sunflower SPL transcription factor gene family and establishes a strong platform for future research on sunflower responses to abiotic stress tolerance.

Genome-Wide Analysis of SPL/miR156 Module and Its Expression Analysis in Vegetative and Reproductive Organs of Oil Palm (Elaeis guineensis)

The SPL (SQUAMOSA-promoter binding protein-like) gene family is one of the largest plant transcription factors and is known to be involved in the regulation of plant growth, development, and stress responses. The genome-wide analysis of SPL gene members in a diverse range of crops has been elucidated. However, none of the genome-wide studies on the SPL gene family have been carried out for oil palm, an important oil-yielding plant. In this research, a total of 24 EgSPL genes were identified via a genome-wide approach. Phylogenetic analysis revealed that most of the EgSPLs are closely related to the Arabidopsis and rice SPL gene members. EgSPL genes were mapped onto the only nine chromosomes of the oil palm genome. Motif analysis revealed conservation of the SBP domain and the occurrence of 1-10 motifs in EgSPL gene members. Gene duplication analysis demonstrated the tandem duplication of SPL members in the oil palm genome. Heatmap analysis indicated the significant expression of SPL genes in shoot and flower organs of oil palm plants. Among the identified EgSPL genes, a total 14 EgSPLs were shown to be targets of miR156. Real-time PCR analysis of 14 SPL genes showed that most of the EgSPL genes were more highly expressed in female and male inflorescences of oil palm plants than in vegetative tissues. Altogether, the present study revealed the significant role of EgSPL genes in inflorescence development.

Heat shock factor binding protein BrHSBP1 regulates seed and pod development in Brassica rapa

Plant heat shock factor binding proteins (HSBPs) are well known for their implication in the negative regulation of heat stress response (HSR) pathways. Herein, we report on the hitherto unknown functions of HSBP1 in Brassica rapa (BrHSBP1). BrHBSP1 was found to be predominant in flower buds and young leaves, while its segmental duplicate, BrHSBP1-like, was abundant in green siliques. Exposure to abiotic stress conditions, such as heat, drought, cold, and H2O2, and to phytohormones was found to differentially regulate BrHSBP1. The activity of BrHSBP1-GFP fusion proteins revealed their cellular localization in nuclei and cytosols. Transgenic overexpression of BrHSBP1 (BrHSBP1OX) improved pod and seed sizes, while CRISPR-Cas BrHSBP1 knock-out mutants (Brhsbp1_KO) were associated with aborted seed and pod development. The transcriptomic signatures of BrHSBP1OX and Brhsbp1_KO lines revealed that 360 and 2381 genes, respectively, were differentially expressed (Log2FC≥2, padj<0.05) expressed relative to control lines. In particular, developmental processes, including plant reproductive structure development (RSD)-related genes, were relatively downregulated in Brhsbp1_KO. Furthermore, yeast two-hybrid assays confirmed that BrHSBP1 can physically bind to RSD and other genes. Taking the findings together, it is clear that BrHSBP1 is involved in seed development via the modulation of RSD genes. Our findings represent the addition of a new regulatory player in seed and pod development in B. rapa.

Genome-Wide Analysis of SQUAMOSA-Promoter-Binding Protein-like Family in Flowering Pleioblastus pygmaeus

SQUAMOSA Promoter-Binding Protein-Like (SPL) family is well-known for playing an important role in plant growth and development, specifically in the reproductive process. Bamboo plants have special reproductive characteristics with a prolonged vegetative phase and uncertain flowering time. However, the underlying functions of SPL genes in reproductive growth are undisclosed in bamboo plants. In the study, a total of 28 SPLs were screened from an ornamental dwarf bamboo species, Pleioblastus pygmaeus. Phylogenetic analysis indicates that 183 SPLs from eight plant species can be classified into nine subfamilies, and the 28 PpSPLs are distributed among eight subfamilies. Homologous analysis shows that as many as 32 pairs of homologous genes were found between P. pygmaeus and rice, and 83 pairs were found between P. pygmaeus and Moso bamboo, whose Ka/Ks values are all <1. MiRNA target prediction reveals that 13 out of the 28 PpSPLs have recognition sites complementary to miRNA156. To screen the SPLs involved in the reproductive growth of bamboo plants, the mRNA abundance of the 28 PpSPLs was profiled in the different tissues of flowering P. pygmaeus and non-flowering plants by RNA-Seq. Moreover, the relative expression level of eight PpSPLs is significantly higher in flowering P. pygmaeus than that in non-flowering plants, which was also validated by RT-qPCR. Combined with phylogenetic analysis and homologous analysis, the eight significant, differentially expressed PpSPLs were identified to be associated with the reproductive process and flower organ development. Among them, there are four potential miRNA156-targeting PpSPLs involved in the flowering process. Of significant interest in the study is the identification of 28 SPLs and the exploration of four key flowering-related SPLs from P. pygmaeus, which provides a theoretic basis for revealing the underlying functions of SPLs in the reproductive growth of bamboo plants.

Selection and validation of reference genes for RT-qPCR analysis of different organs at various development stages in Caragana intermedia

Reverse transcription quantitative PCR (RT-qPCR) is a technique widely used to investigate the expression of genes. An appropriate reference gene (RG) is essential for RT-qPCR analysis to obtain accurate and reliable results. Caragana intermedia plays an important role in afforestation as a bush. However, due to the lack of appropriate RGs, the research on development-related genes is limited. In this study, the selection for suitable RGs of different organs at various development stages to normalize the results of RT-qPCR about development-related genes was performed. To test the expression stability across all samples, we used the software algorithms such as geNorm, NormFinder, BestKeeper, and RefFinder to evaluate all the candidate RGs. Our results showed that CiEF1α was the most stable RG with little fluctuation among all samples. In addition, CiGAPDH in roots, CiSKIP1 in stems and leaves, and CiEF1α in different organs were selected as the most stable RGs. To confirm the applicability of the most stable RGs, the relative expression of CiWRKY17 was normalized using different candidate RGs. Taken together, our research laid a foundation for the study of development-related genes in C. intermedia.

Long-day photoperiod and cool temperature induce flowering in cassava: Expression of signaling genes

Cassava is a staple food crop in the tropics, and is of particular importance in Africa. Recent development of genomic selection technology have improved the speed of cassava breeding; however, cassava flower initiation and development remains a bottleneck. The objectives of the current studies were to elucidate the effect of photoperiod, temperature and their interactions on the time of flowering and flower development in controlled environments, and to use RNA-sequencing to identify transcriptome expression underlying these environmental responses. Compared to a normal tropical day-length of 12 h, increasing the photoperiod by 4 h or decreasing the air temperature from 34/31 to 22°/19°C (day/night) substantially hastened the time to flowering. For both photoperiod and temperature, the environment most favorable for flowering was opposite the one for storage root harvest index. There was a pronounced treatment interaction: at warm day-time temperatures, percent flowering was low, and photoperiod had little effect. In contrast, at cooler temperatures, percent flowering increased, and long-day (LD) photoperiod had a strong effect in hastening flowering. In response to temperature, many differentially expressed genes in the sugar, phase-change, and flowering-time-integrator pathways had expression/flowering patterns in the same direction as in Arabidopsis (positive or negative) even though the effect of temperature on flowering operates in the reverse direction in cassava compared to Arabidopsis. Three trehalose-6-phosphate-synthase-1 (TPS1) genes and four members of the SPL gene family had significantly increased expression at cool temperature, suggesting sugar signaling roles in flower induction. In response to LD photoperiod, regulatory genes were expressed as in Arabidopsis and other LD flowering plants. Several hormone-related genes were expressed in response to both photoperiod and temperature. In summary, these findings provide insight on photoperiod and temperature responses and underlying gene expression that may assist breeding programs to manipulate flowering for more rapid crop improvement.

Genome-Wide Identification and Characterization of SPL Family Genes in Chenopodium quinoa

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode a large family of plant-specific transcription factors that play important roles in plant growth, development, and stress responses. However, there is little information available on SPL genes in Chenopodiaceae. Here, 23 SPL genes were identified and characterized in the highly nutritious crop Chenopodium quinoa. Chromosome localization analysis indicated that the 23 CqSPL genes were unevenly distributed on 12 of 18 chromosomes. Two zinc finger-like structures and a nuclear location signal were present in the SBP domains of all CqSPLs, with the exception of CqSPL21/22. Phylogenetic analysis revealed that these genes were classified into eight groups (group I–VIII). The exon–intron structure and motif composition of the genes in each group were similar. Of the 23 CqSPLs, 13 were potential targets of miR156/7. In addition, 5 putative miR156-encoding loci and 13 putative miR157-encoding loci were predicted in the quinoa genome, and they were unevenly distributed on chromosome 1–4. The expression of several Cqu-MIR156/7 loci was confirmed by reverse transcription polymerase chain reaction in seedlings. Many putative cis-elements associated with light, stress, and phytohormone responses were identified in the promoter regions of CqSPLs, suggesting that CqSPL genes are likely involved in the regulation of key developmental processes and stress responses. Expression analysis revealed highly diverse expression patterns of CqSPLs among tissues. Many CqSPLs were highly expressed in leaves, flowers, and seeds, and their expression levels were low in the roots, suggesting that CqSPLs play distinct roles in the development and growth of quinoa. The expression of 13 of 23 CqSPL genes responded to salt treatment (11 up-regulated and 2 down-regulated). A total of 22 of 23 CqSPL genes responded to drought stress (21 up-regulated and 1 down-regulated). Moreover, the expression of 14 CqSPL genes was significantly altered following cadmium treatment (3 up-regulated and 11 down-regulated). CqSPL genes are thus involved in quinoa responses to salt/drought and cadmium stresses. These findings provide new insights that will aid future studies of the biological functions of CqSPLs in C. quinoa.

Deciphering Pleiotropic Signatures of Regulatory SNPs in Zea mays L. Using Multi-Omics Data and Machine Learning Algorithms

Maize is one of the most widely grown cereals in the world. However, to address the challenges in maize breeding arising from climatic anomalies, there is a need for developing novel strategies to harness the power of multi-omics technologies. In this regard, pleiotropy is an important genetic phenomenon that can be utilized to simultaneously enhance multiple agronomic phenotypes in maize. In addition to pleiotropy, another aspect is the consideration of the regulatory SNPs (rSNPs) that are likely to have causal effects in phenotypic development. By incorporating both aspects in our study, we performed a systematic analysis based on multi-omics data to reveal the novel pleiotropic signatures of rSNPs in a global maize population. For this purpose, we first applied Random Forests and then Markov clustering algorithms to decipher the pleiotropic signatures of rSNPs, based on which hierarchical network models are constructed to elucidate the complex interplay among transcription factors, rSNPs, and phenotypes. The results obtained in our study could help to understand the genetic programs orchestrating multiple phenotypes and thus could provide novel breeding targets for the simultaneous improvement of several agronomic traits.

Genome-wide identification, gene cloning, subcellular location and expression analysis of SPL gene family in P. granatum L

BACKGROUNDS

Pomegranate is an excellent tree species with nutritional, medicinal, ornamental and ecological values. Studies have confirmed that SPL factors play an important role in floral transition and flower development.

RESULTS

Used bioinformatics methods, 15 SPL (SQUAMOSA promoter-binding protein-like) genes were identified and analyzed from the 'Taishanhong' pomegranate (P. granatum L.) genome. Phylogenetic analysis showed that PgSPLs were divided into six subfamilies (G1 ~ G6). PgSPL promoter sequences contained multiple cis-acting elements associated with abiotic stress or hormonal response. Based on the transcriptome data, expression profiles of different tissues and different developmental stages showed that PgSPL genes had distinct temporal and spatial expression characteristics. The expression analysis of miR156 in small RNA sequencing results showed that miR156 negatively regulated the expression of target genes. qRT-PCR analysis showed that the expression levels of PgSPL2, PgSPL3, PgSPL6, PgSPL11 and PgSPL14 in leaves were significantly higher than those in buds and stems (p < 0.05). The expression levels of PgSPL5, PgSPL12 and PgSPL13 in flower buds were significantly higher than that in leaves and stems (p < 0.05). The full-length of coding sequence of PgSPL5 and PgSPL13 were obtained by homologous cloning technology. The full length of PgSPL5 is 1020 bp, and PgSPL13 is 489 bp, which encodes 339 and 162 amino acids, respectively. Further investigation revealed that PgSPL5 and PgSPL13 proteins were located in the nucleus. Exogenous plant growth regulator induction experiments showed that PgSPL5 was up-regulated in leaves and stems. PgSPL13 was up-regulated in leaves and down-regulated in stems. When sprayed with 6-BA, IBA and PP333 respectively, PgSPL5 and PgSPL13 were up-regulated most significantly at P2 (bud vertical diameter was 5.1 ~ 12.0 mm) stage of bisexual and functional male flowers.

CONCLUSIONS

Our findings suggested that PgSPL2, PgSPL3, PgSPL6, PgSPL11 and PgSPL14 played roles in leaves development of pomegranate. PgSPL5, PgSPL12 and PgSPL13 played roles in pomegranate flower development. PgSPL5 and PgSPL13 were involved in the response process of different plant hormone signal transduction in pomegranate development. This study provided a robust basis for further functional analyses of SPL genes in pomegranate.

Genome-Wide Characterization and Expression Analysis of the SBP-Box Gene Family in Sweet Orange (Citrus sinensis)

SBP-box is an important plant-specific transcription factor family and is involved in diverse biological processes. Here, we identified a total of 15 SBP-BOX genes in the important fruit crop sweet orange (Citrus sinensis) and characterized their gene structures, conserved domain and motif, chromosomal location, and cis-acting regulatory elements. SBP genes were classified into four subfamilies based on the amino acid sequence homology, and the classification is equally strongly supported by the gene and protein structures. Our analysis revealed that segmental duplication events were the main driving force in the evolution of CsSBP genes, and gene pairs might undergo extensive purifying selection. Further synteny analysis of the SBP members among sweet orange and other plant species provides valuable information for clarifying the CsSBP family evolutionary relationship. According to publicly available RNA-seq data and qRT-PCR analysis from various sweet orange tissues, CsSBP genes may be expressed in different tissues and developmental stages. Gene expression analysis showed variable expression profiles of CsSBP genes under various abiotic stresses, such as high and low-temperature, salt, and wound treatments, demonstrating the potential role of SBP members in sweet orange response to abiotic stress. Noticeably, all CsSBP genes were also downregulated in sweet orange upon the infection of an important fungal pathogen Diaporthe citri. Our results provide valuable information for exploring the role of SBP-Box in sweet orange.

Genome-wide identification, phylogenetic analysis, and expression analysis of the SPL gene family in orchardgrass (Dactylis glomerata L.)

SPL (SQUAMOSA promoter binding protein-like) is a plant-specific transcription factor family that contains the conserved SBP domain, which plays a vital role in the vegetative-to-reproductive phase transition, flowering development and regulation, tillering/branching, and stress responses. Although the SPL family has been identified and characterized in various plant species, limited information about it has been obtained in orchardgrass, which is a critical forage crop worldwide. In this study, 17 putative DgSPL genes were identified among seven chromosomes, and seven groups that share similar gene structures and conserved motifs were determined by phylogenetic analysis. Of these, eight genes have potential target sites for miR156. cis-Element and gene ontology annotation analysis indicated DgSPLs may be involved in regulating development and abiotic stress responses. The expression patterns of eight DgSPL genes at five developmental stages, in five tissues, and under three stress conditions were determined by RNA-seq and qRT-PCR. These assays indicated DgSPLs are involved in vegetative-to-reproductive phase transition, floral development, and stress responses. The transient expression analysis in tobacco and heterologous expression assays in yeast indicated that miR156-targeted DG1G01828.1 and DG0G01071.1 are nucleus-localized proteins, that may respond to drought, salt, and heat stress. Our study represents the first systematic analysis of the SPL family in orchardgrass. This research provides a comprehensive assessment of the DgSPL family, which lays the foundation for further examination of the role of miR156/DgSPL in regulating development and stress responses in forages grasses.

Comparative analysis of the SPL gene family in five Rosaceae species: Fragaria vesca, Malus domestica, Prunus persica, Rubus occidentalis, and Pyrus pyrifolia

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are very important for the plant growth and development. Here 15 RoSPLs were identified in Rubus occidentalis. The conserved domains and motifs, phylogenetic relationships, posttranscriptional regulation, and physiological function of the 92 SPL family genes in Fragaria vesca, Malus domestica, Prunus persica, R. occidentalis, and Pyrus pyrifolia were analyzed. Sequence alignment and phylogenetic analysis showed the SPL proteins had sequence conservation, some FvSPLs could be lost or developed, and there was a closer relationship between M. domestica and P. pyrifolia, F. vesca and R. occidentalis, respectively. Genes with similar motifs clustering together in the same group had their functional redundancy. Based on the function of SPLs in Arabidopsis thaliana, these SPLs could be involved in vegetative transition from juvenile to adult, morphological change in the reproductive phase, anthocyanin biosynthesis, and defense stress. Forty-eight SPLs had complementary sequences of miR156, of which nine PrpSPLs in P. persica and eight RoSPLs in R. occidentalis as the potential targets of miR156 were reported for the first time, suggesting the conservative regulatory effects of miR156 and indicating the roles of miR156-SPL modules in plant growth, development, and defense response. It provides a basic understanding of SPLs in Rosaceae plants.

Isoprenoid-Derived Metabolites and Sugars in the Regulation of Flowering Time: Does Day Length Matter?

In plants, a diverse set of pathways regulate the transition to flowering, leading to remarkable developmental flexibility. Although the importance of photoperiod in the regulation of flowering time is well known, increasing evidence suggests the existence of crosstalk among the flowering pathways regulated by photoperiod and metabolic pathways. For example, isoprenoid-derived phytohormones (abscisic acid, gibberellins, brassinosteroids, and cytokinins) play important roles in regulating flowering time. Moreover, emerging evidence reveals that other metabolites, such as chlorophylls and carotenoids, as well as sugar metabolism and sugar accumulation, also affect flowering time. In this review, we summarize recent findings on the roles of isoprenoid-derived metabolites and sugars in the regulation of flowering time and how day length affects these factors.

Phylogenetic Analysis of the SQUAMOSA Promoter-Binding Protein-Like Genes in Four Ipomoea Species and Expression Profiling of the IbSPLs During Storage Root Development in Sweet Potato (Ipomoea batatas)

As a major plant-specific transcription factor family, SPL genes play a crucial role in plant growth, development, and stress tolerance. The SPL transcription factor family has been widely studied in various plant species; however, systematic studies on SPL genes in the genus Ipomoea are lacking. Here, we identified a total of 29, 27, 26, and 23 SPLs in Ipomoea batatas, Ipomoea trifida, Ipomoea triloba, and Ipomoea nil, respectively. Based on the phylogenetic analysis of SPL proteins from model plants, the Ipomoea SPLs were classified into eight clades, which included conserved gene structures, domain organizations and motif compositions. Moreover, segmental duplication, which is derived from the Ipomoea lineage-specific whole-genome triplication event, was speculated to have a predominant role in Ipomoea SPL expansion. Particularly, tandem duplication was primarily responsible for the expansion of SPL subclades IV-b and IV-c. Furthermore, 25 interspecific orthologous groups were identified in Ipomoea, rice, Arabidopsis, and tomato. These findings support the expansion of SPLs in Ipomoea genus, with most of the SPLs being evolutionarily conserved. Of the 105 Ipomoea SPLs, 69 were predicted to be the targets of miR156, with seven IbSPLs being further verified as targets using degradome-seq data. Using transcriptomic data from aboveground and underground sweet potato tissues, IbSPLs showed diverse expression patterns, including seven highly expressed IbSPLs in the underground tissues. Furthermore, the expression of 11 IbSPLs was validated using qRT-PCR, and two (IbSPL17/IbSPL28) showed significantly increased expression during root development. Additionally, the qRT-PCR analysis revealed that six IbSPLs were strongly induced in the roots under phytohormone treatments, particularly zeatin and abscisic acid. Finally, the transcriptomic data of storage roots from 88 sweet potato accessions were used for weighted gene co-expression network analysis, which revealed four IbSPLs (IbSPL16/IbSPL17/IbSPL21/IbSPL28) clusters with genes involved in "regulation of root morphogenesis," "cell division," "cytoskeleton organization," and "plant-type cell wall organization or biogenesis," indicating their potential role in storage root development. This study not only provides novel insights into the evolutionary and functional divergence of the SPLs in the genus Ipomoea but also lays a foundation for further elucidation of the potential functional roles of IbSPLs on storage root development.