BACKGROUND: Digenetic trematodes, including blood flukes, intestinal flukes, liver flukes, lung flukes, and pancreatic flukes, are highly diverse and distributed widely. They affect at least 200 million people worldwide, so better understanding of their global distribution and prevalence are crucial for controlling and preventing human trematodiosis. Hence, this scoping review aims to conduct a comprehensive investigation on the spatio-temporal distribution and epidemiology of some important zoonotic digenetic trematodes.
METHODS: We conducted a scoping review by searching PubMed, Web of Science, Google Scholar, China National Knowledge Infrastructure, and Wanfang databases for articles, reviews, and case reports of zoonotic digenetic trematodes, without any restrictions on the year of publication. We followed the inclusion and exclusion criteria to identify relevant studies. And relevant information of the identified studies were collected and summarized.
RESULTS: We identified a total of 470 articles that met the inclusion criteria and were included in the review finally. Our analysis revealed the prevalence and global distribution of species in Schistosoma, Echinostoma, Isthmiophora, Echinochasmus, Paragonimus, Opisthorchiidae, Fasciolidae, Heterophyidae, and Eurytrema. Although some flukes are distributed worldwide, developing countries in Asia and Africa are still the most prevalent areas. Furthermore, there were some overlaps between the distribution of zoonotic digenetic trematodes from the same genus, and the prevalence of some zoonotic digenetic trematodes was not entirely consistent with their global distribution. The temporal disparities in zoonotic digenetic trematodes may attribute to the environmental changes. The gaps in our knowledge of the epidemiology and control of zoonotic digenetic trematodes indicate the need for large cohort studies in most countries.
CONCLUSIONS: This review provides important insights into the prevalence and global distribution of some zoonotic digenetic trematodes, firstly reveals spatio-temporal disparities in these digenetic trematodes. Countries with higher prevalence rate could be potential sources of transmitting diseases to other areas and are threat for possible outbreaks in the future. Therefore, continued global efforts to control and prevent human trematodiosis, and more international collaborations are necessary in the future.

Three new species of Gyrodactylus are described from three species of bitterling in Donghu Lake, China: Gyrodactylus ocellorhodei n. sp. from Rhodeus ocellatus; G. sinenorhodei n. sp. from Rhodeus sinensis; and G. acheilorhodei n. sp. from Acheilognathus macropterus. All the three new species showed similar opisthaptor morphology, especially the marginal hooks: all had a slender and perpendicular sickle shaft, and flat sickle base with distinct heel and inner arch which was different from the G. rhodei-group species parasitic on bitterling. Multivariate analyses based on hamulus and marginal hooks suggested that these three new species cannot be completely distinguished, despite some morphology divergence observed in certain less reliable morphometric features, such as hamulus root length, ventral bar total length and process shape. These three new species shared an identical 18S ribosomal RNA gene sequence, while the variation in the Internal Transcribed Spacers (ITS1-ITS2) sequence among them (8.4-11.2%, K2P) far exceeded the 1% ITS sequence difference that had been suggested as a threshold for species delimitation of Gyrodactylus. Phylogenetic analysis based on ITS1-ITS2 showed that all these sequenced Gyrodactylus spp. parasitic on the subfamily Acheilognathinae host formed a monophyletic group. However, a clear differentiation (18.9-20.9%, K2P of ITS1-ITS2) could be found between the subgroup from China (G. ocellorhodei n. sp., G. sinenorhodei n. sp. and G. acheilorhodei n. sp.) and that from Europe (G. rhodei).

