Skip to main content


Genome-wide characterization of SPL family in Medicago truncatula reveals the novel roles of miR156/SPL module in spiky pod development

Article metrics

  • 740 Accesses



SQUAMOSA Promoter Binding Protein-Likes (SPLs) proteins are plant-specific transcription factors that play many crucial roles in plant growth and development. However, there is little information about SPL family in the model legume Medicago truncatula.


In this study, a total of 23 MtSPL genes were identified in M. truncatula genome, in which 17 of the MtSPLs contained the putative MtmiR156 binding site at the coding or 3′ UTR regions. Tissue-specific expression pattern analysis showed that most MtmiR156-targeted MtSPLs were highly expressed in seed and pod. The observation of MtmiR156B-overexpressing plants reveals that MtmiR156/MtSPL modules are not only involved in the development of leaves and branches, but also in the seed pod development, especially the formation of spine on pod.


The spines on pods are developed in many plant species, which allow pods to adhere to the animals, and then be transported on the outside. This study sheds light on the new function of SPL family in seed dispersal by controlling the formation of spiky pod, and provides insights on understanding evolutionary divergence of the members of SPL gene family among plant species.


Transcription factors (TFs) play an essential role in regulatory networks of many important developmental processes by activating or repressing the transcription of downstream target genes [1]. The SQUAMOSA Promoter-Binding Protein-Likes (SPLs) proteins are plant-specific TFs and all of them have a highly conserved SBP domain with proximately 78 amino acids in length [2,3,4]. The SBP domain contains two zinc-finger like structure (Cys-Cys-Cys-His and Cys-Cys-His-Cys) and one nuclear localization signal (NLS) motif [5, 6]. SPL proteins specifically recognize and bind to the cis-element TNCGTACAA, which is located at the promoters of target genes with GTAC as its core sequence [5, 7]. The initial two SPL genes, AmSBP1 and AmSBP2, are firstly identified from the Antirrhinum majus, and they act as transcriptional activators and regulate the expression of floral meristem identity gene SQUAMOSA [3]. The SPL genes have been reported to play the important role in regulation of multiple aspects of plant growth and development, including leaf and shoot development [8, 9], vegetative phase change [10, 11], flowering [12], sporogenesis [13], male fertility [14], plant hormone signal transduction [15], and copper homeostasis [16]. So far, the SPL families have been identified and characterized in several plant species, such as Arabidopsis thaliana [2], Oryza sativa [17], Solanum lycopersicum [18], Gossypium raimondii [19], Vitis vinifera [20], Brassica napus [21], Glycine max [22], Prunus mume [23], Arachis hypogaea L. [24], Phyllostachys edulis [25] and Capsicum annuum L. [26].

Many microRNAs have been identified, which regulate gene expression by binding and cleaving their target mRNAs [27, 28]. miRNA156 is one of the highly conserved miRNA families in plants [29]. In Arabidopsis, in total 16 SPLs are identified and termed as AtSPL1 to AtSPL16, respectively [5]. Among them, 10 are targets of AtmiR156, which are AtSPL2-AtSPL6, AtSPL9-AtSPL11, AtSPL13, and AtSPL15 [9, 12, 30,31,32]. In rice, 19 OsSPL genes have been identified, and 11 of them are targeted by OsmiR156 [17]. Most Arabidopsis miR156 binding sites of the targeted-SPL genes are located in the downstream of the SBP domain at the coding sequences, while in AtSPL3, AtSPL4, and AtSPL5, they are located in the 3′ UTR of the mRNAs [12].

These miR156-targeted SPL genes play redundant roles in plant morphology and development among different species. Several studies show that SPL genes are involved in the regulation of flower and inflorescence development. In Arabidopsis, AtSPL3, AtSPL4 and AtSPL5 have similar functions and play vital roles in vegetative phase change and floral transition [7, 9, 12]. Furthermore, AtSPL3 can bind directly on the promoter regions of AP1, LFY and FUL to activate their expression in controlling the timing of flower formation [33]. In addition, SOC1 and FT regulate the expression of AtSPL3, AtSPL4 and AtSPL5 in response to photoperiod and gibberellin (GA) signals to promote flowering [34]. In wheat, two paralogous SPL genes, TaSPL20 and TaSPL21, are highly expressed in the lemma and palea, and ectopic expression of them in rice can promote panicle branching [35]. SPL genes also play important roles in regulating lateral organ and shoot development. In Arabidopsis, AtSPL2, AtSPL10 and AtSPL11 redundantly regulate proper development of lateral organs in association with shoot maturation in the reproductive phase [36]. Besides, mutant phenotype analysis shows that AtSPL9 and AtSPL15 act redundantly in regulation of the juvenile-to-adult phase transition [11]. In maize, SPL transcription factor TASSELSHEATH4 plays an essential role in bract development and the establishment of meristem boundaries [37]. Two duplicated loci, UNBRANCHED2 and UNBRANCHED3, affect crop productivity traits by regulating the rate of lateral primordia initiation [38]. In addition, a series of studies show that SPL genes are involved in the regulation of seed and fruit development. In rice, the OsmiR156-targeted SPL gene, OsSPL14/IDEAL PLANT ARCHITECTURE 1 (IPA1), regulates shoot branching in the vegetative stage and enhances rice grain yield at the productive stage [39, 40]. OsSPL16/GRAIN WIDTH 8 (GW8) acts as a positive regulator of cell proliferation and is involved in control of grain size, shape and quality [41]. In tomato, the SPL gene COLORLESS NON-RIPENING (CNR) is critical for normal fruit ripening, and mutation in CNR results in colorless fruits with a substantial loss of cell-to-cell adhesion [42].

