The mitochondrial genome and ribosomal operon of Brachycladium goliath (Digenea: Brachycladiidae) recovered from a stranded minke whale

Pathological, morphological, and molecular characterization of the infection by Brachycladium goliath in humpback whale Megaptera novaeangliae

The humpback whale (Megaptera novaeangliae) is a cosmopolitan migratory, seasonal mysticete that frequents the Brazilian coast. Strands of specimens may occur during the migratory stay in the country. In 2021 and 2022, three live humpback whales stranded on the coast of Rio Grande do Sul and Santa Catarina states in southern Brazil. After euthanasia, specimens were necropsied, and organs were thoroughly examined for lesions. Grossly, in all three cases, the liver exhibited multifocal, irregular, firm, white areas on the hepatic capsule, which extended into the parenchyma. On the cut surface, the livers were yellow to pale brown with orangish to greenish areas, the bile ducts were prominent, thickened, and severely dilated, and leaf-shaped flukes were found inside of them. Additionally, one case showed moderate atrophy of the right hepatic lobe. The histological findings included dilation of bile ducts, hyperplasia of the bile duct epithelium, marked inflammatory infiltration of lymphocytes, plasma cells, and eosinophils, and portal fibrosis. The parasite Brachycladium goliath was both morphologically and molecularly identified based on diagnostic key for trematodes and the original description of the species, and the amplification and sequencing of the ITS-2 region, respectively. Even though hepatic injury was not the primary cause of stranding, it may have contributed to the debilitation of the whales. To the authors' knowledge, this is the first study that reports M. novaeangliae as a definitive host of B. goliath and that describes the lesions caused by the parasite in cetaceans.

The complete mitochondrial genome of 3 species of allocreadiids (Digenea, Allocreadiidae): characterization and phylogenetic position within the order Plagiorchiida

Trematodes of the family Allocreadiidae are primarily found in the intestines of freshwater fishes around the world. The family includes 15 genera and c. 130 species. The last 2 decades have witnessed an increase in the genetic library of its species. Molecular data have been crucial for species delimitation and species description within Allocreadiidae and for understanding their evolutionary and biogeographical history and classification. Here, the mitogenomes of 3 species of allocreadiids were obtained using high throughput sequencing methods. Mitogenomes were compared with other members of the order Plagiorchiida to determine their molecular composition, gene rearrangement and phylogenetic interrelationships. The complete circular mitogenomes of Allocreadium lobatum, Creptotrematina aguirrepequenoi and Wallinia mexicana were 14 424, 13 769 and 13 924 bp long respectively, comprising 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 2 non-coding regions. Gene arrangements were identical to other Xiphidiatan trematodes. Phylogenetic analyses using the mitogenomes revealed Allocreadiidae as a monophyletic group closely related to other members of the suborder Xiphidiata; A. lobatum was yielded as the sister taxon of C. aguirrepequenoi + W. mexicana. Our study increases the complete mitochondrial genome library of trematodes and strengthens our understanding of the phylogenetic relationships and classification of this parasite group.

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.

Characterization of the Complete Mitochondrial Genome and Phylogenetic Analyses of Eurytrema coelomaticum (Trematoda: Dicrocoeliidae)

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.

A first record of digenean parasites of the dwarf sperm whale Kogia sima with morphological and molecular information

Two species of digenean trematodes of the family Brachycladiidae were obtained from two male dwarf sperm whales Kogia sima that stranded along the island of Kyushu, southern Japan in 2017. From the liver of the first animal, a single, large gravid specimen of a digenean species was collected. The morphological features were consistent with those of the genus Brachycladium. The worm had a large body and was characterized by anterior caeca without lateral diverticula, the shape of testes, ovary, and eggs. Molecular analyses using gene sequences of the 28S rRNA and the mitochondrial NADH dehydrogenase subunit 3 also supported the inclusion of this specimen into the genus Brachycladium. The identity of this worm is undetermined due to the lack of information on the genus and is reported as Brachycladium sp. From the cranial sinuses of the second animal, 33 specimens of digeneans were collected that were morphologically identified as Nasitrema gondo. This report documents a new host record for N. gondo, and the sequence information is provided for this digenean for the first time. This is the second record of digenean parasites for the family Kogiidae, and the first record with morphological and molecular information. The possibility of digenean infection in the liver and cranial sinus should be kept in mind during the necropsy of stranded kogiids.

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.

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.

