Complete mitochondrial genome of the giant liver fluke Fascioloides magna (Digenea: Fasciolidae) and its comparison with selected trematodes

Phylogenetic relationship of Prohemistomum vivax to other trematodes based on the internal transcribed spacer region and mitochondrial genes

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 giant liver fluke in Europe: A review of Fascioloides magna within cervids and livestock with considerations on an expanding snail-fluke transmission risk

The giant liver fluke, Fascioloides magna, common in North America and introduced into Europe in the 19th century, is an underappreciated model system for epidemiological studies on biological invasions, interactions with other native parasites and for health impact assessments within both definitive and intermediate snail host populations. Owing to its first contemporary appearance in Europe and then its subsequent spread, fascioloidosis has become a fluke-livestock/wildlife-snail combination of increasing interest for veterinarians, parasitologists and population geneticists. Here, we present a description of its recent epidemiology, biogeography and biology, inclusive of host species lists. Special emphasis is placed upon known definitive hosts of F. magna within Europe and considerations upon this fluke's local intermediate snail hosts. This helps us envisage plausible future epidemiological scenarios for further expansion across Europe, potentially even invasion into the UK. In line with others who draw attention upon needs for better systematic monitoring of putative risk-areas of fluke transmission, we close by highlighting why better surveillance of F. magna across continental Europe, and neighbouring territories, is justified.

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.

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.

Characterization and comparative analysis of the complete mitochondrial genome of Azygia hwangtsiyui Tsin, 1933 (Digenea), the first for a member of the family Azygiidae

Azygia hwangtsiyui (Trematoda, Azygiidae), a neglected parasite of predatory fishes, is little-known in terms of its molecular epidemiology, population ecology and phylogenetic study. In the present study, the complete mitochondrial genome of A. hwangtsiyui was sequenced and characterized: it is a 13,973 bp circular DNA molecule and encodes 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as two non-coding regions. The A+T content of the A. hwangtsiyui mitogenome is 59.6% and displays a remarkable bias in nucleotide composition with a negative AT skew (-0.437) and a positive GC skew (0.408). Phylogenetic analysis based on concatenated amino acid sequences of twelve protein-coding genes reveals that A. hwangtsiyui is placed in a separate clade, suggesting that it has no close relationship with any other trematode family. This is the first characterization of the A. hwangtsiyui mitogenome, and the first reported mitogenome of the family Azygiidae. These novel datasets of the A. hwangtsiyui mt genome represent a meaningful resource for the development of mitochondrial markers for the identification, diagnostics, taxonomy, homology and phylogenetic relationships of trematodes.

Characterization and phylogenetic properties of the complete mitochondrial genome of Fascioloides jacksoni (syn. Fasciola jacksoni) support the suggested intergeneric change from Fasciola to Fascioloides (Platyhelminthes: Trematoda: Plagiorchiida)

