The complete mitochondrial genome of Gastrothylax crumenifer (Gastrothylacidae, Trematoda) and comparative analyses with selected trematodes

Evolutionary affinities and morphological characterization of the enigmatic Zonocotyle bicaecata (Trematoda: Paramphistomoidea: Zonocotylidae) from the Upper Paraná River basin

Neotropical fish amphistomes represent a highly diverse group within the Paramphistomoidea, with wide distribution across major South American hydrological drainages. However, the limited molecular characterization of these taxa has impeded a comprehensive assessment of their evolutionary relationships and the systematic relevance of morphological features in classification schemes. Our study, based on the critical evaluation of the type material of both nominal species of Zonocotyle (type genus of the monotypic Zonocotylidae), and newly collected specimens of Zonocotyle bicaecata from Steindachnerina insculpta (Curimatidae) in the Upper Paraná River basin, Brazil, presents a morphological reappraisal of Z. bicaecata and provides molecular data (28S rDNA, ITS1-5.8S-ITS2 region, and COI mtDNA) to assess its phylogenetic relationships. Our phylogenetic analyses confirm this species belongs to the Paramphistomoidea. The most comprehensive analyses (based on 28S and COI) further indicate a close relationship with other fish amphistomes from the Neotropical region. Additionally, we emphasized the necessity for a new classification within Paramphistomoidea and briefly discussed the host range of Zonocotyle among curimatid fish hosts.

Amphistomes

Amphistomes, commonly referred to as 'stomach' or 'rumen' flukes because of the localization of these flukes in the stomach of ruminants, are digenetic trematodes distinguished by the absence of an oral sucker and the position of the ventral sucker or acetabulum at the posterior end of the body. The body is characterized by a leaf-like fleshy structure, pink or red with a large posterior sucker. Amphistomes are an important group of parasites since they cause 'amphistomiasis' (variously known as paramphistomosis/amphistomosis), a serious disease of great economic importance in ruminants worldwide. These parasites have a broad spectrum of definitive hosts together with a wide geographical distribution. Though they form a continuous evolutional lineage from fishes to mammals, amphistomes mainly inhabit the rumen and reticulum of ruminant mammals, while some species occur in the large intestine or parenteric sites of ruminants, pigs, equines and man.

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.

The complete mitochondrial genome of Ogmocotyle ailuri: gene content, composition and rearrangement and phylogenetic implications

Trematodes of the genus Ogmocotyle are intestinal flukes that can infect a variety of definitive hosts, resulting in significant economic losses worldwide. However, there are few studies on molecular data of these trematodes. In this study, the mitochondrial (mt) genome of Ogmocotyle ailuri isolated from red panda (Ailurus fulgens) was determined and compared with those from Pronocephalata to investigate the mt genome content, genetic distance, gene rearrangements and phylogeny. The complete mt genome of O. ailuri is a typical closed circular molecule of 14 642 base pairs, comprising 12 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and 2 non-coding regions. All genes are transcribed in the same direction. In addition, 23 intergenic spacers and 2 locations with gene overlaps were determined. Sequence identities and sliding window analysis indicated that cox1 is the most conserved gene among 12 PCGs in O. ailuri mt genome. The sequenced mt genomes of the 48 Plagiorchiida trematodes showed 5 types of gene arrangement based on all mt genome genes, with the gene arrangement of O. ailuri being type I. Phylogenetic analysis using concatenated amino acid sequences of 12 PCGs revealed that O. ailuri was closer to Ogmocotyle sikae than to Notocotylus intestinalis. These data enhance the Ogmocotyle mt genome database and provide molecular resources for further studies of Pronocephalata taxonomy, population genetics and systematics.

Differentiating paramphistome species in cattle using DNA barcoding coupled with high-resolution melting analysis (Bar-HRM)

Paramphistomosis is caused by paramphistome or amphistome parasites, including Fischoederius elongatus, Gastrothylax crumenifer, Orthocoelium parvipapillatum, and Paramphistomum epiclitum. The control and prevention of these parasite outbreaks are difficult because of the wide occurrence of these species. Besides, the clinical manifestations and their egg characteristics are similar to those of other intestinal flukes in the paramphistome group, leading to misdiagnosis. Here, we employed DNA barcoding using NADH dehydrogenase (ubiquinone, alpha 1) (ND1) and cytochrome c oxidase subunit I (COI), coupled with high-resolution melting analysis (Bar-HRM), for species differentiation. As a result, ParND1_3 and ParCOI4 resulted in positive amplification in the paramphistomes and Fasciola gigantica, with significantly different melting curves for each species. The melting temperatures of each species obtained clearly differed. Regarding sensitivity, the limit of detection (LoD) for all species of paramphistomes was 1 pg/µl. Our findings suggest that Bar-HRM using ParND1_3 is highly suitable for the differentiation of paramphistome species. This approach can be used in parasite detection and epidemiological studies in cattle.

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.