No Tail No Fail: Life Cycles of the Zoogonidae (Digenea)

The Zoogonidae is the only digenean family where known cercariae lack the tail but actively search for the second intermediate host. However, the data on the zoogonid life cycles are scarce. In the present study, we elucidated and verified life cycles of the Zoogonidae from the White Sea. Using rDNA data, we showed that Pseudozoogonoides subaequiporus utilizes gastropods from the family Buccinidae as the first intermediate host and protobranch bivalves as the second one. This life cycle can be facultatively truncated: some cercariae of P. subaequiporus encyst within the daughter sporocysts. Molecular data also confirmed previous hypotheses on Zoogonoides viviapus life cycle with buccinid gastropods acting as the first intermediate hosts, and annelids and bivalves as the second intermediate hosts. We demonstrated the presence of short tail primordium in the developing cercariae of both species. Based on the reviewed and our own data, we hypothesize that the emergence of tailless cercariae in the evolution of the Zoogonidae is linked to the switch to non-arthropod second intermediate hosts, and that it possibly happened only in the subfamily Zoogoninae. Basally branching zoogonids have retained the ancestral second intermediate host and might have also retained the tail.

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.

A Phylogenetic Re-Evaluation of the Stenakrine Opecoelids (Trematoda, Digenea: Opecoeloidea) with Some Taxonomic Novelties

The Opecoeloidea is a large group of xiphidiate digeneans parasitizing marine and freshwater fishes. According to the current taxonomic model, this superfamily contains only one family with numerous subfamilies. This study is devoted to the members of the Stenakrinae. Based on phylogenetic analysis of concatenated sequences of 18S and 28S rRNA genes of stenakrine opecoelids Caudotestis dobrovolski, C. cf. dobrovolski, Hexagrammia zhukovi, Stenakron vetustum, as well as the deep-sea xiphidiate digenean Zdzitowieckitrema incognitum, which so far has had an ambiguous phylogenetic status, we erect a new opecoeloid family, the Zdzitowieckitrematidae fam. nov. The genera Holsworthotrema and Scorpidotrema are removed from the Stenakrinae to the Scorpidotrematinae subfam. nov. within the Opecoelidae. We also remove the Stenakrinae from the Opecoelidae and recognize it as a separate family within the Opecoeloidea. The Stenakridae stat. nov. is a sister taxon to a well-supported Opecoelidae. The Zdzitowieckitrematidae occupies a sister position relative to the stenakrids and the opecoelids taken together. All three families are clearly phylogenetically distinct, however convincing morphological differences are revealed only between the Zdzitowieckitrematidae and the Stenakridae and between the Opecoelidae and the Stenakridae.

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.

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.

Parasites of Selected Freshwater Snails in the Eastern Murray Darling Basin, Australia

Aquatic snails serve an important role in the ecosystem. They also play an essential role in the life cycle of many parasites as hosts and may pose risks to animal and human health. In Australia, the role of snails in the transmission of parasites of livestock is well studied. However, despite the country’s unique biodiversity and wildlife, little is known about the role of snails in the transmission and survival of parasites in other ecosystems, including aquatic and aquaculture systems. This study aimed to determine the occurrence of parasites in freshwater snails in the eastern Murray Darling Basin. A total of 275 snails were collected from various localities, including aquaculture fishery ponds and natural creeks during the summer and autumn months in the southern hemisphere. Three different species of freshwater snails, all common to the area, were found, including Bullastra lessoni (n = 11), Isidorella hainesii (n = 157), and Haitia acuta (n = 107), of which 9.1%, 1.3%, and 4.7%, respectively, were found to be harboring various developmental stages of Trematoda. No other parasite was found in the examined snails. Parasites were identified as Choanocotyle hobbsi, Plagiorchis sp. and Petasiger sp. based on the sequences of their ITS2, 18S, and 28S ribosomal DNA region. Herein, we report a native parasite Choanocotyle hobbsi in an introduced snail, Haitia acuta, from both natural and aquaculture ponds. As there are no genetic sequences for adult specimens of Petasiger spp. and Plagiorchis spp. collected in Australia for comparison, whether the specimens collected in this study are the larval stage of one of the previously described species or are a new, undescribed species cannot yet be determined. Our results also suggest snails collected from aquaculture ponds may be infected with considerably more parasites.

The biodiversity of marine trematodes: then, now and in the future

