Molecular phylogeny of the higher taxa of Odonata (Insecta) inferred from COI, 16S rRNA, 28S rRNA, and EF1-α sequences

In this study, we sequenced both two mitochondrial genes (COI and 16S rRNA) and nuclear genes (28S rRNA and elongation factor-1α) from 71 species of Odonata that represent 7 superfamilies in 3 suborders. Phylogenetic testing for each two concatenated gene sequences based on function (ribosomal vs protein-coding genes) and origin (mitochondrial vs nuclear genes) proved limited resolution. Thus, four concatenated sequences were utilized to test the previous phylogenetic hypotheses of higher taxa of Odonata via Bayesian inference (BI) and maximum likelihood (ML) algorithms, along with the data partition by the BI method. As a result, three slightly different topologies were obtained, but the BI tree without partition was slightly better supported by the topological test. This topology supported the suborders Anisoptera and Zygoptera each being a monophyly, and the close relationship of Anisozygoptera to Anisoptera. All the families represented by multiple taxa in both Anisoptera and Zygoptera were consistently revealed to each be a monophyly with the highest nodal support. Unlike consistent and robust familial relationships in Zygoptera those of Anisoptera were partially unresolved, presenting the following relationships: ((((Libellulidae + Corduliidae) + Macromiidae) + Gomphidae + Aeshnidae) + Anisozygoptera) + (((Coenagrionidae + Platycnemdidae) + Calopterygidae) + Lestidae). The subfamily Sympetrinae, represented by three genera in the anisopteran family Libellulidae, was not monophyletic, dividing Crocothemis and Deielia in one group together with other subfamilies and Sympetrum in another independent group.