The mitochondrial genome of Gyrodactylus salaris (Platyhelminthes: Monogenea), a pathogen of Atlantic salmon (Salmo salar)

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.

Sequencing and Reconstructing Helminth Mitochondrial Genomes Directly from Genomic Next-Generation Sequencing Data

We present a detailed method for extraction of high-molecular weight genomic DNA suitable for numerous DNA sequencing applications, and a straightforward in silico approach for reconstructing novel mitochondrial (mt) genomes directly from total genomic DNA extracts derived from next-generation sequencing (NGS) data sets. The in silico post-sequencing pipeline described is fast, accurate, and highly efficient, with modest memory requirements that can be performed using a standard desktop computer. The approach is particularly effective for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information currently available and overcomes many of the limitations of traditional strategies. The described methodologies are also applicable for metagenomics sequencing from mixed or pooled samples containing multiple species and subsequent specific assembly of specific mitochondrial genomes.

Eudiplozoon nipponicum (Monogenea, Diplozoidae) and its adaptation to haematophagy as revealed by transcriptome and secretome profiling

Background

Ectoparasites from the family Diplozoidae (Platyhelminthes, Monogenea) belong to obligate haematophagous helminths of cyprinid fish. Current knowledge of these worms is for the most part limited to their morphological, phylogenetic, and population features. Information concerning the biochemical and molecular nature of physiological processes involved in host–parasite interaction, such as evasion of the immune system and its regulation, digestion of macromolecules, suppression of blood coagulation and inflammation, and effect on host tissue and physiology, is lacking. In this study, we report for the first time a comprehensive transcriptomic/secretome description of expressed genes and proteins secreted by the adult stage of Eudiplozoon nipponicum (Goto, 1891) Khotenovsky, 1985, an obligate sanguivorous monogenean which parasitises the gills of the common carp (Cyprinus carpio).

Results

RNA-seq raw reads (324,941 Roche 454 and 149,697,864 Illumina) were generated, de novo assembled, and filtered into 37,062 protein-coding transcripts. For 19,644 (53.0%) of them, we determined their sequential homologues. In silico functional analysis of E. nipponicum RNA-seq data revealed numerous transcripts, pathways, and GO terms responsible for immunomodulation (inhibitors of proteolytic enzymes, CD59-like proteins, fatty acid binding proteins), feeding (proteolytic enzymes cathepsins B, D, L1, and L3), and development (fructose 1,6-bisphosphatase, ferritin, and annexin). LC-MS/MS spectrometry analysis identified 721 proteins secreted by E. nipponicum with predominantly immunomodulatory and anti-inflammatory functions (peptidyl-prolyl cis-trans isomerase, homolog to SmKK7, tetraspanin) and ability to digest host macromolecules (cathepsins B, D, L1).

Conclusions

In this study, we integrated two high-throughput sequencing techniques, mass spectrometry analysis, and comprehensive bioinformatics approach in order to arrive at the first comprehensive description of monogenean transcriptome and secretome. Exploration of E. nipponicum transcriptome-related nucleotide sequences and translated and secreted proteins offer a better understanding of molecular biology and biochemistry of these, often neglected, organisms. It enabled us to report the essential physiological pathways and protein molecules involved in their interactions with the fish hosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07589-z.

Species delimitation of Gyrodactylus (Monogenea: Gyrodactylidae) infecting the southernmost cyprinids (Actinopterygii: Cyprinidae) in the New World

The genus Gyrodactylus von Nordmann, 1832 represents one of the most diverse and widespread taxa within Monogenea, with approximately 500 species described worldwide. Thirty-three species of Gyrodactylus have been recorded in Mexico, and in the last two decades, at least 26 new species have been described mainly from freshwater fish families such as poeciliids, goodeids, profundulids, characids, and cichlids. In this study, we describe two new species of Gyrodactylus infecting freshwater cyprinids based on morphological and molecular characteristics. Gyrodactylus ticuchi n. sp. and Gyrodactylus tobala n. sp. were recovered from Notropis moralesi de Buen and N. imeldae Cortés, respectively, captured in five localities from the State of Oaxaca, Mexico. The new species differ slightly from their congeners in the morphology of the haptoral hard parts and the male copulatory organ. Sequences of the Internal Transcribed Spacers rDNA (ITS1-5.8S-ITS2), cytochrome oxidase subunit I (cox1), and the D2 + D3 domains of the large subunit (28S rDNA) were obtained from multiple specimens and analyzed using Maximum Likelihood (ML) and Bayesian Inference (BI). Phylogenetic hypotheses using ITS rDNA, cox1, and 28S rDNA genes recovered two new species of Gyrodactylus from N. moralesi and N. imeldae; we briefly discuss their phylogenetic relationship with other congeners. These gyrodactylids represent the first species described in species of Notropis from southern Mexico, the cyprinids exhibiting the southernmost distribution in the New World.

