Influence of data conflict and molecular phylogeny of major clades in Cimicomorphan true bugs (Insecta: Hemiptera: Heteroptera)

Cimicomorpha, which consists of 16 families representing more than 19,400 species, is the largest infraorder in Heteroptera, Insecta. We present the first molecular phylogenetic investigation of family relationships of Cimicomorpha, including 46 taxa from 12 of 16 Cimicomorphan families. Three genes, with a total of 3277 bp of sequence data (nuclear 18S rDNA: 2022 bp, 28S rDNA: 755 bp, and mitochondrial 16S rDNA: 498 bp) were analyzed. Data partitions were analyzed separately and in combination, by employing ML (maximum likelihood), MP (maximum parsimony), and Bayesian methods. As saturation was detected in substitutions of 16S rDNA, influence of data conflict in combined analyses was further explored by three methods: the incongruence length difference (ILD) test, the partitioned Bremer support (PBS), and the partition addition bootstrap alteration approach (PABA). PBS and PABA approaches suggested that 16S rDNA was not very suitable for addressing relationships at this level in Cimicomorpha. Our results also supported the nabid–cimicoid lineage for Cimicoidea proposed by Schuh and Štys [Schuh, R.T., Štys, P., 1991. Phylogenetic analysis of Cimicomorphan family relationships (Heteroptera). J. NY Entomol. Soc. 99 (3), 298–350]. Data incongruence and the utility of the three genes were briefly discussed.