Molecular characterization of the SPL gene family in Populus trichocarpa

TaSPL6B, a member of the Squamosa promoter binding protein-like family, regulates shoot branching and florescence in Arabidopsis thaliana

BACKGROUND

Squamosa promoter-binding protein-like (SPL) proteins are essential to plant growth and development as plant-specific transcription factors. However, the functions of SPL proteins in wheat need to be further explored.

RESULTS

We cloned and characterized TaSPL6B of wheat in this study. Analysis of physicochemical properties revealed that it contained 961 amino acids and had a molecular weight of 105 kDa. Full-length TaSPL6B transcription activity was not validated in yeast and subcellular localization analysis revealed that TaSPL6B was distributed in the nucleus. Ectopic expression of TaSPL6B in Arabidopsis led to increasing number of branches and early flowering. TaSPL6B was highly transcribed in internodes of transgenic Arabidopsis. The expression of AtSMXL6/AtSMXL7/AtSMXL8 (homologous genes of TaD53) was markedly increased, whereas the expression of AtSPL2 (homologous genes of TaSPL3) and AtBRC1 (homologous genes of TaTB1) was markedly reduced in the internodes of transgenic Arabidopsis. Besides, TaSPL6B, TaSPL3 and TaD53 interacted with one another, as demonstrated by yeast two-hybrid and bimolecular fluorescence complementation assays. Therefore, we speculated that TaSPL6B brought together TaD53 and TaSPL3 and enhanced the inhibition effect of TaD53 on TaSPL3 through integrating light and strigolactone signaling pathways, followed by suppression of TaTB1, a key repressor of tillering.

CONCLUSIONS

As a whole, our findings contribute to a better understanding of how SPL genes work in wheat and will be useful for further research into how TaSPL6B affects yield-related traits in wheat.

Genome-wide identification and analysis of SPL gene family in chickpea (Cicer arietinum L.)

A transcription factor in plants encodes SQUAMOSA promoter binding protein-like (SPL) serves a broad spectrum of important roles for the plant, like, growth, flowering, and signal transduction. A gene family that encodes SPL proteins is documented in various model plant species, including Arabidopsis thaliana and Oryza sativa. Chickpea (Cicer arietinum), a leguminous crop, has not been thoroughly explored with regard to the SPL protein-encoding gene family. Chickpea SPL family genes were located and characterized computationally using a genomic database. Gene data of chickpea were obtained from the phytozome repository and was examined using bioinformatics methods. For investigating the possible roles of SPLs in chickpea, genome-wide characterization, expression, as well as structural analyses of this SPL gene family were performed. Cicer arietinum genome had 19 SPL genes, whereas, according to phylogenetic analysis, the SPLs in chickpea are segregated among four categories: Group-I has 2 introns, Group-II and IV have 1-2 introns (except CaSPL13 and CaSPL15 having 3 introns), and Group-III has 9 introns (except CaSPL1 and CaSPL11 with 1 and 8 introns, respectively). The SBP domain revealed that SPL proteins featured two zinc-binding sites, i.e., C3H and C2HC and one nuclear localization signal. All CaSPL proteins are found to contain highly conserved motifs, i.e., Motifs 1, 2, and 4, except CaSPL10 in which Motifs 1 and 4 were absent. Following analysis, it was found that Motifs 1 and 2 of the chickpea SBP domain are Zinc finger motifs, and Motif 4 includes a nuclear localization signal. All pairs of CaSPL paralogs developed by purifying selection. The CaSPL promoter investigation discovered cis-elements that are responsive to stress, light, and phytohormones. Examination of their expression patterns highlighted major CaSPLs to be evinced primarily among younger pods and flowers. Indicating their involvement in the plant's growth as well as development, along with their capacity to react as per different situations by handling the regulation of target gene's expression, several CaSPL genes are also expressed under certain stress conditions, namely, cold, salt, and drought. The majority of the CaSPL genes are widely expressed and play crucial roles in terms of the plant's growth, development, and responses to the environmental-stress conditions. Our work provides extensive insight into the gene family CaSPL, which might facilitate further studies related to the evolution and functions of the SPL genes for chickpea and other plant species.

A Comprehensive Analysis of the Peanut SQUAMOSA Promoter Binding Protein-like Gene Family and How AhSPL5 Enhances Salt Tolerance in Transgenic Arabidopsis

SPL (SQUAMOSA promoter binding protein-like), as one family of plant transcription factors, plays an important function in plant growth and development and in response to environmental stresses. Despite SPL gene families having been identified in various plant species, the understanding of this gene family in peanuts remains insufficient. In this study, thirty-eight genes (AhSPL1-AhSPL38) were identified and classified into seven groups based on a phylogenetic analysis. In addition, a thorough analysis indicated that the AhSPL genes experienced segmental duplications. The analysis of the gene structure and protein motif patterns revealed similarities in the structure of exons and introns, as well as the organization of the motifs within the same group, thereby providing additional support to the conclusions drawn from the phylogenetic analysis. The analysis of the regulatory elements and RNA-seq data suggested that the AhSPL genes might be widely involved in peanut growth and development, as well as in response to environmental stresses. Furthermore, the expression of some AhSPL genes, including AhSPL5, AhSPL16, AhSPL25, and AhSPL36, were induced by drought and salt stresses. Notably, the expression of the AhSPL genes might potentially be regulated by regulatory factors with distinct functionalities, such as transcription factors ERF, WRKY, MYB, and Dof, and microRNAs, like ahy-miR156. Notably, the overexpression of AhSPL5 can enhance salt tolerance in transgenic Arabidopsis by enhancing its ROS-scavenging capability and positively regulating the expression of stress-responsive genes. These results provide insight into the evolutionary origin of plant SPL genes and how they enhance plant tolerance to salt stress.

SPL16 and SPL23 mediate photoperiodic control of seasonal growth in Populus trees

Perennial trees in boreal and temperate regions undergo growth cessation and bud set under short photoperiods, which are regulated by phytochrome B (phyB) photoreceptors and PHYTOCHROME INTERACTING FACTOR 8 (PIF8) proteins. However, the direct signaling components downstream of the phyB-PIF8 module remain unclear. We found that short photoperiods suppressed the expression of miR156, while upregulated the expression of miR156-targeted SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE 16 (SPL16) and SPL23 in leaves and shoot apices of Populus trees. Accordingly, either overexpression of MIR156a/c or mutagenesis of SPL16/23 resulted in the attenuation of growth cessation and bud set under short days (SD), whereas overexpression of SPL16 and SPL23 conferred early growth cessation. We further showed that SPL16 and SPL23 directly suppressed FLOWERING LOCUS T2 (FT2) expression while promoted BRANCHED1 (BRC1.1 and BRC1.2) expression. Moreover, we revealed that PIF8.1/8.2, positive regulators of growth cessation, directly bound to promoters of MIR156a and MIR156c and inhibited their expression to modulate downstream pathways. Our results reveal a connection between the phyB-PIF8 module-mediated photoperiod perception and the miR156-SPL16/23-FT2/BRC1 regulatory cascades in SD-induced growth cessation. Our study provides insights into the rewiring of a conserved miR156-SPL module in the regulation of seasonal growth in Populus trees.

Genome-wide investigation of SQUAMOSA promoter binding protein-like genes in Liriodendron and functional characterization of LcSPL2

The plant-specific SQUAMOSA promoter-binding protein-like (SPL) transcription factors play a pivotal role in various developmental processes, including leaf morphogenesis and vegetative to reproductive phase transition. Liriodendron chinense and Liriodendron tulipifera are widely used in landscaping due to their tulip-like flowers and peculiar leaves. However, the SPL gene family in Liriodendron has not been identified and systematically characterized. We systematically identified and characterized the SPL family members in Liriodendron, including phylogeny, gene structure and syntenic analyses. Subsequently, we quantified the expression patterns of LcSPLs across various tissue sites through transcription-quantitative polymerase chain reaction (RT-qPCR) assays and identified the target gene, LcSPL2. Finally, we characterized the functions of LcSPL2 via ectopic transformation. Altogether, 17 LcSPL and 18 LtSPL genes were genome-widely identified in L. chinense and L. tulipifera, respectively. All the 35 SPLs were grouped into 9 clades. Both species had three SPL gene pairs arising from segmental duplication events, and the LcSPLs displayed high collinearity with the L. tulipifera genome. RT-qPCR assays showed that SPL genes were differentially expressed in different tissues, especially. Because LcSPL2 is highly expressed in pistils and leaves, it was selected to describe the SPL gene family of L. chinense by ectopic expression. We showed that overexpression of LcSPL2 in Arabidopsis thaliana resulted in earlier flowering and fewer rosette leaves. Moreover, we observed that overexpression of LcSPL2 in A. thaliana up-regulated the expression levels of four genes related to flower development. This study identified SPL genes in Liriodendron and characterized the function of LcSPL2 in advancing flower development.

