Comparative mitogenomics of the zoonotic parasite Echinostoma revolutum resolves taxonomic relationships within the 'E. revolutum' species group and the Echinostomata (Platyhelminthes: Digenea)

Mitogenomics of the zoonotic parasite Echinostoma miyagawai and insights into the evolution of tandem repeat regions within the mitochondrial non-coding control region

Echinostoma miyagawai is a cosmopolitan parasite within the Echinostomatidae and is a cause of human echinostomiasis. Species within the family have been a challenge to disentangle with E. miyagawai being synonyms of several other Echinostoma species. However, complete mitochondrial genomes have been shown to be vital in distinguishing echinostomatid species, but detailed comparisons of not only gene content but also structural features have been limited. Using long range sequencing techniques, the complete mitochondrial genome of E. miyagawai was sequenced and compared to other members of Echinostomatidae. In total 12 protein coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes were identified, as was an extensive noncoding control region (CR), consisting of 2 types of multiple tandem repeat units. Phylogenetic analyses of complete mitochondrial genomes corresponded to previous studies on single mitochondrial genes and nuclear ribosomal nuclear markers confirmed E. miyagawai to be within in the ‘Echinostoma revolutum’ group. The tandem repeat units found in the CR contained promoter sequences containing domains typical of initiation sites for replication and transcription as well as several palindromic regions which were shared between echinostomatid species. The study illustrates not only the utility complete mitogenomes in disentangling the relationship between these parasite species, but also provides some insight into the potential adaptations and other evolutionary processes that may govern the divergence of mitochondrial genomes for the first time in echinostomatids.

The ribosomal transcription units of five echinostomes and their taxonomic implications for the suborder Echinostomata (Trematoda: Platyhelminthes)

Five newly obtained nuclear ribosomal transcription unit (rTU) sequences from Echinostomatidae and Echinochasmidae are presented. The inter- and intrafamilial relationships of these and other families in the suborder Echinostomata are also analyzed. The sequences obtained are the complete rTU of Artyfechinostomum malayanum (9,499 bp), the near-complete rTU of Hypoderaeum conoideum (8,076 bp), and the coding regions (from 5'-terminus of 18S to 3'-terminus of 28S rRNA gene) in Echinostoma revolutum (6,856 bp), Echinostoma miyagawai (6,854 bp), and Echinochasmus japonicus (7,150 bp). Except for the longer first internal transcribed spacer (ITS1) in Echinochasmus japonicus, all genes and spacers were almost identical in length. Comprehensive maximum-likelihood phylogenies were constructed using the PhyML software package. The datasets were either the concatenated 28S + 18S rDNA sequences (5.7-5.8 kb) from 60 complete rTUs of 19 families or complete 28S sequences only (about 3.8-3.9 kb) from 70 strains or species of 22 families. The phylogenetic trees confirmed Echinostomatoidea as monophyletic. Furthermore, a detailed phylogeny constructed from alignments of 169 28S D1-D3 rDNA sequences (1.1-1.3 kb) from 98 species of 50 genera of 10 families, including 154 echinostomatoid sequences (85 species/42 genera), clearly indicated known generic relationships within Echinostomatidae and Echinochasmidae and relationships of families within Echinostomata and several other suborders. Within Echinostomatidae, Echinostoma, Artyfechinostomum, and Hypoderaeum appeared as monophyletic, while Echinochasmus (Echinochasmidae) was polyphyletic. The Echinochasmidae are a sister group to the Psilostomidae. The datasets provided here will be useful for taxonomic reappraisal as well as studies of evolutionary and population genetics in the superfamily Echinostomatoidea, the sole superfamily in the suborder Echinostomata.

Trematode Genomics and Proteomics

Trematode infections stand out as one of the frequently overlooked tropical diseases, despite their wide global prevalence and remarkable capacity to parasitize diverse host species and tissues. Furthermore, these parasites hold significant socio-economic, medical, veterinary and agricultural implications. Over the past decades, substantial strides have been taken to bridge the information gap concerning various "omic" tools, such as proteomics and genomics, in this field. In this edition of the book, we highlight recent progress in genomics and proteomics concerning trematodes with a particular focus on the advances made in the past 5 years. Additionally, we present insights into cutting-edge technologies employed in studying trematode biology and shed light on the available resources for exploring the molecular facets of this particular group of parasitic helminths.

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

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

