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

Complete mitochondrial genomes of Dactylogyrus crucifer and Dactylogyrus zandti reveal distinct patterns of codon usage within Dactylogyrus

Monogeneans of the genus Dactylogyrus Diesing, 1850, the largest genus in the family Dactylogyridae, mostly parasitize the gills of cyprinoid hosts; however, only 3 Dactylogyrus' mitochondrial genomes (mitogenomes) are studied so far. The aim of this research is to extend our understanding of the mitogenomes of Dactylogyrus. We sequenced the mitogenomes of D. crucifer and D. zandti isolated from Rutilus rutilus and Abramis brama orientalis in northwest China, and then we compared these mitogenomes with other monogeneans. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of D. crucifer and D. zandti and characterized the mitogenomes to understand the gene structure, gene identity, the secondary structures of the 22 tRNA genes, and relative synonymous codon usage. We used the analytic Bayesian Information and Maximum Likelihood methods to determine their associated phylogenetic trees. The mitogenomes of D. crucifer and D. zandti were 14,403 and 18,584 bp, respectively. Organization and positioning of these genes were in accordance with Dactylogyrus lamellatus and Dactylogyrus tuba. The nucleotide composition of Dactylogyridae was different from other families of Monogenea, and the A+T count of genus Dactylogyrus (54 - 58.4 %) was lower than other genus species of the family Dactylogyridea (63.9 - 78.4 %) in protein-coding genes. Dactylogyrus members displayed a codon usage bias. The relative synonymous codon used by Dactylogyrus was not conserved and was lower than other monogeneans. The codon use patterns of closely-related species isolated from closely-related hosts were identical. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrus isolated from host subfamily Leuciscinae formed a sister-group. Our results contributed significantly to an increased database of mitogenomes, more than 50 %, for Dactylogyrus that may help future studies of mitochondrial genes and codon uses for the analysis of monogenean phylogenetics.

The complete mitochondrial genome of Gyrodactylus pseudorasborae (Platyhelminthes: Monogenea) with a phylogeny of Gyrodactylidae parasites

BACKGROUND

Gyrodactylus von Nordmann, 1832, a genus of viviparous parasites within the family Gyrodactylidae, contains one of the largest nominal species in the world. Gyrodactylus pseudorasborae Ondračková, Seifertová & Tkachenko, 2023 widely distributed in Europe and China, although its mitochondrial genome remains unclear. This study aims to sequence the mitogenome of G.pseudorasborae and clarify its phylogenetic relationship within the Gyrodactylidea. The mitochondrial genome of G. pseudorasborae was amplified in six parts from a single parasite, sequenced using primer walking, annotated and analyzed using bioinformatic tools.

RESULTS

The mitochondrial genome of G. pseudorasborae is 14,189 bp in length, containing 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and two major non-coding regions (NCR: NC1 and NC2). The overall A + T content of the mitogenome is 73.1%, a medium content compared with all reported mitochondrial genomes of monogeneans. The mitogenome of G. pseudorasborae presents a clear bias in nucleotide composition with a negative AT skew and a positive GC skew except for NCR. All tRNAs have the typical cloverleaf secondary structure except for tRNACys, tRNASer1, and tRNASer2, which lack the dihydrouridine (DHU) arm. Furthermore, one repetitive non-coding region of 32 bp repeats occurred in the NC1 region with poly-T stretch, stem-loop structure, and TAn motif. The gene order is identical to the mitochondrial genomes reported from other Gyrodactylus species except Gyrodactylus nyanzae Paperna, 1973 and Gyrodactylus sp. FZ-2021. Phylogenetic analyses show that G. pseudorasborae and Gyrodactylus parvae You, Easy & Cone, 2008 cluster together with high nodal support based on 12 PCGs sequences and amino acid sequences, Gyrodactylidae forms independent monophyletic clade within Gyrodactylidea.

CONCLUSION

Both the mitochondrial genome and phylogenetic analyses support G. pseudorasborae is a member of the genus Gyrodactylus and Gyrodactylidae forms an independent monophyletic clade within Gyrodactylidea. Furthermore, the mitochondrial genome of G. pseudorasborae is the shortest in the Gyrodactylidea species compared with size differences in NCR.

First report of the complete mitochondrial genome of Carpomya pardalina (Bigot) (Diptera: Tephritidae) and phylogenetic relationships with other Tephritidae