Fascioloides jacksoni (syn. Fasciola jacksoni, Cobbold, 1869) (Platyhelminthes: Echinostomatoidea), is a liver fluke that causes severe morbidity and mortality of Asian elephants (Elephas maximus maximus). Understandings on molecular diagnosis, epidemiology, genetics and evolution of this flatworm are limited. In this study, we present the complete mitochondrial DNA (mt) sequence of 14,952 bp obtained from an individual fluke and comparative characterization of mitogenomic features with fasciolids, primarily, Fascioloides magna and other taxa in the superfamily Echinostomatoidea. Taxonomic relationship within and between Echinostomatoidea, Opisthorchioidea and Paramphistomoidea in the order Plagiorchiida, are also taxonomically considered. The complete circular mt molecule of Fas. jacksoni contained 12 protein-coding, two ribosomal RNA, 22 transfer RNA genes, and a non-coding region (NCR) rich in tandem repeat units. As common in digenean trematodes, Fas. jacksoni has the usual gene order, the absence of atp8 and the overlapped region by 40 bp between nad4L and nad4 genes. The NCR located between tRNA^Glu (trnE) and cox3 contained nine nearly identical tandem repeat units (TRs of 113 bp each). Special DHU-arm missing tRNAs for Serine were found for both, tRNA^S1(AGN) and tRNA^S2(UCN). Base composition indicated that cox1 of Fas. jacksoni showed the lowest (11.8% to Fas. magna, 12.9 - 13.6% to Fasciola spp. and 18.1% to Fasciolopsis buski) and nad6 the highest divergence rate (19.2%, 23.8-26.5% and 27.2% to each fasciolid group), respectively. A clear bias in nucleotide composition, as of 61.68%, 62.88% and 61.54%, with a negative AT-skew of the corresponding values (-0.523, -0.225 and - 0.426) for PCGs, MRGs and mtDNA for Fas. jacksoni and likewise data for the fasciolids. Phylogenetic analysis confirmed the sister branch of Fas. jacksoni and Fas. magna with the nodal support of 100%, clearly separated from the taxonomically recognized Fasciola spp. With the previous studies, mitogenomic data presented in this study are strongly supportive for Fasciola jacksoni reappraisal as Fascioloides jacksoni in the Fascioloides genus.

Description and phylogenetic analyses of ribosomal transcription units from species of Fasciolidae (Platyhelminthes: Digenea)

Many members of Fasciolidae are common trematodes in cattle, buffaloes, sheep, elephants, pigs, with some capable of infecting humans also. In this study, the complete or near-complete sequences of ribosomal transcription unit (rTU or rDNA), each of Fasciola hepatica (Australia), Fascioloides jacksoni (Sri Lanka), Fasciolopsis buski (Vietnam) and three isolates of F. gigantica (Vietnam), were obtained and characterized. The full length of rDNA for each F. hepatica, 'hybrid' Fasciola sp., Fas. jacksoni and Fa. Buski, was 7657 bp, 7966 bp, 7781 bp and 8361 bp, with the complete intergenic spacer region (IGS) (862 bp, 1170 bp, 987 bp and 561 bp), respectively. The rDNA of two 'pure' F. gigantica isolates from Vietnam was 6794 bp with unsequenced IGS. For 28S rRNA genes the Fasciola spp. are equal, 1958 bp for 18S, 160 bp for 5.8S, 3863 bp and 454 bp for ITS1 but ITS2 differ by one nucleotide (Thymine) (359 or 360 bp). The ITS1 of the sensu lato Fa. buski has some distinguishable features, 286 bp for ITS2, 3862 bp for 28S and four repeat units of 356-361 bp each found in ITS1. The 28S rDNA analysis showed the lowest level of divergence (0-0.57%) between F. hepatica and F. gigantica and higher (2.23-2.62%) and highest (6-6.42%) for Fas. jacksoni and Fasciolopsis, respectively. The tree of 43 strains/species clearly produced a well-supported phylogeny, where 18 fasciolids consistently grouped, forming a discrete Fasciolidae clade, distinct from Philophthalmidae, Echinostomatidae and Echinochasmidae in Echinostomatoidea. Fascioloides jacksoni is outside Fasciola spp.: basal with Fas. magna, as previously demonstrated.

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