BACKGROUND: Gyrodactylus is a lineage of monogenean flatworm ectoparasites exhibiting many features that make them a suitable model to study the host-parasite coevolutionary dynamics. Previous coevolutionary studies of this lineage mainly relied on low-power datasets (a small number of samples and a single molecular marker) and (now) outdated algorithms.
METHODS: To investigate the coevolutionary relationship of gyrodactylids and their fish hosts in high resolution, we used complete mitogenomes (including two newly sequenced Gyrodactylus species), a large number of species in the single-gene dataset, and four different coevolutionary algorithms.
RESULTS: The overall coevolutionary fit between the parasites and hosts was consistently significant. Multiple indicators confirmed that gyrodactylids are generally highly host-specific parasites, but several species could parasitize either multiple (more than 5) or phylogenetically distant fish hosts. The molecular dating results indicated that gyrodactylids tend to evolve towards high host specificity. Speciation by host switch was identified as a more important speciation mode than co-speciation. Assuming that the ancestral host belonged to Cypriniformes, we inferred four major host switch events to non-Cypriniformes hosts (mostly Salmoniformes), all of which occurred deep in the evolutionary history. Despite their relative rarity, these events had strong macroevolutionary consequences for gyrodactylid diversity. For example, in our dataset, 57.28% of all studied gyrodactylids parasitized only non-Cypriniformes hosts, which implies that the evolutionary history of more than half of all included lineages could be traced back to these major host switch events. The geographical co-occurrence of fishes and gyrodactylids determined the host use by these gyrodactylids, and geography accounted for most of the phylogenetic signal in host use.
CONCLUSIONS: Our findings suggest that the coevolution of Gyrodactylus flatworms and their hosts is largely driven by geography, phylogeny, and host switches.

Prohemistomum vivax is a zoonotic small cyathocotylid trematode that inhabits the intestines of fish-eating birds and mammals. Here, we amplified the internal transcribed spacer (ITS) sequence and six mitochondrial protein-coding genes (PCGs) from P. vivax. The ITS region was 1389 base pairs long and had a partial 18S ribosomal RNA gene, a full ITS1, 5.8S rRNA, and ITS2 sequence, and a partial 28S rRNA gene. The ITS region of P. vivax showed a minimum pairwise distance (0.3-0.6%) from the ITS sequences of Cyathocotylidae sp. 1 and 2 metacercariae from Clarias gariepinus. This result suggests that these metacercariae belong to P. vivax metacercariae. We first amplified mitochondrial genes from P. vivax, including cytochrome c oxidase subunit III (cox3) partial sequence; tRNA-His, cytochrome b (cytb), and NADH dehydrogenase subunit 4L (nad4L) complete sequences; and NADH dehydrogenase subunit 4 (nad4), cytochrome c oxidase I (cox1), and NADH dehydrogenase subunit 5 (nad5) partial sequences. P. vivax was most closely related to Cyathocotyle prussica (NC_039780) and Holostephanus sp. (OP082179), with cox1, cox3, and cytb genes conserved among the three trematodes. The ML phylogenetic tree of ITS sequences supports the order Diplostomida, divided into two main clades (the superfamily Diplostomoidea and Schistosomatoidea). The phylogeny of concatenated amino acid sequences of P. vivax six PCGs revealed that diplostomoids and Clinostomum sp. evolved in a clade with Plagiorchiida members, away from Schistosoma species. These results may yield ribosomal and mitochondrial genetic markers for molecular epidemiological investigations of cyathocotylid intestinal flukes.

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3\'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Eurytrema coelomaticum, a pancreatic fluke, is recognized as a causative agent of substantial economic losses in ruminants. This infection, commonly referred to as eurytrematosis, is a significant concern due to its detrimental impact on livestock production. However, there is a paucity of knowledge regarding the mitochondrial genome of E. coelomaticum. In this study, we performed the initial sequencing of the complete mitochondrial genome of E. coelomaticum. Our findings unveiled that the mitochondrial genome of E. coelomaticum spans a length of 15,831 bp and consists of 12 protein-coding genes, 22 tRNA genes, two rRNA genes, and two noncoding regions. The A+T content constituted 62.49% of the genome. Moreover, all 12 protein-coding genes of E. coelomaticum exhibit the same arrangement as those of E. pancreaticum and other published species belonging to the family Dicrocoeliidae. The presence of a short string of additional amino acids (approximately 20~23 aa) at the N-terminal of the cox1 protein in both E. coelomaticum and E. pancreaticum mitochondrial genomes has contributed to the elongation of the cox1 gene in genus Eurytrema, surpassing that of all previously sequenced Dicrocoeliidae. The phylogenetic analysis displayed a close relationship between E. coelomaticum and E. pancreaticum, along with a genus-level association between Eurytrema and Lyperosomum. These findings underscore the importance of mitochondrial genomic data for comparative studies of Dicrocoeliidae and even Digenea, offering valuable DNA markers for future investigations in the systematic, epidemiological, and population genetic studies of this parasite and other digenean trematodes.

