Closing the mitochondrial circle on paraphyly of the Monogenea (Platyhelminthes) infers evolution in the diet of parasitic flatworms

Molecular and SEM studies on Thaparocleidus vistulensis (Siwak, 1932) (Monopisthocotyla, Ancylodiscoididae)

Presenting new molecular and scanning electron microscope (SEM) features, this study gives additional data to the better knowledge of Thaparocleidus vistulensis (Siwak, 1932) (Monopisthocotyla, Ancylodiscoididae), a parasite of the European catfish Silurus glanis Linnaeus, 1758 (Siluriformes, Siluridae) cultured in a commercial fish farm in Hungary. In addition, notes on the early development of sclerotized anchors are also provided. The main morphological difference of T. vistulensis compared to other congeneric species is associated with the male copulatory organ, which exhibits 5-7 loops in the middle of the penis length and a long open V-shaped sclerotized accessory piece, dividing terminally into two parts, securing the terminal part of the penis tube. The present study provides for the first time molecular characterization data based on the 2694 bp long nucleotide sequence of rDNA (ITS1, 5.8S, ITS2, and flanked with partial 18S and partial 28S) submitted in GenBank with the accession number OR916383. A phylogenetic tree based on ITS1 sequences supports a well-defined clade including T. vistulensis, forming a sister group with T. siluri, a species-specific monopisthocotylan parasite to S. glanis. The morphological characterization of T. vistulensis, especially for the male copulatory organ, together with the molecular data in the present study, extends knowledge about this monopisthocotylan species and provides new information for future phylogeny studies.

The evolution of endoparasitism and complex life cycles in parasitic platyhelminths

Within flatworms, the vast majority of parasitism is innate to Neodermata, the most derived and diversified group of the phylum Platyhelminthes.1,2 The four major lineages of Neodermata maintain various combinations of life strategies.3 They include both externally (ecto-) and internally feeding (endo-) parasites. Some lineages complete their life cycles directly by infecting a single host, whereas others succeed only through serial infections of multiple hosts of various vertebrate and invertebrate groups. Food sources and modes of digestion add further combinatorial layers to the often incompletely understood mosaic of neodermatan life histories. Their evolutionary trajectories have remained molecularly unresolved because of conflicting evolutionary inferences and a lack of genomic data.4 Here, we generated transcriptomes for nine early branching neodermatan representatives and performed detailed phylogenomic analyses to address these critical gaps. Polyopisthocotylea, mostly hematophagous ectoparasites, form a group with the mostly hematophagous but endoparasitic trematodes (Trematoda), rather than sharing a common ancestor with Monopisthocotylea, ectoparasitic epithelial feeders. Phylogenetic placement of the highly specialized endoparasitic Cestoda alters depending on the model. Regardless of this uncertainty, this study brings an unconventional perspective on the evolution of platyhelminth parasitism, rejecting a common origin for the endoparasitic lifestyle intrinsic to cestodes and trematodes. Instead, our data indicate that complex life cycles and invasion of vertebrates' gut lumen, the hallmark features of these parasites, evolved independently within Neodermata. We propose the demise of the traditionally recognized class Monogenea and the promotion of its two subclasses to the class level as Monopisthocotyla new class and Polyopisthocotyla new class.

Monogenoidea Communities Differ Between Male Morphotypes of Bluegill Sunfish (Lepomis macrochirus)

PURPOSE

Bluegill sunfish (Lepomis macrochirus) are a North American sport fish that exhibit two different male morphotypes. Alpha-males are large, colorful, territorial, and have significant parental investment, while β-males are small, drab, and have two reproductive phenotypes, neither of which involves parental investment. Due to the two different reproductive strategies, the level of interaction between congenerics varies, which may impact parasites that are spread via close contact, such as Monogenoidea that infect the gills. Monogeneans are ectoparasites on the gills and skin of fish hosts that may cause significant pathology in high numbers and be an indicator of host behavior and interaction between hosts.

METHODS

In this study, 328 L. macrochirus (106 α-males, 92 β-males, and 130 females) were necropsied from 8 lakes and ponds in northwestern Virginia to identify and enumerate monogenean parasites from the gills.

RESULTS

Alpha-males had a significantly greater parasite abundance and species richness compared to β-males. This may be due to the increased size and surface area of gills in α-males, increased interaction with females during mating, and stationary behavior when guarding nests, which enhanced the chances of α-males contracting the parasites. This also led to significant differences in the monogenean communities infecting the two morphotypes, which were also significantly influenced by the size of the hosts.

CONCLUSIONS

It is important in future studies regarding parasitism that behavioral morphotypes within a sex, such as α-male and β-male L. macrochirus in this study, be treated separately as behavioral and morphometric differences between them can potentially lead to differences in parasitism.

Mitochondrial phylogenomics provides conclusive evidence that the family Ancyrocephalidae is deeply paraphyletic

BACKGROUND

Unresolved taxonomic classification and paraphyly pervade the flatworm class Monogenea: the class itself may be paraphyletic and split into Polyopisthocotylea and Monopisthocotylea; there are some indications that the monopisthocotylean order Dactylogyridea may also be paraphyletic; single-gene markers and some morphological traits indicate that the family Ancyrocephalidae is paraphyletic and intertwined with the family Dactylogyridae.

METHODS

To attempt to study the relationships of Ancyrocephalidae and Monopisthocotylea using a phylogenetic marker with high resolution, we sequenced mitochondrial genomes of two fish ectoparasites from the family Dactylogyridae: Dactylogyrus simplex and Dactylogyrus tuba. We conducted phylogenetic analyses using three datasets and three methods. Datasets were ITS1 (nuclear) and nucleotide and amino acid sequences of almost complete mitogenomes of almost all available Monopisthocotylea mitogenomes. Methods were maximum likelihood (IQ-TREE), Bayesian inference (MrBayes) and CAT-GTR (PhyloBayes).

RESULTS

Both mitogenomes exhibited the ancestral gene order for Neodermata, and both were compact, with few and small intergenic regions and many and large overlaps. Gene sequences were remarkably divergent for nominally congeneric species, with only trnI exhibiting an identity value > 80%. Both mitogenomes had exceptionally low A + T base content and AT skews. We found evidence of pervasive compositional heterogeneity in the dataset and indications that base composition biases cause phylogenetic artefacts. All six mitogenomic analyses produced unique topologies, but all nine analyses produced topologies that rendered Ancyrocephalidae deeply paraphyletic. Mitogenomic data consistently resolved the order Capsalidea as nested within the Dactylogyridea.

CONCLUSIONS

The analyses indicate that taxonomic revisions are needed for multiple Polyopisthocotylea lineages, from genera to orders. In combination with previous findings, these results offer conclusive evidence that Ancyrocephalidae is a paraphyletic taxon. The most parsimonious solution to resolve this is to create a catch-all Dactylogyridae sensu lato clade comprising the current Ancyrocephalidae, Ancylodiscoididae, Pseudodactylogyridae and Dactylogyridae families, but the revision needs to be confirmed by another marker with a sufficient resolution.

Exploring Evolutionary Relationships within Neodermata Using Putative Orthologous Groups of Proteins, with Emphasis on Peptidases

The phylogenetic relationships within Neodermata were examined based on putative orthologous groups of proteins (OGPs) from 11 species of Monogenea, Trematoda, and Cestoda. The dataset included OGPs from BUSCO and OMA. Additionally, peptidases were identified and evaluated as phylogenetic markers. Phylogenies were inferred using the maximum likelihood method. A network analysis and a hierarchical grouping analysis of the principal components (HCPC) of orthologous groups of peptidases were performed. The phylogenetic analyses showed the monopisthocotylean monogeneans as the sister-group of cestodes, and the polyopisthocotylean monogeneans as the sister-group of trematodes. However, the sister-group relationship between Monopisthocotylea and Cestoda was not statistically well supported. The network analysis and HCPC also showed a cluster formed by polyopisthocotyleans and trematodes. The present study supports the non-monophyly of Monogenea. An analysis of mutation rates indicated that secreted peptidases and inhibitors, and those with multiple copies, are under positive selection pressure, which could explain the expansion of some families such as C01, C19, I02, and S01. Whilst not definitive, our study presents another point of view in the discussion of the evolution of Neodermata, and we hope that our data drive further discussion and debate on this intriguing topic.