Near-complete phylogeny of extant Crocodylia (Reptilia) using mitogenome-based data

Species of the order Crocodylia are mostly large, predatory and semi-aquatic reptiles. Crocodylia, the closest living relatives of birds, first appeared in the Late Cretaceous period. In the present study, the complete mitochondrial (mt) genomes of 19 Crocodylia species, including two species (Melanosuchus niger and Caiman yacare) that have not been previously sequenced for mitogenomes, were processed through Illumina sequencing to offer genetic resources and compare with the mitogenomes of Crocodylia species reported previously. In addition, a high-resolution phylogenetic tree of nearly all current recognized species of Crocodylia is constructed based on mitogenomic data. Phylogenetic analyses support monophyly of three families: Alligatoridae (four genera: Alligator, Caiman, Melanosuchus and Paleosuchus), Crocodylidae (three genera: Crocodylus, Mecistops and Osteolaemus) and Gavialidae (two genera: Gavialis and Tomistoma). The tree topology is generally similar to previous studies. Molecular dating suggests that the first split within Crocodylia date back to the Upper Cretaceous (approx. 86.75 Mya). The estimated time to the most recent common ancestor (TMRCA) of Alligatoridae is 53.33 Mya and that of Crocodylidae and Gavialidae is 50.13 Mya, which might be closely linked to climate changes during the Late Palaeocene and Early Eocene. Additionally, this study proves that the diversification rate within Crocodylia began to increase from the Late Eocene (about 36 Mya) and two diversification peak periods of Crocodylia (0–10 Mya and 10–20 Mya) are disclosed, which is roughly consistent with the estimated crocodylian species richness through time. Combining all these clues, we can suggest that climate fluctuation may have played a decisive role in the speciation of Crocodylia.

Laser capture microdissection in combination with mass spectrometry: Approach to characterization of tissue-specific proteomes of Eudiplozoon nipponicum (Monogenea, Polyopisthocotylea)

Eudiplozoon nipponicum (Goto, 1891) is a hematophagous monogenean ectoparasite which inhabits the gills of the common carp (Cyprinus carpio). Heavy infestation can lead to anemia and in conjunction with secondary bacterial infections cause poor health and eventual death of the host. This study is based on an innovative approach to protein localization which has never been used in parasitology before. Using laser capture microdissection, we dissected particular areas of the parasite body without contaminating the samples by surrounding tissue and in combination with analysis by mass spectrometry obtained tissue-specific proteomes of tegument, intestine, and parenchyma of our model organism, E. nipponicum. We successfully verified the presence of certain functional proteins (e.g. cathepsin L) in tissues where their presence was expected (intestine) and confirmed that there were no traces of these proteins in other tissues (tegument and parenchyma). Additionally, we identified a total of 2,059 proteins, including 72 peptidases and 33 peptidase inhibitors. As expected, the greatest variety was found in the intestine and the lowest variety in the parenchyma. Our results are significant on two levels. Firstly, we demonstrated that one can localize all proteins in one analysis and without using laboratory animals (antibodies for immunolocalization of single proteins). Secondly, this study offers the first complex proteomic data on not only the E. nipponicum but within the whole class of Monogenea, which was from this point of view until recently neglected.

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.

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.

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.

Revision of Gyrodactylus salaris phylogeny inspired by new evidence for Eemian crossing between lineages living on grayling in Baltic and White sea basins

