Characterization of complete mitochondrial genome and ribosomal operon for Park, 1938 (Digenea: Haploporidae) by means of next-generation sequencing data.

We obtained new data on the complete mitochondrial DNA (mtDNA) and the ribosomal operon of the trematode (Digenea: Haploporata: Haploporidae), an intestinal parasite of , using next-generation sequencing. The mtDNA of . contained 13,965 bp, including 12 protein-coding genes, two ribosomal genes, 22 transport RNA (tRNA) genes and a non-coding region. The ribosomal operon of . was 10,644 bp in length, including ETS1 (1449 bp), 18S ribosomal RNA (rRNA) gene (1988 bp), ITS1 ribosomal DNA (rDNA) (558 bp), 5.8S rRNA gene (157 bp), ITS2 rDNA (274 bp), 28S rRNA gene (4152 bp) and ETS2 (2066 bp). Phylogenetic analysis based on mtDNA protein-coding regions showed that . was closely related to , a species from the same suborder Haploporata. Bayesian phylogenetic tree topology was the most reliable and confirmed the validity of the Haploporata. The results of sequence cluster analysis based on codon usage bias demonstrated some agreement with the results of the phylogenetic analysis. In particular, spp. were differentiated from the other members of Digenea and the members of Pronocephalata were localized within the same cluster. and . fell within different clusters. The grouping of . and . within the same cluster was obtained on the basis of frequencies of 13 specified codons, of which three codon pairs were degenerate. A similarity was found between two haploporid species and two spp. in the presence of TTG start codon of the mitochondrial 5 gene. Our results confirmed the taxonomical status of the Haploporata identified in the previous studies and revealed some characteristic features of the codon usage in its representatives.