Trematodes are the richest class of platyhelminths in the marine environment, infecting all classes of marine vertebrates as sexual adults and many phyla of marine invertebrates as part of their life cycles. Despite the cryptic nature of their existence (almost all marine trematodes are internal parasites), they have been the focus of study for almost 250 years, with the first species described in 1774. Here we review progress in the study of the "biodiversity" of these parasites, contrasting the progress made in the last 50 years (post-1971) to that in the almost 200 years before it (pre-1972). We consider an understanding of biodiversity to require knowledge of the species present in the system, an understanding of their evolutionary relationships (which informs higher classification), and, specifically for trematodes, an understanding of their complex life cycles. The fauna is now large, comprising well over 5,000 species. Although species description continues, we see evidence of a slow-down in all aspects of discovery. There has been only one completely new family identified since 1984 and the proposal of new genera is in decline as is the description of new species, especially for those of tetrapods. However, the extent to which this slow-down reflects an approach to the richness asymptote is made uncertain by changes in the field; reduced effort and difficulty of study may be important components of the effect. Regardless of how close we are to a complete description of the fauna, we infer that the outline is well-understood although the details are not. Adoption of molecular methodologies over the last 40 years have complemented morphometric analyses to facilitate objective recognition of species; however, despite these objective data, there is still inconsistency between authors on species delimitation. Molecular methodologies have also completely revolutionised inference of relationships at all levels, from within genera to between orders, and underpinned elucidation of novel life cycles. We expect the next 50 years to produce further dividends from technological innovations. The backdrop to the field will be global environmental concerns and the growing problem of funding for basic biodiversity studies.

Life cycle truncation in Digenea, a case study of Neophasis spp. (Acanthocolpidae)

Truncated life cycles may emerge in digeneans if the second intermediate host is eliminated, and the first intermediate host, the mollusc, takes up its role. To understand the causes of this type of life cycle truncation, we analyzed closely related species of the genus Neophasis (Acanthocolpidae) with three-host and two-host life cycles. The life cycle of Neophasis anarrhichae involves two hosts: wolffishes of the genus Anarhichas as the definitive host and the common whelk Buccinum undatum as the intermediate host. Neophasis oculata, a closely related species with a three-host life cycle, would be a suitable candidate for the comparison, but some previous data on its life cycle seem to be erroneous. In this study, we aimed to redescribe the life cycle of N. oculata and to verify the life cycle of N. anarrhichae using molecular and morphological methods. Putative life cycle stages of these two species from intermediate hosts were linked with adult worms from definitive hosts using ribosomal molecular data: 18S, ITS1, 5.8S-ITS2, 28S. These markers did not differ within the species and were only slightly different between them. Intra- and interspecific variability was also estimated using mitochondrial COI gene. In the constructed phylogeny Neophasis spp. formed a common clade with two other genera of the Acanthocolpidae, Tormopsolus and Pleorchis. We demonstrated that the first intermediate hosts of N. oculata were gastropods Neptunea despecta and B. undatum (Buccinoidea). Shorthorn sculpins Myoxocephalus scorpius were shown to act as the second intermediate and definitive hosts of N. oculata. The previous reconstruction of the two-host life cycle of N. anarrhichae was reaffirmed. We suggest that life cycle truncation in N. anarrhichae was initiated by an acquisition of continuous morphogenesis in the hermaphroditic generation and supported by a strong prey-predator relationship between A. lupus and B. undatum.

•

Life cycles of digeneans Neophasis oculata and N. anarrhichae were redescribed.

•

Phylogenetic position of genus Neophasis within family Acanthocolpidae was established.

•

Causes of life cycle truncation within Neophasis are discussed.

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.

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.

An opisthorchiid concept of the genus Liliatrema (Trematoda: Plagiorchiida: Opisthorchioidea): an unexpected systematic position

