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

Morpho-Molecular Features and Phylogenetic Relationships of Metorchis butoridi Oschmarin, 1963 (Trematoda: Opisthorchiidae) from East Asia

Adult trematodes of the genus Metorchis were found in the gallbladders of ducklings that had been experimentally fed freshwater fishes of the genera Rhynchocypris and Rhodeus that were naturally infected by Metorchis metacercariae. Some of the trematodes were identified as Metorchis ussuriensis, whose morphology of developmental stages and molecular data had previously been described in detail. The other trematodes were confirmed as species Metorchis butoridi on the basis of morphological features: subterminal oral sucker, vitelline follicles with interrupted bands, and rosette-shaped testes. An analysis of phylogenetic relationships within Opisthorchiidae using nuclear and mitochondrial markers confirmed that the obtained trematodes were actually from the genus Metorchis. The morphological and molecular features indicated that a number of trematodes found in East Asia and described as Metorchis orientalis belong to M. butoridi. Also, the M. orientalis individuals from Europe are, in fact, representatives of another Metorchis species.

The complete mitochondrial genome of Morishitium polonicum (Trematoda, Cyclocoelidae) and its phylogenetic implications

Trematodes can adversely impact the health and survival of wild animals. The trematode family Cyclocoelidae, which includes large digenean bird parasites, lacks molecular analysis, and reclassifications have not been supported. This study produced the first fully assembled and annotated mitochondrial genome sequence for the trematode Morishitium polonicum. The whole length of the M. polonicum (GenBank accession number: OP930879) mitogenome is 14083 bp, containing 22 transfer ribonucleic acids (tRNAs), 2 ribosomal RNAs (rRNAs, rrnL and rrnS), and a noncoding control section (D-loop) 13777 to 13854 bp in length. The 12 PCG areas have 3269 codons and a total length of 10053 bp, which makes up 71.38% of the mitochondrial genome's overall sequence. Most (10/12) of the PCGs that code for proteins begin with ATG, while the nad4L and nad1 genes have a GTG start codon. Phylogenetic analysis using the concatenated nucleotide sequences of 12 PCGs, and the ML tree analysis results showed that M. polonicum is more closely related to with Echinostomatidae and Fasciolidae, which indicates that the family Cyclocoelidae is more closely associated with Echinochasmidae. This study provides mtDNA information, and analysis of mitogenomic structure and evolution. Moreover, we aimed to understand the phylogenetic relationships of this fluke.

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.

The complete mitochondrial genome of Prosthogonimus cuneatus and Prosthogonimus pellucidus (Trematoda: Prosthogonimidae), their features and phylogenetic relationships in the superfamily Microphalloidea

Prosthogonimus cuneatus and Prosthogonimus pellucidus (Trematoda: Prosthogonimidae) are common flukes of poultry and other birds which can cause severe impacts on animal health and losses to the poultry industry. However, there are limited studies on the molecular epidemiology, population genetics, and systematics of Prosthogonimus species. In the present study, the complete mitochondrial (mt) genomes of P. cuneatus and P. pellucidus were determined to be 14,829 bp and 15,013 bp in length, respectively. Both mt genomes contain 12 protein-coding genes (PCGs) (cox1-3, nad1-6, nad4L, cytb, and atp6), 22 transfer RNA genes, two ribosomal RNA genes, and one non-coding region. Our comparative analysis shows that the atp6 genes of P. cuneatus and P. pellucidus are longer than any previously published atp6 genes of other trematodes. The lengths of the atp6 genes of P. cuneatus and P. pellucidus in this study seem unusual, and should therefore be studied further. The mt genes of P. cuneatus and P. pellucidus are transcribed in the same direction, and the gene arrangements are identical to those of Plagiorchis maculosus, Tamerlania zarudnyi, and Tanaisia sp., but different from those of Eurytrema pancreaticum, Dicrocoelium chinensis, and Brachycladium goliath. The mt genome A + T contents of P. cuneatus and P. pellucidus are 64.47% and 65.34%, respectively. In the 12 PCGs, ATG is the most common initiation codon, whereas TAG is the most common termination codon. The sequence identity of the same 12 PCGs among the eight trematodes (P. cuneatus, P. pellucidus, Pl. maculosus, D. chinensis, E. pancreaticum, B. goliath, T. zarudnyi, Tanaisia sp.) of Xiphidiata are 55.5%-81.7% at the nucleotide level and 43.9%-82.5% at the amino acid level. The nucleotide similarities among the complete mt genomes of the eight trematodes range from 54.1%-81.5%. Phylogenetic analysis based on the aligned concatenated amino acid sequences of the 12 PCGs shows that P. cuneatus and P. pellucidus cluster together and are sister to T. zarudnyi and Tanaisia sp., and this clade is more closely related to E. pancreaticum, Dicrocoelium spp. and Lyperosomum longicauda in the family Dicrocoeliidae, than it is to species in the families Plagiorchiidae and Brachycladiidae. These are the first reported complete mt genomes of Prosthogonimidae, and these data will provide additional molecular resources for further studies of Prosthogonimidae taxonomy, population genetics, and systematics.

