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.
Pentatomomorpha is one of the most biodiverse infraorders among the true bugs（Hemiptera： Heteroptera）. Phylogenetic relationships among the superfamilies within this infraorder have been uncertain, especially for the Eutrichophora. The previous studies were based on morphological characters, or just mitochondrial or nuclear genes, or only partial 18 S rDNA and COI. In this study, we used maximum likelihood（ML） and Bayesian inference（BI） based on massive molecular datasets（18S rDNA, 28 S rDNA, Hox and mitochondrial genes totaling 21 loci and 12,538 characters） to infer a robust phylogeny for this terrestrial group. Results strongly support the monophyly of all superfamilies; the superfamily status of Aradoidea and the following relationships：（Aradoidea ＋（Pentatomoidea ＋（Coreoidea ＋（Lygaeoidea ＋ Pyrrhocoroidea）））） in Pentatomomorpha, and（Coreoidea ＋（Lygaeoidea ＋ Pyrrhocoroidea）） in Eutrichophora. Our results suggest that sampling greater numbers of genes is an effective tool for resolving phylogenetic problems.