•

The complete mt genomes of Diplodiscus japonicus and D. mehari were determined first time.

•

There is only one NCR in Diplodiscus japonicus and D. mehari complete mt genomes.

•

Phylogenetic analyses revealed two monophyletic groups for the flukes parasitic in amphibians and ruminants.

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.

First next-generation sequencing data for Haploporidae (Digenea: Haploporata): characterization of complete mitochondrial genome and ribosomal operon for Parasaccocoelium mugili Zhukov, 1971

The first data on a whole mitochondrial genome of Haploporidae, Parasaccocoelium mugili (Digenea: Haploporata: Haploporidae) was generated using the next-generation sequencing (NGS) approach. We sequenced the complete mitochondrial DNA (mtDNA) and ribosomal operon of Parasaccocoelium mugili, intestine parasite of mullet fish. The mtDNA of P. mugili contained 14,021 bp, including 12 protein-coding genes, two ribosomal genes, 22 tRNA genes, and non-coding region. The ribosomal operon of P. mugili was 8308 bp in length, including 18S rRNA gene (1981 bp), ITS1 rDNA (955 bp), 5.8S rRNA gene (157 bp), ITS2 rDNA (268 bp), 28S rRNA gene (4180 bp), and ETS (767 bp). We used the mtDNA protein-coding regions to make phylogenetic reconstructions of Haploporidae. Additionally, we performed the sequence cluster analysis based on codon usage bias of most of currently available mitochondrial genome data for trematodes. The observed gene arrangement in mtDNA sequence of P. mugili is identical to those of Plagiorchis maculosus (Rudolphi, 1802). Results of maximum likelihood (ML) phylogenetic analysis showed that P. mugili was closely related to Paragonimus species from the suborder Xiphidiata. The results of sequence cluster analysis based on codon usage bias showed that P. mugili has the highest similarity with Plagiorchis maculosus (Xiphidiata). Our results do not contradict to proposing a new suborder for Haploporoidea-Haploporata. On the basis of obtained results, the relationship between mitochondrial protein-coding gene rearrangements and synonymous nucleotide substitutions in mitochondrial genomes has been suggested.

Characterization of the Complete Mitochondrial Genome of Fischoederius elongatus Derived from Cows in Shanghai, China

A study was conducted to reveal the characterization of the complete mitochondrial genome of Fischoederius elongatus derived from cows in Shanghai, China. Results indicated that the complete mt genome of F. elongatus was 14,288 bp and contained 12 protein-coding genes (cox1-3, nad1-6, nad4L, atp6, and cytb), 22 transfer RNA genes, and two ribosomal RNA genes (l-rRNA and s-rRNA). The overall A + T content of the mt genome was 63.83%, and the nucleotide composition was A (19.83%), C (9.75%), G (26.43%), and T (44.00%). A total of 3284 amino acids were encoded by current F. elongatus isolate mt genome, TTT (Phe) (9.84%) and TTG (Leu) (7.73%) codon were the most frequent amino acids, whereas the ACC (Thr) (0.06%), GCC (Ala) (0.09%), CTC (Leu) (0.09%), and AAC (Asn) (0.09%) codon were the least frequent ones. At the third codon position of F. elongatus mt protein genes, T (50.82%) was observed most frequently and C (5.85%) was the least one. The current results can contribute to epidemiology diagnosis, molecular identification, taxonomy, genetic, and drug development researches about this parasite species in cattle.

Characterization of the complete mitochondrial genome of the echinostome Echinostoma miyagawai and phylogenetic implications

