Complete mitochondrial genome of Haplorchis taichui and comparative analysis with other trematodes

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 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.

The complete mitogenome of the Asian lung fluke Paragonimus skrjabini miyazakii and its implications for the family Paragonimidae (Trematoda: Platyhelminthes)

The complete circular mitogenome of Paragonimus skrjabini miyazakii (Platyhelminthes: Paragonimidae) from Japan, obtained by PacBio long-read sequencing, was 17 591 bp and contained 12 protein-coding genes (PCGs), 2 mitoribosomal RNA and 22 transfer RNA genes. The atp8 gene was absent, and there was a 40 bp overlap between nad4L and nad4. The long non-coding region (4.3 kb) included distinct types of long and short repeat units. The pattern of base usage for PCGs and the mtDNA coding region overall in Asian and American Paragonimus species (P. s. miyazakii, P. heterotremus, P. ohirai and P. kellicotti) and the Indian form of P. westermani was T > G > A > C. On the other hand, East-Asian P. westermani used T > G > C > A. Five Asian and American Paragonimus species and P. westermani had TTT/Phe, TTG/Leu and GTT/Val as the most frequently used codons, whereas the least-used codons were different in each species and between regional forms of P. westermani. The phylogenetic tree reconstructed from a concatenated alignment of amino acids of 12 PCGs from 36 strains/26 species/5 families of trematodes confirmed that the Paragonimidae is monophyletic, with 100% nodal support. Paragonimus skrjabini miyazakii was resolved as a sister to P. heterotremus. The P. westermani clade was clearly separate from remaining congeners. The latter clade was comprised of 2 subclades, one of the East-Asian and the other of the Indian Type 1 samples. Additional mitogenomes in the Paragonimidae are needed for genomic characterization and are useful for diagnostics, identification and genetic/ phylogenetic/ epidemiological/ evolutionary studies of the Paragonimidae.

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.

Mitophylogenomics of the zoonotic fluke Echinostoma malayanum confirms it as a member of the genus Artyfechinostomum Lane, 1915 and illustrates the complexity of Echinostomatidae systematics

The complete mitochondrial genome (mitogenome or mtDNA) of the trematode Echinostoma malayanum Leiper, 1911 was fully determined and annotated. The circular mtDNA molecule comprised 12 protein-coding genes (PCGs) (cox1 - 3, cob, nad1 - 6, nad4L, atp6), two mitoribosomal RNAs (MRGs) (16S or rrnL and 12S or rrnS), and 22 transfer RNAs (tRNAs or trn), and a non-coding region (NCR) rich in long and short tandem repeats (5.5 LRUs/336 bp/each and 7.5 SRUs/207 bp/each). The atp8 gene is absent and the 3' end of nad4L overlaps the 5' end of nad4 by 40 bp. Special DHU-arm missing tRNAs for Serine were found for both tRNA^Ser1(AGN) and tRNA^Ser2(UCN). Codons of TTT (for phenylalanine), TTG (for leucine), and GTT (for valine) were the most, and CGC (for Arginine) was the least frequently used. A similar usage pattern was seen in base composition, AT and GC skewness for PCGs, MRGs, and mtDNA* (coding cox3 to nad5) in E. malayanum and Echinostomatidae. The nucleotide use is characterized by (T > G > A > C) for PCGs/mtDNA*, and by (T > G ≈ A > C) for MRGs. E. malayanum exhibited the lowest genetic distance (0.53%) to Artyfechinostomum sufrartyfex, relatively high to the Echinostoma congeners (13.20-13.99%), higher to Hypoderaeum conoideum (16.18%), and the highest to interfamilial Echinochasmidae (26.62%); Cyclocoelidae (30.24%); and Himasthlidae (25.36%). Topology indicated the monophyletic position between E. malayanum/A. sufrartyfex and the group of Echinostoma caproni, Echinostoma paraensei, Echinostoma miyagawai, and Echinostoma revolutum, rendering Hypoderaeum conoideum and unidentified Echinostoma species paraphyletic. The strictly closed genomic/taxonomic/phylogenetic features (including base composition, skewness, codon usage/bias, genetic distance, and topo-position) reinforced Echinostoma malayanum to retake its generic validity within the Artyfechinostomum genus.

