Abstract
Evolution of unisexual flowers involves extreme changes in floral development. Spinach is one of the species to discern the formation and evolution of dioecy. MADS-box gene family is involved in regulation of floral organ identity and development and in many other plant developmental processes. However, there is no systematic analysis of MADS-box family genes in spinach. A comprehensive genome-wide analysis and transcriptome profiling of MADS-box genes were undertaken to understand their involvement in unisexual flower development at different stages in spinach. In total, 54 MADS-box genes found to be unevenly located across 6 chromosomes and can be divided into type I and type II genes. Twenty type I MADS-box genes are subdivided into Mα, Mβ and Mγ subgroups. While thirty-four type II SoMADSs consist of 3 MIKC*, and 31 MIKCC -type genes including sixteen floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in spinach. Gene structure, motif distribution, physiochemical properties, gene duplication and collinearity analyses for these genes are performed in detail. Promoters of both types of SoMADS genes contain mainly MeJA and ABA response elements. Expression profiling indicated that MIKCc genes exhibited more dynamic and intricate expression patterns compared to M-type genes and the majority of type-II genes AP1, SVP, and SOC1 sub-groups showed female flower-biased expression profiles, suggesting their role in carpel development, while PI showed male-biased expression throughout flower developmental stages, suggesting their role in stamen development. These results provide genomic resources and insights into spinach dioecious flower development and expedite spinach improvement.