Pseudomonas aeruginosa is a Gram-negative bacterium capable of widespread niches, which is also one of the main bacteria that cause patient infection. The metabolic diversity of Pseudomonas aeruginosa is an essential factor in adapting to a variety of environments. Based on the previous studies, adaptive genetic variation in the glycerol kinase GlpK, the glycerol 3-phosphotransferase, contributes to the fitness of bacteria in human bodies, such as Mycobacterium tuberculosis and Escherichia coli. Thus, this study aimed to explore the molecular evolution and function of glpK in P. aeruginosa. Using extensive population genomic data, we have identified the prevalence of two glpK copies in P. aeruginosa that clustered into distinct branches, which were later known as Clade 1 and 2. The evolution analysis revealed that glpK in Clade 1 derived from an ancestral P. aeruginosa species and the other from an ancient horizontal gene transfer event. In addition, we confirmed that the GlpK in Clade 2 still retained glycerol kinase activity but was much weaker than that of GlpK in Clade 1. We demonstrated the importance of the critical amino acid Q70 in GlpK glycerol kinase activity by point mutation. Furthermore, Co-expression network analysis implied that the two glpK copies of P. aeruginosa regulate separate networks and may be a strategy to improve fitness in P. aeruginosa.