Morphological and molecular identification of Cryptocotyle lingua metacercariae isolated from Atlantic cod (Gadus morhua) from Danish seas and whiting (Merlangius merlangus) from the English Channel

Trematode larvae (metacercariae) causing black spot disease occur frequently in gills, fins, skin and the superficial muscle layers of marine fish. Species within the genus Cryptocotyle Lühe, 1899 are frequently associated with this disease. Descriptions of the metacercarial stage are relatively limited and none has hitherto been reported from fish from the English Channel. The present study reports the morphological and molecular identifications of encysted black spot-inducing parasites from whiting (Merlangius merlangus) and Atlantic cod (Gadus morhua) caught respectively from the north coast of France (English Channel) and from Danish sea waters. Metacercariae were characterised morphologically based on microscopic observations and molecularly using Sanger sequencing of fragments of the mitochondrial cox1 gene and rDNA ITS region. Morphological data were compared with available data in the literature. Phylogenetic trees including reference sequences were built to confirm morphological and molecular identifications. This survey constitutes the first description of C. lingua metacercariae in the English Channel ecosystems.

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.

Gene Duplication and Gain in the Trematode Atriophallophorus winterbourni Contributes to Adaptation to Parasitism

Gene duplications and novel genes have been shown to play a major role in helminth adaptation to a parasitic lifestyle because they provide the novelty necessary for adaptation to a changing environment, such as living in multiple hosts. Here we present the de novo sequenced and annotated genome of the parasitic trematode Atriophallophorus winterbourni and its comparative genomic analysis to other major parasitic trematodes. First, we reconstructed the species phylogeny, and dated the split of A. winterbourni from the Opisthorchiata suborder to approximately 237.4 Ma (±120.4 Myr). We then addressed the question of which expanded gene families and gained genes are potentially involved in adaptation to parasitism. To do this, we used hierarchical orthologous groups to reconstruct three ancestral genomes on the phylogeny leading to A. winterbourni and performed a GO (Gene Ontology) enrichment analysis of the gene composition of each ancestral genome, allowing us to characterize the subsequent genomic changes. Out of the 11,499 genes in the A. winterbourni genome, as much as 24% have arisen through duplication events since the speciation of A. winterbourni from the Opisthorchiata, and as much as 31.9% appear to be novel, that is, newly acquired. We found 13 gene families in A. winterbourni to have had more than ten genes arising through these recent duplications; all of which have functions potentially relating to host behavioral manipulation, host tissue penetration, and hiding from host immunity through antigen presentation. We identified several families with genes evolving under positive selection. Our results provide a valuable resource for future studies on the genomic basis of adaptation to parasitism and point to specific candidate genes putatively involved in antagonistic host-parasite adaptation.

Discovery of Carcinogenic Liver Fluke Metacercariae in Second Intermediate Hosts and Surveillance on Fish-Borne Trematode Metacercariae Infections in Mekong Region of Myanmar

Countries of lower Mekong regions are highly alarmed by the spread of fish-borne trematode infections, i.e., small liver flukes and minute intestinal flukes especially in Thailand, Lao People’s Democratic Republic (Lao PDR), Vietnam, Cambodia and Myanmar. Moreover; the incidence of cholangiocarcinoma has also been increasing in the endemic area of liver fluke infections. Only a few reports have been published concerning the fish-borne trematodes infections in the central region of Myanmar. However; there is still a lack of information regarding the status of trematodes infections in second intermediate hosts in the Mekong region of Myanmar. Therefore, we conducted surveillance on the distribution of trematode metacercariae in small cyprinoid fishes collected from the Mekong region of Myanmar. A total of 689 fishes (12 different species of cyprinoid fishes) have been collected and examined by pepsin digestion methods. We discovered four species of fish-borne trematode metacercariae infections, i.e., carcinogenic liver fluke, Opisthorchis viverrini; minute intestinal flukes, Haplorchis taichui; Haplorchis pumilio and Haplorchoides sp. in Tachileik, the Mekong Region of Myanmar. The outcome of this study could be a useful index for the fish-borne zoonotic trematode epidemiology in the Mekong area. Besides, the results of our study contribute to filling the gap of information necessary for the control and prevention of fish-borne trematode zoonotic infections in the Mekong region.