In this research, grayling-specific Gyrodactylus salarisMalmberg, 1957 isolates from Baltic Sea basin were collected in Sweden for the first time. Samples were obtained in three drainage systems: Kalixälven (River Kaitum), Ljungan (River Sölvbacka strömmar), and Umeälven (River Juktån). Three molecular markers were analysed: nuclear ITS rDNA (Internal Transcribed Spacer) and ADNAM1 (Anonymous DNA Marker 1), and mitochondrial cox1 gene. As a result, four new mitochondrial haplotypes were identified (III-C1tt, III-C1tt_ht, IX-A1tt and X-A1tt). The ADNAM1 analyses resulted in revealing two new alleles (WS4 and BS9) and two new genotypes (T6 and T7). T7 seems to be an indicator of ancient crossing between Baltic and White Sea lineages of the parasite which happened during a first 3000-year period of Eemian interglacial about 130,000 years ago in the connection between Baltic and White Sea. Molecular clock estimates were adjusted, revealing the mean substitution rate and the divergence rate among branches of 3.6% (95% HPD: 2.2%–5.2%) and 7.2% per million years, respectively. As a result, cox1 phylogeny rooted with the introgressed haplotypes has been revised and altered in accordance to new data, revealing fourteen equidistant lineages five of which have been excluded from the study. Based on the new phylogenetic approach, including the molecular clock, this work suggests an overall revision of G. salaris phylogeny and attempts at precisely drawing the division of lineages within this polytypic species as well as proposes unification in nomenclature for its strains.

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.

Comparative analyses within Gyrodactylus (Platyhelminthes: Monogenea) mitochondrial genomes and conserved polymerase chain reaction primers for gyrodactylid mitochondrial DNA

In this study, we describe the complete mitochondrial genomes of Gyrodactylus brachymystacis and Gyrodactylus parvae infecting rainbow trout (Oncorhynchus mykiss) and the invasive topmouth gudgeon (Pseudorasbora parva), respectively. The two circular genomes have a common genome organization found in other Gyrodactylus species. Comparative analyses of mitochondrial genomes from six Gyrodactylus species were carried out to determine base composition, codon usage, transfer RNA and ribosomal RNA genes, major non-coding regions, and nucleotide diversity within the genus. We also provide the first universal models of the secondary structures of rrnS and rrnL for this group thereby promoting utilization of these genetic markers. Universal primers provided herein can be used to obtain more mitochondrial information for pathogen identification and may reveal different levels of molecular phylogenetic inferences for this lineage.

Ancient and modern genome shuffling: Reticulate mito-nuclear phylogeny of four related allopatric species of Gyrodactylus von Nordmann, 1832 (Monogenea: Gyrodactylidae), ectoparasites on the Eurasian minnow Phoxinus phoxinus (L.) (Cyprinidae)

Phylogenetic analyses including four allopatric species of Gyrodactylus von Nordmann, 1832 on the Eurasian minnow Phoxinus phoxinus (L.) (Cyprinidae) revealed incongruence between the nuclear ITS1-5.8S-ITS2 and mitochondrial cox1 phylogenies due to ancient hybridisation. Gyrodactylus pannonicus Molnár, 1968 was sampled close to its type-locality, the upper reaches of River Tisza, tributary of Danube in the Black Sea Basin. Faunistic search detected three new related species with maximum composite likelihood distances in cox1 between 16.8-23.2% (tentatively 1.3 to 1.8 My of divergence). Gyrodactylus albolacustris n. sp. recorded in the White Sea Basin, eastern Baltic Basin and Mongolia was close to G. pannonicus in the nuclear ITS (divergence of 0.9%), but diverged in cox1 by 19.8%. The Mongolian isolate of G. albolacustris n. sp. diverged from the European isolates in cox1 by 8.9%, suggesting 0.7 My of isolation. The two other new species differed from G. pannonicus by >4% in ITS and some large indels in ITS1, and by >20% in cox1. Gyrodactylus danastriae n. sp. was found in River Strwiąż, a tributary of the River Dniester (Black Sea Basin) and was characterised by smaller size of anchors and by 29-41 bp dimorphic insertion in ITS1. Gyrodactylus botnicus n. sp. is considered endemic in the Baltic Basin, but was also found in the White Sea Basin as a postglacial immigrant, where it had hybridised with G. albolacustris n. sp. in spite of the high divergence in ITS (3.9%) and cox1 (22%). The discordant nuclear and mitochondrial phylogenies revealed an ancient mitochondrial introgression: G. albolacustris n. sp. was derived from a hybridisation combining proto-pannonicus ITS with proto-danastriae mitochondria, perhaps 1.3 My ago. The postglacial hybridisation of G. albolacustris n. sp. (as the donor of mtDNA^alb and ITS^alb) and G. botnicus n. sp. (donor of the ITS^bot) offered a model of shuffling of the genomic components: the process of the homogenisation and stabilisation of nuclear ITS (concerted evolution) and the lineage sorting has hardly begun.

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.