Prevalence of fish-borne zoonotic trematode infection in Jilin Province, China

Fish-borne zoonotic trematodes (FZTs) are the most serious food-borne parasites in Asia and have become a burden to public health and a new challenge in food safety. In Jilin Province, China, the prevalence of FZTs in intermediate and definitive hosts has not been extensively explored. In the present study, we investigated the prevalence of FZTs in Jilin Province, China. From July to November 2020, a total of 132 freshwater snails (the first intermediate host of FZTs), 4122 wild freshwater fishes (the second intermediate host of FZTs) and 143 fecal samples from canines, ducks and swine (the definitive host of FZTs) were collected from the Yitong River basin of Jilin Province. FZT species were identified by morphological observation combined with internal transcribed spacer (ITS) sequence analysis. The prevalence of FZTs was then calculated. The results showed that the prevailing species of FZTs in Jilin Province, China, were Clonorchis sinensis, Metorchis orientalis and Echinochasmus japonicus. The total prevalence of FZTs was 29.74% (1226/4122) in fish, the total infection rates were 2.27% (3/132) in snails, 75.00% (21/28) in canines and 37.18% (29/78) in ducks. The coinfection rates of the two trematodes were 13.39% (552/4122) in fish, 35.71% (10/28) in canines and 7.69% (6/78) in ducks. The coinfection rate of the three flukes was 2.60% (107/4122) in fish. Nine of the 12 fish species examined were infected with FZT metacercariae.

• C.

sinensis, M. orientalis and E. japonicus were reported in Jilin province of China.

•

The prevalence of FZTs in freshwater fish were 29.74% (1226/4122).

•

Four fish species were identified as new second intermediate hosts of E. japonicus.

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.

Comparative Analysis of Mitochondrial Genomes among Twelve Sibling Species of the Genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae) and Phylogenetic Analysis

Simple Summary

Atkinsoniella is a large genus of 98 species across the world and 88 species recorded in China within the globally distributed subfamily Cicadellinae, which is phytophagous, and some of which have been reported as important agricultural pests. Some Atkinsoniella species are very similar in morphological characteristics, making accurate identification at species level confusing. To provide further evidence toward understanding the relationships within the genus Atkinsoniella and subfamily Cicadellinae, mitogenomes of 12 Atkinsoniella sibling species were obtained and annotated. Their characteristics were comparatively analyzed. In addition, the comprehensive phylogenetic relationship within the subfamily Cicadellinae was determined based on three mitochondrial datasets using both the maximum-likelihood (ML) and Bayesian inference (BI) methods. The results suggested that the genus Atkinsoniella was recovered as a monophyletic group. The branches of the 12 newly sequenced species were clearly separated, with most nodes receiving strong support in all analyses, indicating that mitogenomics is an effective method for identifying closely related species and understanding their phylogenetic and evolutionary relationships.