Revision of the systematics of the Polystomoidinae (Platyhelminthes, Monogenea, Polystomatidae) with redefinition of Polystomoides Ward, 1917 and Uteropolystomoides Tinsley, 2017

Polystomatids are platyhelminth parasites that infect mainly amphibians and freshwater turtles. For more than seven decades, chelonian polystomes were classified into three genera according to the number of hamuli, i.e. absent for Neopolystoma, one pair for Polystomoidella and two pairs for Polystomoides. Following re-examination of morphological characters, seven new genera were erected the past six years, namely Apaloneotrema, Aussietrema, Fornixtrema, Manotrema, Pleurodirotrema, Uropolystomoides and Uteropolystomoides. However, the polyphyly of Neopolystoma and Polystomoides on the one hand, and the nested position of Uteropolystomoides within a clade encompassing all Neopolystoma and Polystomoides spp. on the other, still raised questions about the validity of these genera. We therefore re-examined several types, paratypes and voucher specimens, and investigated the molecular phylogeny of polystomes sampled from the oral cavity of North American turtles to re-evaluate their systematic status. We show that all Polystomoides Ward, 1917, sensu Du Preez et al., 2022, Neopolystoma Price, 1939, sensu Du Preez et al., 2022 and Uteropolystomoides Tinsley, 2017 species, display vaginae that are peripheral and extend well beyond the intestine. We thus reassign all species of the clade to Polystomoides and propose nine new combinations; however, although Uteropolystomoides is nested within this clade, based on its unique morphological features, we propose to keep it as a valid taxon. Polystomoides as redefined herein groups all polystome species infecting either the oral cavity or the urinary bladder of cryptodires, with peripheral vaginae and with or without two pairs of small hamuli. Uteropolystomoides nelsoni (Du Preez & Van Rooyen 2015), originally described from Pseudemys nelsoni Carr is now regarded as Uteropolystomoides multifalx (Stunkard, 1924) n. comb. infecting three distinct Pseudemys species of North America.

Neuromuscular organization and haptoral armament of Polyclithrum ponticum (Monogenea: Gyrodactylidae)

Most gyrodactylids have a haptor armed with a pair of hamuli, two connecting bars and 16 marginal hooks. In some gyrodactylids, however, the haptor is disc-shaped and reinforced by additional sclerites. The genus Polyclithrum has arguably the most elaborate haptor in this group. This study aimed to gain better understanding of the anatomy of Polyclithrum by examining neuromusculature and haptoral armament of Polyclithrum ponticum, a species parasitizing Mugil cephalus in the Black Sea, with emphasis on haptoral sclerites and musculature in connection with host-attachment mechanisms. Musculature was stained by phalloidin, the nervous system by anti-serotonin and anti-FMRFamide antibodies, and haptoral sclerites were visualized in reflected light. The study provided new information on sclerites: in addition to previously described supplementary sclerites (A1-6), ear-shaped sclerites (ESSs) and two paired groups of ribs, reflected light revealed a rod-shaped process on the ESSs and a pair of small posterior sclerites. The sclerites were shown to be operated by 16 muscles, the most prominent of which were two transverse muscles connecting the hamular roots, three muscles attached to sclerite A2, the muscle fibres of anterior ribs and a set of extrinsic muscles. The nervous system consists of a pair of cerebral ganglia connected by a commissure and three pairs of nerve cords that unite in the haptor to form a loop between the opposite cords. The arrangement of sclerites and muscles suggests that Polyclithrum initiates the attachment by clamping a host's surface with longitudinally folded haptor and then secures its position with marginal hooks.

Evolutionary Adaptations of Parasitic Flatworms to Different Oxygen Tensions

During the evolution of the Earth, the increase in the atmospheric concentration of oxygen gave rise to the development of organisms with aerobic metabolism, which utilized this molecule as the ultimate electron acceptor, whereas other organisms maintained an anaerobic metabolism. Platyhelminthes exhibit both aerobic and anaerobic metabolism depending on the availability of oxygen in their environment and/or due to differential oxygen tensions during certain stages of their life cycle. As these organisms do not have a circulatory system, gas exchange occurs by the passive diffusion through their body wall. Consequently, the flatworms developed several adaptations related to the oxygen gradient that is established between the aerobic tegument and the cellular parenchyma that is mostly anaerobic. Because of the aerobic metabolism, hydrogen peroxide (H₂O₂) is produced in abundance. Catalase usually scavenges H₂O₂ in mammals; however, this enzyme is absent in parasitic platyhelminths. Thus, the architecture of the antioxidant systems is different, depending primarily on the superoxide dismutase, glutathione peroxidase, and peroxiredoxin enzymes represented mainly in the tegument. Here, we discuss the adaptations that parasitic flatworms have developed to be able to transit from the different metabolic conditions to those they are exposed to during their life cycle.

The complete mitochondrial genomes of Paradiplozoon yarkandense and Paradiplozoon homoion confirm that Diplozoidae evolve at an elevated rate

BACKGROUND

Diplozoidae are monogenean (Monogenea: Polyopisthocotylea) fish parasites characterised by a unique life history: two larvae permanently fuse into an X-shaped "Siamese" organism. Taxonomy and phylogeny of Diplozoidae and Polyopisthocotylea remain unresolved due to the unavailability of molecular markers with sufficiently high resolution. Mitogenomes may be a suitable candidate, but there are currently only 12 available for the Polyopisthocotylea (three for Diplozoidae). The only available study of diplozoid mitogenomes found unique base composition patterns and elevated evolution rates in comparison with other Monogenean mitogenomes.

METHODS

To further explore their evolution and generate molecular data for evolutionary studies, we sequenced the complete mitogenomes of two Diplozoidae species, Paradiplozoon homoion and Paradiplozoon yarkandense, and conducted a number of comparative mitogenomic analyses with other polyopisthocotyleans.

RESULTS

We found further evidence that mitogenomes of Diplozoidae evolve at a unique, elevated rate, which was reflected in their exceptionally long branches, large sizes, unique base composition, skews, and very low gene sequence similarity levels between the two newly sequenced species. They also exhibited remarkably large overlaps between some genes. Phylogenetic analysis of Polyopisthocotylea resolved all major taxa as monophyletic, and Mazocraeidea was split into two major clades: (Diplozoidae) + (all four remaining families: Diclidophoridae, Chauhaneidae, Mazocraeidae and Microcotylidae). It also provided further confirmation that the genus Paradiplozoon is paraphyletic and requires a taxonomic revision, so the two species may have to be renamed Indodiplozoon homoion and Diplozoon yarkandense comb. nov.

CONCLUSIONS

Although our findings indicate that mitogenomes may be a promising tool for resolving the phylogeny of Polyopisthocotylea, elevated evolutionary rates of Diplozoidae may cause phylogenetic artefacts, so future studies should pay caution to this problem. Furthermore, as the reason for their elevated evolution remains unknown, Diplozoidae are a remarkably interesting lineage for other types of evolutionary mitogenomic studies.

Novel molecular data for monogenean parasites of sparid fishes in the Mediterranean and a molecular phylogeny of the Microcotylidae Taschenberg, 1879

During a study of the monogeneans of four sparid fishes (Diplodus vulgaris, Pagellus bogaraveo, Pagrus pagrus and Sparus aurata) from the Western Mediterranean off Algeria, a large collection of parasites was characterised molecularly (28S rRNA and cox1 genes). A total of 46 partial sequences (23 for each gene) were generated from 38 isolates of monogeneans which included four species (Atrispinum acarne, Microcotyle erythrini (sensu stricto), Sparicotyle chrysophrii and Prostatomicrocotylinae gen. sp.) of the family Microcotylidae, two putative species of the family Capsalidae (Encotyllabe spp.), and one species (Choricotyle chrysophryi) of the family Diclidophoridae. Our study provides (i) the first molecular data for the Capsalidae in the Mediterranean; (ii) the first record of a member of the Prostatomicrocotylinae in the Mediterranean and in a sparid fish (D. vulgaris); (iii) the first cox1 sequences for A. acarne; (iv) the second record of M. erythrini (s.s.) from P. pagrus; and (v) the second confirmed by molecular data record of S. chrysophrii in wild populations of S. aurata. The first phylogenetic hypotheses for the family Microcotylidae developed here, revealed the monophyly of the subfamily Prostatomicrocotylinae and the genus Microcotyle but the relationships among the subfamilies were still largely unresolved with the best represented subfamily Microcotylinae being polyphyletic. Our results highlight the importance of molecular methods in the assessment of monogenean diversity and the need for a thorough taxon-sampling approach to increase the accuracy of phylogenetic reconstruction of the relationships of the large and taxonomically complex polyopisthocotylean family Microcotylidae.

