Genomic organization, phylogenetic comparison, and expression profiles of the SPL family genes and their regulation in soybean

SQUAMOSA Promoter-Binding Protein-Like (SPL) genes form a major family of plant-specific transcription factors and play an important role in plant growth and development. In this study, we report the identification of 41 SPL genes (GmSPLs) in the soybean genome. Phylogenetic analysis revealed that these genes were divided into five groups (groups 1-5). Further, exon/intron structure and motif composition revealed that the GmSPL genes are conserved within their same group. The N-terminal zinc finger 1 (Zn1) of the SBP domain was a CCCH (Cys3His1) and the C terminus zinc finger 2 (Zn2) was a CCHC (Cys2HisCys) type. The 41 GmSPL genes were distributed unevenly on 17 of the 20 chromosomes, with tandem and segmental duplication events. We found that segmental duplication has made an important contribution to soybean SPL gene family expansion. The Ka/Ks ratios revealed that the duplicated GmSPL genes evolved under the effect of purifying selection. In addition, 17 of the 41 GmSPLs were found as targets of miR156; these might be involved in their posttranscriptional regulation through miR156. Importantly, RLM-RACE analysis confirmed the GmmiR156-mediated cleavage of GmSPL2a transcript in 2-4 mm stage of soybean seed. Alternative splicing events in 9 GmSPLs were detected which produces transcripts and proteins of different lengths that may modulate protein signaling, binding, localization, stability, and other properties. Expression analysis of the soybean SPL genes in various tissues and different developmental stages of seed suggested distinct spatiotemporal patterns. Differences in the expression patterns of miR156-targeted and miR156-non-targeted soybean SPL genes suggest that miR156 plays key functions in soybean development. Our results provide an important foundation for further uncovering the crucial roles of GmSPLs in the development of soybean and other biological processes.