Carpomya pardalina is known for its potential invasiveness, which poses a significant and alarming threat to Cucurbitaceae crops. It is considered a highly perilous pest species that requires immediate attention for quarantine and prevention. Due to the challenges in distinguishing pests of the Tephritidae family based on morphological characteristics, it is imperative to elucidate the mitochondrial genomic information of C. pardalina. In this study, the mitochondrial genome sequence of C. pardalina was determined and analyzed using next-generation sequencing. The results revealed that the mitogenome sequence had a total length of 16,257 bp, representing a typical circular molecule. It consisted of 13 PCGs, two rRNA genes, 22 tRNA genes and a non-coding region. The structure and organization of the mitochondrial genome of C. pardalina were found to be typical and similar to the published homologous sequences of other fruit flies in the Tephritidae family. Phylogenetic analysis confirmed that C. pardalina belongs to the Carpomya genus, which is consistent with traditional morphological taxonomy. Additionally, Carpomya and Rhagoletis were identified as sister groups. This study presents the first report of the complete mitochondrial genome of C. pardalina, which can serve as a valuable resource for future investigations in species diagnosis, evolutionary biology, prevention and control measures.

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.

The complete mitogenome of Paratanakia chii (Cypriniformes; Cyprinidae)

Paratanakia chii is a bitterling fish of the genus Paratanakia, subfamily Acheilognathinae and family Cyprinidae. The mitochondrial DNA sequence of P. chii is reported in this paper. The complete mitochondrial genome of P. chii is 16,575 bp in length, including 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and 1 displacement loop (D-loop). The genome sequence is consistent with those of most other carp. The majority of PCGs have AT- (Met) start codons and TA- end codons. The A + T contents of the genome, PCGs, transfer RNAs (tRNAs), and ribosomal RNAs (rRNAs) are 56.92%, 58.07%, 56.34%, and 54.21%, respectively. Phylogenetic analysis showed that P. chii is most closely related to Tanankia himantegus. These data will benefit relative ecological and phylogenetic studies.

Characterization of the complete mitochondrial genome of Nippotaenia mogurndae Yamaguti and Miyata, 1940 (Cestoda: Nippotaeniidae)

In this study, we report the first complete mitochondrial genome of the tapeworm Nippotaenia mogurndae in the order Nippotaeniidea Yamaguti, 1939. This mitogenome, which is 14,307 base pairs (bp) long with an A + T content of 72.2%, consists of 12 protein-coding genes, 22 transfer RNA (tRNA) genes, two rRNA genes, and two non-coding regions. Most tRNAs have a conventional cloverleaf structure, but trnS1 and trnR lack dihydrouridine arms of tRNA. The two largest non-coding regions, NCR1 (220 bp) and NCR2 (817 bp), are located between trnY and trnS2 and between nad5 and trnG, respectively. Phylogenetic analyses of mitogenomic data indicate that N. mogurndae is closely related to tapeworms in the order Cyclophyllidea.

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.

The Complete Mitochondrial Genome of Pennella sp. Parasitizing Thunnus albacares

In the study, the parasite from the yellowfin tuna (Thunnus albacares) was separated, and morphological observation and molecular identification were carried out. Our results showed that the parasite was similar to Pennella sp. Its cephalothorax was covered by spherical to spherical non-branched nipples of almost the same size, which were very similar in shape and arrangement. A pair of slightly larger, the unbranched antenna was present on the outer margin of the small papillae-covered area. The gene sequence of COX1 with a length of 1,558 bp in the mitochondria of the parasite was 100% similar to Pennella sp. (MZ934363). The mitochondrial genome had a total length of 14,620 bp. It consisted of 36 genes (12 protein-coding, 22 transfer RNAs and 2 ribosomal RNAs) and a dummy control region, but the mitochondrial genome had no ATP8 gene. Morphological observation showed that Pennella sp. was dark red, with a convex cephalothorax, with a total length of 8.42 cm, parasitic on the dorsal side of yellowfin tuna. Pennella sp. included the cephalothorax, neck, trunk, abdomen and egg belt. This study was the first report on the mitochondrial genome of Pennella sp. The results provide basic data for further identifying the parasites of Pennella genus.

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.

Mitochondrial genomes of two Polydora (Spionidae) species provide further evidence that mitochondrial architecture in the Sedentaria (Annelida) is not conserved