Heteromicrocotyla polyorchis Unnithan, 1961 (Monogenea: Heteromicrocotylidae), a gill parasite of the yellow-spotted trevally, Carangoides fulvoguttatus (Carangidae) from Saudi Arabia: Morphology and phylogeny

Little information for parasitic infections of Carangoides fulvoguttatus was recorded. The present study was intended to investigate the gill parasite Heteromicrocotyla polyorchis of this fish and to provide a full morphological description and clarify its taxonomic status through phylogenetic analysis of the 28S rRNA gene region. A total of sixty fish specimens have been collected from the studied area (the Red Sea in Jeddah Province, Saudi Arabia) and gills were isolated and examined for identification of parasites. Using light electron microscopy, the recovered monogenean parasite's morphology was exhaustively characterized and described. Microscope examinations found that this parasite species represent Heteromicrocotyla polyorchis, and it could be distinguished from congeners of the same genus by armed genital atrium and cirrus sac, follicular post-ovarian testes, unique shape and number of clamps on both haptor sides, and the dorsally curved tip of the male copulatory organ. Morphological features were combined with molecular analysis of the 28S rRNA gene region. The selected gene for the isolated Heteromicrocotyla species was analyzed using appropriate primers to assist in phylogeny with those in the GenBank database. The present monogenean species was characterized by unique genetic sequences that were analyzed and deposited in the GenBank for the first time under the accession number MW406473. Phylogenetic analyses reported that the maximum identity between the current Heteromicrocotyla species and taxa of Heteromicrocotylidae was between 91.42 and 92.09% and confirmed its taxonomic status in this family with a well-distinct clade. The present study supported the second report of H. polyorchis as carangid fish ectoparasites and investigates the first appearance in C. fulvoguttatus inhabiting Saudi Arabia.

Life-cycle complexity in helminths: What are the benefits?

Parasitic worms (i.e., helminths) commonly infect multiple hosts in succession. With every transmission step, they risk not infecting the next host and thus dying before reproducing. Given this risk, what are the benefits of complex life cycles? Using a dataset for 973 species of trophically transmitted acanthocephalans, cestodes, and nematodes, we tested whether hosts at the start of a life cycle increase transmission and whether hosts at the end of a life cycle enable growth to larger, more fecund sizes. Helminths with longer life cycles, that is, more successive hosts, infected conspicuously smaller first hosts, slightly larger final hosts, and exploited trophic links with lower predator-prey mass ratios. Smaller first hosts likely facilitate transmission because of their higher abundance and because parasite propagules were the size of their normal food. Bigger definitive hosts likely increase fecundity because parasites grew larger in big hosts, particularly endotherms. Helminths with long life cycles attained larger adult sizes through later maturation, not faster growth. Our results indicate that complex helminth life cycles are ubiquitous because growth and reproduction are highest in large, endothermic hosts that are typically only accessible via small intermediate hosts, that is, the best hosts for growth and transmission are not the same.

The complete mitochondrial genome of the Capsaloides cristatus (Platyhelminthes, Monogenea), a pathogen of the sailfish (Istiophorus platypterus)

This study reports the complete mitochondrial genome of the Capsaloides cristatus (Monogenea: Capsalidae) collected from the gill lamella of Istiophorus platypterus. The total length of the mitogenome was 13,948 bp, containing 12 typical platyhelminthic protein-coding genes, 22 tRNA genes, 2 rRNA genes and a putative non-coding region, with the atp8 gene being absent. The total A + T content was 65.99%, which was significantly higher than that of the C + G content (34.01%). There were two kinds of start codons (ATG and GTG) and three kinds of terminated codons (TAA, TAG and TGA) in the 12 protein-coding genes. Phylogentic analysis revealed close relationships among the genera Capsaloides, Capsala, Benedenia and Neobenedenia with high bootstrap value supported. This study will provide useful molecular data for a better understanding of the species identification and phylogenetic position of C. cristatus.

Ultrastructural patterns of the excretory ducts of basal neodermatan groups (Platyhelminthes) and new protonephridial characters of basal cestodes

BACKGROUND

The flatworms (Lophotrochozoa: Platyhelminthes) are one of the major phyla of invertebrates but their interrelationships are still not well understood including unravelling the most closely related taxon of the Neodermata, which includes exclusively obligate parasites of all main groups of vertebrates with some 60,000 estimated species. Recent phylogenomic studies indicate that the freshwater 'microturbellarian' Bothrioplana semperi may be the closest ancestor to the Neodermata, but this hypothesis receives little morphological support. Therefore, additional morphological and ultrastructural characters that might help understand interrelations within the Neodermata are needed.

METHODS

Ultrastructure of the excretory ducts of representatives of the most basal parasitic flatworms (Neodermata), namely monocotylid (Monopisthocotylea) and chimaericolid (Polyopisthocotylea) monogeneans, aspidogastreans (Trematoda), as well as gyrocotylidean and amphilinidean tapeworms (Cestoda), were studied using transmission electron microscopy (TEM).

RESULTS

The present study revealed the same pattern of the cytoarchitecture of excretory ducts in all studied species of the basal neodermatans. This pattern is characterised by the presence of septate junctions between the adjacent epithelial cells and lateral ciliary flames along different levels of the excretory ducts. Additionally, a new character was observed in the protonephridial terminal cell of Gyrocotyle urna, namely a septate junction between terminal and adjacent duct cells at the level of the distal extremity of the flame tuft. In Amphilina foliacea, a new type of protonephridial cell with multiple flame bulbs and unique character of its weir, which consists of a single row of the ribs, is described. A remarkable difference has been observed between the structure of the luminal surface of the excretory ducts of the studied basal neodermatan groups and B. semperi.

CONCLUSIONS

The present study does not provide ultrastructural support for a close relationship between the Neodermata and B. semperi.

Population genetics and comparative mitogenomic analyses reveal cryptic diversity of Amphioctopus neglectus (Cephalopoda: Octopodidae)

This study presented 96 cox1 and 76 cox3 genes of Amphioctopus neglectus populations. Three distinct lineages were formed from phylogenetic trees and networks constructed using haplotypes. Mitogenomes of A. neglectus-a and A. neglectus-b as the representatives of two lineages separated from population genetics were sequenced to compare with A. neglectus at the genome-level. Amphioctopus neglectus-a showed significant differences with A. neglectus, mainly reflected in gene length, intergenic regions and the secondary structure of tandem repeat motifs. Notably, two sequence deletions in mitogenomes of the two representative species were detected in different positions of major non-coding regions, which were the most distinct differences with A. neglectus. Pairwise genetic distances and the phylogenetic analysis supported the relationship of (A. neglectus-a + (A. neglectus + A. neglectus-b)). This study suggested that A. neglectus-a should be considered as a potential cryptic species of this complex, while A. neglectus-b needed further verification to be defined.

The first mitochondrial genomes of endosymbiotic rhabdocoels illustrate evolutionary relaxation of atp8 and genome plasticity in flatworms

The first three mitochondrial (mt) genomes of endosymbiotic turbellarian flatworms are characterised for the rhabdocoels Graffilla buccinicola, Syndesmis echinorum and S. kurakaikina. Interspecific comparison of the three newly obtained sequences and the only previously characterised rhabdocoel, the free-living species Bothromesostoma personatum, reveals high mt genomic variability, including numerous rearrangements. The first intrageneric comparison within rhabdocoels shows that gene order is not fully conserved even between congeneric species. Atp8, until recently assumed absent in flatworms, was putatively annotated in two sequences. Selection pressure was tested in a phylogenetic framework and is shown to be significantly relaxed in this and another protein-coding gene: cox1. If present, atp8 appears highly derived in platyhelminths and its functionality needs to be addressed in future research. Our findings for the first time allude to a large degree of undiscovered (mt) genomic plasticity in rhabdocoels. It merits further attention whether this variation is correlated with a symbiotic lifestyle. Our results illustrate that this phenomenon is widespread in flatworms as a whole and not exclusive to the better-studied neodermatans.