Comparative analysis of SPL transcription factors from streptophyte algae and embryophytes reveals evolutionary trajectories of SPL family in streptophytes

SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode plant-specific transcription factors which are important regulators of diverse plant developmental processes. We took advantage of available genome sequences of streptophyte algae representatives to investigate the relationships of SPL genes between freshwater green algae and land plants. Our analysis showed that streptophyte algae, hornwort and liverwort genomes encode from one to four SPL genes which is the smallest set, in comparison to other land plants studied to date. Based on the phylogenetic analysis, four major SPL phylogenetic groups were distinguished with Group 3 and 4 being sister to Group 1 and 2. Comparative motif analysis revealed conserved protein motifs within each phylogenetic group and unique bryophyte-specific motifs within Group 1 which suggests lineage-specific protein speciation processes. Moreover, the gene structure analysis also indicated the specificity of each by identifying differences in exon-intron structures between the phylogenetic groups, suggesting their evolutionary divergence. Since current understanding of SPL genes mostly arises from seed plants, the presented comparative and phylogenetic analyzes from freshwater green algae and land plants provide new insights on the evolutionary trajectories of the SPL gene family in different classes of streptophytes.

Genome-Wide Identification and Expression Analysis of the SBP-Box Gene Family in Loquat Fruit Development

The loquat (Eriobotrya japonica L.) is a special evergreen tree, and its fruit is of high medical and health value as well as having stable market demand around the world. In recent years, research on the accumulation of nutrients in loquat fruit, such as carotenoids, flavonoids, and terpenoids, has become a hotspot. The SBP-box gene family encodes transcription factors involved in plant growth and development. However, there has been no report on the SBP-box gene family in the loquat genome and their functions in carotenoid biosynthesis and fruit ripening. In this study, we identified 28 EjSBP genes in the loquat genome, which were unevenly distributed on 12 chromosomes. We also systematically investigated the phylogenetic relationship, collinearity, gene structure, conserved motifs, and cis-elements of EjSBP proteins. Most EjSBP genes showed high expression in the root, stem, leaf, and inflorescence, while only five EjSBP genes were highly expressed in the fruit. Gene expression analysis revealed eight differentially expressed EjSBP genes between yellow- and white-fleshed fruits, suggesting that the EjSBP genes play important roles in loquat fruit development at the breaker stage. Notably, EjSBP01 and EjSBP19 exhibited completely opposite expression patterns between white- and yellow-fleshed fruits during fruit development, and showed a close relationship with SlCnr involved in carotenoid biosynthesis and fruit ripening, indicating that these two genes may participate in the synthesis and accumulation of carotenoids in loquat fruit. In summary, this study provides comprehensive information about the SBP-box gene family in the loquat, and identified two EjSBP genes as candidates involved in carotenoid synthesis and accumulation during loquat fruit development.

Genome-Wide Identification and Expression Analysis of the SQUAMOSA Promoter-Binding Protein-like (SPL) Transcription Factor Family in Catalpabungei

As a plant-specific transcription factor, the SPL gene family plays a critical role in plant growth and development. Although the SPL gene family has been identified in diverse plant species, there have been no genome-wide identification or systematic study reports on the SPL gene family in Catalpa bungei. In this study, we identified 19 putative SPL gene family members in the C. bungei genome. According to the phylogenetic relationship, they can be divided into eight groups, and the genes in the same group have a similar gene structure and conserved motifs. Synteny analysis showed that fragment duplication played an important role in the expansion of the CbuSPL gene family. At the same time, CbuSPL genes have cis-acting elements and functions related to light response, hormone response, growth and development, and stress response. Tissue-specific expression and developmental period-specific expression analysis showed that CbuSPL may be involved in flowering initiation and development, flowering transition, and leaf development. In addition, the ectopic expression of CbuSPL4 in Arabidopsis confirmed that it can promote early flowering and induce the expression of related flowering genes. These systematic research results will lay a foundation for further study on the functional analysis of SPL genes in C. bungei.

Genome-wide identification and analysis of the evolution and expression pattern of the SBP gene family in two Chimonanthus species

BACKGROUND

Chimonanthus praecox and Chimonanthus salicifolius are closely related species that diverged approximately six million years ago. While both C. praecox and C. salicifolius could withstand brief periods of low temperatures of - 15 °C. Their flowering times are different, C. praecox blooms in early spring, whereas C. salicifolius blooms in autumn. The SBP-box (SQUAMOSA promoter-binding protein) is a plant-specific gene family that plays a crucial vital role in regulating plant flowering. Although extensively studied in various plants, the SBP gene family remains uncharacterized in Calycanthaceae.

METHODS AND RESULTS

We conducted genome-wide identification of SBP genes in both C. praecox and C. salicifolius and comprehensively characterized the chromosomal localization, gene structure, conserved motifs, and domains of the identified SBP genes. In total, 15 and 18 SBP genes were identified in C. praecox and C. salicifolius, respectively. According to phylogenetic analysis, the SBP genes from Arabidopsis, C. praecox, and C. salicifolius were clustered into eight groups. Analysis of the gene structure and conserved protein motifs showed that SBP proteins of the same subfamily have similar motif structures. The expression patterns of SBP genes were analyzed using transcriptome data. The results revealed that more than half of the genes exhibited lower expression levels in leaves than in flowers, suggesting their potential involvement in the flower development process and may be linked to the winter and autumn flowering of C. praecox and C. salicifolius.

CONCLUSION

Thirty-three SBPs were identified in C. praecox and C. salicifolius. The evolutionary characteristics and expression patterns were examined in this study. These results provide valuable information to elucidate the evolutionary relationships of the SBP family and help determine the functional characteristics of the SBP genes in subsequent studies.

Comparative analysis of the Squamosa Promoter Binding-Like (SPL) gene family in Nicotiana benthamiana and Nicotiana tabacum

SQUAMOSA PROMOTER BINDING-LIKE (SPL) proteins constitute a large family of transcription factors known to play key roles in growth and developmental processes, including juvenile-to-adult and vegetative-to-reproductive phase transitions. This makes SPLs interesting targets for precision breeding in plants of the Nicotiana genus used as e.g. recombinant biofactories. We report the identification of 49 SPL genes in Nicotiana tabacum cv. K326 and 43 SPL genes in Nicotiana benthamiana LAB strain, which were classified into eight phylogenetic groups according to the SPL classification in Arabidopsis. Exon-intron gene structure and DNA-binding domains were highly conserved between homeologues and orthologues. Thirty of the NbSPL genes and 33 of the NtSPL genes were found to be possible targets of microRNA 156. The expression of SPL genes in leaves was analysed by RNA-seq at three different stages, revealing that genes not under miR156 control were in general constitutively expressed at high levels, whereas miR156-regulated genes showed lower expression, often developmentally regulated. We selected the N. benthamiana SPL13_1a gene as target for a CRISPR/Cas9 knock-out experiment. We show here that a full knock-out in this single gene leads to a significant delay in flowering time, a trait that could be exploited to increase biomass for recombinant protein production.

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.

Genome-Wide Identification and Expression Analysis of the SPL Gene Family in Three Orchids