Contrary to the early evidence, which indicated that the mitochondrial architecture in one of the two major annelida clades, Sedentaria, is relatively conserved, a handful of relatively recent studies found evidence that some species exhibit elevated rates of mitochondrial architecture evolution. We sequenced complete mitogenomes belonging to two congeneric shell-boring Spionidae species that cause considerable economic losses in the commercial marine mollusk aquaculture: Polydora brevipalpa and Polydora websteri. The two mitogenomes exhibited very similar architecture. In comparison to other sedentarians, they exhibited some standard features, including all genes encoded on the same strand, uncommon but not unique duplicated trnM gene, as well as a number of unique features. Their comparatively large size (17,673 bp) can be attributed to four non-coding regions larger than 500 bp. We identified an unusually large (putative) overlap of 14 bases between nad2 and cox1 genes in both species. Importantly, the two species exhibited completely rearranged gene orders in comparison to all other available mitogenomes. Along with Serpulidae and Sabellidae, Polydora is the third identified sedentarian lineage that exhibits disproportionally elevated rates of mitogenomic architecture rearrangements. Selection analyses indicate that these three lineages also exhibited relaxed purifying selection pressures.

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.

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.

Morphology is not a reliable taxonomic tool for the genus Lernaea: molecular data and experimental infection reveal that L. cyprinacea and L. cruciata are conspecific

Background

Species belonging to the genus Lernaea are cosmopolitan parasites that can infect many different freshwater fish hosts. Due to a high degree of morphological intraspecific variability and high levels of interspecific similarities, their classification is extremely difficult and controversial. Although the suitability of the shape of cephalic horns has been questioned decades ago by some experimental infection studies, this character still plays the central role in the identification of Lernaea spp.

Methods

We used the nominal species Lernaea cyprinacea and Lernaea cruciata to test the hypothesis that the shape of the anchor can exhibit host-induced morphological variability, and that the two taxa may be synonymous.

Results

We examined 517 wild or farmed specimens of five host fish species (four cyprinids and a mosquitofish), and found that all 16 parasite specimens collected from mosquitofish could be morphologically identified as L. cruciata, whereas the remaining 25 parasite specimens were all identified as L. cyprinacea. We experimentally infected goldfish and mosquitofish specimens with offspring (copepodids) of a single L. cyprinacea specimen: the adult parasites from goldfish were morphologically identified as L. cyprinacea, and those from mosquitofish as L. cruciata. We then used molecular data to corroborate that all these specimens are conspecific.

Conclusions

Our results suggest that L. cyprinacea and L. cruciata may be synonyms, misidentified as different species as a result of host-induced morphological variation. Given the current shortage of molecular data for the genus Lernaea, in order to resolve the taxonomy of this genus (determine the exact number of species), future studies should aim to sequence as much molecular data as possible, and conduct further experimental infections.

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 of the Caecal Fluke of Poultry, Postharmostomum commutatum, as the First Representative from the Superfamily Brachylaimoidea

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

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.

Characterization of the complete mitochondrial genome of Brentisentisyangtzensis Yu & Wu, 1989 (Acanthocephala, Illiosentidae)

The mitogenome of Brentisentisyangtzensis is 13,864 bp in length and has the circular structure typical of metazoans. It contains 36 genes: 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and 12 protein-encoding genes (PCGs). All genes are transcribed from the same strand. Thirteen overlapping regions were found in the mitochondrial genome. The overall A+T content of B.yangtzensis is 68.3% versus 31.7% of G+C content (A = 27.8%, T = 40.5%, C = 9.0%, G = 22.7%). B.yangtzenensis (Illiosentidae) and Leptorhynchoidesthecatus (Rhadinorhynchidae) form a sister clade, showing the relatively close relationship between the Illiosentidae and the Rhadinorhynchidae. The mitochondrial gene arrangements of acanthocephalan species are relatively conserved, with only a few translocations of tRNAs (trnS1, trnS2, trnV, and trnK) detected. An identical gene order was found both in a sister clade (Centrorhynchusaluconis and Plagiorhynchustransversus) and across different classes (B.yangtzensis (Palaeacanthocephala), Acanthosentischeni (Eoacanthocephala) and Macracanthorhynchushirudinaceus (Archiacanthocephala), Oncicolaluehei and L.thecatus (Palaeacanthocephala)). More studies and more sequences of acanthocephalan species are needed to gain a clear understanding of the phylogenetic relationships.

A novel gene arrangement among the Stylommatophora by the complete mitochondrial genome of the terrestrial slug Meghimatium bilineatum (Gastropoda, Arionoidea)