Mitochondrial genomes and 28S rDNA contradict the proposed obsoletion of the order Tetraonchidea (Platyhelminthes: Monogenea)

Due to the incongruence of morphology-based hypotheses and scarcity of molecular data, validity of the order Tetraonchidea remains contentious. The only complete mitogenome currently available for the entire order is that of Paratetraonchoides inermis (Tetraonchoididae). To study the phylogeny of Tetraonchidea from mitogenomic perspective, we sequenced the first mitogenome for the family Tetraonchidae: Tetraonchus monenteron (Tetraonchidea). To get a nuclear-data perspective, we also sequenced nuclear 28S rDNA gene of both species. The mitogenome of T. monenteron does not have high A + T content, nor tRNA pseudo-genes, both of which were unique features reported in P. inermis. However, T. monenteron exhibits a unique gene order, with a large number of tRNA rearrangements in comparison to P. inermis and other monogeneans. Phylogenetic analyses conducted using Bayesian inference and maximum likelihood methods, complemented with partitioning, consistently support the sister-group relationship of T. monenteron (Tetraonchidae) and P. inermis (Tetraonchoididae). This is also partially supported by the 28S rDNA data and two morphologic apomorphies. This close relationship of Tetraonchidae and Tetraonchoididae challenges the latest major morphology-based classification, which proposed obsoletion of the Tetraonchidea order, and grouped Tetraonchoididae into the Gyrodactylidea clade. The validity of this order shall have to be further confirmed with more data.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Evidence for Adaptive Selection in the Mitogenome of a Mesoparasitic Monogenean Flatworm Enterogyrus malmbergi

Whereas a majority of monogenean flatworms are ectoparasitic, i.e., parasitize on external surfaces (mainly gills) of their fish hosts, Enterogyrus species (subfamily Ancyrocephalinae) are mesoparasitic, i.e., parasitize in the stomach of the host. As there are numerous drastic differences between these two environments (including lower oxygen availability), we hypothesized that this life-history innovation might have produced adaptive pressures on the energy metabolism, which is partially encoded by the mitochondrial genome (OXPHOS). To test this hypothesis, we sequenced mitochondrial genomes of two Ancyrocephalinae species: mesoparasitic E. malmbergi and ectoparasitic Ancyrocephalus mogurndae. The mitogenomic architecture of E. malmbergi is mostly standard for monogeneans, but that of A. mogurndae exhibits some unique features: missing trnL2 gene, very low AT content (60%), a non-canonical start codon of the nad2 gene, and exceptionally long tandem-repeats in the non-coding region (253 bp). Phylogenetic analyses produced paraphyletic Ancyrocephalinae (with embedded Dactylogyrinae), but with low support values. Selective pressure (PAML and HYPHY) and protein structure analyses all found evidence for adaptive evolution in cox2 and cox3 genes of the mesoparasitic E. malmbergi. These findings tentatively support our hypothesis of adaptive evolution driven by life-history innovations in the mitogenome of this species. However, as only one stomach-inhabiting mesoparasitic monogenean was available for this analysis, our findings should be corroborated on a larger number of mesoparasitic monogeneans and by physiological studies.

The complete mitochondrial genome and description of a new cryptic species of Benedenia Diesing, 1858 (Monogenea: Capsalidae), a major pathogen infecting the yellowtail kingfish Seriola lalandi Valenciennes in the South-East Pacific

Background

The monogenean Benedenia seriolae parasitizes fishes belonging to the genus Seriola, represents a species complex, and causes substantial impact on fish welfare in aquaculture systems worldwide. This study reports, for the first time, the complete mitochondrial genome of B. humboldti n. sp., a new cryptic species from the South-East Pacific (SEP).

Methods

The mitogenome of B. humboldti n. sp. was assembled from short Illumina 150 bp pair-end reads. The phylogenetic position of B. humboldti n. sp. among other closely related congeneric and confamiliar capsalids was examined using mitochondrial protein-coding genes (PCGs). Morphology of B. humboldti n. sp. was examined based on fixed and stained specimens.

Results

The AT-rich mitochondrial genome of B. humboldti is 13,455 bp in length and comprises 12 PCGs (atp8 was absent as in other monogenean genomes), 2 ribosomal RNA genes, and 22 transfer RNA genes. All protein-coding, ribosomal RNA, and transfer RNA genes are encoded on the H-strand. The gene order observed in the mitochondrial genome of B. humboldti n. sp. was identical to that of B. seriolae from Japan but different from that of B. seriolae from Australia. The genetic distance between B. humboldti n. sp. and B. seriolae from Japan was high. Minor but reliable differences in the shape of the penis were observed between Benedenia humboldti n. sp. and congeneric species.

Conclusions

Phylogenetic analyses based on PCGs in association with differences in the shape of the penis permitted us to conclude that the material from the South-East Pacific represents a new species of Benedenia infecting S. lalandi off the coast of Chile. The discovery of this parasite represents the first step to improving our understanding of infestation dynamics and to develop control strategies for this pathogen infecting the farmed yellowtail kingfish, Seriola lalandi, in the South-East Pacific.

Mitochondrial Genomes of Two Thaparocleidus Species (Platyhelminthes: Monogenea) Reveal the First rRNA Gene Rearrangement among the Neodermata

Phylogenetic framework for the closely related Ancylodiscoidinae and Ancyrocephalinae subfamilies remains contentious. As this issue was never studied using a large molecular marker, we sequenced the first two Ancylodiscoidinae mitogenomes: Thaparocleidus asoti and Thaparocleidus varicus. Both mitogenomes had two non-coding regions (NCRs) that contained a number of repetitive hairpin-forming elements (RHE). Due to these, the mitogenome of T. asoti (16,074 bp) is the longest among the Monogenea; especially large is its major NCR, with 3500 bp, approximately 1500 bp of which could not be sequenced (thus, the total mitogenome size is ≈ 17,600 bp). Although RHEs have been identified in other monopisthocotyleans, they appear to be independently derived in different taxa. The presence of RHEs may have contributed to the high gene order rearrangement rate observed in the two mitogenomes, including the first report of a transposition of rRNA genes within the Neodermata. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrinae embedded within the Ancyrocephalinae (paraphyly), whereas Ancylodiscoidinae formed a sister-group with them. This was also supported by the gene order analysis. 28S rDNA dataset produced polyphyletic Dactylogyridae and Ancyrocephalinae. The phylogeny of the two subfamilies shall have to be further evaluated with more data.

Homoplasy or plesiomorphy? Reconstruction of the evolutionary history of mitochondrial gene order rearrangements in the subphylum Neodermata

Recent mitogenomic studies have exposed a gene order (GO) shared by two classes, four orders and 31 species ('common GO') within the flatworm subphylum Neodermata. There are two possible hypotheses for this phenomenon: convergent evolution (homoplasy) or shared ancestry (plesiomorphy). To test those, we conducted a meta-analysis on all available mitogenomes to infer the evolutionary history of GO in Neodermata. To improve the resolution, we added a newly sequenced mitogenome that exhibited the common GO, Euryhaliotrema johni (Ancyrocephalinae), to the dataset. Phylogenetic analyses conducted on two datasets (nucleotides of all 36 genes and amino acid sequences of 12 protein coding genes) and four algorithms (MrBayes, RAxML, IQ-TREE and PhyloBayes) produced topology instability towards the tips, so ancestral GO reconstructions were conducted using TreeREx and MLGO programs using all eight obtained topologies, plus three unique topologies from previous studies. The results consistently supported the second hypothesis, resolving the common GO as a plesiomorphic ancestral GO for Neodermata, Cestoda, Monopisthocotylea, Cestoda + Trematoda and Cestoda + Trematoda + Monopisthocotylea. This allowed us to trace the evolutionary GO scenarios from each common ancestor to its descendants amongst the Monogenea and Cestoda classes, and propose that the common GO was most likely retained throughout all of the common ancestors, leading to the extant species possessing the common GO. Neodermatan phylogeny inferred from GOs was largely incongruent with all 11 topologies described above, but it did support the mitogenomic dataset in resolving Polyopisthocotylea as the earliest neodermatan branch. Although highly derived GOs might be of some use in resolving isolated taxonomic and phylogenetic uncertainties, we conclude that, due to the discontinuous nature of their evolution, they tend to produce artefactual phylogenetic relationships, which makes them unsuitable for phylogenetic reconstruction in Neodermata. Wider and denser sampling of neodermatan mitogenomic sequences will be needed to infer the evolutionary pathways leading to the observed diversity of GOs with confidence.