Abstract

The herbivorous leafhopper genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae), a large genus of subfamily Cicadellinae, consists of 98 valid species worldwide and 88 species recorded in China. Some species of the genus are very similar in morphological characteristics, so they are difficult to identify accurately. In this study, 12 mitochondrial genomes of Atkinsoniella species with similar morphological characteristics were first obtained through high-throughput sequencing, which featured a typical circular molecule of 15,034–15,988 bp in length. The arrangement and orientation of 37 genes were identical to those of typical Cicadellidae mitogenomes. The phylogenetic relationship within the subfamily Cicadellinae was reconstructed using maximum-likelihood (ML) and Bayesian inference (BI) methods based on three concatenated datasets. The topological structures of the six obtained phylogenetic trees were highly consistent. The results suggested that Atkinsoniella was recovered as a monophyletic group and emerged as a sister group with the monophyletic clade of Bothrogonia, Paracrocampsa (part), and Draeculacephala (part). The branches of the 12 newly sequenced species were clearly separated, with most nodes receiving strong support in all analyses. In addition, the key to the 12 Atkinsoniella species was provided to identify species according to morphological characteristics. This study further promotes research on the classification, genetics, evolution, and phylogeny of the genus Atkinsoniella and subfamily Cicadellinae.

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.

The potential use of mitochondrial ribosomal genes (12S and 16S) in DNA barcoding and phylogenetic analysis of trematodes

Background

Genetic markers like the nuclear ribosomal RNA (rRNA) genes, internal transcribed spacer regions, mitochondrial protein-coding genes, and genomes have been utilized for molecular identification of parasitic trematodes. However, challenges such as the design of broadly applicable primers for the vast number of species within Digenea and the genetic markers’ ability to provide sufficient species-level resolution limited their utility. This study presented novel and broadly applicable primers using the mitochondrial 12S and 16S rRNA genes for Digenea and aimed to show their suitability as alternative genetic markers for molecular identification of orders Plagiorchiida, Echinostomida, and Strigeida.

Results

Our results revealed that the mitochondrial 12S and 16S rRNA genes are suitable for trematode molecular identification, with sufficient resolution to discriminate closely related species and achieve accurate species identification through phylogenetic placements. Moreover, the robustness of our newly designed primers to amplify medically important parasitic trematodes encompassing three orders was demonstrated through successful amplification. The convenience and applicability of the newly designed primers and adequate genetic variation of the mitochondrial rRNA genes can be useful as complementary markers for trematode molecular-based studies.

Conclusions

We demonstrated that the mitochondrial rRNA genes could be alternative genetic markers robust for trematode molecular identification and potentially helpful for DNA barcoding where our primers can be widely applied across the major Digenea orders. Furthermore, the potential of the mitochondrial rRNA genes for molecular systematics can be explored, enhancing their appeal for trematode molecular-based studies. The novelty of utilizing the mitochondrial rRNA genes and the designed primers in this study can potentially open avenues for species identification, discovery, and systematics in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08302-4.

An annotated checklist of the eukaryotic parasites of humans, exclusive of fungi and algae

The classification of "parasites" in the medical field is a challenging notion, a group which historically has included all eukaryotes exclusive of fungi that invade and derive resources from the human host. Since antiquity, humans have been identifying and documenting parasitic infections, and this collective catalog of parasitic agents has expanded considerably with technology. As our understanding of species boundaries and the use of molecular tools has evolved, so has our concept of the taxonomy of human parasites. Consequently, new species have been recognized while others have been relegated to synonyms. On the other hand, the decline of expertise in classical parasitology and limited curricula have led to a loss of awareness of many rarely encountered species. Here, we provide a comprehensive checklist of all reported eukaryotic organisms (excluding fungi and allied taxa) parasitizing humans resulting in 274 genus-group taxa and 848 species-group taxa. For each species, or genus where indicated, a concise summary of geographic distribution, natural hosts, route of transmission and site within human host, and vectored pathogens are presented. Ubiquitous, human-adapted species as well as very rare, incidental zoonotic organisms are discussed in this annotated checklist. We also provide a list of 79 excluded genera and species that have been previously reported as human parasites but are not believed to be true human parasites or represent misidentifications or taxonomic changes.

