Molecular evolution of the HD-ZIP I gene family in legume genomes

• Strong purifying selections and the slow protein-level evolution were demonstrated.
• Synteny quality indicates complicated evolution processes.
• The blocks identified in this work were formed after the monocot–dicot split.
• A large-scale duplication predated speciation.
• The ancient large-scale duplication was followed by gene loss and rearrangement.

Homeodomain leucine zipper I (HD-ZIP I) genes were used to increase the plasticity of plants by mediating external signals and regulating growth in response to environmental conditions. The way genomic histories drove the evolution of the HD-ZIP I family in legume species was described; HD-ZIP I genes were searched in Lotus japonicus, Medicago truncatula, Cajanus cajan and Phaseolus vulgaris, and then divided into five clades through phylogenetic analysis. Microsynteny analysis was made based on genomic segments containing the HD-ZIP I genes. Some pairs turned out to conform with syntenic genome regions, while others corresponded to those that were inverted, expanded, or contracted after the divergence of legumes. Besides, we dated their duplications by Ks analysis and demonstrated that all the blocks were formed after the monocot–dicot split; we observed Ka/Ks ratios representing strong purifying selections in the four legume species which might have been followed by gene loss and rearrangement.