Identification and Evolutionary Analysis of Cotton (Gossypium hirsutum) WOX Family Genes and Their Potential Function in Somatic Embryogenesis

WUSCHEL-related homeobox (WOX) proteins participate profoundly in plant development and stress responses. As the difficulty of somatic embryogenesis severely constrains cotton genetic modification, in this study, we identified and comprehensively analyzed WOX genes in cotton. As a result, 40 WOX genes were identified in the upland cotton genome. All these cotton WOX genes were classified into three clades, ancient, intermediate, and modern clades, based on the phylogenetic analysis of previous studies. The majority (24) of the cotton WOX genes belonged to the modern clade, in which all gene members contain the vital functional domain WUS-box, which is necessary for plant stem cell regulation and maintenance. Collinearity analysis indicated that the WOX gene family in cotton expanded to some degree compared to Arabidopsis, especially in the modern clade. Genome duplication and segmental duplication may greatly contribute to expansion. Hormone-response- and abiotic-stress-response-related cis-acting regulatory elements were widely distributed in the promoter regions of cotton WOX genes, suggesting that the corresponding functions of stress responses and the participation of development processes were involved in hormone responses. By RNA sequencing, we profiled the expression patterns of cotton WOX genes in somatic embryogenesis. Only about half of cotton WOX genes were actively expressed during somatic embryogenesis; different cotton WOX genes may function in different development stages. The most representative, GhWOX4 and GhWOX13, may function in almost all stages of somatic embryogenesis; GhWOX2 and GhWOX9 function in the late stages of embryo patterning and embryo development during cotton somatic embryogenesis. Co-expression analysis showed that the cotton WOXs co-expressed with genes involved in extensive genetic information processing, including DNA replication, DNA repair, homologous recombination, RNA transport, protein processing, and several signaling and metabolism pathways, in which plant hormones signal transduction, MAPK signaling pathways, phosphatidylinositol signaling systems, and ABC transporters, as well as the metabolism of fatty acid; valine, leucine, and isoleucine biosynthesis; and cutin, suberine, and wax biosynthesis, were most significantly enriched. Taken together, the present study provides useful information and new insights into the functions of cotton WOX genes during somatic embryogenesis. The specific regulatory roles of some WOX genes in somatic embryogenesis are worthy of further functional research.

Automated summary

In this study, we identified and analyzed WOX genes in cotton and found 40 WOX genes in the upland cotton genome. These genes were classified into three clades: ancient, intermediate, and modern clades. The majority of the cotton WOX genes belonged to the modern clade and exhibited expansion compared to Arabidopsis. Promoter analysis revealed the presence of hormone-response- and abiotic-stress-response-related cis-acting regulatory elements in the cotton WOX genes. Expression analysis during somatic embryogenesis showed that different WOX genes functioned in different stages of development. Co-expression analysis revealed the involvement of WOX genes in multiple genetic processes and metabolic pathways. This study provides valuable information and insights into the functions of cotton WOX genes during somatic embryogenesis.