The Complete Mitochondrial Genome of the Caecal Fluke of Poultry, Postharmostomum commutatum, as the First Representative from the Superfamily Brachylaimoidea

Postharmostomum commutatum (Platyhelminthes: Brachylaimoidea), a parasite of the caeca of poultry, has been frequently reported from many countries and regions, including China. However, the molecular epidemiology, population genetics and phylogenetics of this parasite are poorly understood. In the present study, we determined and characterized the complete mitochondrial (mt) genome of P. commutatum, as the first representative from the superfamily Brachylaimoidea. The mt genome of P. commutatum is a circular DNA molecule of 13,799 bp in size and encodes the complete set of 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as a typical control region. The mt genome of P. commutatum presents a clear bias in nucleotide composition with a negative AT-skew on average (-0.306) and a positive GC-skew on average (0.466). Phylogenetic analyses showed that P. commutatum (superfamily Brachylaimoidea) and other ten members of the order Diplostomida were recovered as sister groups of the order Plagiorchiida, indicating that the order Diplostomida is paraphyletic. This is the first mt genome of any member of the superfamily Brachylaimoidea and should represent a rich source of genetic markers for molecular epidemiological, population genetic and phylogenetic studies of parasitic flukes of socio-economic importance in poultry.

Characterization of the mitochondrial genome sequences of the liver fluke Amphimerus sp. (Trematoda: Opisthorchiidae) from Ecuador and phylogenetic implications

Amphimerus Barker, 1911 is a liver fluke infecting several animal species and humans. Being a digenetic trematode of the Opisthorchiidae family, Amphimerus is closely related to the genera Metorchis, Clonorchis and Opisthorchis. Recently, a high prevalence of Amphimerus infection in humans, cats, and dogs had been demonstrated in a tropical Pacific region of Ecuador. Hence, we determined and characterized the entire mt genome sequences of adult liver flukes, morphologically identified as Amphimerus, collected in the endemic region of Ecuador, and examined its phylogenetic relationships with flukes in the Opisthorchiidae family using Bayesian inference (BI) based on the concatenated amino acid sequences and partial cox1 sequences. The complete mt genome sequence (15, 151 bp in length) of the Amphimerus sp. contains 35 genes, including 12 protein-coding genes (PCGs, without atp8), two rRNAs (rrnL and rrnS) and 21 tRNAs, lacking trnG. The gene content and arrangement of the Ecuadorian Amphimerus mt genome was similar to those of other trematodes in the Opisthorchiidae family. All genes in the circular mt genome of Amphimerus sp. are transcribed from the same strand in one direction, with the A + T content of 60.77%. Genetic distances between Amphimerus sp. and other genera in Opisthorchiidae were rather high, ranging from 26.86% to 28.75% at nucleotide level and 29.37%-31.12% at amino acid level. Phylogenetic analysis placed the Ecuadorian Amphimerus within the branch of Opisthorchiidae, but very distinct from Opisthorchis. Our results indicate that the liver fluke Amphimerus from Ecuador does not belong to the genus Opisthorchis, and that it should be assigned under the genus Amphimerus. The determination of the mt genome of the Ecuadorian Amphimerus provides a new genetic resource for future studies on taxonomy and molecular epidemiology of Opisthorchiidae trematodes.

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.

Updated molecular phylogenetic data for Opisthorchis spp. (Trematoda: Opisthorchioidea) from ducks in Vietnam

BACKGROUND