Cytochrome b as a more promising marker for analysing the distribution vector for Metagonimus suifunensis (Trematoda: Heterophyidae)

In this study of Metagonimus suifunensis (M. suifunensis) in the Russian Southern Far East, the variability of the full-length sequences of the cytochrome b (cytb) mtDNA gene was assessed for the first time. In addition, the cox1 mtDNA gene sequences were also obtained for this species from new localities. In total, 87 and 81 sequences of the cytb and cox1 genes, respectively, were used in the current study. The cytb gene proved more promising and revealed two haplogroups that are associated with the spatial distribution of the species: geographical isolation caused the fixation of differences between northern and southern populations. In addition, the results obtained for the cytb gene opened up new perspectives in the analysis of sequences of the cox1 gene, which was not sufficiently effective as a sole marker. Based on data for both mitochondrial genes, molecular processes influencing the formation of the modern population were analysed for M. suifunensis. The new data confirmed the previously expressed opinion that this species colonized the study territory from north to south and will form the basis for determining possible ways of its further expansion, which is important for predicting the emergence of new foci of metagonimosis.

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 Notocotylus sp. (Trematoda, Notocotylidae) and its phylogenetic implications

The parasite genus Notocotylus comprises at least 50 species colonizing mainly aquatic birds and to a lesser extent some mammals, particularly rodents. Here trematode specimens isolated from a wild black swan were characterized and identified to belong to the genus Notocotylus via morphological and molecular analyses. Phylogenetic position of the isolate among other trematodes was determined based on the ribosomal internal transcribed spacer (ITS) 1 and 2. The complete mitochondrial (mt) genome of the isolate was amplified, sequenced, assembled, analyzed, and annotated. The isolate has an AT-rich mt genome (14,317 bp in length) that comprises 12 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes. The Notocotylus isolate identified in this study has relatively high mt genome sequence identity and identical gene content and arrangement to a known Notocotylidae species, Ogmocotyle sikae. The isolate formed a genetic clade with O. sikae in phylogenetic analysis of the concatenated PCG amino acid sequences. Compared to the ITS, the trematode mt genome appears more informative for resolving high-level phylogenies. To the best of our knowledge, this is the first study exploring the complete mt genome for the genus Notocotylus, and it offers a novel genomic resource that has important implications for trematode phylogenetics.

Trematode diversity reflecting the community structure of Danish freshwater systems: molecular clues

Background

Digenean trematodes are parasitic platyhelminths that use several hosts in their life cycles and are thereby embedded in various ecosystems affected by local environmental conditions. Their presence in a habitat will reflect the presence of different host species and, as such, they can serve as ecological indicators. Only limited information on the occurrence of trematodes and their link to other trophic levels in the Danish freshwater ecosystems is currently available.Therefore, the main aim of the present study was to increase our knowledge in this field.

Methods

Snails were sampled from 21 freshwater lakes in Denmark, following which shedding procedures were performed, cercariae were recoved and the released parasites were identified using molecular tools (PCR and sequencing).

Results

A total of 5657 snail hosts belonging to ten species were identified, revealing a highly diverse parasite fauna comprising 22 trematode species. The overall trematode prevalence was 12.6%, but large variations occurred between host species. The snail host Lymnaea stagnalis showed the highest prevalence and also exhibited the highest diversity, accounting for 47.6% of the species richness.

Conclusions

This survey contributes updated information on parasite–host relations and compatibility and may assist in describing the ecological structure of the investigated Danish freshwater ecosystems.

First case report of Metorchis orientalis from Black Swan

