The phylogenetic relationships among the Nepomorpha (Hemiptera: Heteroptera) were analyzed based on nuclear protein coding genes-partial sequences of four Hox genes: abd-A, Dfd, Ubx and pb. Fifteen taxa, of which eleven were ingroups representing ten of the eleven families, were selected to produce a phylogenetic inference of nepomorphan superfamilial or familial relationships via parsimony and Bayesian techniques. Our results supported the monophylies of Nepomorpha, Naucoroidea (Aphelocheiridae + Naucoridae), Nepoidea (Belostomatidae + Nepidae), Ochteroidea (Ochteridae + Gelastocoridae) and Pieoidea (Pleidae + Helotrephidae); the most basal lineage of Ochteroidea; Notonectoidea contained Notonectidae only and formed a new sister relationship with (Pleoidea + Naucoroidea); and the sister relationship of (Nepoidea + Corixoidea). The Hox genes are suitable molecular markers to resolve phylogenetic relationships at superfamily or family levels in Heteroptera.