Medicago truncatula is a fast-emerging model legume for functional genomics study. However, the information and function of SPL gene family in M. truncatula are largely unknown. In this study, we reported the genome-wide identification and characterization of SPL genes in M. truncatula. Totally, 23 MtSPL genes were identified, and their phylogenetic relationship, protein motifs, gene structures, and chromosomal locations were analyzed. Furthermore, we found that most MtmiR156-targeted MtSPL genes were highly expressed in pod and seed. Overexpression of MtmiR156B transgenic plants displayed the small seeds and pods, especially the loss of pod spine. These findings demonstrate that MtmiR156-targeted MtSPL genes play the novel roles in pod and seed development in M. truncatula.


Genome-wide identification of MtSPL genes

In order to identify SPL genes in M. truncatula genome, we executed a genome-wide search of MtSPLs by protein BLAST using the 16 AtSPLs and 19 OsSPLs sequences against the Medicago truncatula Genome Database. Initially, a total of 24 putative MtSPL sequences were obtained. Medtr8g066240 was excluded from the MtSPL gene family due to the absence of a complete SBP domain. Therefore, totally 23 MtSPLs with the conserved SBP domain were characterized in genome of M. truncatula (Additional file 1). The MtSPL genes were named according to their closest Arabidopsis orthologs (Fig. 1). The protein lengths of MtSPLs ranged from 143 to 1025 amino acids, and the gene locus, isoelectric point, intron number, and chromosome location of the MtSPL genes were shown in Table 1.

Fig. 1

Phylogenetic analysis of SPL proteins from M. truncatula, Arabidopsis and rice. The Neighbor-Joining (NJ) phylogenetic tree was constructed using full length SPL protein sequences in MEGA 7.0 with 1000 replicates

Table 1 The SPL gene family in M. truncatula

Phylogenetic analysis and chromosomal localization of MtSPL genes

To further achieve the evolutionary relationship between MtSPL genes and other SPLs, a phylogenetic tree containing 16 AtSPLs, 19 OsSPLs, and 23 MtSPLs was constructed using MEGA7.0 with Neighbor-Joining method (Fig. 1). According to the phylogenetic analyses, the 58 SPL proteins were clustered into seven distinct groups (named G1-G7), and each group contained at least one SPL from A. thaliana, rice, and M. truncatula (Fig. 1). To determine the chromosome distribution of MtSPL genes in M. truncatula, the 23 MtSPL genes were located on each chromosome based on the M. truncatula genome data. These MtSPL genes showed uneven distribution on the M. truncatula chromosomes (Fig. 2a). Chromosome 8, 2, 7, and 1 contained 6, 5, 4, and 3 MtSPL genes, respectively. Both chromosome 3 and 4 contained 2 MtSPL genes, chromosome 5 contained only one MtSPL gene, and no MtSPL gene located on chromosome 6.

Fig. 2

Chromosomal distribution, phylogenetic relationship and gene structure of MtSPL genes. a The distribution of MtSPL genes on M. truncatula chromosomes based on the genome annotation. The chromosomes are indicated with different colors. The scale (Mb) represents the lengths of the chromosomes. b Phylogenetic analysis and exon/intron structures of MtSPL genes. Exons and introns are indicated by black wedges and black lines. The scale (Kb) represents the lengths of the MtSPL genes

Conserved motifs and gene structure analysis of MtSPL genes

To further understand the structural diversity of the MtSPL genes, gene exon/intron structure analysis was carried out using online Gene Structure Display Server tool. The exon/intron structures of the 23 MtSPL genes were generated by alignment of MtSPL gene coding sequences with their corresponding genomic sequences. Gene structure illustrations displayed the high variation in the number of introns, from 0 to 10 (Fig. 2b). To gain a better understanding of the protein sequence characteristics of the MtSPLs, online MEME search was performed to analyze the conserved domains and motifs. Besides the conserved SBP domain (motif 4), in total 20 conserved motifs were identified in MtSPLs (Fig. 3a, Additional file 2). The conserved SBP domain consisted of three motifs: zinc finger 1 (Zn-1, C3H), zinc finger 2 (Zn-2, C2HC), and a nuclear localization signal (Fig. 3b, c).

Fig. 3

Conserved domains and motifs in MtSPL proteins. a The full length MtSPL protein sequences were used to execute the motif search on MEME tool. Domains and motifs were represented by the boxes with different numbers and colors. b Alignment of the conserved SBP domain in MtSPL proteins. Multiple SBP domain sequences alignment was performed using Jalview software. Two Zn-finger structures (Zn-1, Cys3His and Zn-2, Cys2HisCys) and one NLS are indicated. c Sequence logo of the conserved SBP domain of MtSPLs. Sequence logo was obtained from Weblogo online software. The overall height of each stack indicates the sequence conservation at each position, while the height of the letters within each stack indicates the relative frequency of the corresponding amino acid

Analysis of MtmiR156 and its target sequences in MtSPLs

In order to understand the function of miR156 and miR156-targeted MtSPL genes, we searched the miRBase Database and found 10 MtmiR156 genes, MtmiR156A-MtmiR156J, in M. truncatula genome (Fig. 4a, Additional file 3a). Based on the MtmiR156 precursor sequences, the stem-loop structures of MtmiR156 were found (Additional file 3a). Previous studies showed that miR156 could bind to the coding or 3′ UTR region of the SPL genes and reduce gene activity. Then, we used the online psRNATarget tool to find MtmiR156 complementary sequences in the MtSPL transcripts. By comparison of the MtmiR156 mature sequences and the MtSPL transcript sequences, we found that total 17 MtSPL genes have the MtmiR156 binding sites, 13 of which were located in coding regions and 4 in 3′ UTR regions, respectively (Fig. 4b, Additional file 3b).

Fig. 4