Mitochondrial genomes of two diplectanids (Platyhelminthes: Monogenea) expose paraphyly of the order Dactylogyridea and extensive tRNA gene rearrangements

BACKGROUND

Recent mitochondrial phylogenomics studies have reported a sister-group relationship of the orders Capsalidea and Dactylogyridea, which is inconsistent with previous morphology- and molecular-based phylogenies. As Dactylogyridea mitochondrial genomes (mitogenomes) are currently represented by only one family, to improve the phylogenetic resolution, we sequenced and characterized two dactylogyridean parasites, Lamellodiscus spari and Lepidotrema longipenis, belonging to a non-represented family Diplectanidae.

RESULTS

The L. longipenis mitogenome (15,433 bp) contains the standard 36 flatworm mitochondrial genes (atp8 is absent), whereas we failed to detect trnS1, trnC and trnG in L. spari (14,614 bp). Both mitogenomes exhibit unique gene orders (among the Monogenea), with a number of tRNA rearrangements. Both long non-coding regions contain a number of different (partially overlapping) repeat sequences. Intriguingly, these include putative tRNA pseudogenes in a tandem array (17 trnV pseudogenes in L. longipenis, 13 trnY pseudogenes in L. spari). Combined nucleotide diversity, non-synonymous/synonymous substitutions ratio and average sequence identity analyses consistently showed that nad2, nad5 and nad4 were the most variable PCGs, whereas cox1, cox2 and cytb were the most conserved. Phylogenomic analysis showed that the newly sequenced species of the family Diplectanidae formed a sister-group with the Dactylogyridae + Capsalidae clade. Thus Dactylogyridea (represented by the Diplectanidae and Dactylogyridae) was rendered paraphyletic (with high statistical support) by the nested Capsalidea (represented by the Capsalidae) clade.

CONCLUSIONS

Our results show that nad2, nad5 and nad4 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in the Diplectanidae. The unique gene order pattern further suggests discontinuous evolution of mitogenomic gene order arrangement in the Class Monogenea. This first report of paraphyly of the Dactylogyridea highlights the need to generate more molecular data for monogenean parasites, in order to be able to clarify their relationships using large datasets, as single-gene markers appear to provide a phylogenetic resolution which is too low for the task.

Three new Diplozoidae mitogenomes expose unusual compositional biases within the Monogenea class: implications for phylogenetic studies

BACKGROUND

As the topologies produced by previous molecular and morphological studies were contradictory and unstable (polytomy), evolutionary relationships within the Diplozoidae family and the Monogenea class (controversial relationships among the Discocotylinea, Microcotylinea and Gastrocotylinea suborders) remain unresolved. Complete mitogenomes carry a relatively large amount of information, sufficient to provide a much higher phylogenetic resolution than traditionally used morphological traits and/or single molecular markers. However, their implementation is hampered by the scarcity of available monogenean mitogenomes. Therefore, we sequenced and characterized mitogenomes belonging to three Diplozoidae family species, and conducted comparative genomic and phylogenomic analyses for the entire Monogenea class.

RESULTS

Taxonomic identification was inconclusive, so two of the species were identified merely to the genus level. The complete mitogenomes of Sindiplozoon sp. and Eudiplozoon sp. are 14,334 bp and 15,239 bp in size, respectively. Paradiplozoon opsariichthydis (15,385 bp) is incomplete: an approximately 2000 bp-long gap within a non-coding region could not be sequenced. Each genome contains the standard 36 genes (atp8 is missing). G + T content and the degree of GC- and AT-skews of these three mitogenome (and their individual elements) were higher than in other monogeneans. nad2, atp6 and nad6 were the most variable PCGs, whereas cox1, nad1 and cytb were the most conserved. Mitochondrial phylogenomics analysis, conducted using concatenated amino acid sequences of all PCGs, indicates that evolutionary relationships of the three genera are: (Eudiplozoon, (Paradiplozoon, Sindiplozoon)); and of the three suborders: (Discocotylinea, (Microcotylinea, Gastrocotylinea)). These intergeneric relationships were also supported by the skewness and principal component analyses.

CONCLUSIONS

Our results show that nad2, atp6 and nad6 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in Diplozoidae. Nucleotide bias and codon and amino acid usage patterns of the three diplozoid mitogenomes are more similar to cestodes and trematodes than to other monogenean flatworms. This unusual mutational bias was reflected in disproportionately long branches in the phylogram. Our study offsets the scarcity of molecular data for the subclass Polyopisthocotylea to some extent, and might provide important new insights into the evolutionary history of the three genera and three suborders.

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.

Spermiogenesis and spermatozoon ultrastructure in basal polyopisthocotylean monogeneans, Hexabothriidae and Chimaericolidae, and their significance for the phylogeny of the Monogenea

Sperm ultrastructure provides morphological characters useful for understanding phylogeny; no study was available for two basal branches of the Polyopisthocotylea, the Chimaericolidea and Diclybothriidea. We describe here spermiogenesis and sperm in Chimaericola leptogaster (Chimaericolidae) and Rajonchocotyle emarginata (Hexabothriidae), and sperm in Callorhynchocotyle callorhynchi (Hexabothriidae). Spermiogenesis in C. leptogaster and R. emarginata shows the usual pattern of most Polyopisthocotylea with typical zones of differentiation and proximo-distal fusion of the flagella. In all three species, the structure of the spermatozoon is biflagellate, with two incorporated trepaxonematan 9 + "1" axonemes and a posterior nucleus. However, unexpected structures were also seen. An alleged synapomorphy of the Polyopisthocotylea is the presence of a continuous row of longitudinal microtubules in the nuclear region. The sperm of C. leptogaster has a posterior part with a single axoneme, and the part with the nucleus is devoid of the continuous row of microtubules. The spermatozoon of R. emarginata has an anterior region with membrane ornamentation, and posterior lateral microtubules are absent. The spermatozoon of C. callorhynchi has transverse sections with only dorsal and ventral microtubules, and its posterior part shows flat sections containing a single axoneme and the nucleus. These findings have important implications for phylogeny and for the definition of synapomorphies in the Neodermata. We point out a series of discrepancies between actual data and interpretation of character states in the matrix of a phylogeny of the Monogenea. Our main conclusion is that the synapomorphy "lateral microtubules in the principal region of the spermatozoon" does not define the Polyopisthocotylea but is restricted to the Mazocraeidea.

Sequencing of the complete mitochondrial genome of a fish-parasitic flatworm Paratetraonchoides inermis (Platyhelminthes: Monogenea): tRNA gene arrangement reshuffling and implications for phylogeny

BACKGROUND

Paratetraonchoides inermis (Monogenea: Tetraonchoididae) is a flatworm parasitising the gills of uranoscopid fishes. Its morphological characteristics are ambiguous, and molecular data have never been used to study its phylogenetic relationships, which makes its taxonomic classification controversial. Also, several decades of unsuccessful attempts to resolve the relationships within the Monogenea present a strong indication that morphological datasets may not be robust enough to be used to infer evolutionary histories. As the use of molecular data is currently severely limited by their scarcity, we have sequenced and characterized the complete mitochondrial (mt) genome of P. inermis. To investigate its phylogenetic position, we performed phylogenetic analyses using Bayesian inference and maximum likelihood approaches using concatenated amino acid sequences of all 12 protein-coding genes on a dataset containing all available monogenean mt genomes.

RESULTS

The circular mt genome of P. inermis (14,654 bp) contains the standard 36 genes: 22 tRNAs, two rRNAs, 12 protein-encoding genes (PCGs; Atp8 is missing) and a major non-coding region (mNCR). All genes are transcribed from the same strand. The A + T content of the whole genome (82.6%), as well as its elements, is the highest reported among the monogeneans thus far. Three tRNA-like cloverleaf structures were found in mNCR. Several results of the phylogenomic analysis are in disagreement with previously proposed relationships: instead of being closely related to the Gyrodactylidea, Tetraonchidea exhibit a phylogenetic affinity with the Dactylogyridea + Capsalidea clade; and the order Capsalidea is neither basal within the subclass Monopisthocotylea, nor groups with the Gyrodactylidea, but instead forms a sister clade with the Dactylogyridea. The mt genome of P. inermis exhibits a unique gene order, with an extensive reorganization of tRNAs. Monogenea exhibit exceptional gene order plasticity within the Neodermata.