Liliatrema is a small genus of trematodes consisting of two species. Its systematic position has long been debated, partly because of the confusing reports about the structure of male terminal genitalia. Here we test the phylogenetic position of the genus Liliatrema using data on complete 18S rRNA and partial 28S rRNA gene sequences obtained for Liliatrema skrjabini. We also provide a detailed description of terminal genitalia in adult specimens of L. sobolevi and metacercariae of both Liliatrema species. The results of the 28S rDNA-based phylogenetic analysis indicate that Liliatrema falls within a well-supported clade, which also includes Apophallus and traditional opisthorchiids. This clade, in turn, is nested within a well-supported clade, containing Euryhelmis, Cryptocotyle and Scaphanocephalus. In the 18S+8S rDNA analysis, Liliatrema appears as a sister-taxon to the Cryptocotyle + Euryhelmis group. The Liliatrema + (Cryptocotyle + Euryhelmis) clade is a well-supported sister-group to the traditional opisthorchiids. The morphology of the terminal genitalia of the liliatrematids also corresponds to that of the opisthorchioids. Thus, the results of our morphological and phylogenetic analyses favour an unexpected conclusion that the genus Liliatrema belongs to the Opisthorchioidea. We propose that the genera Liliatrema, Apophallus, Euryhelmis, Cryptocotyle and Scaphanocephalus belong, respectively, within the subfamilies Liliatrematinae, Apophallinae, Euryhelminthinae and Cryptocotylinae of the family Opisthorchiidae.

A new species of deropristid trematode from the sterlet Acipenser ruthenus (Actinopterygii: Acipenseridae) and revision of superfamily affiliation of the family Deropristidae

A new species, Skrjabinopsolus nudidorsalis sp. nov. is described from the sterlet Acipenser ruthenus, caught in the River Volga basin (Russia). This species differs from previously described congeners by the absence of vitelline follicles on the dorsal side of the body. The complete 18S rRNA and partial 28S rRNA gene sequences obtained for S. nudidorsalis are the first molecular data for the family Deropristidae. The results of phylogenetic analysis indicate that Deropristidae is sister to the Monorchiidae + Lissorchiidae group. The results of the phylogenetic study contradict the current taxonomic hypothesis that Deropristidae belongs to the superfamily Lepocreadioidea and allow inclusion of this family in Monorchioidea. The morphological similarity of deropristids to other monorchioids is recognizable from the presence of a bipartite internal seminal vesicle, spinous cirrus and a voluminous, armed metraterm.

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.

Molecular phylogenetics and mitogenomics of three avian dicrocoeliids (Digenea: Dicrocoeliidae) and comparison with mammalian dicrocoeliids

Background

The Dicrocoeliidae are digenetic trematodes mostly parasitic in the bile ducts and gall bladder of various avian and mammalian hosts. Until recently their systematics was based on morphological data only. Due to the high morphological uniformity across multiple dicrocoeliid taxa and insufficient knowledge of relative systematic value of traditionally used morphological characters, their taxonomy has always been unstable. Therefore, DNA sequence data provide a critical independent source of characters for phylogenetic inference and improvement of the system.

Methods

We examined the phylogenetic affinities of three avian dicrocoeliids representing the genera Brachylecithum, Brachydistomum and Lyperosomum, using partial sequences of the nuclear large ribosomal subunit (28S) RNA gene. We also sequenced the complete or nearly complete mitogenomes of these three isolates and conducted a comparative mitogenomic analysis with the previously available mitogenomes from three mammalian dicrocoeliids (from 2 different genera) and examined the phylogenetic position of the family Dicrocoeliidae within the order Plagiorchiida based on concatenated nucleotide sequences of all mitochondrial genes (except trnG and trnE).

Results

Combined nucleotide diversity, Kimura-2-parameter distance, non-synonymous/synonymous substitutions ratio and average sequence identity analyses consistently demonstrated that cox1, cytb, nad1 and two rRNAs were the most conserved and atp6, nad5, nad3 and nad2 were the most variable genes across dicrocoeliid mitogenomes. Phylogenetic analyses based on mtDNA sequences did not support the close relatedness of the Paragonimidae and Dicrocoeliidae and suggested non-monophyly of the Gorgoderoidea as currently recognized.

Conclusions

Our results show that fast-evolving mitochondrial genes atp6, nad5 and nad3 would be better markers than slow-evolving genes cox1 and nad1 for species discrimination and population level studies in the Dicrocoeliidae. Furthermore, the Dicrocoeliidae being outside of the clade containing other xiphidiatan trematodes suggests a need for the re-evaluation of the taxonomic content of the Xiphidiata.

New species and new genus of Pseudohaploporinae (Digenea): Pseudohaploporus pusitestis sp. n. and Parahaploporus elegantus n. g., sp. n. (Digenea: Pseudohaploporinae) from Vietnamese mullet fish