SPL transcription factors regulate important processes such as plant growth and development, metabolic regulation, and abiotic stress. They play crucial roles in the development of flower organs. However, little is known about the characteristics and functions of the SPLs in the Orchidaceae. In this study, Cymbidium goeringii Rchb. f., Dendrobium chrysotoxum Lindl., and Gastrodia elata BI. were used as research objects. The SPL gene family of these orchids was analyzed on a genome-wide scale, and their physicochemical properties, phylogenetic relationships, gene structures, and expression patterns were studied. Transcriptome and qRT-PCR methods were combined to investigate the regulatory effect of SPLs on the development of flower organs during the flowering process (bud, initial bloom, and full bloom). This study identifies a total of 43 SPLs from C. goeringii (16), D. chrysotoxum (17), and G. elata (10) and divides them into eight subfamilies according to the phylogenetic tree. Most SPL proteins contained conserved SBP domains and complex gene structures; half of the genes had introns longer than 10 kb. The largest number and variety of cis-acting elements associated with light reactions were enriched, accounting for about 45% of the total (444/985); 13/43 SPLs contain response elements of miRNA156. GO enrichment analysis showed that the functions of most SPLs were mainly enriched in the development of plant flower organs and stems. In addition, expression patterns and qRT-PCR analysis suggested the involvement of SPL genes in the regulation of flower organ development in orchids. There was little change in the expression of the CgoSPL in C. goeringii, but DchSPL9 and GelSPL2 showed significant expression during the flowering process of D. chrysotoxum and G. elata, respectively. In summary, this paper provides a reference for exploring the regulation of the SPL gene family in orchids.

Genome-Wide Identification and Characterization of the SBP Gene Family in Passion Fruit (Passiflora edulis Sims)

The SQUAMOSA promoter binding proteins (SBPs) gene family plays important roles in plant growth and development. The SBP gene family has been identified and reported in many species, but it has not been well studied in passion fruit. In this study, a total of 14 SBP genes were identified in passion fruit and named from PeSBP1 to PeSBP14 based on their chromosomal distribution. The phylogenetic tree, gene structure, conserved motifs, collinearity analysis, and expression patterns of the identified SBP members were analyzed. We classified the PeSBP genes into eight groups (I to VIII) according to the phylogenetic tree, gene structure, and conserved motifs. Synteny analysis found that 5 homologous gene pairs existed in PeSBP genes and 11 orthologous gene pairs existed between passion fruit and Arabidopsis. Synonymous nucleotide substitution analysis showed that the PeSBP genes were under strong negative selection. The expression pattern of PeSBP genes in seed, root, leaf, and flower showed that nine of the PeSBP genes displayed high expression in the leaf and the flower. The expression patterns of PeSBP3/6/8/9/10 were further detected by qRT-PCR. In addition, differences in the expression levels occurred for each gene in the different flower organs and at the different developmental stages. There were large differences among SBPs based on transcriptional levels under cold, heat, salt, and osmotic stress conditions. Altogether, this study provides an overview of SBP genes in passion fruit and lays the foundation for further functional analysis.

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.

Genome-wide analysis and molecular dissection of the SPL gene family in Fraxinus mandshurica

BACKGROUND

SQUAMOSA promoter binding protein-like (SPL) is a unique family of transcription factors in plants, which is engaged in regulating plant growth and development, physiological and biochemical processes. Fraxinus mandshurica is an excellent timber species with a wide range of uses in northeastern China and enjoys a high reputation in the international market. SPL family analysis has been reported in some plants while SPL family analysis of Fraxinus mandshurica has not been reported.

RESULTS

We used phylogeny, conserved motifs, gene structure, secondary structure prediction, miR156 binding sites, promoter cis elements and GO annotation to systematically analyze the FmSPLs family. This was followed by expression analysis by subcellular localization, expression patterns at various tissue sites, abiotic stress and hormone induction. Because FmSPL2 is highly expressed in flowers it was selected to describe the SPL gene family of Fraxinus mandshurica by ectopic expression. Among them, 10 FmSPL genes that were highly expressed at different loci were selected for expression analysis under abiotic stress (NaCl and Cold) and hormone induction (IAA and ABA). These 10 FmSPL genes showed corresponding trends in response to both abiotic stress and hormone induction. We showed that overexpression of FmSPL2 in transgenic Nicotiana tabacum L. resulted in taller plants, shorter root length, increased root number, rounded leaves, and earlier flowering time.

CONCLUSIONS

We identified 36 SPL genes, which were classified into seven subfamilies based on sequence analysis. FmSPL2 was selected for subsequent heterologous expression by analysis of expression patterns in various tissues and under abiotic stress and hormone induction, and significant phenotypic changes were observed in the transgenic Nicotiana tabacum L. These results provide insight into the evolutionary origin and biological significance of plant SPL. The aim of this study was to lay the foundation for the genetic improvement of Fraxinus mandshurica and the subsequent functional analysis of FmSPL2.

Genome-Wide Analysis of the SPL Gene Family and Expression Analysis during Flowering Induction in Prunus × yedoensis 'Somei-yoshino'

SQUAMOSA Promoter-Binding Protein-Like (SPL) genes encode plant-specific transcription factors which bind to the SQUAMOSA promoter of the MADS-box genes to regulate its expression. It plays important regulatory roles in floral induction and development, fertility, light signals and hormonal transduction, and stress response in plants. In this study, 32 PySPL genes with complete SBP (squamosa promoter binding protein) conserved domain were identified from the genome of Prunus × yedoensis ‘Somei-yoshino’ and analyzed by bioinformatics. 32 PySPLs were distributed on 13 chromosomes, encoding 32 PySPL proteins with different physical and chemical properties. The phylogenetic tree constructed with Arabidopsis thaliana and Oryza sativa can be divided into 10 subtribes, indicating PySPLs of different clusters have different biological functions. The conserved motif prediction showed that the number and distribution of motifs on each PySPL is varied. The gene structure analysis revealed that PySPLs harbored exons ranging from 2 to 10. The predictive analysis of acting elements showed that the promoter of PySPLs contain a large number of light-responsive elements, as well as response elements related to hormone response, growth and development and stress response. The analysis of the PySPLs expressions in flower induction and flower organs based on qRT-PCR showed that PySPL06/22 may be the key genes of flower development, PySPL01/06 and PySPL22 may play a role in the development of sepal and pistil, respectively. The results provide a foundation for the study of SPL transcription factors of Prunus × yedoensis ‘Somei-yoshino’ and provide more reference information of the function of SPL gene in flowering.

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.

Nitric oxide enhances resistance of Pleurotus eryngii to cadmium stress by alleviating oxidative damage and regulating of short-chain dehydrogenase/reductase family

The function and mechanism of nitric oxide (NO) in regulating Pleurotus eryngii biological response to cadmium (Cd) stress was evaluated by using anti-oxidation and short-chain dehydrogenase/reductase (SDR) family analysis. The fresh biomass of P. eryngii mycelia sharply decreased after treatment with 50 µM Cd; the lipid peroxidation and H₂O₂ accumulation in P. eryngii were found responsible for it. Proper exogenous supply of NO (150 µM SNP) alleviated the oxidative damage induced by Cd stress in P. eryngii, which reduced the accumulation of thiobarbituric acid reactive substances (TBARS) and H₂O₂. The activities of antioxidant enzymes (superoxide dismutase, peroxidase) were significantly increased to deal with Cd stress when treated with SNP (150 µM), and the content of proline was also closely related to NO-mediated reduction of Cd toxicity. Moreover, SDR family members were widely involved in the response to Cd stress, especially PleSCH70 gene was observed for the first time in participating in NO-mediated enhancement of Cd tolerance in P. eryngii. Taken together, this study provides new insights in understanding the tolerance mechanisms of P. eryngii to heavy metal and lays a foundation for molecular breeding of P. eryngii to improve its tolerance to environmental stress.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family and its important role in salt stress in Medicago sativa L

BACKGROUND

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are widely present in plants and are involved in signal transduction, the stress response and development. The SPL gene family has been characterized in several model species, such as A. thaliana and G. max. However, there is no in-depth analysis of the SPL gene family in forage, especially alfalfa (Medicago sativa L.), one of the most important forage crops worldwide.

RESULT

