Characterization of the complete mitochondrial genomes of Diplodiscus japonicus and Diplodiscus mehari (Trematoda: Diplodiscidae): Comparison with the members of the superfamily Paramphistomoidea and phylogenetic implication

Diplodiscus japonicus and Diplodiscus mehari (Trematoda: Diplodiscidae) are two important parasites in wood frogs, which have large infection rates and essential importance of ecology, economy and society. In this study, the complete mitochondrial (mt) genomes of D. japonicus and D. mehari were sequenced, then compared with other related trematodes in the superfamily Paramphistomoidea. The complete circular mt sequence of D. japonicus and D. mehari were 14,210 bp and 14,179 bp in length, respectively. Both mt genomes comprised 36 functional subunits, consisting of 12 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and one non-coding region. The mt genes of D. japonicus and D. mehari were transcribed in the same direction, and the gene arrangements were identical to those of Paramphistomoidea trematodes. In the 12 PCGs, GTG was the most common initiation codon, whereas TAG was the most common termination codon. All tRNAs had a typical cloverleaf structure except tRNA Ser1. A comparison with related Paramphistomoidea trematode mt genomes suggested that the cox1 gene of D. mehari was the longest in these trematodes. Phylogenetic analyses revealed that Paramphistomoidea trematodes formed a monophyletic branch, Paramphistomidae and Gastrothylacidae were more closely related than Diplodiscidae. And the further analysis with Pronocephalata branch found that the flukes parasitic in amphibians (frogs) formed one group, and the flukes from ruminants (cattle, sheep, ect) formed another group. Our study demonstrated the importance of sequencing mt genomes of D. japonicus and D. mehari, which will provide significant molecular resources for further studies of Paramphistomoidea taxonomy, population genetics and systematics.