Stylommatophora is a main clade of Gastropoda that encompasses approximately 112 gastropod families and may exceed a total of 30,000 species. Twenty-four complete stylommatophoran mitogenomes have been sequenced to date, yet our understanding of mitochondrial evolution in stylommatophorans is still in its infancy. To further expand the set of available mitogenomes, we sequenced the mitogenome of Meghimatium bilineatum (Arionoidea: Philomycidae), a widespread land slug in East Asia. This is the first report on a mitogenome of the superfamily Arionoidea, and indeed on a terrestrial slug. The mitogenome of Meghimatium bilineatum comprises 13,972 bp and exhibits a novel, highly distinctive gene arrangement among the Stylommatophora. Phylogenetic reconstructions based on the sequences of all protein-coding genes consistently recovered Meghimatium bilineatum as sister-group of the Succineidae. A phylogenetic reconstruction based on gene order, however, suggested a highly divergent tree topology, which is less credible when taking into account prior knowledge of stylommatophoran relationships. Our CREx (Common interval Rearrangement Explorer) analysis suggested that three successive events of tandem duplication random loss (TDRL) best explain the evolutionary process of gene order rearrangement in Meghimatium bilineatum from an ancestral stylommatophoran mitogenome. The present example offers new insights into the mechanisms of mitogenome rearrangements in gastropods at large and into the usefulness of mitogenomic gene order as a phylogenetic marker.

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.

Characterization of the complete mitochondrial genome of Parabreviscolexniepini Xi et al., 2018 (Cestoda, Caryophyllidea)

Parabreviscolexniepini is a recently described caryophyllidean monozoic tapeworm from schizothoracine fish on the Tibetan Plateau. In the present study, the complete mitochondrial genome of P.niepini is determined for the first time. The mitogenome is 15,034 bp in length with an A+T content of 59.6%, and consists of 12 protein-encoding genes, 22 tRNA genes, two rRNA genes, and two non-coding regions. The secondary structure of tRNAs exhibit the conventional cloverleaf structure, except for trnS1^(AGN) and trnR, which lack DHU arms. The anti-codon of trnS1^(AGN) in the mitogenome of P.niepini is TCT. The two major non-coding regions, 567 bp and 1428 bp in size, are located between trnL2 and cox2, trnG and cox3, respectively. The gene order of P.niepini shows a consistent pattern with other caryophyllideans. Phylogenetic analysis based on mitogenomic data indicates that P.niepini has a close evolutionary relationship with tapeworms Breviscolexorientalis and Atractolytocestushuronensis.

The complete mitochondrial genome of Cymothoa indica has a highly rearranged gene order and clusters at the very base of the Isopoda clade

As a result of great diversity in life histories and a large number of described species, taxonomic and phylogenetic uncertainty permeates the entire crustacean order of Isopoda. Large molecular datasets capable of providing sufficiently high phylogenetic resolution, such as mitochondrial genomes (mitogenomes), are needed to infer their evolutionary history with confidence, but isopod mitogenomes remain remarkably poorly represented in public databases. We sequenced the complete mitogenome of Cymothoa indica, a species belonging to a family from which no mitochondrial genome was sequenced yet, Cymothoidae. The mitogenome (circular, 14484 bp, A+T = 63.8%) is highly compact, appears to be missing two tRNA genes (trnI and trnE), and exhibits a unique gene order with a large number of rearrangements. High compactness and the existence of palindromes indicate that the mechanism behind these rearrangements might be associated with linearization events in its evolutionary history, similar to those proposed for isopods from the Armadillidium genus (Oniscidea). Isopods might present an important model system to study the proposed discontinuity in the dynamics of mitochondrial genomic architecture evolution. Phylogenetic analyses (Bayesian Inference and Maximum Likelihood) conducted using nucleotide sequences of all mitochondrial genes resolved Oniscidea and Cymothoida suborders as paraphyletic. Cymothoa indica was resolved as a sister group (basal) to all remaining isopods, which challenges the accepted isopod phylogeny, where Cymothoida are the most derived, and Phreatoicidea the most basal isopod group. There is growing evidence that Cymothoida suborder might be split into two evolutionary distant clades, with parasitic species being the most basal split in the Isopoda clade, but a much larger amount of molecular resources carrying a high phylogenetic resolution will be needed to infer the remarkably complex evolutionary history of this group of animals with confidence.

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.