Two new species of Pseudohaploporinae, Pseudohaploporus pusitestis sp. n. and Parahaploporus elegantus n. g., sp. n., are described from intestines of the Vietnamese mullet fish Moolgarda seheli and Osteomugil cunnesius, respectively. Pseudohaploporus pusitestis sp. n. differs from two known Pseudohaploporus species, P. vietnamensis and P. planiliza, by the absence of a diverticulate hermaphroditic duct and muscular sphincters at the proximal end of the hermaphroditic sac. Metrically, P. pusitestis sp. n. is close to P. vietnamensis and differs from this species and from P. planilizum by lower maximum sizes of most parameters. Parahaploporus elegantus n. g., sp. differs from representatives of Pseudohaploporus by the presence of a single testis and the armament of hermaphroditic duct and is morphologically close to trematodes of the genus Haploporus. However, P. elegantus n. g differs from all known Haploporus species from mugilids of the Indo-West Pacific by the structure of the armament of the hermaphroditic duct and also by size of body, organs and eggs. The validity of designating two new species and a new genus of trematodes is supported by ITS and 28S rDNA sequence data. Phylogenetic reconstructions showed that the new trematodes belong to the Pseudohaploporinae, which formed a well-supported cluster within the monophyletic Haploporidae.

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

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

Testing the higher-level phylogenetic classification of Digenea (Platyhelminthes, Trematoda) based on nuclear rDNA sequences before entering the age of the 'next-generation' Tree of Life

Digenea Carus, 1863 represent a highly diverse group of parasitic platyhelminths that infect all major vertebrate groups as definitive hosts. Morphology is the cornerstone of digenean systematics, but molecular markers have been instrumental in searching for a stable classification system of the subclass and in establishing more accurate species limits. The first comprehensive molecular phylogenetic tree of Digenea published in 2003 used two nuclear rRNA genes (ssrDNA = 18S rDNA and lsrDNA = 28S rDNA) and was based on 163 taxa representing 77 nominal families, resulting in a widely accepted phylogenetic classification. The genetic library for the 28S rRNA gene has increased steadily over the last 15 years because this marker possesses a strong phylogenetic signal to resolve sister-group relationships among species and to infer phylogenetic relationships at higher levels of the taxonomic hierarchy. Here, we have updated the database of 18S and 28S rRNA genes until December 2017, we have added newly generated 28S rDNA sequences and we have reassessed phylogenetic relationships to test the current higher-level classification of digeneans (at the subordinal and subfamilial levels). The new dataset consisted of 1077 digenean taxa allocated to 106 nominal families for 28S and 419 taxa in 98 families for 18S. Overall, the results were consistent with previous higher-level classification schemes, and most superfamilies and suborders were recovered as monophyletic assemblages. With the advancement of next-generation sequencing (NGS) technologies, new phylogenetic hypotheses from complete mitochondrial genomes have been proposed, although their power to resolve deep levels of trees remains controversial. Since data from NGS methods are replacing other widely used markers for phylogenetic analyses, it is timely to reassess the phylogenetic relationships of digeneans with conventional nuclear rRNA genes, and to use the new analysis to test the performance of genomic information gathered from NGS, e.g. mitogenomes, to infer higher-level relationships of this group of parasitic platyhelminths.

Revision of the taxonomic status of Synthesium elongatum (Ozaki, 1935) (Brachycladiidae), an intestinal digenean of narrow-ridged finless porpoise (Neophocaena asiaeorientalis)

Synthesium elongatum (Brachycladiidae) is an intestinal digenean described from the finless porpoise (Neophocaena asiaeorientalis) in Japan. Few records of this species exist and there is a remarkable morphological similarity between S. elongatum and S. tursionis, such that a synonymy between the species has been suggested previously. However, no morphological and/or molecular analysis has been carried out to clarify the taxonomic status of S. elongatum. In this study, we collected specimens of Synthesium sp. from N. asiaeorientalis in western Japan. The specimens possess lobed testes within the third quarter of the body, a round ovary, and vitellaria extending to level of uterine field, which are diagnostic characters for both S. elongatum and S. tursionis. They were morphologically identified to S. elongatum or S. tursionis due to the fact that the available morphometric data for both species overlap remarkably. A molecular analysis of the mitochondrial ND3 gene showed that the pairwise nucleotide distances between these specimens and S. tursionis were small, and phylogenetic analysis showed that these specimens and S. tursionis are in the same clade. Therefore, it was indicated that S. elongatum and S. tursionis are the same species and, consequently, S. elongatum is a synonym of S. tursionis.