Mitogenomic and nuclear ribosomal transcription unit datasets support the synonymy of Paragonimus iloktsuenensis and P. ohirai (Paragonimidae: Platyhelminthes)

The complete mitogenome (mtDNA) of nominal Paragonimus iloktsuenensis (Paragonimidae: Trematoda) and the nuclear ribosomal transcription unit (rTU) coding region (rTU*: from 5'-terminus of 18S to 3'-terminus of 28S rRNA gene, excluding the external spacer region) of this species and of P. ohirai were obtained and used to further support the previously suggested synonymy of these taxa in the P. ohirai complex. The complete mitogenome of P. iloktsuenensis was 14,827 bp long (GenBank: ON961029) and nearly identical to that of P. ohirai (14,818 bp; KX765277), with a 99.12% nucleotide identity. The rTU* was 7543 bp and 6932 bp in these two taxa, respectively. All genes and spacers in the rTU were identical in length, with exception of the first internal transcribed spacer, which contained multiple tandem repeat units (6.7 for P. iloktsuenensis and 5.7 for P. ohirai). There was near 100% identity for the rTU genes. The phylogenetic topology inferred from the mtDNA and from individual gene regions (partial cox1 of 387 bp and the ITS-2 of 282 bp - 285 bp) indicated a very close relationship consistent with synonymy of P. iloktsuenensis and P. ohirai. The datasets provided here will be useful for taxonomic reappraisal as well as studies of evolutionary and population genetics of the genus Paragonimus and family Paragonimidae.

Molecular epidemiological analyses reveal extensive connectivity between Echinostoma revolutum (sensu stricto) populations across Eurasia and species richness of zoonotic echinostomatids in England

Echinostoma revolutum (sensu stricto) is a widely distributed member of the Echinostomatidae, a cosmopolitan family of digenetic trematodes with complex life cycles involving a wide range of definitive hosts, particularly aquatic birds. Integrative taxonomic studies, notably those utilising nad1 barcoding, have been essential in discrimination of E. revolutum (s.s.) within the 'Echinostoma revolutum' species complex and investigation of its molecular diversity. No studies, however, have focussed on factors affecting population genetic structure and connectivity of E. revolutum (s.s.) in Eurasia. Here, we used morphology combined with nad1 and cox1 barcoding to determine the occurrence of E. revolutum (s.s.) and its lymnaeid hosts in England for the first time, in addition to other echinostomatid species Echinoparyphium aconiatum, Echinoparyphium recurvatum and Hypoderaeum conoideum. Analysis of genetic diversity in E. revolutum (s.s.) populations across Eurasia demonstrated haplotype sharing and gene flow, probably facilitated by migratory bird hosts. Neutrality and mismatch distribution analyses support possible recent demographic expansion of the Asian population of E. revolutum (s.s.) (nad1 sequences from Bangladesh and Thailand) and stability in European (nad1 sequences from this study, Iceland and continental Europe) and Eurasian (combined data sets from Europe and Asia) populations with evidence of sub-population structure and selection processes. This study provides new molecular evidence for a panmictic population of E. revolutum (s.s.) in Eurasia and phylogeographically expands the nad1 database for identification of echinostomatids.

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.

A new species of Echinostoma (Trematoda: Echinostomatidae) from the 'revolutum' group found in Brazil: refuting the occurrence of Echinostoma miyagawai (=E. robustum) in the Americas

Although Echinostoma robustum (currently a synonym of E. miyagawai) was reported in the Americas based on molecular data, morphological support on adult parasites is still required. Herein, a new species of Echinostoma is described based on worms found in a chicken from Brazil. Molecular phylogenetic analyses based on 28S (1063 bp), ITS (947 bp) and Nad-1 (442 bp) datasets reveal the inclusion of the new species within Echinostoma ‘revolutum’ species complex. Moreover, it was verified the conspecificity between cercariae previously identified as E. robustum in Brazil [identical ITS and only 0.3% of divergence (1 nucleotide) in Nad-1]. Species discovery analyses show that these two isolates form an independent lineage (species) among Echinostoma spp. Compared to E. miyagawai, the new species presents relatively high divergence in Nad-1 (7.88–9.09%). Morphologically, the specimens are distinguished from all nominal species from the ‘revolutum’ species complex by the more posterior position of the testes (length of post-testicular field as a proportion of body length about 20%). They further differ from E. miyagawai and South American Echinostoma spp. by the higher proportion of forebody to the body length. Therefore, combined molecular and morphological evidence supports the proposal of the species named here as Echinostoma pseudorobustum sp. nov.

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.