In total, 76 putative MsSPL genes were identified in the alfalfa genome with an uneven distribution. Based on their identity and gene structure, these MsSPLs were divided into eight phylogenetic groups. Seventy-three MsSPL gene pairs arose from segmental duplication events, and the MsSPLs on the four subgenomes of individual chromosomes displayed high collinearity with the corresponding M. truncatula genome. The prediction of the cis-elements in the promoter regions of the MsSPLs detected two copies of ABA (abscisic acid)-responsive elements (ABREs) on average, implying their potential involvement in alfalfa adaptation to adverse environments. The transcriptome sequencing of MsSPLs in roots and leaves revealed that 54 MsSPLs were expressed in both tissues. Upon salt treatment, three MsSPLs (MsSPL17, MsSPL23 and MsSPL36) were significantly regulated, and the transcription level of MsSPL36 in leaves was repressed to 46.6% of the control level.

CONCLUSION

In this study, based on sequence homology, we identified 76 SPL genes in the alfalfa. The SPLs with high identity shared similar gene structures and motifs. In total, 71.1% (54 of 76) of the MsSPLs were expressed in both roots and leaves, and the majority (74.1%) preferred underground tissues to aerial tissues. MsSPL36 in leaves was significantly repressed under salt stress. These findings provide comprehensive information regarding the SPB-box gene family for improve alfalfa tolerance to high salinity.

Identification and regulatory network analysis of SPL family transcription factors in Populus euphratica Oliv. heteromorphic leaves

The SQUAMOSA promoter-binding protein-like (SPL) family play a key role in guiding the switch of plant growth from juvenile to adult phases. Populus euphratica Oliv. exhibit typical heterophylly, and is therefore an ideal model for studying leaf shape development. To investigate the role and regulated networks of SPLs in the morphogenesis of P. euphratica heteromorphic leaves. In this study, 33 P. euphratica SPL (PeuSPL) genes were identified from P. euphratica genome and transcriptome data. Phylogenetic analysis depicted the classification of these SPL genes into two subgroups. The expression profiles and regulatory networks of P. euphratica SPL genes analysis displayed that major P. euphratica SPL family members gradually increases from linear to broad-ovate leaves, and they were involved in the morphogenesis regulation, stress response, transition from vegetative to reproductive growth, photoperiod, and photosynthesis etc. 14 circRNAs, and 33 lncRNAs can promote the expression of 12 of the P. euphratica SPLs by co-decoying miR156 in heteromorphic leaf morphogenesis. However, it was found that the effect of PeuSPL2-4 and PeuSPL9 in leaf shape development was contrasting to their homologous genes of Arabidopsis. Therefore, it was suggested that the SPL family were evolutionarily conserved for regulation growth, but were varies in different plant for regulation of the organ development.

The Mulberry SPL Gene Family and the Response of MnSPL7 to Silkworm Herbivory through Activating the Transcription of MnTT2L2 in the Catechin Biosynthesis Pathway

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes, as unique plant transcription factors, play important roles in plant developmental regulation and stress response adaptation. Although mulberry is a commercially valuable tree species, there have been few systematic studies on SPL genes. In this work, we identified 15 full-length SPL genes in the mulberry genome, which were distributed on 4 Morus notabilis chromosomes. Phylogenetic analysis clustered the SPL genes from five plants (Malus × domestica Borkh, Populus trichocarpa, M. notabilis, Arabidopsis thaliana, and Oryza sativa) into five groups. Two zinc fingers (Zn1 and Zn2) were found in the conserved SBP domain in all of the MnSPLs. Comparative analyses of gene structures and conserved motifs revealed the conservation of MnSPLs within a group, whereas there were significant structure differences among groups. Gene quantitative analysis showed that the expression of MnSPLs had tissue specificity, and MnSPLs had much higher expression levels in older mulberry leaves. Furthermore, transcriptome data showed that the expression levels of MnSPL7 and MnSPL14 were significantly increased under silkworm herbivory. Molecular experiments revealed that MnSPL7 responded to herbivory treatment through promoting the transcription of MnTT2L2 and further upregulating the expression levels of catechin synthesis genes (F3′H, DFR, and LAR).

Genome-Wide Discovery of miRNAs with Differential Expression Patterns in Responses to Salinity in the Two Contrasting Wheat Cultivars

Salinity is a serious environmental issue. It has a substantial effect on crop yield, as many crop species are sensitive to salinity due to climate change, and it impact is continuing to increase. Plant microRNAs (miRNAs) contribute to salinity stress response in bread wheat. However, the underlying molecular mechanisms by which miRNAs confer salt tolerance in wheat are unclear. We conducted a genome-wide discovery study using Illumina high throughput sequencing and comprehensive in silico analysis to obtain insight into the underlying mechanisms by which small RNAs confer tolerance to salinity in roots of two contrasting wheat cvv., namely Suntop (salt-tolerant) and Sunmate (salt-sensitive). A total of 191 microRNAs were identified in both cultivars, consisting of 110 known miRNAs and 81 novel miRNAs; 181 miRNAs were shared between the two cultivars. The known miRNAs belonged to 35 families consisted of 23 conserved and 12 unique families. Salinity stress induced 43 and 75 miRNAs in Suntop and Sunmate, respectively. Among them, 14 and 29 known and novel miRNAs were expressed in Suntop and 37 and 38 in Sunmate. In silico analysis revealed 861 putative target mRNAs for the 75 known miRNAs and 52 putative target mRNAs for the 15 candidate novel miRNAs. Furthermore, seven miRNAs including tae-miR156, tae-miR160, tae-miR171a-b, tae-miR319, tae-miR159a-b, tae-miR9657 and novel-mir59 that regulate auxin responsive-factor, SPL, SCL6, PCF5, R2R3 MYB, and CBL-CIPK, respectively, were predicted to contribute to salt tolerance in Suntop. This information helps further our understanding of how the molecular mechanisms of salt tolerance are mediated by miRNAs and may facilitate the genetic improvement of wheat cultivars.

Genome-Wide Characterization of SPL Gene Family in Codonopsis pilosula Reveals the Functions of CpSPL2 and CpSPL10 in Promoting the Accumulation of Secondary Metabolites and Growth of C. pilosula Hairy Root

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors play critical roles in regulating diverse aspects of plant growth and development, including vegetative phase change, plant architecture, anthocyanin accumulation, lateral root growth, etc. In the present study, 15 SPL genes were identified based on the genome data of Codonopsis pilosula, a well-known medicinal plant. Phylogenetic analysis clustered CpSPLs into eight groups (G1-G8) along with SPLs from Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Physcomitrella patens. CpSPLs in the same group share similar gene structure and conserved motif composition. Cis-acting elements responding to light, stress and phytohormone widely exist in their promoter regions. Our qRT-PCR results indicated that 15 CpSPLs were differentially expressed in different tissues (root, stem, leaf, flower and calyx), different developmental periods (1, 2 and 3 months after germination) and various conditions (NaCl, MeJA and ABA treatment). Compared with the control, overexpression of CpSPL2 or CpSPL10 significantly promoted not only the growth of hairy roots, but also the accumulation of total saponins and lobetyolin. Our results established a foundation for further investigation of CpSPLs and provided novel insights into their biological functions. As far as we know, this is the first experimental research on gene function in C. pilosula.

Transcription factor SmSPL7 promotes anthocyanin accumulation and negatively regulates phenolic acid biosynthesis in Salvia miltiorrhiza

Plant-specific SQUAMOSA promoter-binding protein-like (SPL) transcription factors play critical regulatory roles during plant growth and development. However, the functions of SPLs in Salvia miltiorrhiza (SmSPLs; a model medicinal plant) have not been reported. Here, the expression patterns and functions of SmSPL7 were characterized in S. miltiorrhiza. SmSPL7 was expressed in all parts of S. miltiorrhiza, with the highest expression level in the leaves, and could be inhibited by multiple hormones, including methyl jasmonate, auxin, abscisic acid, and gibberellin. SmSPL7 is localized within the nucleus and exhibits robust transcriptional activation activity. Transgenic lines overexpressing SmSPL7 demonstrated pronounced growth inhibition, accompanied by increased anthocyanin accumulation via the genetic activation of the anthocyanin biosynthesis pathway. However, SmSPL7 overexpression significantly decreased salvianolic acid B (SalB) production by inhibiting the transcripts of genes implicated in its biosynthesis pathway. Further analysis indicated that SmSPL7 directly binds to SmTAT1 and Sm4CL9 promoters and blocks their expression to inhibit the biosynthesis of SalB. Taken together, these results indicate that SmSPL7 is a negative regulator of SalB biosynthesis but positively regulates anthocyanin accumulation in S. miltiorrhiza. These findings provide new insights into the functionality of the SPL family while establishing an important foundation for further uncovering the crucial roles of SmSPL7 in the growth of S. miltiorrhiza.

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.

