A report on the complete mitochondrial genome of the trematode Azygia robusta Odhner, 1911, its new definitive host from the Russian Far East, and unexpected phylogeny of Azygiidae within Digenea, as inferred from mitogenome sequences.

New data on the complete mitochondrial genome of Azygia robusta (Azygiidae) were obtained by the next-generation sequencing (NGS) approach. The mitochondrial DNA (mtDNA) of A. robusta had a length of 13 857 bp and included 12 protein-coding genes, two ribosomal genes, 22 transfer RNA genes, and two non-coding regions. The nucleotide sequences of the complete mitochondrial genomes of two A. robusta specimens differed from each other by 0.12 ± 0.03%. Six of 12 protein-coding genes demonstrated intraspecific variation. The difference between the nucleotide sequences of the complete mitochondrial genomes of A. robusta and Azygia hwangtsiyui was 26.95 ± 0.35%; the interspecific variation of protein-coding genes between A. robusta and A. hwangtsiyui ranged from 20.5 ± 0.9% (cox1) to 30.7 ± 1.2% (nad5). The observed gene arrangement in the mtDNA sequence of A. robusta was identical to that of A. hwangtsiyui. Codon usage and amino acid frequencies were highly similar between A. robusta and A. hwangtsiyui. The results of phylogenetic analyses based on mtDNA protein-coding regions showed that A. robusta is closely related to A. hwangtsiyui (belonging to the same suborder, Azygiida) that formed a distinct early-diverging branch relative to all other Digenea. A preliminary morphological analysis of paratypes of the two azygiid specimens studied showed visible morphological differences between them. The specimen extracted from Sakhalin taimen (Parahucho perryi) was most similar to A. robusta. Thus, we here provide the first record of a new definitive host, P. perryi, for A. robusta and also molecular characteristics of the trematode specimens.

Characterization of complete mitochondrial genome and ribosomal operon for Carassotrema koreanum 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 Carassotrema koreanum (Digenea: Haploporata: Haploporidae), an intestinal parasite of Carassius auratus, using next-generation sequencing. The mtDNA of C. koreanum 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 C. koreanum 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 C. koreanum was closely related to Parasaccocoelium mugili, 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, Schistosoma spp. were differentiated from the other members of Digenea and the members of Pronocephalata were localized within the same cluster. Carassotrema koreanum and P. mugili fell within different clusters. The grouping of C. koreanum and P. mugili 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 Dicrocoelium spp. in the presence of TTG start codon of the mitochondrial nad5 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.