Metorchis orientalis belongs to the genus Metorchis of Opisthorchiidae, which mainly parasitizes in liver and bile ducts of waterfowl, causing liver dysfunction of the host. It has been reported that M. orientalis also infects humans. As a natural species in Australia and a popular ornamental animal, Black Swan (Cygnus atratus) has been imported into many countries. At present there has been no report of M. orientalis infection in Black Swan. In the present study M. orientalis infection in Black Swan was identified by a combination of different techniques, including morphological observation and molecular analysis. M. orientalis adults were found in the gallbladder and bile duct of a three-year-old female Black Swan, which was further confirmed by internal transcribed spacer (ITS) sequence analysis. In addition, the intermediate and definitive hosts of M. orientalis from the ‘Qing’ lake (a man-made lake in Changchun, China) that Black Swan lived were investigated and the infection route was preliminarily determined. Parafossarulus striatulus functioned as the first intermediate host which contained M. orientalis DNA, and fishes such as Pseudorasbora parva and Rhodeinae served as the second intermediate hosts with M. orientalis metacercariae in the fish flesh. M. orientalis eggs were found in the feces of three other Swans and six ducks that lived in the ‘Qing’ lake. This was the first reported case about M. orientalis infection of Black Swan. Our study described the course of the infection and provided new information about potential carriers and disseminators of M. orientalis.

•

M. orientalis infection was firstly reported in Black Swan in China.

•

P. striatulus served as the first intermediate host, P. parva and Rhodeinae as the second intermediate host. •Swans and ducks served as the definitive hosts.

Mitochondrial Genome Sequence of Echinostoma revolutum from Red-Crowned Crane (Grus japonensis)

Echinostoma revolutum is a zoonotic food-borne intestinal trematode that can cause intestinal bleeding, enteritis, and diarrhea in human and birds. To identify a suspected E. revolutum trematode from a red-crowned crane (Grus japonensis) and to reveal the genetic characteristics of its mitochondrial (mt) genome, the internal transcribed spacer (ITS) and complete mt genome sequence of this trematode were amplified. The results identified the trematode as E. revolutum. Its entire mt genome sequence was 15,714 bp in length, including 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and one non-coding region (NCR), with 61.73% A+T base content and a significant AT preference. The length of the 22 tRNA genes ranged from 59 bp to 70 bp, and their secondary structure showed the typical cloverleaf and D-loop structure. The length of the large subunit of rRNA (rrnL) and the small subunit of rRNA (rrnS) gene was 1,011 bp and 742 bp, respectively. Phylogenetic trees showed that E. revolutum and E. miyagawai clustered together, belonging to Echinostomatidae with Hypoderaeum conoideum. This study may enrich the mitochondrial gene database of Echinostoma trematodes and provide valuable data for studying the molecular identification and phylogeny of some digenean trematodes.

Characterization of the complete mitochondrial genome sequence of Tracheophilus cymbius (Digenea), the first representative from the family Cyclocoelidae

Tracheophilus cymbius (Trematoda: Cyclocoelidae) is a common tracheal fluke of waterfowl, causing serious loss in the poultry industry. However, taxonomic identification of T. cymbius remains controversial and confused. Mitochondrial (mt) genomes can provide genetic markers for the identification of closely related species. We determined the mt genome of T. cymbius and reconstructed phylogenies with other trematodes. The T. cymbius mt genome is 13,760 bp in size, and contains 12 protein-coding genes (cox 1-3, nad 1-6, nad 4L, cyt b and atp 6), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes and one non-coding region. All are transcribed in the same direction. The A + T content is 62.82%. ATG and TAG are the most common initiation and termination codons, respectively. Phylogenetic analyses of concatenated nucleotide sequences show T. cymbius grouping in suborder Echinostomata, and clustering together, with high statistical support, as a sister taxon with Echinochasmus japonicus (Echinochasmidae), the two forming a distinct branch rooted to the ancestor of all Echinostomatidae and Fasciolidae species. This is the first report of the T. cymbius mt genome, and the first reported mt genome within the family Cyclocoelidae. These data will provide a significant resource of molecular markers for studying the taxonomy, population genetics and systematics of trematodes.

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 Echinostoma miyagawai: Comparisons with closely related species and phylogenetic implications