Comprehensive Analysis of the SBP Family in Blueberry and Their Regulatory Mechanism Controlling Chlorophyll Accumulation

SQUAMOSA Promoter Binding Protein (SBP) family genes act as central players to regulate plant growth and development with functional redundancy and specificity. Addressing the diversity of the SBP family in crops is of great significance to precisely utilize them to improve agronomic traits. Blueberry is an important economic berry crop. However, the SBP family has not been described in blueberry. In the present study, twenty VcSBP genes were identified through data mining against blueberry transcriptome databases. These VcSBPs could be clustered into eight groups, and the gene structures and motif compositions are divergent among the groups and similar within each group. The VcSBPs were differentially expressed in various tissues. Intriguingly, 10 VcSBPs were highly expressed at green fruit stages and dramatically decreased at the onset of fruit ripening, implying that they are important regulators during early fruit development. Computational analysis showed that 10 VcSBPs were targeted by miR156, and four of them were further verified by degradome sequencing. Moreover, their functional diversity was studied in Arabidopsis. Noticeably, three VcSBPs significantly increased chlorophyll accumulation, and qRT-PCR analysis indicated that VcSBP13a in Arabidopsis enhanced the expression of chlorophyll biosynthetic genes such as AtDVR, AtPORA, AtPORB, AtPORC, and AtCAO. Finally, the targets of VcSBPs were computationally identified in blueberry, and the Y1H assay showed that VcSBP13a could physically bind to the promoter region of the chlorophyll-associated gene VcLHCB1. Our findings provided an overall framework for individually understanding the characteristics and functions of the SBP family in blueberry.

Hierarchical Action of Mulberry miR156 in the Vegetative Phase Transition

The vegetative phase transition is a prerequisite for flowering in angiosperm plants. Mulberry miR156 has been confirmed to be a crucial factor in the vegetative phase transition in Arabidopsis thaliana. The over-expression of miR156 in transgenic Populus × canadensis dramatically prolongs the juvenile phase. Here, we find that the expression of mno-miR156 decreases with age in all tissues in mulberry, which led us to study the hierarchical action of miR156 in mulberry. Utilizing degradome sequencing and dual-luciferase reporter assays, nine MnSPLs were shown to be directly regulated by miR156. The results of yeast one-hybrid and dual-luciferase reporter assays also revealed that six MnSPLs could recognize the promoter sequences of mno-miR172 and activate its expression. Our results demonstrate that mno-miR156 performs its role by repressing MnSPL/mno-miR172 pathway expression in mulberry. This work uncovered a miR156/SPLs/miR172 regulation pathway in the development of mulberry and fills a gap in our knowledge about the molecular mechanism of vegetative phase transition in perennial woody plants.

Pleiotropic function of the SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE gene TaSPL14 in wheat plant architecture

MAIN CONCLUSION

The function of SQUAMOSA PROMOTER-BINDING PROTEIN-BOX gene TaSPL14 in wheat is similar to that of OsSPL14 in rice in regulating plant height, panicle length, spikelet number, and thousand-grain weight of wheat, but differs during tiller development. TaSPL14 may regulate spike development via ethylene-response gene EIN3-LIKE 1 (TaEIL1), ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR 2.11 (TaRAP2.11), and ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR 1 (TaERF1), but not DENSE AND ERECT PANICLE 1 (TaDEP1) in wheat. The SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE gene OsSPL14 from rice is considered to be a major determinant of ideal plant architecture consisting of few unproductive tillers, more grains per spike, and high resistance of stems to lodging. However, the function of its orthologous gene, TaSPL14, in wheat is unknown. Here, we reported the functional similarities and differences between TaSPL14 and OsSPL14. Similar to OsSPL14 knock-outs in rice, wheat TaSPL14 knock-out plants exhibited decreased plant height, panicle length, spikelet number, and thousand-grain weight. In contrast to OsSPL14, however, TaSPL14 did not affect tiller number. Transcriptome analysis revealed that the expression of genes related to ethylene response was significantly decreased in young spikes of TaSPL14 knock-out lines as compared with wild type. TaSPL14 directly binds to the promoters of the ethylene-response genes TaEIL1, TaRAP2.11, and TaERF1, and promotes their expression, suggesting that TaSPL14 might regulate wheat spike development via the ethylene-response pathway. The elucidation of TaSPL14 will contribute to understanding of the molecular mechanisms that underlie wheat plant architecture.

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.

CRISPR-based assessment of genomic structure in the conserved SQUAMOSA promoter-binding-like gene clusters in rice

Although SQUAMOSA promoter-binding-like (SPL) transcription factors are important regulators of development in rice (Oryza sativa), prior assessments of the SPL family have been limited to single genes. A functional comparison across the full gene family in standardized genetic backgrounds has not been reported previously. Here, we demonstrate that the SPL gene family in rice is enriched due to the most recent whole genome duplication (WGD). Notably, 10 of 19 rice SPL genes (52%) cluster in four units that have persisted for at least 50 million years. We show that SPL gene grouping and retention following WGD is widespread in angiosperms, suggesting the conservatism and importance of this gene arrangement. We used Cas9 editing to generate transformation lines for all 19 SPL genes in a common set of backgrounds, and found that knockouts of 14 SPL genes exhibited defects in plant height, 10 exhibited defects in panicle size, and nine had altered grain lengths. We observed subfunctionalization of genes in the paleoduplicated pairs, but little evidence of neofunctionalization. Expression of OsSPL3 was negatively correlated with that of its closest neighbor in its synteny group, OsSPL4, and its sister paired gene, OsSPL12, in the opposing group. Nucleotide diversity was lower in eight of the nine singleton genes in domesticated rice, relative to wild rice, whereas the reverse was true for the paired genes. Together, these results provide functional information on eight previously unexamined OsSPL family members and suggest that paleoduplicate pair redundancy benefits plant survival and innovation.

Molecular characterization and expression analysis of the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE gene family in Paeonia suffruticosa

A total of 16 PsSPL genes were identified in tree peony. PsSPLs potentially regulated flowering time, lateral bud and seed development, and the juvenile-to-adult phase transition. SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors are important for plant growth and development. Here, we report the identification of 16 full-length PsSPLs in tree peony (Peaonia suffruticosa Andr.) and 9 PsSPLs that have miR156 target sites. Phylogenetic analysis of the relationship of SPLs in P. suffruticosa and Arabidopsis suggested that they can be classified into six groups, and PsSPLs were highly correlated with Arabidopsis SPLs counterparts in the same group. Cis-element of promoter region analysis suggested that PsSPL genes play roles in physiological processes and developmental events. Expression analysis indicated that most PsSPL genes exhibited high expression levels in the tissues and organs examined here. The increasing expression levels of PsSPL1, PsSPL2, PsSPL8, PsSPL9, PsSPL12, and PsSPL16, and decreasing expression levels of PsSPL1A and PsSPL1B in buds over time suggested that they were probably regulated by the juvenile-to-adult phase transition. In addition, the expression profiles of PsSPL genes in different developmental buds and seeds suggested that PsSPL2, PsSPL3, PsSPL9, PsSPL10, PsSPL13, and PsSPL13A were important genes for regulating the flowering time of the tree peony; PsSPL2 and PsSPL8 might play a role in suppressing lateral bud development, and PsSPL2, PsSPL13, and PsSPL14 positively controlled grain size and number, and pod branching. These results provide a foundation for future functional analysis of PsSPL genes in tree peony growth and development.

Global analysis of SBP gene family in Brachypodium distachyon reveals its association with spike development

