In silico genome-wide identification and comprehensive characterization of the BES1 gene family in soybean

The BES1 transcription factor family play a central role in brassinosteroid signaling pathway that regulates a wide range of plant growth and developmental processes, as well as resistances to various stresses. However, no comprehensive study of the BES1 gene family in soybean has been reported. In this work, 16 GmBES1-like genes were identified in soybean, which could be divided into two clades based on their phylogenetic relationships, gene structures and motif compositions. We then examined their duplication status and evolutionary models. The result showed that most of the GmBES1-like genes have duplicated counterparts generated from the recent Glycine WGD event, and these genes are originated from 6 distinct ancestors before the Gamma WGT event. We further studied the expression profiles of GmBES1-like genes, and found their spatio-temporal and stressed expression patterns varied tremendously. For example, GmBES1-5 and GmBES1-6 were highly expressed in almost every sample, whereas GmBES1-7 and GmBES1-8 were not expressed. Additionally, interaction network analysis revealed the presence of 3 clusters between GmBES1-like genes and other associated genes, implying that they have both the conserved and divergent functions. Lastly, we analyzed the genetic diversity of GmBES1-like genes in 302 resequenced wild, landrace and improved soybean accessions. It showed that most of these genes are well conserved, and they are not changed during domestication and improvement. These results provide insights into the characterization of GmBES1 family and lay the foundation for further functional study of such genes.