miR156 sequences and binding sites of MtmiR156 in MtSPL genes. a Mature miR156 sequences in A. thaliana and M. truncatula. Shaded zones indicate the conserved sequences. b The MtSPLs regulated by MtmiR156. The gray boxes represent the CDS of MtSPL genes. The gray lines represent 3′ UTR of the MtSPL genes. The light blue boxes represent the conserved SBP domain. The red lines represent the MtmiR156 target site

Expression profile of MtSPL genes in different organs

The expression pattern of a gene is often correlated with its function. In order to understand the developmental functions of the MtSPL genes, we investigated the expression profiles of the 23 MtSPL genes by quantitative real-time PCR (qRT-PCR) in six different organs, including roots, stems, leaves, flowers, pods, and seeds. qRT-PCR results showed that the relative expression levels and patterns of the MtSPL genes were varied in these organs (Fig. 5). For example, the non-MtmiR156-targeted MtSPLs (MtSPL1, MtSPL7, MtSPL12, MtSPL14 and MtSPL16) were expressed in all of the organs tested, while MtSPL8 was highly expressed in flower and pod (Fig. 5a). The MtmiR156-targeted MtSPL genes also showed differential expression profiles (Fig. 5b). Most MtmiR156-targeted MtSPLs, such as MtSPL4, MtSPL5A, MtSPL5B, MtSPL6B, MtSPL6C, MtSPL11, MtSPL13B and MtSPL15A, were highly detected in pod, implying their specific roles in pod development. MtSPL2, MtSPL5B, MtSPL10A and MtSPL13A were expressed at high levels in seed, indicating the possible involvement in seed development.

Fig. 5

Expression patterns of MtSPL genes in six different organs. a Expression patterns of non-MtmiR156-targeted MtSPL genes. b Expression patterns of MtmiR156-targeted MtSPL genes

MtmiR156-targeted MtSPLs play important roles in pod and seed development

To investigate the possible roles of the MtmiR156-targeted MtSPLs in growth and development of M. truncatula, the genomic sequence of MtmiR156B was cloned and introduced into wild type plants under the regulation of the cauliflower mosaic virus 35S promoter (Fig. 6a). Seven positive transgenic lines were obtained based on the PCR results (Additional file 4). The MtmiR156B was highly expressed in two transgenic plants (Fig. 6b). The MtmiR156B-overexpressing plants exhibited small leaves, increased lateral branches and reduced plant height, indicating that MtmiR156B-targeted MtSPLs play important roles in leaf morphogenesis, branching and stem elongation (Fig. 6c-e). Moreover, reduced function of the MtmiR156B-targeted MtSPLs also led to the defects in reproductive development. Compared with wild type, the spikes of the flowers in MtmiR156B-overexpressing plants were absent (Fig. 6f), however, the development of stamen and carpel was normal (Fig. 6g).

Fig. 6

Developmental phenotypes of the MtmiR156B-overexpressing plants. a Schematic illustration of vector used for MtmiR156B overexpression. b Transcript levels of MtmiR156B in wild type and 35S:MtmiR156B transgenic plants. c Two-month-old wild type and the 35S:MtmiR156B transgenic plants. d The number of primary branch in the wild type and the 35S:MtmiR156B transgenic plants. P < 0.001. e Adult leaves of the wild type and 35S:MtmiR156B transgenic plants. f Flower phenotype in the wild type and 35S:MtmiR156B transgenic plants. Arrowheads indicated the developmental defects of spike in MtmiR156B-overexpressing plants. g The side view of the central carpel of wild type and 35S:MtmiR156B transgenic plants. Bars: 5 cm (c), 1 cm (e), 2 mm (f), 2 mm (g)

Most MtmiR156-targeted MtSPLs were highly expressed in pod, implying that they may play redundant roles in pod development. Compared with the long slender spines developed on the pod surface in wild type, overexpression of MtmiR156B led to the reduction in pod size and the loss of spines on pod surface (Fig. 7a-d). The seed number in each pod, seed size and weight were significantly reduced in the MtmiR156B-overexpressing plants, compared with those in wild type (Fig. 7e-g). These results demonstrate that the MtmiR156/MtSPL regulation module is critical for the pod and seed development.

Fig. 7

MtmiR156-targeted MtSPL genes are involved in pod and seed development. a and b The vertical (a) and side (b) view of pods of wild type and 35S:MtmiR156B transgenic plants. c and d Close view of the pods in (a) and (b). e Seeds of wild type and 35S:MtmiR156B transgenic plants. f Seeds number in each pod in wild type and 35S:MtmiR156B transgenic plants. g Weight of 100 seeds of wild type and 35S:MtmiR156B transgenic plants. P < 0.01, P < 0.001. Bars: 1 cm (a, b), 2 mm (c, d), 5 mm (e)

To further determine which MtSPL genes are involved in pod wall and seed development, the expression levels of all the MtmiR156-targeted SPL genes were analyzed in the pod wall and seed of the wild type and MtmiR156B-overexpressing plants. qRT-PCR results demonstrated that the expression of nine MtSPLs was significantly reduced in the pod wall of the MtmiR156B-overexpressing plants (Fig. 8a). The expression of MtSPL5A, MtSPL5B, MtSPL15A, and MtSPL15B showed over two-fold decrease in the transgenic plants, indicating that those MtSPL genes play crucial roles in spiky pod wall development (Fig. 8a). While, eleven MtSPLs, especially MtSPL5A, MtSPL6B, MtSPL6C, MtSPL10A and MtSPL13B, were downregulated in seed of the MtmiR156B-overexpressing plants, indicating the functional redundancy among those MtSPL genes during seed development (Fig. 8b).

Fig. 8

Expression levels of MtmiR156-targeted SPL genes in pod wall and seed of the wild type and 35S:MtmiR156B-overexpressing plants. a Expression levels of the MtmiR156-targeted SPLs in pod wall. b Expression levels of the MtmiR156-targeted SPLs in seed. P < 0.05, P < 0.01, P < 0.001


