The first complete mitochondrial genome of a parasitic isopod supports Epicaridea Latreille, 1825 as a suborder and reveals the less conservative genome of isopods.

The complete mitochondrial genome sequence of the holoparasitic isopod Gyge ovalis (Shiino, 1939) has been determined. The mitogenome is 14,268 bp in length and contains 34 genes: 13 protein-coding genes, two ribosomal RNA, 19 tRNA and a control region. Three tRNA genes (trnE, trnI and trnS1) are missing. Most of the tRNA genes show secondary structures which derive from the usual cloverleaf pattern except for trnC which is characterised by the loss of the DHU-arm. Compared to the isopod ground pattern and Eurydice pulchra Leach, 1815 (suborder Cymothoida Wägele, 1989), the genome of G. ovalis shows few differences, with changes only around the control region. However, the genome of G. ovalis is very different from that of non-cymothoidan isopods and reveals that the gene order evolution in isopods is less conservative compared to other crustaceans. Phylogenic trees were constructed using maxiumum likelihood and Bayesian inference analyses based on 13 protein-coding genes. The results do not support the placement of G. ovalis with E. pulchra and Bathynomus sp. in the same suborder; rather, G. ovalis appears to have a closer relationship to Ligia oceanica (Linnaeus, 1767), but this result suggests a need for more data and further analysis. Nevertheless, these results cast doubt that Epicaridea Latreille, 1825 can be placed as an infraorder within the suborder Cymothoida, and Epicaridea appears to also deserve subordinal rank. Further development of robust phylogenetic relationships across Isopoda Latreille, 1817 will require more genetic data from a greater diversity of taxa belonging to all isopod suborders.