An opisthorchiid liver fluke was recently reported from ducks (Anas platyrhynchos) in Binh Dinh Province of Central Vietnam, and referred to as "Opisthorchis viverrini-like". This species uses common cyprinoid fishes as second intermediate hosts as does Opisthorchis viverrini, with which it is sympatric in this province. In this study, we refer to the liver fluke from ducks as "Opisthorchis sp. BD2013", and provide new sequence data from the mitochondrial (mt) genome and the nuclear ribosomal transcription unit. A phylogenetic analysis was conducted to clarify the basal taxonomic position of this species from ducks within the genus Opisthorchis (Digenea: Opisthorchiidae).

METHODS

Adults and eggs of liver flukes were collected from ducks, metacercariae from fishes (Puntius brevis, Rasbora aurotaenia, Esomus metallicus) and cercariae from snails (Bithynia funiculata) in different localities in Binh Dinh Province. From four developmental life stage samples (adults, eggs, metacercariae and cercariae), the complete cytochrome b (cob), nicotinamide dehydrogenase subunit 1 (nad1) and cytochrome c oxidase subunit 1 (cox1) genes, and near-complete 18S and partial 28S ribosomal DNA (rDNA) sequences were obtained by PCR-coupled sequencing. The alignments of nucleotide sequences of concatenated cob + nad1 + cox1, and of concatenated 18S + 28S were separately subjected to phylogenetic analyses. Homologous sequences from other trematode species were included in each alignment.

RESULTS

Phylogenetic trees were inferred from concatenated (cob + nad1 + cox1) nucleotide sequences and combined 18S + 28S nucleotide sequences of five Opisthorchis sp. BD2013 samples and additional reference taxa. Both trees demonstrated the anticipated clustering of taxa within the superfamily Opisthorchioidea, the paraphyly of the genus Opisthorchis and the sister-species relationship of Opisthorchis sp. BD2013 with O. viverrini.

CONCLUSIONS

While it is likely that Opisthorchis sp. BD2013 is distinct from O. viverrini, it is clearly a sister taxon of O. viverrini within the limited number of Opisthorchis species for which appropriate sequence data are available. The new sequences provided here will assist the diagnosis and the taxonomic clarification of the opisthorchiid species.

Fishborne zoonotic heterophyid infections: An update

Fishborne heterophyid trematodes infecting humans are at least 29 species worldwide and belong to 13 genera. Its global burden is much more than 7 million infected people. They include Metagonimus (M. yokogawai, M. takahashii, M. miyatai, M. minutus, and M. katsuradai), Heterophyes (H. heterophyes, H. nocens, H. dispar, and H. aequalis), Haplorchis (H. taichui, H. pumilio, H. yokogawai, and H. vanissimus), Pygidiopsis (P. summa and P. genata), Heterophyopsis (H. continua), Stellantchasmus (S. falcatus), Centrocestus (C. formosanus, C. armatus, C. cuspidatus, and C. kurokawai), Stictodora (S. fuscata and S. lari), Procerovum (P. varium and P. calderoni), Acanthotrema (A. felis), Apophallus (A. donicus), Ascocotyle (A. longa), and Cryptocotyle (C. lingua). Human infections are scattered around the world but the major endemic areas are located in Southeast Asia. The source of human infection is ingestion of raw or improperly cooked fish. The pathogenicity, host-parasite relationships, and clinical manifestations in each species infection are poorly understood; these should be elucidated particularly in immunocompromised hosts. Problems exist in the differential diagnosis of these parasitic infections because of close morphological similarity of eggs in feces and unavailability of alternative methods such as serology. Molecular diagnostic techniques are promising but they are still at an infant stage. Praziquantel has been proved to be highly effective against most of the patients infected with heterophyid flukes. Epidemiological surveys and detection of human infections are required for better understanding of the geographical distribution and global burden of each heterophyid species. In this review, the most updated knowledge on the morphology, biology, epidemiology, pathogenesis and pathology, immunology, clinical manifestations, diagnosis and treatment, and prevention and control of fishborne zoonotic heterophyid infections is provided.

Molecular characterization of the unique Mesostephanus appendiculatus (Trematoda: Cyathocotylidae) by small ribosomal RNA from Egypt