Transcription factors play important roles during the processes of plant growth and development. The SPL genes encode a family of plant-specific transcription factors that contain the conserved SBP domain [2, 43]. In this study, through a genome-wide identification, we obtained 23 MtSPL genes from M. truncatula genome. Phylogenetic analysis showed that the MtSPL genes were more closely related to Arabidopsis than rice SPL genes, indicating that eudicots SPL genes may diverge from a common ancestor [44]. However, the number of MtSPL genes in M. truncatula is greater than that in Arabidopsis and rice, which contain 16 and 19 SPLs, respectively. Sequence homologous analysis suggested that faster gene duplication rates or species-specific gene duplication manners might play important roles in SPL gene expansion in M. truncatula. Gene structure and motif composition analyses showed that MtSPL genes within the same group shared similar motifs and exon/intron organization, suggesting that the functional evolution of SPL genes may be tightly correlated with the diversification of gene structure and conservation of motifs [6, 45].

Some SPL genes were posttranscriptionally regulated by miR156 and involved in multiple developmental processes [4, 9, 11, 32, 36, 39,40,41, 45,46,47,48,49]. Based on the miRNA database information, ten MtmiR156 genes were found in M. truncatula genome. The mature sequences of the miR156 genes between A. thaliana and M. truncatula are similar, indicating the functional conservation of the miR156 in different plant species [29]. It has been reported that 11 of 16 SPLs in Arabidopsis and 11 of 19 SPLs in rice contained the putative miR156 binding sites [2, 7, 11, 12, 17, 36, 46]. By searching the MtmiR156 complementary sequence in MtSPLs mRNAs, we found that 17 MtSPLs out of 23 contained the putative MtmiR156 binding sites, suggesting the conservation of miR156-mediated posttranscriptional regulation in plants.

In Arabidopsis, loss-of-function of multiple SPL genes or overexpression of AtmiR156A resulted in the decreased rosette leaf area [46]. In tomato and petunia, overexpressing AtMIR156B or PhSBP1code-silenced plants produced the higher number of small leaves and branches [47, 50]. Furthermore, overexpression of OsmiR156B in switchgrass and loss-of-function in IPA1 in rice and down-regulation of MsSPL13 in alfalfa also led to the increased number of branches [39, 40, 51, 52]. In accordance with these reports, the similar phenotypes, including more branches and small leaves, were also displayed in MtmiR156B-overexpressing plants. These observations indicate that MtmiR156/MtSPL regulation module plays the conserved roles in vegetative development.

The seeds of legume are developed within an ovary-derived pod, whose walls provide numerous functions, such as the protection of the seeds and the production of photosynthates [53,54,55]. Seed dispersal is the transport of pods/seeds away from the parent plant and has been implicated in rapid plant migration and the spread of invasive [56,57,58,59]. The pods in many plant species develop the spines or stiff hairs, such as Trifolium angustifolium, M. polymorpha and M. truncatula [60, 61]. Such structures are very important for seed migration, because seed pods can adhere to animals by means of spines or hairs, and be transported on the outside [60, 62]. So, the proper development of spines/hairs on pod is critical for seed dispersal, along with species diversity or ecological invasion. The pod and pod wall of M. truncatula are helical and thick with spines [63,64,65]. In this study, we found that most MtmiR156-targeted MtSPL genes were highly expressed in pod and seed. Moreover, seed size and number in MtmiR156B-overexpressing plants were decreased. Importantly, the development of spines on pod was also defective, due to the downregulation of several MtmiR156-targeted MtSPLs. These observations indicate that the MtmiR156/MtSPL regulation module may contribute to the genetic variability through the regulation of pod morphology.

In Arabidopsis, ectopic expression of the TaSPL16 results in early flowering and increase of organ size and yield [66]. This finding implies that SPL is possible for the improvement of seed production in legume species. In this study, 35S:MtmiR156B plants show defective spines on pod in M. trunctula. However, the number of lateral branches is increased dramatically. The biomass is a critical trait in evaluation of forage grass. Therefore, downregulation of targeted SPL genes by overexpression of miR156 in legume forage, such as alfalfa, may provide a helpful tool to improve forage production.


In this study, we performed genome-wide analyses and identified SPL genes in M. truncatula. The genetic redundancy of MtSPL genes hinders the discoveries of their potential functions. However, the phenotypes of MtmiR156B-overexpressing plants reveal that MtmiR156/MtSPL modules are not only involved in the development of leaves and branches, but also indirectly contribute to seed dispersal by controlling the formation of spine on pods. Characterization of the loss-of-function MtSPL mutants may help to provide insight into the roles of MtmiR156/MtSPL module in the development of spine of pod, and shed light on the new function of SPL family among plant species.


Plant materials and growth conditions

Medicago truncatula R108 ecotype was used as the wild type, which is obtained from the Noble Research Institution, USA. The seeds of 35S:MtmiR156B and wild type were scarified with sandpaper and treated at 4 °C for 5 days. The geminated seeds were transferred to nursery seedling plate (4 × 4 × 5 cm Length, Width, Height) for 2 weeks. Then, the seedlings were transferred to Luqing soil mix (soil: vermiculite = 3:1) and grown in the greenhouse at 22 °C (day) / 22 °C (night) with 16 h (day) / 8 h (night) photoperiod, and relative humidity at 70 to 80%.

Identification and phylogenetic analysis of SPL genes in M. truncatula