Echinostoma miyagawai (Trematoda: Echinostomatidae) is a common parasite of poultry that also infects humans. Es. miyagawai belongs to the "37 collar-spined" or "revolutum" group, which is very difficult to identify and classify based only on morphological characters. Molecular techniques can resolve this problem. The present study, for the first time, determined, and presented the complete Es. miyagawai mitochondrial genome. A comparative analysis of closely related species, and a reconstruction of Echinostomatidae phylogeny among the trematodes, is also presented. The Es. miyagawai mitochondrial genome is 14,416 bp in size, and contains 12 protein-coding genes (cox1-3, nad1-6, nad4L, cytb, and atp6), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one non-coding region (NCR). All Es. miyagawai genes are transcribed in the same direction, and gene arrangement in Es. miyagawai is identical to six other Echinostomatidae and Echinochasmidae species. The complete Es. miyagawai mitochondrial genome A + T content is 65.3%, and full-length, pair-wise nucleotide sequence identity between the six species within the two families range from 64.2-84.6%. The Es. miyagawai sequences is most similar to Echinostoma caproni. Sequence difference are 15.0-33.5% at the nucleotide level, and 8.6-44.2% at the amino acid level, among the six species, for the 12 protein-coding genes. ATG and TAG are the most common initiation and termination codons, respectively. Twenty of the Es. miyagawai transfer RNA genes transcribe products of the conventional cloverleaf structure, while two of the transfer RNA genes, namely trnS1^(AGC) and trnS2^(UGA), have unpaired D-arms. Phylogenetic analyses using our mitochondrial data indicate that Es. miyagawai is closely related to other Echinostomatidae species, except for Echinostoma hortense, which forms a distinct paraphyletic branch, and Echinochasmus japonicus, which is outside the clade containing all other Echinostomatidae species. These phylogenetic results support the elevation of subfamily Echinostomatidae. Our dataset also provides a significant resource of molecular markers to study the taxonomy, population genetics, and systematics of the echinostomatids.

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.

De novo assembly and functional annotations of the transcriptome of Metorchis orientalis (trematoda: Opisthorchiidae)

Metorchis orientalis is a neglected zoonotic parasite, living in the gallbladder and bile duct of poultry and some mammals as well as humans. In spite of its economic and medical importance, the information known about the transcriptome and genome of M. orientalis is limited. In this study, we performed de novo sequencing, transcriptome assembly and functional annotations of the adult M. orientalis, obtained about 77.4 million high-quality clean reads, among which the length of the transcript contigs ranged from 100 to 11,249 nt with mean length of 373 nt and N50 length of 919 nt. We then assembled 31,943 unigenes, of which 20,009 (62.6%) were annotated by BLASTn and BLASTx searches against the available database. Among these unigenes, 19,795 (62.0%), 3407 (10.7%), 10,620 (33.2%) of them had significant similarity in the NR, NT and Swiss-Prot databases, respectively; 5744 (18.0%) and 4678 (14.6%) unigenes were assigned to GO and COG, respectively; and 9099 (28.5%) unigenes were identified and mapped onto 256 pathways in the KEGG Pathway database. Furthermore, we found that 98 (1.08%) unigenes were related to bile secretion and 5 (0.05%) to primary bile acid biosynthesis pathways category. The characterization of these transcriptomic data has implications for the better understanding of the biology of M. orientalis, and will facilitate the development of intervention agents for this and other pathogenic flukes of human and animal health significance.

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.

Metacercaria Infection Status of Fishborne Zoonotic Trematodes, Except for Clonorchis sinensis in Fish from the Heilongjiang Province, China