Comparative mitogenomics supports synonymy of the genera Ligula and Digramma (Cestoda: Diphyllobothriidae)

Background

After observing differences in the number of reproductive complexes per proglottid within the genus Ligula, the genus Digramma was erected. However, the validity of Digramma has been previously questioned due to a low variability in the cox1, nad1 and ITS rDNA sequences between the two genera. We undertook a study to greatly increase the amount of sequence data available for resolution of this question by sequencing and characterizing the complete mitogenomes of Digramma interrupta and Ligula intestinalis.

Results

The circular mtDNA molecules of Digramma interrupta and Ligula intestinalis are 13,685 bp and 13,621 bp in size, respectively, both comprising 12 PCGs, 22 tRNA genes, two rRNA genes, and two mNCRs. Both mitogenomes exhibit the same gene order and share 92.7% nucleotide identity, compared with 85.8–86.5% to the most closely related genus Dibothriocephalus. Each gene from D. interrupta and L. intestinalis is almost of the same size, and the sequence identity ranges from 87.5% (trnD) to 100% (trnH, trnQ and trnV). NCR2 sequences of D. interrupta and L. intestinalis are 249 bp and 183 bp in length, respectively, which contributes to the main difference in length between their complete mitogenomes. A sliding window analysis of the 12 PCGs and two rRNAs indicated nucleotide diversity to be higher in nad5, nad6, nad2, nad4 and cox3, whereas the most conserved genes were rrnL and rrnS. Lower sequence identity was also found in nad2, nad4, nad5, nad6 and cox3 genes between the two diphyllobothriids. Within the Diphyllobothriidae, phylogenetic analysis indicated Ligula and Digramma to be most closely related to one another, forming a sister group with Dibothriocephalus.

Conclusions

Owing to higher nucleotide diversity, the genes nad2, nad4, nad5, nad6 and cox3 should be considered optimal candidates to use as molecular markers for population genetics and species identification between the two closely related species. The phylogenetic results in combination with the comparative analysis of the two mitogenomes, consistently support the congeneric status of L. intestinalis and D. interrupta.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2910-9) contains supplementary material, which is available to authorized users.

Sequencing of complete mitochondrial genomes confirms synonymization of Hyalomma asiaticum asiaticum and kozlovi, and advances phylogenetic hypotheses for the Ixodidae

Phylogeny of hard ticks (Ixodidae) remains unresolved. Mitochondrial genomes (mitogenomes) are increasingly used to resolve phylogenetic controversies, but remain unavailable for the entire large Hyalomma genus. Hyalomma asiaticum is a parasitic tick distributed throughout the Asia. As a result of great morphological variability, two subspecies have been recognised historically; until a morphological data-based synonymization was proposed. However, this hypothesis was never tested using molecular data. Therefore, objectives of this study were to: 1. sequence the first Hyalomma mitogenome; 2. scrutinise the proposed synonymization using molecular data, i.e. complete mitogenomes of both subspecies: H. a. asiaticum and kozlovi; 3. conduct phylogenomic and comparative analyses of all available Ixodidae mitogenomes. Results corroborate the proposed synonymization: the two mitogenomes are almost identical (99.6%). Genomic features of both mitogenomes are standard for Metastriata; which includes the presence of two control regions and all three "Tick-Box" motifs. Gene order and strand distribution are perfectly conserved for the entire Metastriata group. Suspecting compositional biases, we conducted phylogenetic analyses (29 almost complete mitogenomes) using homogeneous and heterogeneous (CAT) models of substitution. The results were congruent, apart from the deep-level topology of prostriate ticks (Ixodes): the homogeneous model produced a monophyletic Ixodes, but the CAT model produced a paraphyletic Ixodes (and thereby Prostriata), divided into Australasian and non-Australasian clades. This topology implies that all metastriate ticks have evolved from the ancestor of the non-Australian branch of prostriate ticks. Metastriata was divided into three clades: 1. Amblyomminae and Rhipicephalinae (Rhipicephalus, Hyalomma, Dermacentor); 2. Haemaphysalinae and Bothriocrotoninae, plus Amblyomma sphenodonti; 3. Amblyomma elaphense, basal to all Metastriata. We conclude that mitogenomes have the potential to resolve the long-standing debate about the evolutionary history of ticks, but heterogeneous evolutionary models should be used to alleviate the effects of compositional heterogeneity on deep-level relationships.

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.