The complete mitochondrial sequence of 17,030 bp was obtained from Echinostoma revolutum and characterized with those of previously reported members of the superfamily Echinostomatoidea, i.e. six echinostomatids, one echinochasmid, five fasciolids, one himasthlid, and two cyclocoelids. Relationship within suborders and between superfamilies, such as Echinostomata, Pronocephalata, Troglotremata, Opisthorchiata, and Xiphiditata, are also considered. It contained 12 protein-coding, two ribosomal RNA, 22 transfer RNA genes and a tandem repetitive consisting non-coding region (NCR). The gene order, one way-positive transcription, the absence of atp8 and the overlapped region by 40 bp between nad4L and nad4 genes were similar as in common trematodes. The NCR located between tRNAGlu (trnE) and cox3 contained 11 long (LRUs) and short repeat units (SRUs) (seven LRUs of 317 bp, four SRUs of 207 bp each), and an internal spacer sequence between LRU7 and SRU4 specifying high-level polymorphism. Special DHU-arm missing tRNAs for Serine were found for both tRNAS1(AGN) and tRNAS2(UCN). Echinostoma revolutum indicated the lowest divergence rate to E. miyagawai and the highest to Tracheophilus cymbius and Echinochasmus japonicus. The usage of ATG/GTG start and TAG/TAA stop codons, the AT composition bias, the negative AT-skewness, and the most for Phe/Leu/Val and the least for Arg/Asn/Asp codons were noted. Topology indicated the monophyletic position of E. revolutum to E. miyagawai. Monophyly of Echinostomatidae and Fasciolidae was clearly solved with respect to Echinochasmidae, Himasthlidae, and Cyclocoelidae which were rendered paraphyletic in the suborder Echinostomata.

A new species of Clinostomum Leidy, 1856 in East Asia based on genomic and morphological data

Metacercariae of Clinostomum Leidy, 1856 are frequently encountered in freshwater fish. In 2015, a provisional species of Clinostomum in People's Republic of China (PRC) was distinguished from C. complanatum (Rudolphi, 1819) in Europe based on divergent cytochrome c oxidase I (CO1). However, in subsequent studies in East Asia, the same divergent CO1 genotype was identified as C. complanatum. These matching sequences suggest that either the provisional East Asian species was incorrectly distinguished from C. complanatum in 2015 or that C. complanatum in East Asia was misidentified in later studies. We tested these alternatives by sequencing the mitochondrial genome of C. complanatum in Italy, which was 5.7% divergent from a previously published sequence from Clinostomum in PRC, including differences in 80 of 3390 (2.4%) translated amino acids. Partial CO1 sequences of specimens from PRC and those from Italy, Romania, and Turkey also each formed reciprocally monophyletic clades. Partial CO1 from the East Asian clade varied by mean 3.6% (range 2.4-4.8%) from C. complanatum from Italy, Romania, and Turkey; mean intra-clade CO1 variation was 0.3% (range 0-1.9%). Metacercariae from Europe and East Asia display significant morphometric variation, and data from the literature suggest morphological differences in the genital complex of adults. Although sequences of nuclear rDNA did not differ between isolates from the west and East Asia, taken together, these results lead us to describe a new species of Clinostomum.

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.

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.

Characterization of the complete mitochondrial genome of Uvitellina sp., representative of the family Cyclocoelidae and phylogenetic implications

Mitochondrial (mt) DNA has been useful in revealing the phylogenetic relationship of eukaryotic organisms including flatworms. Therefore, the use of mitogenomic data for the comparative and phylogenetic purposes is needed for those families of digenetic trematodes for which the mitogenomic data are still missing. Molecular data with sufficiently rich informative characters that can better resolve species identification, discrimination, and membership in different genera is also required for members of some morphologically difficult families of trematodes bearing few autapomorphic characters among its members. Here, the internal transcribed spacer (ITS) region of nuclear ribosomal DNA (rDNA) and the complete mt genome of the trematode Uvitellina sp. (Cyclocoelidae: Haematotrephinae) was determined and annotated. The mt genome of this avian trematode is 14,217 bp in length, containing 36 genes plus a single non-coding region. The ITS rDNA sequences were used for the pairwise sequence comparison of Uvitellina sp. with European cyclocoelid species, and the mitochondrial 12 protein-coding genes (PCGs) and two ribosomal RNA genes were used to evaluate the position of the family within selected trematodes. The ITS rDNA analysis of Uvitellina sp. showed less nucleotide differences with Hyptiasmus oculeus (16.77%) than with other European cyclocoelids (18.63-23.58%). The Bayesian inference (BI) analysis using the 12 mt PCGs and two rRNA genes supported the placement of the family Cyclocoelidae within the superfamily Echinostomatoidea (Plagiorchiida: Echinostmata). The availability of the mt genome sequences of Uvitellina sp. provides a novel resource of molecular markers for phylogenetic studies of Cyclocoelidae and other trematodes.

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.