A new subfamily, Pseudohaploporinae subfam. n. (Digenea: Haploporidae), with morphometric and molecular analyses of two new species: Pseudohaploporus vietnamensis n. g., sp. n. and Pseudohaploporus planilizum n. g., sp. n. from Vietnamese mullet

Adult specimens representing Pseudohaploporus vietnamensis n. g., sp. n., Pseudohaploporus planilizum n. g., sp. n., and Pseudohaploporus sp. were collected from the intestines of Osteomugil engeli (Bleeker, 1858), Moolgarda seheli (Forsskål, 1775), and Planiliza subviridis (Valenciennes, 1836) from Vietnam. The new genus Pseudohaploporus n. g. is similar to the genus Haploporus Looss, 1902 in morphology and host-specificity. Species of Pseudohaploporus and Haploporus both possess an elongated body, a long oesophagus, an intestinal bifurcation behind the ventral sucker, a large hermaphroditic bursa, vitellarium consisting of two compact accumulations, concretions in the excretory vesicle, and numerous small eggs that contain miracidia with a pigment eye, and both infect mullet. In contrast to Haploporus, both species of the genus Pseudohaploporus possess two testes. These species differ from each other morphometrically. Pseudohaploporus planilizum sp. n., in comparison with P. vietnamensis sp. n., possess a larger maximum and mean size for the body length, oral and ventral suckers, pharynx, oesophagus length, hermaphrodite bursa, and egg width. Molecular data (ITS and 28S rDNA sequences) support the validity of the new genus and the two new species and indicate the presence of another representative of this genus, Pseudohaploporus sp. Phylogenetic reconstructions show that Pseudohaploporus formed a well-supported clade within a large clade that represents the Haploporidae Nicoll, 1914. Based on morphological and molecular data for P. vietnamensis sp. n. and P. planilizum sp. n., we erected a new subfamily, Pseudohaploporinae subfam. n., which appears to be a sister taxon to Haploporinae Nicoll, 1914. Molecular data show that Pseudohaploporus differ from other haploporids at a level similar to that between other subfamilies.

Validity of the Diplostomoidea and Diplostomida (Digenea, Platyhelminthes) upheld in phylogenomic analysis

Higher systematics within the Digenea, Carus 1863 have been relatively stable since a phylogenetic analysis of partial nuclear ribosomal markers (rDNA) led to the erection of the Diplostomida Olson, Cribb, Tkach, Bray, and Littlewood, 2003. However, recent mitochondrial (mt) genome phylogenies suggest this order might be paraphyletic. These analyses show members of two diplostomidan superfamilies are more closely related to the Plagiorchiida La Rue, 1957 than to other members of the Diplostomida. A recent phylogeny based on partial cytochrome c oxidase I also indicates one of the groups implicated, the Diplostomoidea Poirier, 1886, is non-monophyletic. To determine if these results were robust to additional taxon sampling, we analyzed mt genomes from seven diplostomoids in three families. To choose between phylogenetic alternatives based on mt genomes and the prior rDNA-based topology, we analyzed hundreds of ultra-conserved genomic elements assembled from shotgun sequencing. The Diplostomida was paraphyletic in the mt genome phylogeny but supported in the ultra-conserved genomic element phylogeny. We speculate this mitonuclear discordance is related to ancient, rapid radiation in the Digenea. Both ultra-conserved genomic elements and mt genomes support the monophyly of the Diplostomoidea and show congruent relationships within it. The Cyathocotylidae Mühling, 1898 are early diverging descendants of a paraphyletic clade of Diplostomidae Poirier, 1886, in which are nested members of the Strigeidae Railliet, 1919; the results support prior suggestions that the Crassiphialinae Sudarikov, 1960 will rise to the family level. Morphological traits of diplostomoid metacercariae appear to be more useful for differentiating clades than those of adults. We describe a new species of Cotylurus Szidat, 1928, resurrect a species of Hysteromorpha Lutz, 1931, and find support for a species of Alaria Schrank, 1788 of contested validity. Complete rDNA operons from seven diplostomoid species are provided as a resource for future studies.

The first next-generation sequencing approach to the mitochondrial phylogeny of African monogenean parasites (Platyhelminthes: Gyrodactylidae and Dactylogyridae)

BACKGROUND