The complete mitochondrial genome of Metorchis orientalis (Trematoda: Opisthorchiidae): Comparison with other closely related species and phylogenetic implications

Metorchis orientalis (Trematoda: Opisthorchiidae) is an important trematode infecting many animals and humans, causing metorchiasis. In the present study, the complete mitochondrial (mt) genome of M. orientalis was sequenced. The complete mt genome of M. orientalis is 13,834 bp circular DNA molecule and contains 12 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. The gene content and arrangement of M. orientalis is the same as those of Opisthorchiidae trematodes (Opisthorchis viverrini, Opisthorchis felineus and Clonorchis sinensis), but distinct from Schistosoma spp. Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes with three different computational algorithms (Bayesian inference, maximum likelihood and maximum parsimony) revealed that M. orientalis and O. viverrini represent sister taxa. The mt genome provides a novel genetic marker for further studies of the identification, classification and molecular epidemiology of Opisthorchiidae trematodes, and should have implications for the diagnosis, prevention and control of metorchiasis in animals and humans.

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.

Neotropical Monogenoidea. 58. Three new species of Gyrodactylus (Gyrodactylidae) from Scleromystax spp. (Callichthyidae) and the proposal of COII gene as an additional fragment for barcoding gyrodactylids

Based on molecular markers (COII and ITS1-ITS2) and morphological data, we describe three new Neotropical species of Gyrodactylus von Nordmann, 1832 from Scleromystax barbatus (Quoy et Gaimard) and Scleromystax macropterus (Regan) from southern Brazil. The three new species can be distinguished from each other by sequences of both molecular markers and morphology of hooks and anchors. Gyrodactylus bueni sp. n. is characterised by having hook with shaft curved, heel straight, shelf straight, toe pointed, anchor with superficial root slender, elongate and male copulatory organ armed with two rows of spinelets. Gyrodactylus major sp. n. presents hook with shaft, point curved, proximal shaft straight, heel convex, shelf convex, toe concave, anchor with superficial root robust and male copulatory organ armed with two rows of spinelets. Gyrodactylus scleromystaci sp. n. presents hook with shaft, point recurved, heel convex, shelf convex, toe pointed, anchor with superficial root curved and male copulatory organ armed with two rows of spinelets. These species appear to be closely related to other species of Gyrodactylus known from other species of Callichthyidae. These new species, however, differ by the comparative morphology of the haptoral hard structures and molecular data. Comparative analysis of sequences from these species of Gyrodactylus suggests that the COII gene may represent an important marker for the taxonomy of species of Gyrodactylidae and, perhaps, for species of other lineages of Monogenoidea.

The mitochondrial genome of Paragyrodactylus variegatus (Platyhelminthes: Monogenea): differences in major non-coding region and gene order compared to Gyrodactylus

BACKGROUND

Paragyrodactylus Gvosdev and Martechov, 1953, a viviparous genus of ectoparasite within the Gyrodactylidae, contains three nominal species all of which infect Asian river loaches. The group is suspected to be a basal lineage within Gyrodactylus Nordmann, 1832 sensu lato although this remains unclear. Further molecular study, beyond characterization of the standard Internal Transcribed Spacer region, is needed to clarify the evolutionary relationships within the family and the placement of this genus.

METHODS

The mitochondrial genome of Paragyrodactylus variegatus You, King, Ye and Cone, 2014 was amplified in six parts from a single worm, sequenced using primer walking, annotated and analyzed using bioinformatic tools.

RESULTS

The mitochondrial genome of P. variegatus is 14,517 bp, containing 12 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a major non-coding region (NCR). The overall A + T content of the mitochondrial genome is 76.3%, which is higher than all reported mitochondrial genomes of monogeneans. All of the 22 tRNAs have the typical cloverleaf secondary structure, except tRNACys, tRNASer1 and tRNASer2 that lack the dihydrouridine (DHU) arm. There are six domains (domain III is absent) and three domains in the inferred secondary structures of the large ribosomal subunit (rrnL) and small ribosomal subunit (rrnS), respectively. The NCR includes six 40 bp tandem repeat units and has the double identical poly-T stretches, stem-loop structure and some surrounding structure elements. The gene order (tRNAGln, tRNAMet and NCR) differs in arrangement compared to the mitochondrial genomes reported from Gyrodactylus spp.

CONCLUSION