SQUAMOSA-promoter binding like proteins (SBPs/SPLs) are plant specific transcription factors targeted by miR156 and involved in various biological pathways, playing multi-faceted developmental roles. This gene family is not well characterized in Brachypodium. We identified a total of 18 SBP genes in B. distachyon genome. Phylogenetic analysis revealed that SBP gene family in Brachypodium expanded through large scale duplication. A total of 10 BdSBP genes were identified as targets of miR156. Transcript cleavage analysis of selected BdSBPs by miR156 confirmed their antagonistic connection. Alternative splicing was observed playing an important role in BdSBPs and miR156 interaction. Characterization of T-DNA Bdsbp9 mutant showed reduced plant growth and spike length, reflecting its involvement in the spike development. Expression of a majority of BdSBPs elevated during spikelet initiation. Specifically, BdSBP1 and BdSBP3 differentially expressed in response to vernalization. Differential transcript abundance of BdSBP1, BdSBP3, BdSBP8, BdSBP9, BdSBP14, BdSBP18 and BdSBP23 genes was observed during the spike development under high temperature. Co-expression network, protein-protein interaction and biological pathway analysis indicate that BdSBP genes mainly regulate transcription, hormone, RNA and transport pathways. Our work reveals the multi-layered control of SBP genes and demonstrates their association with spike development and temperature sensitivity in Brachypodium.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family in wheat

BACKGROUND

Members of the plant-specific SPL gene family (squamosa promoter-binding protein -like) contain the SBP conserved domain and are involved in the regulation of plant growth and development, including the development of plant flowers and plant epidermal hair, the plant stress response, and the synthesis of secondary metabolites. This family has been identified in various plants. However, there is no systematic analysis of the SPL gene family at the genome-wide level of wheat.

RESULTS

In this study, 56 putative TaSPL genes were identified using the comparative genomics method; we renamed them TaSPL001 - TaSPL056 on their chromosomal distribution. According to the un-rooted neighbor joining phylogenetic tree, gene structure and motif analyses, the 56 TaSPL genes were divided into 8 subgroups. A total of 81 TaSPL gene pairs were designated as arising from duplication events and 64 interacting protein branches were identified as involve in the protein interaction network. The expression patterns of 21 randomly selected TaSPL genes in different tissues (roots, stems, leaves and inflorescence) and under 4 treatments (abscisic acid, gibberellin, drought and salt) were detected by quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The wheat genome contains 56 TaSPL genes and those in same subfamily share similar gene structure and motifs. TaSPL gene expansion occurred through segmental duplication events. Combining the results of transcriptional and qRT-PCR analyses, most of these TaSPL genes were found to regulate inflorescence and spike development. Additionally, we found that 13 TaSPLs were upregulated by abscisic acid, indicating that TaSPL genes play a positive role in the abscisic acid-mediated pathway of the seedling stage. This study provides comprehensive information on the SPL gene family of wheat and lays a solid foundation for elucidating the biological functions of TaSPLs and improvement of wheat yield.

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

BACKGROUND

SPL (SQUAMOSA-promoter binding protein-like) proteins form a large family of plant-specific transcription factors that play essential roles in various aspects of plant growth and development. They are potentially important candidates for genetic improvement of agronomic traits. However, there were limited information about the SPL genes in Jatropha curcas, an important biofuel plant.

RESULTS

In Jatropha, 15 JcSPL genes were identified. Phylogenetic analysis revealed that most of the JcSPLs were closely related to SPLs from woody plant rather than herbaceous plant and distantly related to monocotyledon SPLs. Gene structure, conserved motif and repetitive sequence analysis indicated diverse and specific functions of some JcSPL genes. By combination of target prediction and degradome sequencing analysis, 10 of the 15 JcSPLs were shown to be targets of JcmiR156. Quantitative PCR analysis showed diversified spatial-temporal expression patterns of JcSPLs. It is interesting that the expression levels of JcSPL3 were the highest in all tissues examined in 7- or 10-year-old plants and exhibited increasing trend with plant age, suggesting its important role in the regulation of age development in Jatropha. Overexpression of JcSPL3 in Arabidopsis resulted in earlier flowering time, shorter silique length and reduced biomass of roots.

CONCLUSIONS

Through comprehensive and systematic analysis of phylogenetic relationships, conserved motifs, gene structures, chromosomal locations, repetitive sequence and expression patterns, 15 JcSPL genes were identified in Jatropha and characterized in great detail. These results provide deep insight into the evolutionary origin and biological significance of plant SPLs and lay the foundation for further functional characterization of JcSPLs with the purpose of genetic improvement in Jatropha.

Roles of transcription factor SQUAMOSA promoter binding protein-like gene family in papaya (Carica papaya) development and ripening

SQUAMOSA promoter binding protein-like (SPL) family plays vital regulatory roles in plant growth and development. The SPL family in climacteric fruit Carica papaya has not been reported. This study identified 14 papaya SPLs (CpSPL) from papaya genome and analyzed their sequence features, phylogeny, intron/exon structure, conserved motif, miR156-mediated posttranscriptional regulation, and expression patterns. 14 CpSPLs were clustered into 8 groups, and two distinct expression patterns were revealed for miR156-targeted and nontargeted CpSPLs in different tissues and fruit development stages. The expression changes of CpSPLs in ethephon and 1-MCP treated fruit during ripening suggested that the CpSPLs guided by CpmiR156 play crucial roles in ethylene signaling pathway. This study sheds light on the new function of SPL family in fruit development and ripening, providing insights on understanding evolutionary divergence of the members of SPL family among plant species.

Potential function of CbuSPL and gene encoding its interacting protein during flowering in Catalpa bungei

BACKGROUND

"Bairihua", a variety of the Catalpa bungei, has a large amount of flowers and a long flowering period which make it an excellent material for flowering researches in trees. SPL is one of the hub genes that regulate both flowering transition and development.

RESULTS

SPL homologues CbuSPL9 was cloned using degenerate primers with RACE. Expression studies during flowering transition in "Bairihua" and ectopic expression in Arabidopsis showed that CbuSPL9 was functional similarly with its Arabidopsis homologues. In the next step, we used Y2H to identify the proteins that could interact with CbuSPL9. HMGA, an architectural transcriptional factor, was identified and cloned for further research. BiFC and BLI showed that CbuSPL9 could form a heterodimer with CbuHMGA in the nucleus. The expression analysis showed that CbuHMGA had a similar expression trend to that of CbuSPL9 during flowering in "Bairihua". Intriguingly, ectopic expression of CbuHMGA in Arabidopsis would lead to aberrant flowers, but did not effect flowering time.

CONCLUSIONS

Our results implied a novel pathway that CbuSPL9 regulated flowering development, but not flowering transition, with the participation of CbuHMGA. Further investments need to be done to verify the details of this pathway.

Genome-wide identification and characterization, phylogenetic comparison and expression profiles of SPL transcription factor family in B. juncea (Cruciferae)

SQUAMOSA promoter-binding protein-like (SPL), as plant specific transcription factors, is involved in many plant growth and development processes. However, there is less systematical study for SPL transcription factor in B. juncea (Cruciferae). Here, a total of 59 SPL genes classified into eight phylogenetic groups were identified in B. juncea, highly conserved within each ortholog were also found based on gene structure, conserved motif, as well as clustering level. In addition, clustering of SPL domain showed that two zinc finger-like structures and NLS segments were identified in almost of BjuSPLs. Analyzed of putative cis-elements for BjuSPLs demonstrated that SPL transcription factors were involved in adverse environmental changes, such as light, plant stresses and phytohormones response. Expression analysis showed that differentially expressed SPL genes were identified in flower and stem development of Cruciferae; such as BjuSPL3a-B, BjuSPL2b_B and BjuSPL2c_A were significantly expressed in flower; BjuSPL 3b_B and BjuSPL10a_A were significantly expressed in stem node (VP: vegetative period). Moreover, 28 of the 59 BjuSPLs were found involved in their posttranscriptional regulation targeted by miR156. We demonstrated that miR156 negatively regulated BjuSPL10a_A and BjuSPL3b_B to act for stem development in B. juncea.

Roles of the SPL gene family and miR156 in the salt stress responses of tamarisk (Tamarix chinensis)

BACKGROUND

Accumulating evidences show that SPLs are crucial regulators of plant abiotic stress tolerance and the highly conserved module miR156/SPL appears to balance plant growth and stress responses. The halophyte Tamarix chinensis is highly resistant to salt tress. SPLs of T. chinensis (TcSPLs) and theirs roles in salt stress responses remain elusive.

RESULTS