To investigate the metacercarial infections of fishborne zoonotic trematodes (FZT), a total of 6815 freshwater fish (in representing 13 species of 5 families) were collected from Songhua river (n = 2636), Nenjiang river (n = 1935), Mudanjiang river (n = 301), and other lakes or ponds (n = 1943) in 36 representative regions in Heilongjiang Province, China, from August 2012 to December 2015. Metacercariae of four FZT species, that is, Clonorchis sinensis, Metorchis orientalis, Isthmiophora hortensis, and Metagonimus yokogawai, metacercariae were detected in the examination by the artificial digestion method. As the partial data for C. sinensis were previously reported, the remaining three FZT species are to be treated in this study. The overall prevalence of M. orientalis, I. hortensis, and M. yokogawai, metacercariae was 10.54%, 0.28%, and 1.35%, respectively. Metacercariae of M. orientalis were detected in seven fish species, that is, Pseudorasbora parva, Hemiculter leucisculus, Saurogobio dabryi, Rhynchocypris lagowskii, Carassius auratus, Rhodeus ocellatus and Perccottus glehnii. Their prevalences were the highest in false dace, P. parva (26.81%), and in fish from Songhua river (17.94%). Metacercariae of I. hortensis were detected in only one fish species, Misgurnus anguillicaudatus, from Nenjiang river only. Metacercariae of M. yokogawai were detected in three fish species, that is, P. parva, H. leucisculus and S. dabryi. Their prevalences were the highest in sharpbelly, H. leucisculus (6.05%), and in fish from Mudanjiang river (5.65%). This study first demonstrated the existence of M. orientalis, I. hortensis, and M. yokogawai in freshwater fish from Heilongjiang Province, posing a major public health concern. Eight fish species, namely M. anguillicaudatus, P. parva, H. leucisculus, S. dabryi, R. lagowskii, C. auratus, R. ocellatus, and P. glehnii, cannot be eaten raw. Moreover, the findings of this study not only extended the second intermediate host range of FZT, but also improve the information of the distribution of FZT in China.

First report and molecular identification of Opisthorchis viverrini infection in human communities from Lower Myanmar

Opisthorchis viverrini is endemic in the South East Asian region, especially in Cambodia, Lao People's Democratic Republic, Vietnam and Thailand, but there have been no previous records from Myanmar. During stool surveys of rural populations in three regions of Lower Myanmar, Opisthorchis-like eggs were found in 34 out of 364 (9.3%) participants by stool microscopy after using the modified formalin-ether concentration technique. DNA was extracted from these positive stool samples and a portion of the mitochondrial cytochrome c oxidase subunit I (cox1) gene was amplified using the polymerase chain reaction and then sequenced. DNA sequences, successfully obtained from 18 of 34 positive samples (Bago Region, n = 13; Mon State, n = 3; Yangon Region, n = 2), confirmed that the eggs were of O. viverrini. Sequences showed 99.7% identity with O. viverrini mitochondrial cox1 (GenBank accession no. JF739555) but 95%, 88.7%, 82.6% and 81.4% identities with those of Opisthorchis lobatus from Lao People's Democratic Republic (GenBank accession nos. HQ328539-HQ328541), Metorchis orientalis from China (KT239342), Clonorchis sinensis from China (JF729303) and Opisthorchis felineus from Russia (EU921260), respectively. When alignement with other Opisthorchiidae trematodes, 81% similarity with Metorchis bilis from Czech Republic (GenBank accession nos. KT740966, KT740969, KT740970) and Slovakia (GenBank accession nos. KT740971-KT740973), 84.6% similarity with Metorchis xanthosomus from Czech Republic (GenBank accession no. KT740974), 78.6% similarity with M. xanthosomus from Poland (GenBank accession no. KT740968) and 82.2% similarity with Euamphimerus pancreaticus from Czech Republic (GenBank accession no. KT740975) were revealed. This study demonstrated, for the first time, O. viverrini from rural people in Myanmar using molecular methods and is an urgent call for surveillance and control activities against opisthorchiasis in Myanmar.

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.