CONCLUSIONS

This study shows that gene order within monopisthocotylid mt genomes is evolving at uneven rates, which creates misleading evolutionary signals. Furthermore, our results indicate that all previous attempts to resolve the evolutionary history of the Monogenea may have produced at least partially erroneous relationships. This further corroborates the necessity to generate more molecular data for this group of parasitic animals.

Sequencing, characterization and phylogenomics of the complete mitochondrial genome of Dactylogyrus lamellatus (Monogenea: Dactylogyridae)

Despite the worldwide distribution and pathogenicity of monogenean parasites belonging to the largest helminth genus, Dactylogyrus, there are no complete Dactylogyrinae (subfamily) mitogenomes published to date. In order to fill this knowledge gap, we have sequenced and characterized the complete mitogenome of Dactylogyrus lamellatus, a common parasite on the gills of grass carp (Ctenopharyngodon idella). The circular mitogenome is 15,187 bp in size, containing the standard 22 tRNA genes, 2 rRNA genes, 12 protein-encoding genes and a long non-coding region (NCR). There are two highly repetitive regions in the NCR. We have used concatenated nucleotide sequences of all 36 genes to perform the phylogenetic analysis using Bayesian inference and maximum likelihood approaches. As expected, the two dactylogyrids, D. lamellatus (Dactylogyrinae) and Tetrancistrum nebulosi (Ancyrocephalinae), were closely related to each other. These two formed a sister group with Capsalidae, and this cluster finally formed a further sister group with Gyrodactylidae. Phylogenetic affinity between Dactylogyrinae and Ancyrocephalinae was further confirmed by the similarity in their gene arrangement. The sequencing of the first Dactylogyrinae, along with a more suitable selection of outgroups, has enabled us to infer a much better phylogenetic resolution than recent mitogenomic studies. However, as many lineages of the class Monogenea remain underrepresented or not represented at all, a much larger number of mitogenome sequences will have to be available in order to infer the evolutionary relationships among the monogeneans fully, and with certainty.

Autonomy and integration in complex parasite life cycles

Complex life cycles are common in free-living and parasitic organisms alike. The adaptive decoupling hypothesis postulates that separate life cycle stages have a degree of developmental and genetic autonomy, allowing them to be independently optimized for dissimilar, competing tasks. That is, complex life cycles evolved to facilitate functional specialization. Here, I review the connections between the different stages in parasite life cycles. I first examine evolutionary connections between life stages, such as the genetic coupling of parasite performance in consecutive hosts, the interspecific correlations between traits expressed in different hosts, and the developmental and functional obstacles to stage loss. Then, I evaluate how environmental factors link life stages through carryover effects, where stressful larval conditions impact parasites even after transmission to a new host. There is evidence for both autonomy and integration across stages, so the relevant question becomes how integrated are parasite life cycles and through what mechanisms? By highlighting how genetics, development, selection and the environment can lead to interdependencies among successive life stages, I wish to promote a holistic approach to studying complex life cycle parasites and emphasize that what happens in one stage is potentially highly relevant for later stages.

The complete mitochondrial genome of Gyrodactylus gurleyi (Platyhelminthes: Monogenea)

Gyrodactylus gurleyi, was inhabited on the fins and gills of goldfish (Carassius auratus), which belonged to the family Gyrodactylidae. In this study, we sequenced the complete mitochondrial genome of G. gurleyi with the total length of 14 771 bp. The mitogenome contained 12 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and two major non-coding regions (NC1 and NC2). The overall AT content was 72.1%. In phylogenetic analysis, G. gurleyi and G. kobayashii clustered together and then united with the clade of other three Gyrodactylus species (G. salaris, G. thymalli and G. derjavinoides) with high nodal support.

The mitochondrial genome of the egg-laying flatworm Aglaiogyrodactylus forficulatus (Platyhelminthes: Monogenoidea)

BACKGROUND

The rather species-poor oviparous gyrodactylids are restricted to South America. It was suggested that they have a basal position within the otherwise viviparous Gyrodactylidae. Accordingly, it was proposed that the species-rich viviparous gyrodactylids diversified and dispersed from there.

METHODS

The mitochondrial genome of Aglaiogyrodactylus forficulatus was bioinformatically assembled from next-generation illumina MiSeq sequencing reads, annotated, and compared to previously published mitochondrial genomes of other monogenoidean flatworm species.

RESULTS

The mitochondrial genome of A. forficulatus consists of 14,371 bp with an average A + T content of 75.12 %. All expected 12 protein coding, 22 tRNA, and 2 rRNA genes were identified. Furthermore, there were two repetitive non-coding regions essentially consisting of 88 bp and 233 bp repeats, respectively. Maximum Likelihood analyses placed the mitochondrial genome of A. forficulatus in a well-supported clade together with the viviparous Gyrodactylidae species. The gene order differs in comparison to that of other monogenoidean species, with rearrangements mainly affecting tRNA genes. In comparison to Paragyrodactylus variegatus, four gene order rearrangements, i.e. three transpositions and one complex tandem-duplication-random-loss event, were detected.

CONCLUSION

Mitochondrial genome sequence analyses support a basal position of the oviparous A. forficulatus within Gyrodactylidae, and a sister group relationship of the oviparous and viviparous forms.

Comparative genomics of flatworms (platyhelminthes) reveals shared genomic features of ecto- and endoparastic neodermata

The ectoparasitic Monogenea comprise a major part of the obligate parasitic flatworm diversity. Although genomic adaptations to parasitism have been studied in the endoparasitic tapeworms (Cestoda) and flukes (Trematoda), no representative of the Monogenea has been investigated yet. We present the high-quality draft genome of Gyrodactylus salaris, an economically important monogenean ectoparasite of wild Atlantic salmon (Salmo salar). A total of 15,488 gene models were identified, of which 7,102 were functionally annotated. The controversial phylogenetic relationships within the obligate parasitic Neodermata were resolved in a phylogenomic analysis using 1,719 gene models (alignment length of >500,000 amino acids) for a set of 16 metazoan taxa. The Monogenea were found basal to the Cestoda and Trematoda, which implies ectoparasitism being plesiomorphic within the Neodermata and strongly supports a common origin of complex life cycles. Comparative analysis of seven parasitic flatworm genomes identified shared genomic features for the ecto- and endoparasitic lineages, such as a substantial reduction of the core bilaterian gene complement, including the homeodomain-containing genes, and a loss of the piwi and vasa genes, which are considered essential for animal development. Furthermore, the shared loss of functional fatty acid biosynthesis pathways and the absence of peroxisomes, the latter organelles presumed ubiquitous in eukaryotes except for parasitic protozoans, were inferred. The draft genome of G. salaris opens for future in-depth analyses of pathogenicity and host specificity of poorly characterized G. salaris strains, and will enhance studies addressing the genomics of host–parasite interactions and speciation in the highly diverse monogenean flatworms.

The complete mitochondrial genome of Gyrodactylus kobayashii (Platyhelminthes: Monogenea)

The complete mitochondrial genome of Gyrodactylus kobayashii was 14 786 bp in length, containing 12 protein-coding genes (lacking Atp8), 22 tRNA genes, two rRNA genes and two major non-coding regions (NC1 and NC2). The overall A + T content of mitochondrial genome was 71.6%. A close relationship between G. kobayashii and the three Gyrodactylus species (G. salaris, G. thymalli and G. derjavinoides) was uncovered in the phylogenetic tree based on amino acid sequences.

Molecular characterization of Urocleidoides cuiabai and U. malabaricusi (Monogenea: Dactylogyridae) from the trahira fish Hoplias aff. malabaricus in the Paraná River, Brazil

Urocleidoides ectoparasites are mainly found on fish of the neotropical regions. Although molecular research on monogeneans is available, no genetic data exist characterizing species in the Urocleidoides genus. Some DNA sequences have been efficacious in systematic studies and in the reconstruction of phylogenies of fish parasites. Relevant roles have been given to the sequence of the mitochondrial gene of cytochrome c oxidase I (COI). This study characterized COI sequences of the parasites Urocleidoides malabaricusi and U. cuiabai in trahira fish Hoplias aff. malabaricus of the flood plain of the Upper River Paraná, Brazil. The two species under analysis were distinguished by sequencing and analysing a 420-bp fragment of the COI gene, which suggested the existence of the cryptic species U. malabaricusi.