Echinostomes are important intestinal foodborne parasites. Despite their significance as pathogens, characterization of the molecular biology and phylogenetics of these parasites are limited. In the present study, we determined the entire mitochondrial (mt) genome of the echinostome Echinostoma miyagawai (Hunan isolate) and examined the phylogenetic relationship with selected members of the suborder Echinostomata. The complete mt genome of E. miyagawai (Hunan isolate) was 14,468 bp in size. This circular mt genome contained 12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one non-coding region. The gene order and genomic content were identical with its congeners. Phylogenetic analyses (maximum parsimony, maximum likelihood, and Bayesian inference) based on the concatenated amino acid sequences of 12 protein-coding genes strongly supported monophyly for the genus Echinostoma; however, they rejected monophyly for the family Echinostomatidae and the genus Fasciola. The mt genomic data described in this study provides useful genetic markers for studying the population genetics, molecular biology, and phylogenetics of these echinostomes.

The complete mitochondrial genome of Paragonimus ohirai (Paragonimidae: Trematoda: Platyhelminthes) and its comparison with P. westermani congeners and other trematodes

We present the complete mitochondrial genome of Paragonimus ohirai Miyazaki, 1939 and compare its features with those of previously reported mitochondrial genomes of the pathogenic lung-fluke, Paragonimus westermani, and other members of the genus. The circular mitochondrial DNA molecule of the single fully sequenced individual of P. ohirai was 14,818 bp in length, containing 12 protein-coding, two ribosomal RNA and 22 transfer RNA genes. As is common among trematodes, an atp8 gene was absent from the mitogenome of P. ohirai and the 5' end of nad4 overlapped with the 3' end of nad4L by 40 bp. Paragonimusohirai and four forms/strains of P. westermani from South Korea and India, exhibited remarkably different base compositions and hence codon usage in protein-coding genes. In the fully sequenced P. ohirai individual, the non-coding region started with two long identical repeats (292 bp each), separated by tRNAGlu . These were followed by an array of six short tandem repeats (STR), 117 bp each. Numbers of the short tandem repeats varied among P. ohirai individuals. A phylogenetic tree inferred from concatenated mitochondrial protein sequences of 50 strains encompassing 42 species of trematodes belonging to 14 families identified a monophyletic Paragonimidae in the class Trematoda. Characterization of additional mitogenomes in the genus Paragonimus will be useful for biomedical studies and development of molecular tools and mitochondrial markers for diagnostic, identification, hybridization and phylogenetic/epidemiological/evolutionary studies.

New approach to molecular characterization of Paramphistomum cervi and Carmyerius gregarius and comparative analyses with selected trematodes

The two ruminant parasites, Paramphistomum cervi and Carmyerius gregarius, were collected from fresh-slaughtered native cattle at local abattoirs in Sadat district, Menoufia province and identified morphologically, then molecularly by sequencing the nucleotides of 18S ribosomal RNA gene (18S rRNA). The nucleotide sequences of the two isolates were 456 (P. cervi) and 401 bases for (C. gregarius). The data were used along with those of several other helminth species from the GenBank to identify these two species genetically. The nucleotide sequences were aligned using multiple sequence alignments of nucleotides by Clustal W 12.1 V and construct their relationship. Neighbor-joining analytical method was used showing sister relationship between C. gregarius from Sadat district and Gastrodiscoides hominis (EF027096) with relative identity of (98%) due to the presence of single nucleotides polymorphisms (SNPs) in the form of indels as nine nucleotides positions. But when clustering of P. cervi Sadat isolate with Paramphistomoidea sp. S4 isolate P5 (GU735643), this relationship shows complete identity (99%) between them. The homology and diversity was done using Bayesian analyses in MrBayes v3.1. This work will give a useful guide for other researchers for the molecular taxonomic position of Paramphistomatidae spp. in Sadat district among the different species around the world.

The complete mitochondrial genomes of Gnathostoma doloresi from China and Japan

Gnathostoma doloresi is one of the neglected pathogens causing gnathostomiasis. Although this zoonotic parasite leads to significant socioeconomic concerns globally, little is known of its genetics and systematics. In the present study, we sequenced and characterized the complete mitochondrial (mt) genomes of G. doloresi isolates from China and Japan. The lengths of the mt genomes of the G. doloresi China and Japan isolates are 13,809 and 13,812 bp, respectively. Both mt genomes encode 36 genes, including 12 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. The gene order, transcription direction, and genome content are identical with its congener G. spinigerum. Phylogenetic analyses based on concatenated amino acid sequences of 12 PCGs by Bayesian inference (BI) indicated that G. doloresi are closely related to G. spinigerum. Our data provide an invaluable resource for studying the molecular epidemiology, phylogenetics, and population genetics of Gnathostoma spp. and should have implications for further studies of the diagnosis, prevention, and control of gnathostomiasis in humans and animals.