Heteroptera are among the most diverse hemimetabolous insects. Seven infraorders have been recognized within this suborder of Hemiptera. Apart from the well-established sister-group relationship between Cimicomorpha and Pentatomomorpha (= Terheteroptera), the two terminal lineages, the relationships among the other five infraorders are still controversial, of which three (Gerromorpha, Nepomorpha and Leptopodomorpha) are intimately connected to aquatic environments. However, the various and often conflicting available phylogeny hypotheses do not offer a clear background for a connection between diversification and palaeoenvironments. In this study, a molecular data set representing 79 taxa and 10 149 homologous sites is used to infer the phylogenetic relationships within Heteroptera. Bayesian inference, maximum-likelihood and maximum parsimony analyses were employed. The results of phylogenetic inferences largely confirm the widely accepted phylogenetic context. Estimation of the divergence time based on the phylogenetic results revealed that Gerromorpha, Nepomorpha and Leptopodomorpha originated successively during the period from the Late Permian to Early Triassic (269–246 Ma). This timescale is consistent with the origin and radiation time of various aquatic holometabolans. Our results indicate that the aquatic and semi-aquatic true bugs evolved under environmental conditions of high air temperature and humidity in an evolutionary scenario similar to that of the aquatic holometabolans.

The sequences of 18S and 28S rDNAs have been used as molecular markers to resolve phylogenetic relationships of Heteroptera for two decades. The complete sequences of 18S rDNAs have been used in many studies, while in most studies only partial sequences of 28S rDNAs have been used due to technical difficulties of amplifying the complete lengths. In this study, we amplified the complete 18S and 28S rDNA sequences of Eurydema maracandica Oshanin, 1871, and reconstructed the secondary structure models of the corresponding rRNAs. In addition, and more importantly, all of the length variable regions of 18S rRNA were compared among 37 families of Heteroptera based on 140 sequences, and the D3 region of 28S rRNA was compared among 51 families based on 84 sequences. It was found that 8 length variable regions could potentially serve as molecular synapomorphies for some monophyletic groups. Therefore discoveries of more molecular synapomorphies for specific clades can be anticipated from amplification of complete 18S and 28S rDNAs of more representatives of Heteroptera.