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

Morphological and molecular characterisation of Tristoma integrum Diesing, 1850 (Monogenea, Capsalidae), including its complete mitogenome

Capsalids are monopisthocotylean monogenean parasites found on the skin and gills of fish. Capsalines (subfamily Capsalinae) are large-sized capsalids, parasitic on highly prized gamefish, and species of Tristoma parasitise only the gills of swordfish (Xiphias gladius). We obtained specimens of Tristoma integrum Diesing, 1850 from swordfish collected off Algeria in the Mediterranean Sea. Here, we describe the specimens, including the key systematics characters of dorsolateral body sclerites. One specimen was used for a next generation sequencing analysis but a part of it, including the sclerites, was mounted on a permanent slide, drawn, and deposited in a curated collection. We characterised the complete mitogenome, the ribosomal cluster (including 18S and 28S) and additional genes such as Elongation factor 1 alpha (EF1α) and Histone 3. We also retrieved molecular information from the host tissue present in the gut of the monogenean and provide the sequence of the complete rRNA cluster of the host, X. gladius. The mitogenome of T. integrum is 13 968 bp in length and codes for 12 protein, 2 rRNA and 22 tRNA. Phylogenies of capsalids were generated from 28S sequences and concatenated mitochondrial protein-coding genes, respectively. In the 28S phylogeny, most subfamilies based on morphology were not found to be monophyletic, but the Capsalinae were monophyletic. In both phylogenies, the closest member to Tristoma spp. was a member of the Capsaloides. In an Appendix, we report the complex nomenclatural history of Tristoma Cuvier, 1817 and its species.

Complete Mitochondrial Genome of Two Ectoparasitic Capsalids (Platyhelminthes: Monogenea: Monopisthocotylea): Gene Content, Composition, and Rearrangement

The capsalid monogeneans are important pathogens that generally infect marine fishes and have a substantial impact on fish welfare in aquaculture systems worldwide. However, the current mitogenome information on capsalids has received little attention, limiting the understanding of their evolution and phylogenetic relationships with other monogeneans. This paper reports the complete mitochondrial genomes of Capsala katsuwoni and Capsala martinieri for the first time, which we obtained using a next-generation sequencing method. The mitogenomes of C. katsuwoni and C. martinieri are 13,265 and 13,984 bp in length, respectively. Both species contain the typical 12 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The genome compositions show a moderate A+T bias (66.5% and 63.9% for C. katsuwoni and C. martinieri, respectively) and exhibit a negative AT skew but a positive GC skew in both species. One gene block rearrangement was found in C. katsuwoni in comparison with other capsalid species. Instead of being basal to the Gyrodactylidea and Dactylogyridea or being clustered with Dactylogyridea, all species of Capsalidea are grouped into a monophyletic clade. Our results clarify the gene rearrangement process and evolutionary status of Capsalidae and lay a foundation for further phylogenetic studies of monogeneans.

Capsalids (Platyhelminthes: Monogenea) from marine fishes off Okinawa in Japan with the proposal of two new genera

Molecular studies of the Capsalidae suggested that the genus Benedenia is polyphyletic, but a taxonomic organization of the genus that reflects molecular data has not yet been proposed. As a result of molecular analysis (28S rDNA, ITS1-5.8S-ITS2, and cox1 data) including specimens of Benedeniinae newly obtained from Okinawa-jima Island in Japan, two new genera and the revival of Tareenia independent to the genus Benedenia are proposed. Gracilobenedenia n. gen. is distinguished from the other genera of Benedeniinae based on morphological characteristics. This new genus comprises 6 species: G. lutjani n. comb. (type species), G. anticavaginata n. comb., G. rohdei n. comb., G. beverleyburtonae n. comb., G. kuremibai n. gen., n. sp., and G. hichi n. gen., n. sp. Armatobenedenia n. gen. is a monotypic genus for A. armata n. comb. The present molecular phylogenetic analysis showed the independence of Tareenia, and it can be morphologically separated from the other benedeniines. Four species including two new species obtained from Okinawa-jima Island are reported: G. kuremibai n. sp., G. hichi n. sp., G. lutjani n. comb., and Metabenedeniella parva. Furthermore, in the species identification and phylogenetic analysis of capsalids, the usefulness of not only the 28S rDNA but also ITS and the cox1 regions was suggested. These genes were evaluated the efficacy of those regions in DNA barcoding, and the ITS and cox1 regions shown to be useful for DNA barcoding in capsalids compared to the 28S rDNA sequence.

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.

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.

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.

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 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.

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.

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 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.