Monogenean flatworms are the main ectoparasites of fishes. Representatives of the species-rich families Gyrodactylidae and Dactylogyridae, especially those infecting cichlid fishes and clariid catfishes, are important parasites in African aquaculture, even more so due to the massive anthropogenic translocation of their hosts worldwide. Several questions on their evolution, such as the phylogenetic position of Macrogyrodactylus and the highly speciose Gyrodactylus, remain unresolved with available molecular markers. Also, diagnostics and population-level research would benefit from the development of higher-resolution genetic markers. We aim to offer genetic resources for work on African monogeneans by providing mitogenomic data of four species (two belonging to Gyrodactylidae, two to Dactylogyridae), and analysing their gene sequences and gene order from a phylogenetic perspective.

RESULTS

Using Illumina technology, the first four mitochondrial genomes of African monogeneans were assembled and annotated for the cichlid parasites Gyrodactylus nyanzae, Cichlidogyrus halli, Cichlidogyrus mbirizei (near-complete mitogenome) and the catfish parasite Macrogyrodactylus karibae (near-complete mitogenome). Complete nuclear ribosomal operons were also retrieved, as molecular vouchers. The start codon TTG is new for Gyrodactylus and for Dactylogyridae, as is the incomplete stop codon TA for Dactylogyridae. Especially the nad2 gene is promising for primer development. Gene order was identical for protein-coding genes and differed between the African representatives of these families only in a tRNA gene transposition. A mitochondrial phylogeny based on an alignment of nearly 12,500 bp including 12 protein-coding and two ribosomal RNA genes confirms that the Neotropical oviparous Aglaiogyrodactylus forficulatus takes a sister group position with respect to the other gyrodactylids, instead of the supposedly 'primitive' African Macrogyrodactylus. Inclusion of the African Gyrodactylus nyanzae confirms the paraphyly of Gyrodactylus. The position of the African dactylogyrid Cichlidogyrus is unresolved, although gene order suggests it is closely related to marine ancyrocephalines.

CONCLUSIONS

The amount of mitogenomic data available for gyrodactylids and dactylogyrids is increased by roughly one-third. Our study underscores the potential of mitochondrial genes and gene order in flatworm phylogenetics, and of next-generation sequencing for marker development for these non-model helminths for which few primers are available.

Our love-hate relationship with DNA barcodes, the Y2K problem, and the search for next generation barcodes

DNA barcodes are very useful for species identification especially when identification by traditional morphological characters is difficult. However, the short mitochondrial and chloroplast barcodes currently in use often fail to distinguish between closely related species, are prone to lateral transfer, and provide inadequate phylogenetic resolution, particularly at deeper nodes. The deficiencies of short barcode identifiers are similar to the deficiencies of the short year identifiers that caused the Y2K problem in computer science. The resolution of the Y2K problem was to increase the size of the year identifiers. The performance of conventional mitochondrial COI barcodes for phylogenetics was compared with the performance of complete mitochondrial genomes and nuclear ribosomal RNA repeats obtained by genome skimming for a set of caddisfly taxa (Insect Order Trichoptera). The analysis focused on Trichoptera Family Hydropsychidae, the net-spinning caddisflies, which demonstrates many of the frustrating limitations of current barcodes. To conduct phylogenetic comparisons, complete mitochondrial genomes (15 kb each) and nuclear ribosomal repeats (9 kb each) from six caddisfly species were sequenced, assembled, and are reported for the first time. These sequences were analyzed in comparison with eight previously published trichopteran mitochondrial genomes and two triochopteran rRNA repeats, plus outgroup sequences from sister clade Lepidoptera (butterflies and moths). COI trees were not well-resolved, had low bootstrap support, and differed in topology from prior phylogenetic analyses of the Trichoptera. Phylogenetic trees based on mitochondrial genomes or rRNA repeats were well-resolved with high bootstrap support and were largely congruent with each other. Because they are easily sequenced by genome skimming, provide robust phylogenetic resolution at various phylogenetic depths, can better distinguish between closely related species, and (in the case of mitochondrial genomes), are backwards compatible with existing mitochondrial barcodes, it is proposed that mitochondrial genomes and rRNA repeats be used as next generation DNA barcodes.

Filling the gaps in the classification of the Digenea Carus, 1863: systematic position of the Proterodiplostomidae Dubois, 1936 within the superfamily Diplostomoidea Poirier, 1886, inferred from nuclear and mitochondrial DNA sequences