The Duplication and Random Loss Model and Recombination Model together are the most plausible explanations for the variation in gene order. Both morphological characters and characteristics of the mitochondrial genome support Paragyrodactylus as a distinct genus from Gyrodactylus. Considering their specific distribution and known hosts, we believe that Paragyrodactylus is a relict freshwater lineage of viviparous monogenean isolated in the high plateaus of central Asia on closely related river loaches.

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.

MicroRNA loci support conspecificity of Gyrodactylus salaris and Gyrodactylus thymalli (Platyhelminthes: Monogenea)

The monogenean flatworm Gyrodactylus salaris is a serious threat to wild and farmed Atlantic salmon stocks in Norway. Morphologically, the closely related but harmless Gyrodactylus thymalli on grayling can hardly be distinguished from G. salaris. Until now, molecular approaches could not resolve unambiguously whether G. salaris and G. thymalli represent just one polytypic species, two polytypic species or a complex of more than two species. In the first known genome-wide analysis utilizing 37 conserved microRNA loci, the genetic differentiation of seven populations of G. salaris and G. thymalli was assessed. The concatenated alignment spanned 21,742bp including 62 variable positions. A neighbor-joining cluster analysis did not support any host-based or mitochondrial haplotype-based grouping of strains. We conclude that a two species concept for G. salaris and G. thymalli does not reflect meaningful biological entities. Instead, G. salaris and G. thymalli are just one species comprising several pathogenic and non-pathogenic strains on various primary hosts. Following the International Code for Zoological Nomenclature, G. salaris Malmberg, 1957 is the valid species name with G. thymalli Žitňan, 1960 becoming the junior synonym. Accordingly, the range of G. salaris is significantly increased, given that formerly G. salaris-free countries such as e.g., Great Britain are now within the species' natural range. The synonymization of G. salaris and G. thymalli implies severe challenges to current disease management routines, which assume that G. salaris and G. thymalli are readily distinguishable. Protocols for reliable identification of pathogenic and non-pathogenic strains of G. salaris need to be developed.

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.

Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads--a baiting and iterative mapping approach

We present an in silico approach for the reconstruction of complete mitochondrial genomes of non-model organisms directly from next-generation sequencing (NGS) data—mitochondrial baiting and iterative mapping (MITObim). The method is straightforward even if only (i) distantly related mitochondrial genomes or (ii) mitochondrial barcode sequences are available as starting-reference sequences or seeds, respectively. We demonstrate the efficiency of the approach in case studies using real NGS data sets of the two monogenean ectoparasites species Gyrodactylus thymalli and Gyrodactylus derjavinoides including their respective teleost hosts European grayling (Thymallus thymallus) and Rainbow trout (Oncorhynchus mykiss). MITObim appeared superior to existing tools in terms of accuracy, runtime and memory requirements and fully automatically recovered mitochondrial genomes exceeding 99.5% accuracy from total genomic DNA derived NGS data sets in <24 h using a standard desktop computer. The approach overcomes the limitations of traditional strategies for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information at hand and represents a fast and highly efficient in silico alternative to laborious conventional strategies relying on initial long-range PCR. We furthermore demonstrate the applicability of MITObim for metagenomic/pooled data sets using simulated data. MITObim is an easy to use tool even for biologists with modest bioinformatics experience. The software is made available as open source pipeline under the MIT license at https://github.com/chrishah/MITObim.

The complete mitochondrial genome of Dugesia japonica (Platyhelminthes; order Tricladida)

We used two sequencing methods, namely long polymerase chain reaction (PCR) and primer walking, to determine the complete mitochondrial DNA (mtDNA) sequence of Dugesia japonica and most of the mtDNA sequence of Dugesia ryukyuensis. The genome of D. japonica contained 36 genes including 12 of the 13 protein-coding genes characteristic of metazoan mitochondrial genomes, two ribosomal RNA genes, and 22 transfer RNA genes. The genome of D. ryukyuensis contained 33 genes, including 12 protein-coding genes, two ribosomal RNA genes, and 19 transfer RNA genes. The gene order of the mitochondrial genome from the Dugesia species showed no clear homology with either the Neodermata or other free-living Rhabditophora. This indicates that the platyhelminths exhibit great variability in mitochondrial gene order. This is the first complete sequence analysis of the mitochondrial genome of a free-living member of Rhabditophora, which will facilitate further studies on the population genetics and genomic evolution of the Platyhelminthes.