M esostephanus appendiculatus (Family: Cyathocotylidae) is one of the unique trematodes that complete their cycles in human and animal intestines in many countries of the world. The main source of its transmission is eating raw or undercooked infected fish muscle. Earliest analyses of genes to different parasites supported the analysis of helminthes either biological or morphological. This paper detected M. appendiculatus sequence with GenBank accession number gb (KY026782). Comparison of M. appendiculatus with other helminthes using BioEdit 7 and MEGA7 program shows some similarity in different points along its sequence. The phylogenetic analysis clarifies that it was closely related to both trematodes (Clinostomum complanatum and Echinochasmus japonicus) and some cestodes of fish origin such as Polyonchobothrium polypteri, Bothriocephalus sp., and Haplobothrium globuliforme. The obtained results provide a good source for genome analysis of M. appendiculatus in relation to other Platyhelminthes.

Molecular characterization and confocal laser scanning microscopic study of Pygidiopsis macrostomum (Trematoda: Heterophyidae) parasites of guppies Poecilia vivipara

Pygidiopsis macrostomum and Ascocotyle (Phagicola) pindoramensis (Digenea: Heterophyidae) parasitize guppies as intermediate hosts and, respectively, fish-eating mammals or birds as definitive hosts. Heterophyids have zoonotic potential, and molecular studies associated with morphological and ecological aspects have helped to clarify their taxonomy and phylogeny. Poecilia vivipara naturally parasitized by metacercariae of both species (100% prevalence) exhibit no external signs of parasitism. In this work, four new sequences of P. macrostomum (18S rDNA, 28S rDNA and ITS2 rDNA) and one new sequence of A. (P.) pindoramensis (mtDNA cox-1) are presented. Phylogeny reconstructions linked P. macrostomum to other heterophyids, but the separation of the Heterophyidae and Opisthorchiidae remains unclear. Additionally, we used indirect immunocytochemistry and the phalloidin-fluorescence techniques allied with confocal laser scanning microscopy to describe muscular and neuronal structures of P. macrostomum. A complex arrangement of muscular fibres is associated with the tegument, suckers, gut and reproductive system. Radial fibres around the ventral sucker are thick, branched and extend to the body wall. High-resolution confocal imaging revealed a typical digenean muscular arrangement and important heterophyid morphological traits. These data will support future control measures to reduce the parasitism in guppies reared in fish farming systems, especially for aquarium and experimental purposes.

Characterization of the complete mitochondrial genome sequence of Homalogaster paloniae (Gastrodiscidae, Trematoda) and comparative analyses with selected digeneans

Gastrodiscidae species are neglected but significant paramphistomes in small ruminants, which can lead to considerable economic losses to the breeding industry of livestock. However, knowledge about molecular ecology, population genetics, and phylogenetic analysis is still limited. In the present study, we firstly sequenced and analyzed the full mitochondrial (mt) genome of Homalogaster paloniae (14,490 bp). The gene contents and organization of the H. paloniae mt genome is the same as that of other digeneans, such as Fasciola hepatica and Paramphistomum cervi. It is interesting that unlike other paramphistomes, H. paloniae is flat in shape which is similar with Fasciola, such as F. hepatica. Phylogenetic analysis of H. paloniae and other 17 selected digeneans using concatenated amino acid sequences of the 12 protein-coding genes showed that Gastrodiscidae is closely related to Paramphistomidae and Gastrothylacidae. The availability of the mt genome sequence of H. paloniae should provide an important foundation for further molecular study of Gastrodiscidae and other digeneans.

Epidemiological and molecular data on heterophyid trematode metacercariae found in the muscle of grey mullets (Osteichthyes: Mugilidae) from Sardinia (western Mediterranean Sea)