In this study, we conducted a systematic analysis of the TcSPLs gene family including 12 members belonging to 7 groups. The physicochemical properties and conserved motifs showed divergence among groups and similarity in each group. The microRNA response elements (MREs) are conserved in location and sequence, with the exception of first MRE within TcSPL5. The miR156-targeted SPLs are identified by dual-luciferase reporter assay of MRE-miR156 interaction. The digital expression gene profiles cluster suggested potential different functions of miR156-targeted SPLs vs non-targeted SPLs in response to salt stress. The expression patterns analysis of miR156-targeted SPLs with a reverse expression trend to TcmiR156 suggested 1 h (salt stress time) could be a critical time point of post-transcription regulation in salt stress responses.

CONCLUSIONS

Our work demonstrated the post-transcription regulation of miR156-targeted TcSPLs and transcription regulation of non-targeted TcSPLs in salt stress responses, and would be helpful to expound the miR156/SPL-mediated molecular mechanisms underlying T. chinensis salt stress tolerance.

Comparative genome analysis of the SPL gene family reveals novel evolutionary features in maize

SPLs are plant-specific transcription factors that play important regulatory roles in plant growth and development. Systematic analysis of the SPL family has been performed in numerous plants, such as Arabidopsis, rice, and Populus. However, no comparative analysis has been performed across different species to examine evolutionary features. In this study, we present a comparative analysis of SPLs in different species. The results showed that 84 SPLs of the four species can be divided into six groups according to phylogeny. We found that most of the SPL-containing regions in maize showed extensive conservation with duplicated regions of rice and sorghum. A gene duplication analysis in maize indicated that ZmSPLs showed a significant excess of segmental duplication. The Ka/Ks analysis indicated that 9 out of 18 duplicated pairs in maize experienced positive selection, while SPL gene pairs of rice and sorghum mainly evolved under purifying selection, suggesting novel evolutionary features for ZmSPLs. The 31 ZmSPLs were further analyzed by describing their gene structure, phylogenetic relationships, chromosomal location, and expression, Among the ZmSPLs, 13 were predicated to be targeted by miR156s and involved in drought stress response. These results provide the foundation for future functional analyses of ZmSPLs.

Roles of the GA-mediated SPL Gene Family and miR156 in the Floral Development of Chinese Chestnut (Castanea mollissima)

Chestnut (Castanea mollissima) is a deciduous tree species with major economic and ecological value that is widely used in the study of floral development in woody plants due its monoecious and out-of-proportion characteristics. Squamosa promoter-binding protein-like (SPL) is a plant-specific transcription factor that plays an important role in floral development. In this study, a total of 18 SPL genes were identified in the chestnut genome, of which 10 SPL genes have complementary regions of CmmiR156. An analysis of the phylogenetic tree of the squamosa promoter-binding protein (SBP) domains of the SPL genes of Arabidopsis thaliana, Populus trichocarpa, and C. mollissima divided these SPL genes into eight groups. The evolutionary relationship between poplar and chestnut in the same group was similar. A structural analysis of the protein-coding regions (CDSs) showed that the domains have the main function of SBP domains and that other domains also play an important role in determining gene function. The expression patterns of CmmiR156 and CmSPLs in different floral organs of chestnut were analyzed by real-time quantitative PCR. Some CmSPLs with similar structural patterns showed similar expression patterns, indicating that the gene structures determine the synergy of the gene functions. The application of gibberellin (GA) and its inhibitor (Paclobutrazol, PP₃₃₃) to chestnut trees revealed that these exert a significant effect on the number and length of the male and female chestnut flowers. GA treatment significantly increased CmmiR156 expression and thus significantly decreased the expression of its target gene, CmSPL6/CmSPL9/CmSPL16, during floral bud development. This finding indicates that GA might indirectly affect the expression of some of the SPL target genes through miR156. In addition, RNA ligase-mediated rapid amplification of the 5' cDNA ends (RLM-RACE) experiments revealed that CmmiR156 cleaves CmSPL9 and CmSPL16 at the 10th and 12th bases of the complementary region. These results laid an important foundation for further study of the biological function of CmSPLs in the floral development of C. mollissima.

Genome-Wide Identification and Characterization of SQUAMOSA-Promoter-Binding Protein (SBP) Genes Involved in the Flowering Development of Citrus Clementina

SQUAMOSA-promoter binding protein (SBP)-box genes encode a family of plant-specific transcription factors that play vital roles in plant growth and development. In this study, 15 SBP-box genes were identified and isolated from Citrus clementina (CclSBPs), where 10 of these genes were predicted to be putative targets of Citrus clementina microRNA156 (CclmiR156). The 15 CclSBP genes could be classified into six groups based on phylogenetic analysis, diverse intron–exon structure, and motif prediction, similar to the SQUAMOSA promoter binding protein-like (SPL) gene family of Populus trichocarpa and Arabidopsis thaliana. Furthermore, CclSBPs classified into a group/subgroup have similar gene structures and conserved motifs, implying their functional redundancy. Tissue-specific expression analysis of CclSBPs demonstrated their diversified expression patterns. To further explore the potential role of CclSBPs during floral inductive water deficits, the dynamic changes of the 15 CclSBPs were investigated during floral inductive water deficits, and the results showed that some CclSBPs were associated with floral induction. Among these genes, CclSBP6 was not homologous to the Arabidopsis SBP-box gene family, and CclSBP7 was regulated by being alternatively spliced. Therefore, CclSBP6 and CclSBP7 were genetically transformed in Arabidopsis. Overexpression of the two genes changed the flowering time of Arabidopsis.

Genome-wide identification and expression analysis of the SPL gene family in woodland strawberry Fragaria vesca

SQUAMOSA promoter-binding protein-like (SPL) is a class of plant-specific transcription factors that play critical roles in regulating plant growth and development. However, little systematic research on SPL genes has been conducted in strawberry. In this study, 14 SPL genes were identified in the genome of woodland strawberry (Fragaria vesca), one of the model plants of the family Rosaceae. Chromosome localization analysis indicated that the 14 FvSPL genes were unevenly distributed on six chromosomes. Phylogenetic analysis indicated that the FvSPL proteins could be clustered into six groups (G1 to G6). Genes with similar structure were classified into the same group, implying their functional redundancy. In addition, nine out of the 14 FvSPL genes, belonging to G1, G2, and G5, were found to be the putative targets of FvmiR156 genes. Expression analysis indicated FvSPL genes exhibited highly diverse expression patterns in the tissues and organs examined. The transcript levels of most FvmiR156-targeted FvSPL genes in fruit were lower than those non-miR156-targeted genes. In addition, the expression of the FvmiR156-targeted FvSPL genes decreased during fruit ripening, whereas the expression of FvmiR156 genes increased in fruit during this process. The results provide a foundation for future functional analysis of FvSPL genes in strawberry growth and development.

Genome-wide analysis of the SPL/miR156 module and its interaction with the AP2/miR172 unit in barley

The SQUAMOSA-promoter binding like (SPL) gene family encodes transcription factors that have been shown in many species to influence plant growth and development, but information about these genes in barley (Hordeum vulgare L.) is limited. This study identified 17 barley SPL genes, within eight distinct groups, that are orthologs of SPL genes described in Arabidopsis, wheat, and rice. Sixteen barley SPLs undergo alternative splicing. Seven SPLs contain a putative miR156 target site and the transcript levels of the miR156-targeted HvSPLs (HvSPL3, 13 and 23) were lower in vegetative than in reproductive phase but this was true also for some SPL genes such as HvSPL6 that were not regulated by miR156. Because SPL gene products regulate miR172, which is also involved in floral development, the expression of miR172 was studied. An antagonistic expression pattern of miR156 and miR172b during the vegetative and the reproductive phases signifies their apparent function in barley growth phase transition. Characterization of a barley mir172 mutant having an abnormal, indeterminate spikelet phenotype suggests the possible feedback role of AP2/miR172 module on HvSPL genes. This is the first comprehensive analysis of the miR156/SPL/miR172 axis in barley that provides a basis to elucidate their roles in various biological processes.

Genome-wide identification and characterization of the SBP-box gene family in Petunia

Background