The complete mitochondrial genome of an ectoparasitic monopisthocotylean fluke Benedenia hoshinai (Monogenea: Platyhelminthes)

An exponential growth of mitochondrial genome information has brought significant progress in understanding the organismal phylogeny and mitochondrial genome evolution for many metazoans including platyhelminth groups. In this study, we determined the complete mitochondrial genome sequence for Benedenia hoshinai, an ectoparasitic monogenean species, and compared it with its congener Benedenia seriolae. The complete mitochondrial genome is 13,554 bp in length and contains 12 protein-coding genes (lacking the atp8 gene), 2 rRNA genes, and 22 tRNA genes, all encoded in the same direction as found in all other platyhelminth species sequenced to date. The gene arrangement of B. hoshinai mtDNA is almost identical to B. seriolae, differing only by the translocation of trnT between cox1 and rrnL. It is unclear whether the shared position of trnT between B. hoshinai and Gyrodactylus represents evidence for their phylogenetic affinity; testing this hypothesis requires further mitogenomic evidence.

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

The mitochondrial genome of Gyrodactylus derjavinoides (Platyhelminthes: Monogenea)--a mitogenomic approach for Gyrodactylus species and strain identification

Systematists and evolutionary biologists are constantly on the lookout for new sources of characters to discriminate amongst taxa and estimate interrelationships within and between taxa. Entire mitochondrial genomes provide a wealth of data, both at the nucleotide and amino acid level. Molecular markers are of particular utility when applied to small, morphologically conserved taxa, as is the case for many monogenean ectoparasites of fish. Gyrodactylus species display a considerable degree of anatomical conservatism, complicating diagnostics based solely on morphology, and some are significant pests of wild and cultured fish. Here we sequenced the complete mitochondrial genome of Gyrodactylus derjavinoides Malmberg, Collins, Cunningham & Behiar 2007, one of the most frequently found gyrodactylid species on salmonids in Scandinavia, and compared it with the recently published genomes of Gyrodactylus salaris Malmberg, 1957 and Gyrodactylus thymalli Zitnan 1960. Through comparative sliding window analysis we identified regions of high sequence variability and designed new primer sequences. In total, 6 new primer pairs have been developed, amplifying fragments of cox1, cox3, nad1, nad2, nad4, nad5 and atp6. Together, they amplify regions capturing almost half the nucleotide variability present in the complete mitochondrial genome. These degenerate primers should also work for other Gyrodactylus species parasitizing salmonids. In addition, we developed a multiplex assay that simultaneously amplifies four fragments in a single PCR reaction. Besides the diagnostic value, these fragments can be used for studying the transmission dynamics of Gyrodactylus, providing crucial information for an improved understanding of the spread and epidemiology of these important fish pathogens.

Hybrid origin of Baltic salmon-specific parasite Gyrodactylus salaris: a model for speciation by host switch for hemiclonal organisms

Host switching explains the high species number of ectoparasitic, viviparous, mainly parthenogenetic but potentially hermaphroditic flatworms of the genus Gyrodactylus. The starlike mitochondrial phylogeny of Gyrodactylus salaris suggested parallel divergence of several clades on grayling (also named as Gyrodactylus thymalli) and an embedded sister clade on Baltic salmon. The hypothesis that the parasite switched from grayling to salmon during the glacial diaspora was tested using a 493-bp nuclear DNA marker ADNAM1. The parasites on salmon in lakes Onega and Ladoga were heterozygous for divergent ADNAM1 alleles WS1 and BS1, found as nearly fixed in grayling parasites in the White Sea and Baltic Sea basins, respectively. In the Baltic salmon-specific mtDNA clade, the WS/BS heterozygosity was maintained in 23 out of the 24 local clones. The permanently heterozygous clade was endemic in the Baltic Sea basin, and it had accumulated variation in mtDNA (31 variable sites on 1600 bp) and in the alleles of the nuclear locus (two point mutations and three nucleotide conversions along 493 bp). Mendelian shuffling of the nuclear alleles between the local clones indicated rare sex within the clade, but the WS/BS heterozygosity was lost in only one salmon hatchery clone, which was heterozygous WS1/WS3. The Baltic salmon-specific G. salaris lineage was monophyletic, descending from a single historical hybridization and consequential host switch, frozen by permanent heterozygosity. A possible time for the hybridization of grayling parasite strains from the White Sea and Baltic Sea basins was during the Eemian interglacial 132 000 years bp. Strains having a separate divergent mtDNA observed on farmed rainbow trout, and on salmon in Russian lake Kuito were suggested to be clones derived from secondary and tertiary recombination events.

