WUSCHEL-related homeobox (WOX) is a family of transcription factors that are unique to plants and is characterized by the presence of a homeodomain. The WOX transcription factor plays an important role in regulating plant growth and development and the response to abiotic stress. Soybean is one of the most important oil crops worldwide. In this study, based on the available genome data of soybean, the WOX gene family was identified by bioinformatics analysis. The chromosome distribution, gene and protein structures, phylogenetic relationship and gene expression patterns of this family were comprehensively compared. The results showed that a total of 33 putative WOX genes in the soybean genome were found and then designated as GmWOX1- GmWOX33, which were distributed across 19 chromosomes except chromosome 16. Multiple sequence analysis of the GmWOX gene family revealed a highly conserved homeodomain. Phylogenetic tree analysis showed that 33 WOX genes could be divided into three major clades (modern/WUS, intermediate and ancient) in soybean. Of these 33 WOX genes, some showed differential expression patterns in the tested tissues (leaves, pods, unopen and open flowers, nodules, seed, roots, root hairs, stems, shoot apical meristems and shoot tips). In addition, the expression profile and qRT-PCR analysis showed that most of the GmWOX genes responded to different abiotic stress treatments (cold and drought). According to the expression pattern of GmWOX genes in the high regeneration capacity soybean material P3, overexpression of GmWOX18 was selected for function analysis. The overexpression of GmWOX18 increased the regeneration ability of clustered buds. The results will provide valuable information for further studies on the roles of WOX genes in regulating soybean growth, development and responses to abiotic stress, as well as a basis for the functional identification and analysis of WOX genes in soybean.