The Diplostomida Olson, Cribb, Tkach, Bray & Littlewood, 2003 is the less diverse order of the two orders within the subclass Digenea Carus, 1863 and is currently classified into three superfamilies, i.e. Brachylaimoidea Joyeux & Foley, 1930, Diplostomoidea Poirier, 1886, and Schistosomatoidea Stiles & Hassall, 1898. Although the suprageneric-level relationships have been elucidated with the use of molecular markers, the lack of representation of some groups obscure the phylogenetic relationships among families, rendering the classification unstable. Here, we tested the phylogenetic position of the family Proterodiplostomidae Dubois, 1936 based on partial 28S rDNA and complete 18S rDNA sequences for Crocodilicola pseudostoma (Willemoes-Suhm, 1870), a crocodile parasite that has been found as a progenetic metacercaria parasitising the pale catfish Rhamdia guatemalensis (Günther) in Mexico and in other siluruforms in the Neotropics. We augmented the representation of the species, genera and families within the Diplostomida, including mostly representatives of the superfamily Diplostomoidea, and assembled a dataset that contains 49 species for the 28S rRNA gene, and 45 species for the 18S rRNA gene. Additionally, we explored the phylogenetic signal of the mitochondrial gene cox1 in reconstructing the phylogenetic relationships of selected members of the superfamily. Our analyses showed that the family Proterodiplostomidae is the sister taxon to the paraphyletic Diplostomidae Poirier, 1886 and Strigeidae Railliet, 1919, with Cyathocotylidae Mühling, 1898 + Brauninidae Wolf, 1903 as their sister group. Analysis of concatenated 18S + 28S sequences revealed the Liolopidae Odhner, 1912 as the basal group of the superfamily Diplostomoidea, although analyses of independent datasets showed that the position of this family remains uncertain. Analysis based on cox1 unequivocally resolved the Proterodiplostomidae as the sister taxon to the Diplostomidae and Strigeidae, although the Cyathocotylidae was nested in a different clade, along with brachylaimoids and schistosomatoids.

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.

The ribosomal transcription units of Haplorchis pumilio and H. taichui and the use of 28S rDNA sequences for phylogenetic identification of common heterophyids in Vietnam

Background

Heterophyidiasis is now a major public health threat in many tropical countries. Species in the trematode family Heterophyidae infecting humans include Centrocestus formosanus, Haplorchis pumilio, H. taichui, H. yokogawai, Procerovum varium and Stellantchasmus falcatus. For molecular phylogenetic and systematic studies on trematodes, we need more prospective markers for taxonomic identification and classification. This study provides near-complete ribosomal transcription units (rTU) from Haplorchis pumilio and H. taichui and demonstrates the use of 28S rDNA sequences for identification and phylogenetic analysis.

Results

The near-complete ribosomal transcription units (rTU), consisting of 18S, ITS1, 5.8S, ITS2 and 28S rRNA genes and spacers, from H. pumilio and H. taichui from human hosts in Vietnam, were determined and annotated. Sequence analysis revealed tandem repetitive elements in ITS1 in H. pumilio and in ITS2 in H. taichui. A phylogenetic tree inferred from 28S rDNA sequences of 40 trematode strains/species, including 14 Vietnamese heterophyid individuals, clearly confirmed the status of each of the Vietnamese species: Centrocestus formosanus, Haplorchis pumilio, H. taichui, H. yokogawai, Procerovum varium and Stellantchasmus falcatus. However, the family Heterophyidae was clearly not monophyletic, with some genera apparently allied with other families within the superfamily Opisthorchioidea (i.e. Cryptogonimidae and Opisthorchiidae). These families and their constituent genera require substantial re-evaluation using a combination of morphological and molecular data. Our new molecular data will assist in such studies.

Conclusions

The 28S rDNA sequences are conserved among individuals within a species but varied between genera. Based on analysis of 40 28S rDNA sequences representing 19 species in the superfamily Opisthorchioidea and an outgroup taxon (Alaria alata, family Diplostomidae), six common human pathogenic heterophyids were identified and clearly resolved. The phylogenetic tree inferred from these sequences again confirmed anomalies in molecular placement of some members of the family Heterophyidae and demonstrates the need for reappraisal of the entire superfamily Opisthorchioidea. The new sequences provided here supplement those already available in public databases and add to the array of molecular tools that can be used for the diagnosis of heterophyid species in human and animal infections.