This study is a contribution to the molecular taxonomy and epidemiology of heterophyid (Trematoda: Heterophyidae) metacercariae found in the muscle of Mugilidae (Osteichthyes) from Sardinia (western Mediterranean Sea). Sixty specimens of mugilids (13 Chelon labrosus, 18 Liza aurata, 6 Liza ramada, 8 Liza saliens, 15 Mugil cephalus) were examined and 17,899 metacercariae isolated in 95 % of the hosts. Four types of metacercariae were identified: Heterophyes sp. (n = 14,113), Heterophyes sp. -small (1225), Stictodora sp. (1606), and Ascocotyle (Phagicola) sp. (955). The experimental infection of a hamster with Heterophyes sp. metacercariae produced six adults identified as Heterophyes heterophyes and two as Heterophyes cf. nocens. The morphology of Heterophyes sp. -small metacercariae matched with that of Heterophyes dispar. The sequence analysis of the ITS2 and 28S portions of rDNA confirmed the morphological identification of metacercariae, showing four clusters. All adults grouped together with the Heterophyes sp. metacercariae, whereas adults of Heterophyes nocens from Korea clustered separately, showing that this species is distinguished from H. heterophyes, and suggesting caution in the exclusive use of the number of rodlets of the genital sucker to separate the two species. The presence of metacercariae was high in all hosts; the highest prevalence is of Heterophyes sp. (prevalence ≥78 %; mean intensity ≥135 metacercariae/100 g muscle), and the most heavily infected host is M. cephalus (prevalence = 100 %; mean intensity = 841 metacercariae/100 g muscle).

Complete Mitochondrial Genome of Echinostoma hortense (Digenea: Echinostomatidae)

Echinostoma hortense (Digenea: Echinostomatidae) is one of the intestinal flukes with medical importance in humans. However, the mitochondrial (mt) genome of this fluke has not been known yet. The present study has determined the complete mt genome sequences of E. hortense and assessed the phylogenetic relationships with other digenean species for which the complete mt genome sequences are available in GenBank using concatenated amino acid sequences inferred from 12 protein-coding genes. The mt genome of E. hortense contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 non-coding region. The length of the mt genome of E. hortense was 14,994 bp, which was somewhat smaller than those of other trematode species. Phylogenetic analyses based on concatenated nucleotide sequence datasets for all 12 protein-coding genes using maximum parsimony (MP) method showed that E. hortense and Hypoderaeum conoideum gathered together, and they were closer to each other than to Fasciolidae and other echinostomatid trematodes. The availability of the complete mt genome sequences of E. hortense provides important genetic markers for diagnostics, population genetics, and evolutionary studies of digeneans.

The complete mitochondrial genome of Orthocoelium streptocoelium (Digenea: Paramphistomidae) for comparison with other digeneans

Orthocoelium streptocoelium is a common paramphistome species parasitizing the rumen and/or reticulum of small ruminants, leading to significant losses. This study first determined the complete mitochondrial (mt) genome of O. streptocoelium. The complete mt genome of O. streptocoelium was amplified, sequenced, assembled, analysed and then compared with those of other digeneans. The entire mt genome of O. streptocoelium is 13,800 bp in length, which is smaller than those of other digeneans except for Opisthorchis viverrini. This mt genome contains 12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and two non-coding regions. The arrangement of the O. streptocoelium mt genome is the same as those of other digeneans except for Schistosoma haematobium and Schistosoma spindale. Phylogenetic analyses based on concatenated amino acid sequences of the 12 protein-coding genes representing 16 digeneans were conducted to assess the relationship of O. streptocoelium with other digeneans. The result indicated that O. streptocoelium is closely related to Paramphistomum cervi and Fischoederius elongates, which is in accordance with their relationships by taxonomy. This complete mt genome of O. streptocoelium enriched the mitochondrial genome data of paramphistomes and provided important molecular markers for diagnostics and studies of population variation, epidemiology, ecology and evolution of O. streptocoelium and other digeneans.

The complete mitochondrial genome of Metorchis orientalis (Trematoda: Opisthorchiidae): Comparison with other closely related species and phylogenetic implications

Metorchis orientalis (Trematoda: Opisthorchiidae) is an important trematode infecting many animals and humans, causing metorchiasis. In the present study, the complete mitochondrial (mt) genome of M. orientalis was sequenced. The complete mt genome of M. orientalis is 13,834 bp circular DNA molecule and contains 12 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. The gene content and arrangement of M. orientalis is the same as those of Opisthorchiidae trematodes (Opisthorchis viverrini, Opisthorchis felineus and Clonorchis sinensis), but distinct from Schistosoma spp. Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes with three different computational algorithms (Bayesian inference, maximum likelihood and maximum parsimony) revealed that M. orientalis and O. viverrini represent sister taxa. The mt genome provides a novel genetic marker for further studies of the identification, classification and molecular epidemiology of Opisthorchiidae trematodes, and should have implications for the diagnosis, prevention and control of metorchiasis in animals and humans.