Fast evolutionary rates associated with functional loss in class I glucose transporters of Schistosoma mansoni

Background

The trematode parasite, Schistosoma mansoni, has evolved to switch from oxidative phosphorylation to glycolysis in the presence of glucose immediately after invading the human host. This metabolic switch is dependent on extracellular glucose concentration. Four glucose transporters are encoded in the genome of S. mansoni, however, only two were shown to facilitate glucose diffusion.

Results

By modeling the phase of human host infection, we showed that transporter transcript expression profiles of recently transformed schistosomula have two opposing responses to increased glucose concentrations. Concurring with the transcription profiles, our phylogenetic analyses revealed that S. mansoni glucose transporters belong to two separate clusters, one associated with class I glucose transporters from vertebrates and insects, and the other specific to parasitic Platyhelminthes. To study the evolutionary paths of both groups and their functional implications, we determined evolutionary rates, relative divergence times, genomic organization and performed structural analyses with the protein sequences. We finally used the modelled structures of the S. mansoni glucose transporters to biophysically (i) analyze the dynamics of key residues during glucose binding, (ii) test glucose stability within the active site, and (iii) demonstrate glucose diffusion. The two S. mansoni Platyhelminthes-specific glucose transporters, which seem to be younger than the other two, exhibit slower rates of molecular evolution, are encoded by intron-poor genes, and transport glucose. Interestingly, our molecular dynamic analyses suggest that S. mansoni class I glucose transporters are not able to transport glucose.

Conclusions

The glucose transporter family in S. mansoni exhibit different evolutionary histories. Our results suggested that S. mansoni class I glucose transporters lost their capacity to transport glucose and that this function evolved independently in the Platyhelminthes-specific glucose transporters. Finally, taking into account the differences in the dynamics of glucose transport of the Platyhelminthes-specific transporters of S. mansoni compared to that of humans, we conclude that S. mansoni glucose transporters may be targets for rationally designed drugs against schistosomiasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2144-6) contains supplementary material, which is available to authorized users.

Constraining the Deep Origin of Parasitic Flatworms and Host-Interactions with Fossil Evidence

Novel fossil discoveries have contributed to our understanding of the evolutionary appearance of parasitism in flatworms. Furthermore, genetic analyses with greater coverage have shifted our views on the coevolution of parasitic flatworms and their hosts. The putative record of parasitic flatworms is consistent with extant host associations and so can be used to put constraints on the evolutionary origin of the parasites themselves. The future lies in new molecular clock analyses combined with additional discoveries of exceptionally preserved flatworms associated with hosts and coprolites. Besides direct evidence, the host fossil record and biogeography have the potential to constrain their evolutionary history, albeit with caution needed to avoid circularity, and a need for calibrations to be implemented in the most conservative way. This might result in imprecise, but accurate divergence estimates for the evolution of parasitic flatworms.

A transcriptomic-phylogenomic analysis of the evolutionary relationships of flatworms

The interrelationships of the flatworms (phylum Platyhelminthes) are poorly resolved despite decades of morphological and molecular phylogenetic studies [1, 2]. The earliest-branching clades (Catenulida, Macrostomorpha, and Polycladida) share spiral cleavage and entolecithal eggs with other lophotrochozoans. Lecithoepitheliata have primitive spiral cleavage but derived ectolecithal eggs. Other orders (Rhabdocoela, Proseriata, Tricladida and relatives, and Bothrioplanida) all have derived ectolecithal eggs but have uncertain affinities to one another. The orders of parasitic Neodermata emerge from an uncertain position from within these ectolecithal classes. To tackle these problems, we have sequenced transcriptomes from 18 flatworms and 5 other metazoan groups. The addition of published data produces an alignment of >107,000 amino acids with less than 28% missing data from 27 flatworm taxa in 11 orders covering all major clades. Our phylogenetic analyses show that Platyhelminthes consist of the two clades Catenulida and Rhabditophora. Within Rhabditophora, we show the earliest-emerging branch is Macrostomorpha, not Polycladida. We show Lecithoepitheliata are not members of Neoophora but are sister group of Polycladida, implying independent origins of the ectolecithal eggs found in Lecithoepitheliata and Neoophora. We resolve Rhabdocoela as the most basally branching euneoophoran taxon. Tricladida, Bothrioplanida, and Neodermata constitute a group that appears to have lost both spiral cleavage and centrosomes. We identify Bothrioplanida as the long-sought closest free-living sister group of the parasitic Neodermata. Among parasitic orders, we show that Cestoda are closer to Trematoda than to Monogenea, rejecting the concept of the Cercomeromorpha. Our results have important implications for understanding the evolution of this major phylum.

•

Phylogenomics provide insights into the interrelationships of Platyhelminthes

•

Macrostomorpha are the basalmost rhabditophorans

•

Polycladida are sister group of Lecithoepitheliata/Prorhynchida

•

Bothrioplanida are the free-living sister group of Neodermata

Nuclear genomic signals of the 'microturbellarian' roots of platyhelminth evolutionary innovation

Flatworms number among the most diverse invertebrate phyla and represent the most biomedically significant branch of the major bilaterian clade Spiralia, but to date, deep evolutionary relationships within this group have been studied using only a single locus (the rRNA operon), leaving the origins of many key clades unclear. In this study, using a survey of genomes and transcriptomes representing all free-living flatworm orders, we provide resolution of platyhelminth interrelationships based on hundreds of nuclear protein-coding genes, exploring phylogenetic signal through concatenation as well as recently developed consensus approaches. These analyses robustly support a modern hypothesis of flatworm phylogeny, one which emphasizes the primacy of the often-overlooked 'microturbellarian' groups in understanding the major evolutionary transitions within Platyhelminthes: perhaps most notably, we propose a novel scenario for the interrelationships between free-living and vertebrate-parasitic flatworms, providing new opportunities to shed light on the origins and biological consequences of parasitism in these iconic invertebrates.

Evolution of complex life cycles in trophically transmitted helminths. I. Host incorporation and trophic ascent

Links between parasites and food webs are evolutionarily ancient but dynamic: life history theory provides insights into helminth complex life cycle origins. Most adult helminths benefit by sexual reproduction in vertebrates, often high up food chains, but direct infection is commonly constrained by a trophic vacuum between free-living propagules and definitive hosts. Intermediate hosts fill this vacuum, facilitating transmission to definitive hosts. The central question concerns why sexual reproduction, and sometimes even larval growth, is suppressed in intermediate hosts, favouring growth arrest at larval maturity in intermediate hosts and reproductive suppression until transmission to definitive hosts? Increased longevity and higher growth in definitive hosts can generate selection for larger parasite body size and higher fecundity at sexual maturity. Life cycle length is increased by two evolutionary mechanisms, upward and downward incorporation, allowing simple (one-host) cycles to become complex (multihost). In downward incorporation, an intermediate host is added below the definitive host: models suggest that downward incorporation probably evolves only after ecological or evolutionary perturbations create a trophic vacuum. In upward incorporation, a new definitive host is added above the original definitive host, which subsequently becomes an intermediate host, again maintained by the trophic vacuum: theory suggests that this is plausible even under constant ecological/evolutionary conditions. The final cycle is similar irrespective of its origin (upward or downward). Insights about host incorporation are best gained by linking comparative phylogenetic analyses (describing evolutionary history) with evolutionary models (examining selective forces). Ascent of host trophic levels and evolution of optimal host taxa ranges are discussed.

Phylogeography of Tetrancistrum nebulosi (Monogenea, Dactylogyridae) on the host of mottled spinefoot (Siganus fuscescens) in the South China Sea, inferred from mitochondrial COI and ND2 genes