The biology of gyrodactylid monogeneans: the "Russian-doll killers"

This article reviews the history of gyrodactylid research focussing on the unique anatomy, behaviour, ecology and evolution of the viviparous forms while identifying gaps in our knowledge and directions for future research. We provide the first summary of research on the oviparous gyrodactylids from South American catfish, and highlight the plesiomorphic characters shared by gyrodactylids and other primitive monogeneans. Of these, the most important are the crawling, unciliated larva and the spike sensilla of the cephalic lobes. These characters allow gyrodactylids to transfer between hosts at any stage of the life cycle, without a specific transmission stage. We emphasise the importance of progenesis in shaping the evolution of the viviparous genera and discuss the relative extent of progenesis in the different genera. The validity of the familial classification is discussed and we conclude that the most significant division within the family is between the oviparous and the viviparous genera. The older divisions into Isancistrinae and Polyclithrinae should be allowed to lapse. We discuss approaches to the taxonomy of gyrodactylids, and we emphasise the importance of adequate morphological and molecular data in new descriptions. Host specificity patterns in gyrodactylids are discussed extensively and we note the importance of host shifts, revealed by molecular data, in the evolution of gyrodactylids. To date, the most closely related gyrodactylids have not been found on closely related hosts, demonstrating the importance of host shifts in their evolution. The most closely related species pair is that of G. salaris and G. thymalli, and we provide an account of the patterns of evolution taking place in different mitochondrial clades of this species complex. The host specificity of these clades is reviewed, demonstrating that, although each clade has its preferred host, there is a range of specificity to different salmonids, providing opportunities for complex patterns of survival and interbreeding in Scandinavia. At the same time, we identify trends in systematics and phylogeny relevant to the G. salaris epidemics on Atlantic salmon in Norway, which can be applied more generally to parasite epidemiology and evolution. Although much of gyrodactylid research in the last 30 years has been directed towards salmonid parasites, there is great potential in using other experimental systems, such as the gyrodactylids of poeciliids and sticklebacks. We also highlight the role of glacial lakes and modified river systems during the ice ages in gyrodactylid speciation, and suggest that salmon infecting clades of G. salaris first arose from G. thymalli in such lakes, but failed to spread fully across Scandinavia before further dispersal was ended by rising sea levels. This dispersal has been continued by human activity, leading to the appearance of G. salaris as a pathogen in Norway. We review the history and current status of the epidemic, and current strategies for elimination of the parasite from Norway. Finally, we consider opportunities for further spread of the parasite within and beyond Europe.

Mitochondrial DNA variation of a natural population of Gyrodactylus thymalli (Monogenea) from the type locality River Hnilec, Slovakia

The monogenean flatworm Gyrodactylus thymalli (Zitnan, Helminthologia, 2:266-269, 1960) is considered a harmless ectoparasite on grayling (Thymallus thymallus). The species is closely related to G. salaris Malmberg, 1957 that causes severe gyrodactylosis on Atlantic salmon (Salmo salar) in many Norwegian rivers. In this paper, we study the mitochondrial diversity of a G. thymalli population from one of the type localities Hrable on River Hnilec, Slovakia. By sequencing parts of the mitochondrial NADH dehydrogenase subunit 5 gene, we detected three haplotypes that differ from each other by 2.1-4.1%. The haplotype HnilecI was found most common. Our data suggest that River Hnilec has been colonized independently at least three times with G. thymalli.

The complete mitochondrial DNA sequence of the monogenean Gyrodactylus thymalli (Platyhelminthes: Monogenea), a parasite of grayling (Thymallus thymallus)

We present the complete mitochondrial (mt) genome of Gyrodactylus thymalli, a monogenean ectoparasite on grayling (Thymallus thymallus). The circular genome is 14788 bp in size and includes all 35 genes recognized from other flatworm mt genomes. The overall A+T content of the mt genome is 62.8%. Twenty regions of non-coding DNA ranging from 1 to 111 bp in length were identified in addition to 2 highly conserved large non-coding regions 799 bp and 767 bp in size. Compared to the recently described mt DNA of the closely related G. salaris from Atlantic salmon from Signaldalselva, Norway, the mitochondrial genome of G. thymalli from Hnilec, Slovakia, differs on average by 2.2%.