To identify the SPL genes in M. truncatula genome, firstly, we used 16 AtSPL and 19 OsSPL protein sequences to execute BLASTP search the Medicago truncatula Genome Database ( The AtSPL protein sequences were obtained from The Arabidopsis Information Resource (TAIR, The OsSPL protein sequences were obtained from the Rice Genome Annotation Project ( Totally, 24 putative MtSPL genes were identified in M. truncatula genome using blast with a cut-off E-value >1e− 3. Secondly, we searched the Plant Transcription Factor Database ( and confirmed the blast search result. Thirdly, the 24 MtSPL protein sequences were further analyzed on the NCBI Conserved Domain Search website ( and found that Medtr8g066240 lost the conserved SBP domain. Therefore, Medtr8g066240 was excluded from the putative MtSPL genes, and total 23 MtSPL genes were used for study.

Multiple protein sequence alignment was performed using Jalview software ( The phylogenetic tree for Arabidopsis, rice, and Medicago SPL gene family members was constructed using MEGA 7.0.

Chromosome location and gene structure of MtSPL genes

The informations of MtSPL genes on chromosome location, including chromosome length, gene direction, and gene start and stop position, were obtained from the M. truncatula genome database. Exon / intron structures of MtSPL genes were determined by aligning the coding sequences and their corresponding genomic sequences using the online Gene Structure Display Server (GSDS, website.

The identification of conserved domain and the prediction of MtmiR156-targeted MtSPLs prediction

The online MEME tool ( was used to predict both conserved domains and potential motifs in the 23 MtSPL proteins with the following parameters: maximum number of motifs to find, 20; minimum width of motif, 6; maximum width of motif, 100; minimum number of sites for each motif, 2. The mature sequences of M. truncatula MtmiR156A to MtmiR156J were obtained from miRBase database ( The MtmiR156-targeted MtSPL genes and their binding sites were obtained by searching the gene coding and UTR regions on the psRNATarget ( website.

RNA extraction and gene expression analysis

The samples in 60-day old wild-type plants were used for RNA extraction. For gene expression pattern analysis, the roots, leaves, flowers, pods and seeds samples were harvested from primary roots, fully expanded leaves, fully opened flowers, 20-day post-pollination pods and seeds. To analyze the relative expression levels of MtmiR156-targeted MtSPL genes in the MtmiR156B-overexpressing plants, 20-day post-pollination seeds and pod walls were collected from wild type and 35S:MtmiR156B transgenic plants. To analyze the relative expression levels of MtmiR156B in the MtmiR156B-overexpressing plants, 60-day old fully expanded leaves were collected from wild type and 35S:MtmiR156B transgenic plants.

Total RNA of these organs was extracted using the Trizol-RT Reagent (Molecular Research Center, INC) following the manufacturer’s instructions. The quality and quantity of the extracted RNA were analyzed using Nanodrop 2000 Spectrophotometer (Thermo Scientific, USA). Reverse transcription PCR was performed with 2.5 μg total RNA using Roche First Strand cDNA Synthesis Kit (Roche, USA). Then, the cDNA was diluted to 20 ng/μl with DEPC treated H2O. For quantitative real-time PCR (qRT-PCR) analysis, the primers of the 23 MtSPL genes were designed by Primer Express 3.0 software (Additional file 5). qRT-PCR was executed in triplicate for each organ on CFX Connect™ (Bio-Rad, USA) using FastStart Essential DNA Green Master Kit (Roche, USA). The MtUBI gene was selected as internal control, and the relative expression levels of different MtSPL genes were calculated using 2-CT method [67].

Plasmid construction and plant transformation

To obtain the MtmiR156B overexpression construction, the 830 bp DNA sequence containing the MtmiR156B stem-loop structure was amplified using primer pair MtmiR156B-F/MtmiR156B-R. The MtmiR156B sequence was transferred into the pEarleyGate 100 vector by Gateway LR reaction (Invitrogen, USA) [68]. Then 35S:MtmiR156B destination construct was introduced into Agrobacterium strain EHA105. For stable transformation, leaves of wild type were transformed with EHA105 strain containing the 35S:MtmiR156B vector [69].

Availability of data and materials

Accession numbers: MtSPL1: Medtr1g086250; MtSPL2: Medtr3g085180; MtSPL3: Medtr4g088555; MtSPL4:Medtr2g014200; MtSPL5A:Medtr2g078770; MtSPL5B:Medtr8g463140; MtSPL6A:Medtr5g046670; MtSPL6B:Medtr2g461920; MtSPL6C:Medtr4g109770; MtSPL7:Medtr2g020620; MtSPL8:Medtr8g005960; MtSPL9:Medtr7g444860; MtSPL10A:Medtr8g080680; MtSPL10B:Medtr8g080670; MtSPL11:Medtr8g080690; MtSPL12:Medtr7g110320; MtSPL13A:Medtr8g096780; MtSPL13B:Medtr3g099080; MtSPL13C:Medtr7g028740; MtSPL14:Medtr1g035010; MtSPL15A:Medtr7g092930; MtSPL15B:Medtr1g053715; MtSPL16:Medtr2g046550; MtmiR156B: MIMAT0011057. The sequences of all genes involved in this study are listed in Additional file 1, and also can be found in NCBI database. The plant materials in this study are available from the corresponding author on request.


AP1 :


FT :











Nuclear localization signal


Quantitative real-time PCR



SOC1 :





Transcription factors


Untranslated region


Zinc finger 1


Zinc finger 2


  1. 1.

    Yang ZF, Wang XF, Gu SL, Hu ZQ, Xu H, Xu CW. Comparative study of SBP-box gene family in Arabidopsis and rice. Gene. 2008;407(1–2):1–11.

  2. 2.

    Cardon G, Hohmann S, Klein J, Nettesheim K, Saedler H, Huijser P. Molecular characterisation of the Arabidopsis SBP-box genes. Gene. 1999;237(1):91–104.

  3. 3.

    Klein J, Saedler H, Huijser P. A new family of DNA binding proteins includes putative transcriptional regulators of the Antirrhinum majus floral meristem identity gene SQUAMOSA. Mol Gen Genet. 1996;250(1):7–16.

  4. 4.

    Chen X, Zhang Z, Liu D, Zhang K, Li A, Mao L. SQUAMOSA promoter-binding protein-like transcription factors: star players for plant growth and development. J Integr Plant Biol. 2010;52(11):946–51.

  5. 5.

    Birkenbihl RP, Jach G, Saedler H, Huijser P. Functional dissection of the plant-specific SBP-domain: overlap of the DNA-binding and nuclear localization domains. J Mol Biol. 2005;352(3):585–96.

  6. 6.

    Yamasaki K, Kigawa T, Inoue M, Tateno M, Yamasaki T, Yabuki T, Aoki M, Seki E, Matsuda T, Nunokawa E, et al. A novel zinc-binding motif revealed by solution structures of DNA-binding domains of Arabidopsis SBP-family transcription factors. J Mol Biol. 2004;337(1):49–63.

  7. 7.

    Cardon GH, Hohmann S, Nettesheim K, Saedler H, Huijser P. Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3: a novel gene involved in the floral transition. Plant J. 1997;12(2):367–77.

  8. 8.

    Usami T, Horiguchi G, Yano S, Tsukaya H. The more and smaller cells mutants of Arabidopsis thaliana identify novel roles for SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes in the control of heteroblasty. Development. 2009;136(6):955–64.

  9. 9.

    Wu G, Poethig RS. Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3. Development. 2006;133(18):3539–47.

  10. 10.

    Jung JH, Seo PJ, Kang SK, Park CM. miR172 signals are incorporated into the miR156 signaling pathway at the SPL3/4/5 genes in Arabidopsis developmental transitions. Plant Mol Biol. 2011;76(1–2):35–45.

  11. 11.

    Schwarz S, Grande AV, Bujdoso N, Saedler H, Huijser P. The microRNA regulated SBP-box genes SPL9 and SPL15 control shoot maturation in Arabidopsis. Plant Mol Biol. 2008;67(1–2):183–95.

  12. 12.

    Gandikota M, Birkenbihl RP, Hohmann S, Cardon GH, Saedler H, Huijser P. The miRNA156/157 recognition element in the 3′ UTR of the Arabidopsis SBP box gene SPL3 prevents early flowering by translational inhibition in seedlings. Plant J. 2007;49(4):683–93.

  13. 13.

    Unte US, Sorensen AM, Pesaresi P, Gandikota M, Leister D, Saedler H, Huijser P. SPL8, an SBP-box gene that affects pollen sac development in Arabidopsis. Plant Cell. 2003;15(4):1009–19.

  14. 14.

    Xing SP, Salinas M, Hohmann S, Berndtgen R, Huijser P. miR156-targeted and nontargeted SBP-box transcription factors act in concert to secure male fertility in Arabidopsis. Plant Cell. 2010;22(12):3935–50.

  15. 15.

    Zhang Y, Schwarz S, Saedler H, Huijser P. SPL8, a local regulator in a subset of gibberellin-mediated developmental processes in Arabidopsis. Plant Mol Biol. 2007;63(3):429–39.

  16. 16.

    Yamasaki H, Hayashi M, Fukazawa M, Kobayashi Y, Shikanai T. SQUAMOSA promoter binding protein-Like7 is a central regulator for copper homeostasis in Arabidopsis. Plant Cell. 2009;21(1):347–61.

  17. 17.

    Xie KB, Wu CQ, Xiong LZ. Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice. Plant Physiol. 2006;142(1):280–93.

  18. 18.

    Salinas M, Xing SP, Hohmann S, Berndtgen R, Huijser P. Genomic organization, phylogenetic comparison and differential expression of the SBP-box family of transcription factors in tomato. Planta. 2012;235(6):1171–84.

  19. 19.

    Cai C, Guo W, Zhang B. Genome-wide identification and characterization of SPL transcription factor family and their evolution and expression profiling analysis in cotton. Sci Rep. 2018;8(1):762.

  20. 20.

    Hou H, Li J, Gao M, Singer SD, Wang H, Mao L, Fei Z, Wang X. Genomic organization, phylogenetic comparison and differential expression of the SBP-box family genes in grape. PLoS One. 2013;8(3):e59358.

  21. 21.

    Cheng H, Hao M, Wang W, Mei D, Tong C, Wang H, Liu J, Fu L, Hu Q. Genomic identification, characterization and differential expression analysis of SBP-box gene family in Brassica napus. BMC Plant Biol. 2016;16(1):196.

  22. 22.

    Tripathi RK, Goel R, Kumari S, Dahuja A. Genomic organization, phylogenetic comparison, and expression profiles of the SPL family genes and their regulation in soybean. Dev Genes Evol. 2017;227(2):101–19.

  23. 23.

    Xu Z, Sun L, Zhou Y, Yang W, Cheng T, Wang J, Zhang Q. Identification and expression analysis of the SQUAMOSA promoter-binding protein (SBP)-box gene family in Prunus mume. Mol Gen Genomics. 2015;290(5):1701–15.

  24. 24.

    Li M, Zhao SZ, Zhao CZ, Zhang Y, Xia H, Lopez-Baltazar J, Wan SB, Wang XJ. Cloning and characterization of SPL-family genes in the peanut (Arachis hypogaea L.). Genet Mol Res. 2016;15(1):gmr.15017344.

  25. 25.

    Pan F, Wang Y, Liu H, Wu M, Chu W, Chen D, Xiang Y. Genome-wide identification and expression analysis of SBP-like transcription factor genes in Moso bamboo (Phyllostachys edulis). BMC Genomics. 2017;18(1):486.

  26. 26.

    Zhang HX, Jin JH, He YM, Lu BY, Li DW, Chai WG, Khan A, Gong ZH. Genome-wide identification and analysis of the SBP-box family genes under Phytophthora capsici stress in pepper (Capsicum annuum L.). Front Plant Sci. 2016;7:504.

  27. 27.

    Voinnet O. Origin, biogenesis, and activity of plant MicroRNAs. Cell. 2009;136(4):669–87.

  28. 28.

    Rogers K, Chen XM. Biogenesis, turnover, and mode of action of plant MicroRNAs. Plant Cell. 2013;25(7):2383–99.

  29. 29.

    Axtell MJ, Bowman JL. Evolution of plant microRNAs and their targets. Trends Plant Sci. 2008;13(7):343–9.

  30. 30.

    Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ. Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008;18(10):758–62.

  31. 31.

    Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D. Specific effects of MicroRNAs on the plant transcriptome. Dev Cell. 2005;8(4):517–27.

  32. 32.

    Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP. Prediction of plant microRNA targets. Cell. 2002;110(4):513–20.

  33. 33.

    Yamaguchi A, Wu MF, Yang L, Wu G, Poethig RS, Wagner D. The microRNA-regulated SBP-box transcription factor SPL3 is a direct upstream activator of LEAFY, FRUITFULL, and APETALA1. Dev Cell. 2009;17(2):268–78.

  34. 34.

    Jung JH, Ju Y, Seo PJ, Lee JH, Park CM. The SOC1-SPL module integrates photoperiod and gibberellic acid signals to control flowering time in Arabidopsis. Plant J. 2012;69(4):577–88.

  35. 35.

    Zhang B, Xu W, Liu X, Mao X, Li A, Wang J, Chang X, Zhang X, Jing R. Functional conservation and divergence among Homoeologs of TaSPL20 and TaSPL21, two SBP-box genes governing yield-related traits in Hexaploid wheat. Plant Physiol. 2017;174(2):1177–91.

  36. 36.

    Shikata M, Koyama T, Mitsuda N, Ohme-Takagi M. Arabidopsis SBP-box genes SPL10, SPL11 and SPL2 control morphological change in association with shoot maturation in the reproductive phase. Plant Cell Physiol. 2009;50(12):2133–45.

  37. 37.

    Chuck G, Whipple C, Jackson D, Hake S. The maize SBP-box transcription factor encoded by tasselsheath4 regulates bract development and the establishment of meristem boundaries. Development. 2010;137(8):1243–50.

  38. 38.

    Chuck GS, Brown PJ, Meeley R, Hake S. Maize SBP-box transcription factors unbranched2 and unbranched3 affect yield traits by regulating the rate of lateral primordia initiation. Proc Natl Acad Sci U S A. 2014;111(52):18775–80.

  39. 39.

    Jiao Y, Wang Y, Xue D, Wang J, Yan M, Liu G, Dong G, Zeng D, Lu Z, Zhu X, et al. Regulation of OsSPL14 by OsmiR156 defines ideal plant architecture in rice. Nat Genet. 2010;42(6):541–4.

  40. 40.

    Miura K, Ikeda M, Matsubara A, Song XJ, Ito M, Asano K, Matsuoka M, Kitano H, Ashikari M. OsSPL14 promotes panicle branching and higher grain productivity in rice. Nat Genet. 2010;42(6):545–9.

  41. 41.

    Wang S, Wu K, Yuan Q, Liu X, Liu Z, Lin X, Zeng R, Zhu H, Dong G, Qian Q, et al. Control of grain size, shape and quality by OsSPL16 in rice. Nat Genet. 2012;44(8):950–4.

  42. 42.

    Manning K, Tor M, Poole M, Hong Y, Thompson AJ, King GJ, Giovannoni JJ, Seymour GB. A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet. 2006;38(8):948–52.

  43. 43.

    Guo AY, Zhu QH, Gu XC, Ge S, Yang J, Luo JC. Genome-wide identification and evolutionary analysis of the plant specific SBP-box transcription factor family. Gene. 2008;418(1–2):1–8.

  44. 44.

    Preston JC, Hileman LC. Functional evolution in the plant SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family. Front Plant Sci. 2013;4:80.

  45. 45.

    Zhang SD, Ling LZ, Yi TS. Evolution and divergence of SBP-box genes in land plants. BMC Genomics. 2015;16:787.

  46. 46.

    Xu M, Hu T, Zhao J, Park MY, Earley KW, Wu G, Yang L, Poethig RS. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana. PLoS Genet. 2016;12(8):e1006263.

  47. 47.

    Ferreira e Silva GF, Silva EM, Azevedo Mda S, Guivin MA, Ramiro DA, Figueiredo CR, Carrer H, Peres LE, Nogueira FT. microRNA156-targeted SPL/SBP box transcription factors regulate tomato ovary and fruit development. Plant J. 2014;78(4):604–18.

  48. 48.

    Gou J, Fu C, Liu S, Tang C, Debnath S, Flanagan A, Ge Y, Tang Y, Jiang Q, Larson PR, et al. The miR156-SPL4 module predominantly regulates aerial axillary bud formation and controls shoot architecture. New Phytol. 2017;216(3):829–40.

  49. 49.

    Wang H, Wang H. The miR156/SPL module, a regulatory hub and versatile toolbox, gears up crops for enhanced agronomic traits. Mol Plant. 2015;8(5):677–88.

  50. 50.

    Preston JC, Jorgensen SA, Orozco R, Hileman LC. Paralogous SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes differentially regulate leaf initiation and reproductive phase change in petunia. Planta. 2016;243(2):429–40.

  51. 51.

    Gao R, Gruber MY, Amyot L, Hannoufa A. SPL13 regulates shoot branching and flowering time in Medicago sativa. Plant Mol Biol. 2018;96(1–2):119–33.

  52. 52.

    Fu C, Sunkar R, Zhou C, Shen H, Zhang JY, Matts J, Wolf J, Mann DG, Stewart CN Jr, Tang Y, et al. Overexpression of miR156 in switchgrass (Panicum virgatum L.) results in various morphological alterations and leads to improved biomass production. Plant Biotechnol J. 2012;10(4):443–52.

  53. 53.

    Harvey DM. The translocation of 14C-photosynthate in Pisum sativum L. Ann Bot. 1973;37(4):787–94.

  54. 54.

    Setia RC, Setia N, Malik CP. The podwall structure and function in relation to seed development in some legumes. Phyton Annales Rei Botanicae. 1987;27(2):205–20.

  55. 55.

    Atkins CA, Kuo J, Pate JS, Flinn AM, Steele TW. Photosynthetic pod wall of pea (Pisum sativum L.): distribution of carbon dioxide-fixing enzymes in relation to pod structure. Plant Physiol. 1977;60(5):779–86.

  56. 56.

    Vittoz P, Engler R. Seed dispersal distances: a typology based on dispersal modes and plant traits. Bot Helv. 2007;117(2):109–24.

  57. 57.

    Howe HF, Miriti MN. When seed dispersal matters. Bioscience. 2004;54(7):651–60.

  58. 58.

    Sarmento R, Alves-Costa CP, Ayub A, Mello MAR. Partitioning of seed dispersal services between birds and bats in a fragment of the Brazilian Atlantic Forest. Zoologia-Curitiba. 2014;31(3):245–55.

  59. 59.

    Sorensen AE. Seed dispersal by adhesion. Annu Rev Ecol Syst. 1986;17(1):443–63.

  60. 60.

    Manzano P, Malo JE. Extreme long-distance seed dispersal via sheep. Front Ecol Environ. 2006;4(5):244–8.

  61. 61.

    Vitale M, Pupilli F, Labombarda P, Arcioni S. RAPD analysis reveals a low rate of outcrossing in burr medic (Medicago polymorpha L.). Genet Resour Crop Evol. 1998;45(4):337–42.

  62. 62.

    Bonnin I, Huguet T, Gherardi M, Prosperi J-M, Olivieri I. High level of polymorphism and spatial structure in a selfing plant species, Medicago truncatula (Leguminosae), shown using RAPD markers. Am J Bot. 1996;83(7):843–55.

  63. 63.

    Small E, Brookes BS. Coiling of alfalfa pods in relation to resistance against seed chalcids: additional observations. Can J Plant Sci. 1984;64(3):659–65.

  64. 64.

    Small E, Jomphe M, Brookes B. Medicago-Truncatula F Laxicycla (Leguminosae), a new taxon with loose fruit Coiling promoting rapid germination of the fruit-retained seeds. Plant Syst Evol. 1991;174(1–2):37–46.

  65. 65.

    Lesins KA, Lesins I. Genus Medicago (Leguminosae). A taxogenetic study; 1979.

  66. 66.

    Cao R, Guo L, Ma M, Zhang W, Liu X, Zhao H. Identification and functional characterization of Squamosa promoter binding protein-like gene TaSPL16 in wheat (Triticum aestivum L.). Front Plant Sci. 2019;10:212.

  67. 67.

    Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C-T method. Nat Protoc. 2008;3(6):1101–8.

  68. 68.

    Earley KW, Haag JR, Pontes O, Opper K, Juehne T, Song KM, Pikaard CS. Gateway-compatible vectors for plant functional genomics and proteomics. Plant J. 2006;45(4):616–29.

  69. 69.

    Crane C, Wright E, Dixon RA, Wang ZY. Transgenic Medicago truncatula plants obtained from Agrobacterium tumefaciens-transformed roots and Agrobacterium rhizogenes-transformed hairy roots. Planta. 2006;223(6):1344–54.

Download references


We thank all colleagues in our laboratory for helpful discussions and technical assistance.


This work was supported by grants from Ministry of Science and Technology of China (2016YFD0100500 and 2015CB943500), the National Natural Science Foundation of Shandong Province (ZR2018ZC0334 and ZR2019MC013) and of China (31671507 and 31371235) and the 1000-Talents Plan from China for Young Researchers. All the funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

HW, CZ, and LH conceived the study and designed the experiments. HW, ZL, YX, LK, JS, YL, CF, XW, Z-YW, CZ, and LH performed the experiments and analyzed the data. HW, CZ, and LH wrote the article. All authors have read and approved the final version of the manuscript.

Correspondence to Chuanen Zhou or Lu Han.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

The sequences of all genes involved in this study. (DOCX 31 kb)

Additional file 2:

Multiple amino acid sequences alignment of MtSPL proteins using full-length amino acid sequences. Sequences were aligned using Jalview software. (DOCX 575 kb)

Additional file 3:

Secondary structure of MtmiR156 and regulation of MtSPLs by MtmiR156. a RNA secondary structures of the MtmiR156A-MtmiR156A J were predicted by the online mfold Web Server. The nucleotides with light green color in stem-loop structures indicate the mature MtmiR156 sequences. b Multiple MtSPL genes were regulated by MtmiR156. The deoxyribonucleotide with shaded color indicates the conserved sequences targeted by MtmiR156. (DOCX 548 kb)

Additional file 4:

Molecular characterization of MtmiR156B overexpression lines. PCR analysis was performed using primer pair 35S-F/MtmiR156B-R for regenerated transgenic plants together with the positive control (35S:MtmiR156B), and negative control (Wild-type). (DOCX 168 kb)

Additional file 5:

Primers used in this study. (DOCX 23 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, H., Lu, Z., Xu, Y. et al. Genome-wide characterization of SPL family in Medicago truncatula reveals the novel roles of miR156/SPL module in spiky pod development. BMC Genomics 20, 552 (2019) doi:10.1186/s12864-019-5937-1

Download citation


  • SPL genes
  • miR156
  • Seed development
  • Spiky pod formation
  • Medicago truncatula
  • Legume