BACKGROUND: Schistosoma mekongi is a human blood fluke causing schistosomiasis that threatens approximately 1.5 million humans in the world. Nonetheless, the limited available S. mekongi genomic resources have hindered understanding of its biology and parasite-host interactions for disease management and pathogen control. The aim of our study was to integrate multiple technologies to construct a high-quality chromosome-level assembly of the S. mekongi genome.
METHODS: The reference genome for S. mekongi was generated through integrating Illumina, PacBio sequencing, 10 × Genomics linked-read sequencing, and high-throughput chromosome conformation capture (Hi-C) methods. In this study, we conducted de novo assembly, alignment, and gene prediction to assemble and annotate the genome. Comparative genomics allowed us to compare genomes across different species, shedding light on conserved regions and evolutionary relationships. Additionally, our transcriptomic analysis focused on genes associated with parasite-snail interactions in S. mekongi infection. We employed gene ontology (GO) enrichment analysis for functional annotation of these genes.
RESULTS: In the present study, the S. mekongi genome was both assembled into 8 pseudochromosomes with a length of 404 Mb, with contig N50 and scaffold N50 lengths of 1168 kb and 46,759 kb, respectively. We detected that 43% of the genome consists of repeat sequences and predicted 9103 protein-coding genes. We also focused on proteases, particularly leishmanolysin-like metalloproteases (M8), which are crucial in the invasion of hosts by 12 flatworm species. Through phylogenetic analysis, it was discovered that the M8 gene exhibits lineage-specific amplification among the genus Schistosoma. Lineage-specific expansion of M8 was observed in blood flukes. Additionally, the results of the RNA-seq revealed that a mass of genes related to metabolic and biosynthetic processes were up-regulated, which might be beneficial for cercaria production.
CONCLUSIONS: This study delivers a high-quality, chromosome-scale reference genome of S. mekongi, enhancing our understanding of the divergence and evolution of Schistosoma. The molecular research conducted here also plays a pivotal role in drug discovery and vaccine development. Furthermore, our work greatly advances the understanding of host-parasite interactions, providing crucial insights for schistosomiasis intervention strategies.

Based on morphology and ITS sequence data, we identify and supplementally describe Gyrodactylus pseudorasborae Ondračková, Seifertová & Tkachenko, 2023 on the fins of topmouth gudgeon (Pseudoraspora parva) from freshwaters of southern China. The highest similarity (99.57% and 99.47%) to G. pseudorasborae suggested they were the same species. Prevalence and mean intensity were 45% and 2.3, respectively. The gyrodactylid species morphologically resembled G. pseudorasborae recorded from the same host species P. parva in Czech Republic, Ukraine, and Central China. But there were slight morphological differences in the shape and size of the marginal hook. Comparisons of marginal hook sickles of various Gyrodactylus species suggested that G. pseudorasborae and G. parvae were members of the G. wageneri-group. A molecular phylogeny of G. pseudorasborae with related species is presented and discussed within the context of the mechanism of local evolution of these sister species.

Prosthogonimiasis poses a threat to the reproductive system of poultry and wild birds, which are the definitive hosts of the parasite causing this disease. However, the parasite infection of the second intermediate host (dragonfly), the primary vector of this pathogen, is rarely reported. In this study, the prevalence of Prosthogonimus infection in dragonflies was investigated from June 2019 to October 2022 in Heilongjiang Province, northeast China. The species of metacercariae isolated from dragonfly were identified by morphological characteristics, molecular biology techniques, and animal infection experiments. The results showed that 11 species of dragonflies and one damselfly were identified and among six of the dragonflies infected by Prosthogonimus metacercariae, Sympetrum depressiusculum (28.53%) had the highest infection rate among all positive dragonflies, followed by Sympetrum vulgatum (27.86%) and Sympetrum frequens (20.99%), which are preferred hosts, and the total prevalence was 20.39% (2061/10,110) in Heilongjiang Province. Three species of Prosthogoniumus metacercariae were isolated, including Prosthogonimus cuneatus, Prosthogonimus pullucidus, and Prosthogonimus sp., among which P. cuneatus was the dominant species in dragonflies in Heilongjiang Province. This is the first report on the prevalence of Prosthogonimus in dragonflies in China, which provides baseline data for the control of prosthogonimiasis in Heilongjiang Province and a reference for the prevention of prosthogonimiasis in other areas of China.

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.