Mitochondrial genome of Ogmocotyle sikae and implications for phylogenetic studies of the Notocotylidae trematodes

Ogmocotyle spp. (Trematoda: Digenea: Notocotylidae) are neglected but important trematodes that can infect numerous mammal species, causing significant economic losses to livestock industries. However, there have been few studies on the molecular ecology of these trematodes. We amplified and sequenced the complete mitochondrial (mt) genome of Ogmocotyle sikae (14,307 bp). The gene content and gene arrangement of O. sikae mt genome was similar to that of Dicrocoelium chinensis, except that trnE and trnG were reversed. Phylogenetic analysis of O. sikae and selected parasites using Bayesian inference was performed based on concatenated amino acid sequence datasets conceptually translated from the 12 protein-coding genes. The results indicated that the family Notocotylidae is related to the family Paramphistomatidae. Our description of O. sikae mt genome provides a significant resource of molecular markers for future comparative studies of the Notocotylidae and other trematodes.

Analysis of the complete Fischoederius elongatus (Paramphistomidae, Trematoda) mitochondrial genome

Background

Fischoederius elongates is an important trematode of Paramphistomes in ruminants. Animals infected with F. elongates often don’t show obvious symptoms, so it is easy to be ignored. However it can cause severe economic losses to the breeding industry. Knowledge of the mitochondrial genome of F. elongates can be used for phylogenetic and epidemiological studies.

Findings

The complete mt genome sequence of F. elongates is 14,120 bp in length and contains 12 protein-coding genes, 22 tRNA genes, two rRNA genes and two non-coding regions (LNR and SNR). The gene arrangement of F. elongates is the same as other trematodes, such as Fasciola hepatica and Paramphistomum cervi. Phylogenetic analyses using concatenated amino acid sequences of the 12 protein-coding genes by Maximum-likelihood and Neighbor-joining analysis method showed that F. elongates was closely related to P. cervi.

Conclusion

The complete mt genome sequence of F. elongates should provide information for phylogenetic and epidemiological studies for F. elongates and the family Paramphistomidae.

Mitochondrial and nuclear ribosomal DNA dataset supports that Paramphistomum leydeni (Trematoda: Digenea) is a distinct rumen fluke species

BACKGROUND

Rumen flukes parasitize the rumen and reticulum of ruminants, causing paramphistomiasis. Over the years, there has been considerable debate as to whether Paramphistomum leydeni and Paramphistomum cervi are the same or distant species.

METHODS

In the present study, the complete mitochondrial (mt) genome of P. leydeni was amplified using PCR-based sequencing and compared with that of P. cervi. The second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA (rDNA) of P. leydeni specimens (n = 6) and P. cervi specimens (n = 8) was amplified and then sequenced. Phylogenetic relationship of the concatenated amino acid sequence data for 12 protein-coding genes of the two rumen flukes and selected members of Trematoda was evaluated using Bayesian inference (BI).

RESULTS

The complete mt genome of P. leydeni was 14,050 bp in size. Significant nucleotide difference between the P. leydeni mt genome and that of P. cervi (14.7%) was observed. For genetic divergence in ITS-2, sequence difference between P. leydeni and P. cervi was 3.1%, while no sequence variation was detected within each of them. Phylogenetic analysis indicated that P. leydeni and P. cervi are closely-related but distinct rumen flukes.

CONCLUSIONS

Results of the present study support the proposal that P. leydeni and P. cervi represent two distinct valid species. The mt genome sequences of P. leydeni provide plentiful resources of mitochondrial markers, which can be combined with nuclear markers, for further comparative studies of the biology of P. leydeni and its congeners from China and other countries.