Hexokinases (HXKs) are key enzymes in glucose metabolism and play essential roles in energy homeostasis. In vertebrates, the HXK family has undergone extensive expansion and diversification, giving rise to multiple paralogs with distinct structural and regulatory features. Here, we present a comprehensive evolutionary analysis of the HXK gene family across Metazoa, integrating phylogenetic reconstruction, molecular dating, synteny analysis, structural comparisons, and selection pressure inference. Our results indicate that the 100 kDa hexokinase paralogs originated from domain duplication and fusion events in early vertebrates, followed by lineage-specific duplications giving rise to HXK1, HXK2, HXK3, HXK4, and HXK5. Structural analysis revealed that only HXK2 retains catalytic activity in both N- and C-terminal domains, while the other paralogs exhibit divergence and loss of N-terminal catalytic function. Selection analyses across 470 mammalian genomes revealed paralog- and lineage-specific signatures of episodic positive selection, particularly in HXK3, suggesting adaptive evolution in response to distinct metabolic demands. These findings provide new insights into the functional evolution and specialization of HXKs in vertebrate metabolism and highlight the interplay between gene duplication, structural adaptation, and selective pressures in shaping metabolic diversity.
Jardim‐Messeder et al. (Sat,) studied this question.