Mitochondrial DNA Evidence Supports the Hypothesis that Triodontophorus Species Belong to Cyathostominae

Equine strongyles, the significant nematode pathogens of horses, are characterized by high quantities and species abundance, but classification of this group of parasitic nematodes is debated. Mitochondrial (mt) genome DNA data are often used to address classification controversies. Thus, the objectives of this study were to determine the complete mt genomes of three Cyathostominae nematode species (Cyathostomum catinatum, Cylicostephanus minutus, and Poteriostomum imparidentatum) of horses and reconstruct the phylogenetic relationship of Strongylidae with other nematodes in Strongyloidea to test the hypothesis that Triodontophorus spp. belong to Cyathostominae using the mt genomes. The mt genomes of Cy. catinatum, Cs. minutus, and P. imparidentatum were 13,838, 13,826, and 13,817 bp in length, respectively. Complete mt nucleotide sequence comparison of all Strongylidae nematodes revealed that sequence identity ranged from 77.8 to 91.6%. The mt genome sequences of Triodontophorus species had relatively high identity with Cyathostominae nematodes, rather than Strongylus species of the same subfamily (Strongylinae). Comparative analyses of mt genome organization for Strongyloidea nematodes sequenced to date revealed that members of this superfamily possess identical gene arrangements. Phylogenetic analyses using mtDNA data indicated that the Triodontophorus species clustered with Cyathostominae species instead of Strongylus species. The present study first determined the complete mt genome sequences of Cy. catinatum, Cs. minutus, and P. imparidentatum, which will provide novel genetic markers for further studies of Strongylidae taxonomy, population genetics, and systematics. Importantly, sequence comparison and phylogenetic analyses based on mtDNA sequences supported the hypothesis that Triodontophorus belongs to Cyathostominae.

Fasciolopsis buski (Digenea: Fasciolidae) from China and India may represent distinct taxa based on mitochondrial and nuclear ribosomal DNA sequences

Background

Fasciolopsis buski is a zoonotic intestinal fluke infecting humans and pigs, but it has been seriously neglected. It is yet to know whether there is any genetic diversity among F. buski from different geographical locations, particularly in sequences of nuclear ribosomal DNA (rDNA) and mitochondrial (mt) DNA. Therefore, we determined the sequences of partial 18S, the complete internal transcribed spacer (ITS) rDNA and the complete mt genome of F. buski from China, compared the rDNA and mtDNA sequences with those of isolates from India and Vietnam, and assessed the phylogenetic relationships of this fluke and related fasciolid trematodes based on the mtDNA dataset.

Results

The complete mt genome sequence of F. buski from China is 14,833 bp, with 36 genes, including 12 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes (rrnL and rrnS). The AT content of F. buski from China is 65.12%. The gene content and arrangement of the F. buski mt genome is similar to that of Fascioloides magna. Genetic distances between isolates of F. buski from China and India were high (28.2% in mtDNA, 13.2% in ITS-1 and 9.8% in ITS-2) and distinctly higher than the interspecific differences between Fasciola hepatica and Fasciola gigantica. The rDNA and mtDNA datasets for F. buski from China (isolate from pigs) and Vietnam (isolates from humans) were identical. The intergeneric differences in amino acid and nucleotide sequences among the genera Fasciolopsis, Fascioloides and Fasciola ranged between 24.64–25.56% and 26.35–28.46%, respectively.

Conclusions

Our results indicate that F. buski from China and India may represent distinct taxa, while F. buski in Vietnam and China represent the same species. These findings might have implications for the implementation of appropriate control strategies in different regions. Further studies are needed to decode mtDNA and rDNA sequences of F. buski from various geographical isolates for the better understanding of the species complex of F. buski.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2039-2) contains supplementary material, which is available to authorized users.