Mitochondrial genomes of two eucotylids as the first representatives from the superfamily Microphalloidea (Trematoda) and phylogenetic implications

Background

The Eucotylidae Cohn, 1904 (Superfamily: Microphalloidea), is a family of digeneans parasitic in kidneys of birds as adults. The group is characterized by the high level of morphological similarities among genera and unclear systematic value of morphological characters traditionally used for their differentiation. In the present study, we sequenced the complete or nearly complete mitogenomes (mt genome) of two eucotylids representing the genera Tamerlania (T. zarudnyi) and Tanaisia (Tanaisia sp.). They represent the first sequenced mt genomes of any member of the superfamily Microphalloidea.

Methods

A comparative mitogenomic analysis of the two newly sequenced eucotylids was conducted for the investigation of mitochondrial gene arrangement, contents and genetic distance. Phylogenetic position of the family Eucotylidae within the order Plagiorchiida was examined using nucleotide sequences of mitochondrial protein-coding genes (PCGs) plus RNAs using maximum likelihood (ML) and Bayesian inference (BI) methods. BI phylogeny based on concatenated amino acids sequences of PCGs was also conducted to determine possible effects of silent mutations.

Results

The complete mt genome of T. zarudnyi was 16,188 bp and the nearly complete mt genome of Tanaisia sp. was 13,953 bp in length. A long string of additional amino acids (about 123 aa) at the 5′ end of the cox1 gene in both studied eucotylid mt genomes has resulted in the cox1 gene of eucotylids being longer than in all previously sequenced digeneans. The rrnL gene was also longer than previously reported in any digenean mitogenome sequenced so far. The TΨC and DHU loops of the tRNAs varied greatly between the two eucotylids while the anticodon loop was highly conserved. Phylogenetic analyses based on mtDNA nucleotide and amino acids sequences (as a separate set) positioned eucotylids as a sister group to all remaining members of the order Plagiorchiida. Both ML and BI phylogenies revealed the paraphyletic nature of the superfamily Gorgoderoidea and the suborder Xiphidiata.

Conclusions

The average sequence identity, combined nucleotide diversity and Kimura-2 parameter distances between the two eucotylid mitogenomes demonstrated that atp6, nad5, nad4L and nad6 genes are better markers than the traditionally used cox1 or nad1 for the species differentiation and population-level studies of eucotylids because of their higher variability. The position of the Dicrocoeliidae and Eucotylidae outside the clade uniting other xiphidiatan trematodes strengthened the argument for the need for re-evaluation of the taxonomic content of the Xiphidiata.

Taxonomy of Echinostoma revolutum and 37-Collar-Spined Echinostoma spp.: A Historical Review

Echinostoma flukes armed with 37 collar spines on their head collar are called as 37-collar-spined Echinostoma spp. (group) or 'Echinostoma revolutum group'. At least 56 nominal species have been described in this group. However, many of them were morphologically close to and difficult to distinguish from the other, thus synonymized with the others. However, some of the synonymies were disagreed by other researchers, and taxonomic debates have been continued. Fortunately, recent development of molecular techniques, in particular, sequencing of the mitochondrial (nad1 and cox1) and nuclear genes (ITS region; ITS1-5.8S-ITS2), has enabled us to obtain highly useful data on phylogenetic relationships of these 37-collar-spined Echinostoma spp. Thus, 16 different species are currently acknowledged to be valid worldwide, which include E. revolutum, E. bolschewense, E. caproni, E. cinetorchis, E. deserticum, E. lindoense, E. luisreyi, E. mekongi, E. miyagawai, E. nasincovae, E. novaezealandense, E. paraensei, E. paraulum, E. robustum, E. trivolvis, and Echinostoma sp. IG of Georgieva et al., 2013. The validity of the other 10 species is retained until further evaluation, including molecular analyses; E. acuticauda, E. barbosai, E. chloephagae, E. echinatum, E. jurini, E. nudicaudatum, E. parvocirrus, E. pinnicaudatum, E. ralli, and E. rodriguesi. In this review, the history of discovery and taxonomic debates on these 26 valid or validity-retained species are briefly reviewed.