SQUAMOSA PROMOTER BINDING PROTEIN (SBP)-box genes encode a family of plant-specific transcription factors (TFs) that play important roles in many growth and development processes including phase transition, leaf initiation, shoot and inflorescence branching, fruit development and ripening etc. The SBP-box gene family has been identified and characterized in many species, but has not been well studied in Petunia, an important ornamental genus.

Results

We identified 21 putative SPL genes of Petunia axillaris and P. inflata from the reference genome of P. axillaris N and P. inflata S6, respectively, which were supported by the transcriptome data. For further confirmation, all the 21 genes were also cloned from P. hybrida line W115 (Mitchel diploid). Phylogenetic analysis based on the highly conserved SBP domains arranged PhSPLs in eight groups, analogous to those from Arabidopsis and tomato. Furthermore, the Petunia SPL genes had similar exon-intron structure and the deduced proteins contained very similar conserved motifs within the same subgroup. Out of 21 PhSPL genes, fourteen were predicted to be potential targets of PhmiR156/157, and the putative miR156/157 response elements (MREs) were located in the coding region of group IV, V, VII and VIII genes, but in the 3’-UTR regions of group VI genes. SPL genes were also identified from another two wild Petunia species, P. integrifolia and P. exserta, based on their transcriptome databases to investigate the origin of PhSPLs. Phylogenetic analysis and multiple alignments of the coding sequences of PhSPLs and their orthologs from wild species indicated that PhSPLs were originated mainly from P. axillaris. qRT-PCR analysis demonstrated differential spatiotemperal expression patterns of PhSPL genes in petunia and many were expressed predominantly in the axillary buds and/or inflorescences. In addition, overexpression of PhSPL9a and PhSPL9b in Arabidopsis suggested that these genes play a conserved role in promoting the vegetative-to-reproductive phase transition.

Conclusion

Petunia genome contains at least 21 SPL genes, and most of the genes are expressed in different tissues. The PhSPL genes may play conserved and diverse roles in plant growth and development, including flowering regulation, leaf initiation, axillary bud and inflorescence development. This work provides a comprehensive understanding of the SBP-box gene family in Petunia and lays a significant foundation for future studies on the function and evolution of SPL genes in petunia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4537-9) contains supplementary material, which is available to authorized users.

Genome-wide identification and characterization of SPL transcription factor family and their evolution and expression profiling analysis in cotton

Plant specific transcription factors, SQUAMOSA promoter-binding protein-like (SPL), are involved in many biological processes. However, no systematical study has been reported in cotton. In this study, a total of 177 SPL genes were identified, including 29, 30, 59 and 59 SPLs in Gossypium arboreum, G. raimondii, G. barbadense, and G. hirsutum, respectively. These SPL genes were classified into eight phylogenetical groups. The gene structure, conserved motif, and clustering were highly conserved within each orthologs. Two zinc finger-like structures (Cys3His and Cys2HisCys) and NLS segments were existed in all GrSPLs. Segmental duplications play important roles in SPL family expansion, with 20 genes involved in segmental duplications and 2 in tandem duplications, and ten ortholog pairs in syntenic regions between G. raimondii and A. thaliana. Several putative cis-elements, involved in light, stresses and phytohormones response, were found in the promoter regions of GhSPLs, suggesting that plant responses to those environmental changes may be induced through targeting SPL transcription factors. RNA-seq analysis shows that SPL genes were differentially expressed in cotton; some were highly expressed during fiber initiation and early development. Comparing with other plants, SPL genes show subfunctionalization, lost and/or gain functions in cotton during long-term domestication and evolution.

A Global View of Transcriptome Dynamics During Male Floral Bud Development in Populus tomentosa

To obtain a comprehensive overview of the dynamic transcriptome during male floral bud development in Populus tomentosa, high-throughput RNA-seq was conducted during eight flowering-related stages. Among the 109,212 de novo assembled unigenes, 6,959 were differentially expressed during the eight stages. The overrepresented classed of genes identified by Gene Ontology (GO) enrichment included 'response to environmental stimuli' and 'plant-type spore development'. One-third of the differentially expressed genes were transcription factors (TFs). Several genes and gene families were analyzed in depth, including MADS-box TFs, Squamosa promoter binding protein-like family, receptor-like kinases, FLOWERING LOCUS T/TERMINAL-FLOWER-LIKE 1 family, key genes involved in anther and tapetum development, as well as LEAFY, WUSCHEL and CONSTANS. The results provided new insights into the roles of these and other well known gene families during the annual flowering cycle. To explore the mechanisms regulating poplar flowering, a weighted gene co-expression network was constructed using 98 floral-related genes involved in flower meristem identity and flower development. Many modules of co-expressed genes and hub genes were identified, such as APETALA1 and HUA1. This work provides many new insights on the annual flowering cycle in a perennial plant, and a major new resource for plant biology and biotechnology.

Molecular Characterization of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family in Betula luminifera

As a major family of plant-specific transcription factors, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes play vital regulatory roles in plant growth, development and stress responses. In this study, 18 SPL genes were identified and cloned from Betula luminifera. Two zinc finger-like structures and a nuclear location signal (NLS) segments were existed in the SBP domains of all BlSPLs. Phylogenetic analysis showed that these genes were clustered into nine groups (group I-IX). The intron/exon structure and motif composition were highly conserved within the same group. 12 of the 18 BlSPLs were experimentally verified as the targets of miR156, and two cleavage sites were detected in these miR156-targeted BlSPL genes. Many putative cis-elements, associated with light, stresses and phytohormones response, were identified in the promoter regions of BlSPLs, suggesting that BlSPL genes are probably involved in important physiological processes and developmental events. Tissue-specific expression analysis showed that miR156-targeted BlSPLs exhibited a more differential expression pattern, while most miR156-nontargeted BlSPLs tended to be constitutively expressed, suggesting the distinct roles of miR156-targeted and nontargeted BlSPLs in development and growth of B. luminifera. Further expression analysis revealed that miR156-targeted BlSPLs were dramatically up-regulated with age, whereas mature BlmiR156 level was apparently declined with age, indicating that miR156/SPL module plays important roles in vegetative phase change of B. luminifera. Moreover, yeast two-hybrid assay indicated that several miR156-targeted and nontargeted BlSPLs could interact with two DELLA proteins (BlRGA and BlRGL), which suggests that certain BlSPLs take part in the GA regulated processes through protein interaction with DELLA proteins. All these results provide an important basis for further exploring the biological functions of BlSPLs in B. luminifera.

Identification of MicroRNAs and Their Target Genes Related to the Accumulation of Anthocyanins in Litchi chinensis by High-Throughput Sequencing and Degradome Analysis

Litchi (Litchi chinensis Sonn.) is an important subtropical fruit in southern China and the fruit pericarp has attractive red skin at maturity, which is provided by anthocyanins accumulation. To understand the anthocyanin biosynthesis at post-transcriptional level, we investigated the roles of microRNAs (miRNAs) during fruit coloring. In the present study, four small RNA libraries and a mixed degradome library from pericarps of 'Feizixiao' litchi at different developmental phases were constructed and sequenced by Solexa technology. A total of 78 conserved miRNAs belonging to 35 miRNA families and 41 novel miRNAs were identified via high-throughput sequencing, and 129 genes were identified as their targets by the recently developed degradome sequencing. miR156a and a novel microRNA (NEW41) were found to be differentially expressed during fruit coloring, indicating they might affect anthocyanin biosynthesis through their target genes in litchi. qRT-PCR analysis confirmed the expression changes of miR156a and the novel microRNA (NEW41) were inversely correlated with the expression profiles of their target genes LcSPL1/2 and LcCHI, respectively, suggesting regulatory roles of these miRNAs during anthocyanin biosynthesis. The target genes of miR156a, LcSPL1/2, encode transcription factors, as evidenced by a localization in the nucleus, that might play roles in the regulation of transcription. To further explore the relationship of LcSPL1/2 with the anthocyanin regulatory genes, yeast two-hybrid and BiFC analyses showed that LcSPL1 proteins could interact with LcMYB1, which is the key regulatory gene in anthocyanin biosynthesis in litchi. This study represents a comprehensive expression profiling of miRNAs in anthocyanin biosynthesis during litchi fruit maturity and confirmed that the miR156- SPLs module was conserved in anthocyanin biosynthesis in litchi.