To examine the phylogeographical pattern of Tetrancistrum nebulosi (Monogenea, Dactylogyridae) in the South China Sea, fragments of mitochondrial cytochrome c oxidase subunit I and NADH dehydrogenase subunit 2 genes were obtained for 220 individuals collected from 8 localities along the southeast coast of China and 1 locality in Terengganu, Malaysia. Based on these two genes, two and three distinct clades with geographic signals were revealed on the phylogenetic trees respectively. The divergence between these clades was estimated to occur in the late Pleistocene. Analysis of molecular variance and pairwise F_ST suggested a high rate of gene flow among individuals sampled from the Chinese coast, but with obvious genetic differentiation from the Malaysian population. Mismatch distribution and neutrality tests indicated that the T. nebulosi population experienced expansion in Pleistocene low sea level periods. Vicariance was considered to account for the genetic divergence between Chinese and Malaysian populations, while sea level fluctuations and mainland-island connections during glacial cycles were associated with the slight genetic divergence between the populations along the mainland coast of China and those off Sanya. On the contrary, oceanographic circulations and host migration could lead to genetic homogeneity of populations distributed along the mainland coast of China.

Fauna europaea: helminths (animal parasitic)

Fauna Europaea provides a public web-service with an index of scientific names (including important synonyms) of all living European land and freshwater animals, their geographical distribution at country level (up to the Urals, excluding the Caucasus region), and some additional information. The Fauna Europaea project covers about 230,000 taxonomic names, including 130,000 accepted species and 14,000 accepted subspecies, which is much more than the originally projected number of 100,000 species. This represents a huge effort by more than 400 contributing specialists throughout Europe and is a unique (standard) reference suitable for many users in science, government, industry, nature conservation and education. Helminths parasitic in animals represent a large assemblage of worms, representing three phyla, with more than 200 families and almost 4,000 species of parasites from all major vertebrate and many invertebrate groups. A general introduction is given for each of the major groups of parasitic worms, i.e. the Acanthocephala, Monogenea, Trematoda (Aspidogastrea and Digenea), Cestoda and Nematoda. Basic information for each group includes its size, host-range, distribution, morphological features, life-cycle, classification, identification and recent key-works. Tabulations include a complete list of families dealt with, the number of species in each and the name of the specialist responsible for data acquisition, a list of additional specialists who helped with particular groups, and a list of higher taxa dealt with down to the family level. A compilation of useful references is appended.

The complete mitochondrial genome of Neobenedenia melleni (Platyhelminthes: Monogenea): mitochondrial gene content, arrangement and composition compared with two Benedenia species

The complete mitochondrial (mt) genome sequences of Neobenedenia melleni were determined and compared with those of Benedenia seriolae and B. hoshinai. This circular genome comprises 13,270 bp and includes all 36 typical mt genes found in flatworms. Total AT content of N. melleni is 75.9 %. ATG is the most common start codon, while nad4L is initiated by GTG. All protein-coding genes are predicted to terminate with TAG and TAA. N. melleni has the trnR with a TCG anticodon, which is the same to B. seriolae but different from B. hoshinai (ACG). The mt gene arrangement of N. melleni is similar to that of B. seriolae and B. hoshinai with the exception of three translocations (trnF, trnT and trnG). The overlapped region between nad4L and nad4 was found in the N. melleni mt genome, which was also reported for the published Gyrodactylus species, but it was not found in those of B. seriolae and B. hoshinai, which are non-coding regions instead. The present study provides useful molecular characters for species or strain identification and systematic studies of this parasite.

Transcriptome data reveal Syndermatan relationships and suggest the evolution of endoparasitism in Acanthocephala via an epizoic stage

The taxon Syndermata comprises the biologically interesting wheel animals ("Rotifera": Bdelloidea + Monogononta + Seisonidea) and thorny-headed worms (Acanthocephala), and is central for testing superordinate phylogenetic hypotheses (Platyzoa, Gnathifera) in the metazoan tree of life. Recent analyses of syndermatan phylogeny suggested paraphyly of Eurotatoria (free-living bdelloids and monogononts) with respect to endoparasitic acanthocephalans. Data of epizoic seisonids, however, were absent, which may have affected the branching order within the syndermatan clade. Moreover, the position of Seisonidea within Syndermata should help in understanding the evolution of acanthocephalan endoparasitism. Here, we report the first phylogenomic analysis that includes all four higher-ranked groups of Syndermata. The analyzed data sets comprise new transcriptome data for Seison spec. (Seisonidea), Brachionus manjavacas (Monogononta), Adineta vaga (Bdelloidea), and Paratenuisentis ambiguus (Acanthocephala). Maximum likelihood and Bayesian trees for a total of 19 metazoan species were reconstructed from up to 410 functionally diverse proteins. The results unanimously place Monogononta basally within Syndermata, and Bdelloidea appear as the sister group to a clade comprising epizoic Seisonidea and endoparasitic Acanthocephala. Our results support monophyly of Syndermata, Hemirotifera (Bdelloidea + Seisonidea + Acanthocephala), and Pararotatoria (Seisonidea + Acanthocephala), rejecting monophyly of traditional Rotifera and Eurotatoria. This serves as an indication that early acanthocephalans lived epizoically or as ectoparasites on arthropods, before their complex lifecycle with arthropod intermediate and vertebrate definite hosts evolved.

The complete mitochondrial genome of Tetrancistrum nebulosi (Monogenea: Ancyrocephalidae)

The first complete mitochondrial (mt) genome of Ancyrocephalidae is reported herein. The mt genome of Tetrancistrum nebulosi was 13,392 bp in length containing 12 protein-coding genes (lacking atp8), 22 tRNA genes and 2 rRNA genes. The longest non-coding region was located between nad5 and trnG, and the A + T content was 72.4%. All tRNAs had the typical clover-leaf secondary structure except for trnS1((AGN)), trnR, trnF and trnQ. The rrnL and rrnS subunits were separated by trnC, as documented in the monopisthocotylean groups (Benedenia and Gyrodactylus species), while they were adjacent to each other in the polyopisthocotylean species (Microcotyle sebastis, Polylabris halichoeres and Pseudochauhanea macrorchis).

Problematic barcoding in flatworms: A case-study on monogeneans and rhabdocoels (Platyhelminthes)

Some taxonomic groups are less amenable to mitochondrial DNA barcoding than others. Due to the paucity of molecular information of understudied groups and the huge molecular diversity within flatworms, primer design has been hampered. Indeed, all attempts to develop universal flatworm-specific COI markers have failed so far. We demonstrate how high molecular variability and contamination problems limit the possibilities for barcoding using standard COI-based protocols in flatworms. As a consequence, molecular identification methods often rely on other widely applicable markers. In the case of Monogenea, a very diverse group of platyhelminth parasites, and Rhabdocoela, representing one-fourth of all free-living flatworm taxa, this has led to a relatively high availability of nuclear ITS and 18S/28S rDNA sequences on GenBank. In a comparison of the effectiveness in species assignment we conclude that mitochondrial and nuclear ribosomal markers perform equally well. In case intraspecific information is needed, rDNA sequences can guide the selection of the appropriate (i.e. taxon-specific) COI primers if available.

The monogenean which lost its clamps

Ectoparasites face a daily challenge: to remain attached to their hosts. Polyopisthocotylean monogeneans usually attach to the surface of fish gills using highly specialized structures, the sclerotized clamps. In the original description of the protomicrocotylid species Lethacotyle fijiensis, described 60 years ago, the clamps were considered to be absent but few specimens were available and this observation was later questioned. In addition, genera within the family Protomicrocotylidae have either clamps of the “gastrocotylid” or the “microcotylid” types; this puzzled systematists because these clamp types are characteristic of distinct, major groups. Discovery of another, new, species of the genus Lethacotyle, has allowed us to explore the nature of the attachment structures in protomicrocotylids. Lethacotyle vera n. sp. is described from the gills of the carangid Caranx papuensis off New Caledonia. It is distinguished from Lethacotyle fijiensis, the only other species of the genus, by the length of the male copulatory spines. Sequences of 28S rDNA were used to build a tree, in which Lethacotyle vera grouped with other protomicrocotylids. The identity of the host fish was confirmed with COI barcodes. We observed that protomicrocotylids have specialized structures associated with their attachment organ, such as lateral flaps and transverse striations, which are not known in other monogeneans. We thus hypothesized that the clamps in protomicrocotylids were sequentially lost during evolution, coinciding with the development of other attachment structures. To test the hypothesis, we calculated the surfaces of clamps and body in 120 species of gastrocotylinean monogeneans, based on published descriptions. The ratio of clamp surface: body surface was the lowest in protomicrocotylids. We conclude that clamps in protomicrocotylids are vestigial organs, and that occurrence of “gastrocotylid” and simpler “microcotylid” clamps within the same family are steps in an evolutionary sequence, leading to the absence of these attributes in species of Lethacotyle.