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.
