Genome-wide identification of WOX transcription factors and functional characterization of WOX4 and WOX13 involved in cold stress response in Malus baccata

Introduction Cold stress is a major abiotic threat to apple production. Malus baccata has exceptional cold hardiness and is widely used as a superior cold-resistant rootstock. The WUSCHEL-related homeobox (WOX) transcription factor family regulates plant growth, development and stress adaptation, whereas the functions of WOX genes in cold tolerance of M. baccata remain elusive. Methods In the present work, 19 MbWOX family members were identified and characterized at the genome-wide level. Evolutionary analysis, cis-element prediction, transcriptome profiling and real-time quantitative PCR (RT-qPCR) were performed to screen core cold-responsive genes. Overexpression vectors were constructed and transformed into Arabidopsis seedlings for functional verification. Results Evolutionary analysis revealed that segmental duplication drove the expansion of the MbWOX family, and these genes contained a variety of stress-responsive cis-elements. Combined transcriptome and RT-qPCR analyses confirmed that MbWOX4 and MbWOX13 were core cold-responsive genes with distinct expression patterns. The two genes participated in cold signal transduction by interacting with different transcription factor networks. Functional tests revealed that MbWOX4 and MbWOX13 isoforms differentially modulated seedling cold tolerance under low-temperature stress.