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

Sequence and Haplotype Analyses of Ligula intestinalis in Acanthobrama marmid (Cyprinidae) in Turkey

PURPOSE

Ligulosis caused by Ligula intestinalis adversely affects the fisheries carried out in the lakes and ponds, causing economic losses in the fish industry. In this study, it was aimed to reveal the molecular characterization of L. intestinalis isolates obtained from woodfish (Acanthobrama marmid) in Keban Dam Lake in Elazig province of Turkey by using mt-CO1 gene sequences and to determine the genetic differences and haplotypes between the isolates.

METHODS

In the examination made in terms of L. intestinalis, the intestine of the fish was opened with the help of fine-tipped scissors, the contents were allowed to come out, and the parasites were taken into a petri dish containing phosphate buffered saline (PBS). Then, L. intestinalis plerocercoids were taken into 15 ml falcon tubes containing 70% ethanol and stored at - 20 °C until further analysis. From each isolate, total gDNA was extracted from the plerocercoids. A partial (480 bp) mt-CO1 gene was amplified by PCR and sequenced unidirectionally. The final size of the trimmed sequences was 392 bp for 43 sequences. Sequence and haplotype analyses were performed, followed by phylogenetic analyses.

RESULTS

All isolates were confirmed as L. intestinalis by BLAST analysis. In addition, 87 nucleotide mutation positions were determined among 43 CO1 gene sequences. As a result of the haplotype network performed for the mt-CO1 gene region of L. intestinalis isolates; arranged in a star-like configuration with the main haplotype (Hap05), separated from other haplotypes by 1-6 mutation steps, and 29 haplotypes were identified, covering 13.9% (6/43) of the total isolates. Also, 75 variable (polymorphic) sites were determined, 52 of which were parsimony informative sites.

CONCLUSIONS

The molecular characterization of L. intestinalis in woodfish (A. marmid) was identified for the first time in Turkey.

Historical dispersal and host-switching formed the evolutionary history of a globally distributed multi-host parasite - The Ligula intestinalis species complex

Studies on parasite biogeography and host spectrum provide insights into the processes driving parasite diversification. Global geographical distribution and a multi-host spectrum make the tapeworm Ligula intestinalis a promising model for studying both the vicariant and ecological modes of speciation in parasites. To understand the relative importance of host association and biogeography in the evolutionary history of this tapeworm, we analysed mtDNA and reduced-represented genomic SNP data for a total of 139 specimens collected from 18 fish-host genera across a distribution range representing 21 countries. Our results strongly supported the existence of at least 10 evolutionary lineages and estimated the deepest divergence at approximately 4.99-5.05 Mya, which is much younger than the diversification of the fish host genera and orders. Historical biogeography analyses revealed that the ancestor of the parasite diversified following multiple vicariance events and was widespread throughout the Palearctic, Afrotropical, and Nearctic between the late Miocene and early Pliocene. Cyprinoids were inferred as the ancestral hosts for the parasite. Later, from the late Pliocene to Pleistocene, new lineages emerged following a series of biogeographic dispersal and host-switching events. Although only a few of the current Ligula lineages show narrow host-specificity (to a single host genus), almost no host genera, even those that live in sympatry, overlapped between different Ligula lineages. Our analyses uncovered the impact of historical distribution shifts on host switching and the evolution of host specificity without parallel host-parasite co-speciation. Historical biogeography reconstructions also found that the parasite colonized several areas (Afrotropical and Australasian) much earlier than was suggested by only recent faunistic data.

Comparative Mitochondrial Genomic Analysis Robustly Supported That Cat Tapeworm Hydatigera taeniaeformis (Platyhelminthes: Cestoda) Represents a Species Complex

Hydatigera taeniaeformis is one of the most common intestinal tapeworms that has a worldwide distribution. In this study, we sequenced the complete mitochondrial (mt) genome of H. taeniaeformis from the leopard cat (designated HTLC) and compared it with those of H. taeniaeformis from the cat in China (designated HTCC) and Germany (designated HTCG). The complete mt genome sequence of HTLC is 13,814 bp in size, which is 167 bp longer than that of HTCC and is 74 bp longer than that of HTCG. Across the entire mt genome (except for the two non-coding regions), the sequence difference was 3.3% between HTLC and HTCC, 12.0% between HTLC and HTCG, and 12.1% between HTCC and HTCG. The difference across both nucleotide and amino acid sequences of the 12 protein-coding genes was 4.1 and 2.3% between the HTLC and HTCC, 13.3 and 10.0% between the HTLC and HTCG, and 13.8 and 10.6% between the HTCC and HTCG, respectively. Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes showed the separation of H. taeniaeformis from different hosts and geographical regions into two distinct clades. Our analysis showed that the cat tapeworm H. taeniaeformis represents a species complex. The novel mt genomic datasets provide useful markers for further studies of the taxonomy and systematics of cat tapeworm H. taeniaeformis.

Comparative Characterization of Mitogenomes From Five Orders of Cestodes (Eucestoda: Tapeworms)

The recognized potential of using mitogenomics in phylogenetics and the more accessible use of high-throughput sequencing (HTS) offer an opportunity to investigate groups of neglected organisms. Here, we leveraged HTS to execute the most comprehensive documentation of mitogenomes for cestodes based on the number of terminals sequenced. We adopted modern approaches to obtain the complete mitogenome sequences of 86 specimens representing five orders of cestodes (three reported for the first time: Phyllobothriidea, “Tetraphyllidea” and Trypanorhyncha). These complete mitogenomes represent an increase of 41% of the mitogenomes available for cestodes (61–147) and an addition of 33% in the representativeness of the cestode orders. The complete mitochondrial genomes are conserved, circular, encoded in the same strand, and transcribed in the same direction, following the pattern observed previously for tapeworms. Their length varies from 13,369 to 13,795 bp, containing 36 genes in total. Except for the Trypanorhyncha specimen, the gene order of the other four cestode orders sequenced here suggests that it could be a synapomorphy for the acetabulate group (with a reversion for taenids). Our results also suggest that no single gene can tell all the evolutionary history contained in the mitogenome. Therefore, cestodes phylogenies based on a single mitochondrial marker may fail to capture their evolutionary history. We predict that such phylogenies would be improved if conducted under a total evidence framework. The characterization of the new mitochondrial genomes is the first step to provide a valuable resource for future studies on the evolutionary relationships of these groups of parasites.

Evolutionary transitions in broad tapeworms (Cestoda: Diphyllobothriidea) revealed by mitogenome and nuclear ribosomal operon phylogenetics

Broad tapeworms (Diphyllobothriidea) are parasites whose adults are capable of infecting a wide range of freshwater, marine and terrestrial tetrapods including humans. Previous works examining the evolution of habitat and host use in this group have been hampered by the lack of a well-resolved phylogeny. In order to produce a robust phylogenetic framework for diphyllobothriideans, we sequenced the complete mitochondrial genome of 13 representatives, carefully chosen to cover the major clades, and two outgroup species representing the Spathebothriidea and Haplobothriidea. In addition, complementary data from the nuclear ribosomal operon was sequenced for 10 representative taxa. Mitogenomes and ssrDNA and lsrDNA were used towards elucidating the phylogenetic framework for the Diphyllobothriidea. The Cephalochlamydidae is confirmed as the earliest diverging diphyllobothriidean lineage, and Solenophoridae and Diphyllobothriidae are sister groups. We infer a probable freshwater origin of the diphyllobothriideans. The ancestral condition for life cycle complexity could not be unambiguously resolved. However, we infer exclusive use of a three-host life cycle following the origin of the Solenophoridae + Diphyllobothriidae. Regarding definitive host use, although we infer reptiles as the most likely ancestral condition, this result should be revisited with a more densely sampled phylogeny in future studies. Freshwater habitat is used by the early diverging lineages within the Solenophoridae + Diphyllobothriidae clade. For the latter, habitat use shifts between freshwater and marine environments, and definitive host use includes marine and terrestrial mammals and birds. We use mitochondrial genomes to distinguish Schistocephalus species occurring in different species of sticklebacks and demonstrate conspecificity of Ligula cf. intestinalis specimens collected from two Fennoscandian ringed seal subspecies.

Mitochondrial genomes of two eucotylids as the first representatives from the superfamily Microphalloidea (Trematoda) and phylogenetic implications

Background

The Eucotylidae Cohn, 1904 (Superfamily: Microphalloidea), is a family of digeneans parasitic in kidneys of birds as adults. The group is characterized by the high level of morphological similarities among genera and unclear systematic value of morphological characters traditionally used for their differentiation. In the present study, we sequenced the complete or nearly complete mitogenomes (mt genome) of two eucotylids representing the genera Tamerlania (T. zarudnyi) and Tanaisia (Tanaisia sp.). They represent the first sequenced mt genomes of any member of the superfamily Microphalloidea.

Methods

A comparative mitogenomic analysis of the two newly sequenced eucotylids was conducted for the investigation of mitochondrial gene arrangement, contents and genetic distance. Phylogenetic position of the family Eucotylidae within the order Plagiorchiida was examined using nucleotide sequences of mitochondrial protein-coding genes (PCGs) plus RNAs using maximum likelihood (ML) and Bayesian inference (BI) methods. BI phylogeny based on concatenated amino acids sequences of PCGs was also conducted to determine possible effects of silent mutations.

Results

The complete mt genome of T. zarudnyi was 16,188 bp and the nearly complete mt genome of Tanaisia sp. was 13,953 bp in length. A long string of additional amino acids (about 123 aa) at the 5′ end of the cox1 gene in both studied eucotylid mt genomes has resulted in the cox1 gene of eucotylids being longer than in all previously sequenced digeneans. The rrnL gene was also longer than previously reported in any digenean mitogenome sequenced so far. The TΨC and DHU loops of the tRNAs varied greatly between the two eucotylids while the anticodon loop was highly conserved. Phylogenetic analyses based on mtDNA nucleotide and amino acids sequences (as a separate set) positioned eucotylids as a sister group to all remaining members of the order Plagiorchiida. Both ML and BI phylogenies revealed the paraphyletic nature of the superfamily Gorgoderoidea and the suborder Xiphidiata.

Conclusions

The average sequence identity, combined nucleotide diversity and Kimura-2 parameter distances between the two eucotylid mitogenomes demonstrated that atp6, nad5, nad4L and nad6 genes are better markers than the traditionally used cox1 or nad1 for the species differentiation and population-level studies of eucotylids because of their higher variability. The position of the Dicrocoeliidae and Eucotylidae outside the clade uniting other xiphidiatan trematodes strengthened the argument for the need for re-evaluation of the taxonomic content of the Xiphidiata.

Hybrid de novo whole-genome assembly and annotation of the model tapeworm Hymenolepis diminuta

Despite the use of Hymenolepis diminuta as a model organism in experimental parasitology, a full genome description has not yet been published. Here we present a hybrid de novo genome assembly based on complementary sequencing technologies and methods. The combination of Illumina paired-end, Illumina mate-pair and Oxford Nanopore Technology reads greatly improved the assembly of the H. diminuta genome. Our results indicate that the hybrid sequencing approach is the method of choice for obtaining high-quality data. The final genome assembly is 177 Mbp with contig N50 size of 75 kbp and a scaffold N50 size of 2.3 Mbp. We obtained one of the most complete cestode genome assemblies and annotated 15,169 potential protein-coding genes. The obtained data may help explain cestode gene function and better clarify the evolution of its gene families, and thus the adaptive features evolved during millennia of co-evolution with their hosts.

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 and phylogenomics of the complete mitochondrial genome of the polyzoic cestode Gangesia oligonchis (Platyhelminthes: Onchoproteocephalidea)

The order Onchoproteocephalidea (Eucestoda) was recently erected to accommodate the hook-bearing tetraphyllideans and the proteocephalideans, which are characterized by internal proglottization and a tetra-acetabulate scolex. The recognized subfamilies in the Proteocephalidae appeared to be non-monophyletic based on 28S recombinant DNA (rDNA) sequence data. Other molecular markers with higher phylogenetic resolution, such as large mitochondrial DNA fragments and multiple genes, are obviously needed. Thus the mitochondrial genome of Gangesia oligonchis, belonging to the putative earliest diverging group of the Proteocephalidae, was sequenced. The circular mitogenome of G. oligonchis was 13,958 bp in size, and contained the standard 36 genes: 22 transfer RNA genes, two rRNA genes and 12 protein-coding genes, as well as two major non-coding regions. A short NCR and a large NCR (lNCR) region were 216 bp and 419 bp in size, respectively. Highly repetitive regions in the lNCR region were detected with that of 11 repeat units. The mitogenome of G. oligonchis shared 71.1% nucleotide identity with Testudotaenia sp. WL-2016. Phylogenetic analyses of the complete mitochondrial genomes with Bayesian inference and maximum likelihood methods indicated that G. oligonchis formed a sister clade with Testudotaenia sp. WL-2016 with maximum support. The ordinal topology is (Caryophyllidea, (Diphyllobothriidea, (Bothriocephalidea, (Onchoproteocephalidea, Cyclophyllidea)))). The mitogenomic gene arrangement of G. oligonchis was identical to that of Testudotaenia sp. WL-2016. Both mitogenomic and nuclear sequence data for many more taxa are required to effectively explore the inter-relationships among the Onchoproteocephalidea.

Long-read sequencing reveals a 4.4 kb tandem repeat region in the mitogenome of Echinococcus granulosus (sensu stricto) genotype G1

Background

Echinococcus tapeworms cause a severe helminthic zoonosis called echinococcosis. The genus comprises various species and genotypes, of which E. granulosus (sensu stricto) represents a significant global public health and socioeconomic burden. Mitochondrial (mt) genomes have provided useful genetic markers to explore the nature and extent of genetic diversity within Echinococcus and have underpinned phylogenetic and population structure analyses of this genus. Our recent work indicated a sequence gap (> 1 kb) in the mt genomes of E. granulosus genotype G1, which could not be determined by PCR-based Sanger sequencing. The aim of the present study was to define the complete mt genome, irrespective of structural complexities, using a long-read sequencing method.

Methods

We extracted high molecular weight genomic DNA from protoscoleces from a single cyst of E. granulosus genotype G1 from a sheep from Australia using a conventional method and sequenced it using PacBio Sequel (long-read) technology, complemented by BGISEQ-500 short-read sequencing. Sequence data obtained were assembled using a recently-developed workflow.

Results

We assembled a complete mt genome sequence of 17,675 bp, which is > 4 kb larger than the complete mt genomes known for E. granulosus genotype G1. This assembly includes a previously-elusive tandem repeat region, which is 4417 bp long and consists of ten near-identical 441–445 bp repeat units, each harbouring a 184 bp non-coding region and adjacent regions. We also identified a short non-coding region of 183 bp, which includes an inverted repeat.

Conclusions

We report what we consider to be the first complete mt genome of E. granulosus genotype G1 and characterise all repeat regions in this genome. The numbers, sizes, sequences and functions of tandem repeat regions remain to be studied in different isolates of genotype G1 and in other genotypes and species. The discovery of such ‘new’ repeat elements in the mt genome of genotype G1 by PacBio sequencing raises a question about the completeness of some published genomes of taeniid cestodes assembled from conventional or short-read sequence datasets. This study shows that long-read sequencing readily overcomes the challenges of assembling repeat elements to achieve improved genomes.

How to become a successful invasive tapeworm: a case study of abandoned sexuality and exceptional chromosome diversification in the triploid carp parasite Atractolytocestus huronensis Anthony, 1958 (Caryophyllidea: Lytocestidae)

Background

A cytogenetic analysis of the new local triploid population of the caryophyllidean tapeworm Atractolytocestus huronensis, a unique parthenogenetic species with the ability to colonise new regions, was performed to understand the inner structure of its chromosome complement.

Methods

A karyotype analysis was carried out using classical Giemsa staining and C-banding combined with fluorescent DAPI staining. A hypothesis that triplets are composed from three homologue chromosomes of approximately the same length and same centromere position was tested statistically for multiple dependent variables using a non-parametric Friedman’s ANOVA. The chromosomal location of ribosomal DNA clusters within the nucleolar organization region (NORs) and telomeric (TTAGGG)_n sequences were detected by fluorescent in situ hybridization (FISH). Chromosomes were subjected to AgNO₃ staining in order to determine whether the rDNA sites represent active NORs.

Results

The cytogenetic analysis confirmed the karyotype composed from eight chromosome triplets (3n = 24) as well as the existence of a pair of NORs located on each chromosome of the second triplet. Six NORs varied their activity from cell to cell, and it was reflected in the numbers of nucleoli (from 1 to 5). A huge morphological diversification of homologue chromosomes was originally detected in six out of eight triplets; the homologue elements differed significantly either in length and/or morphology, and some of them carried discernible interstitial telomeric sequences (ITSs), while the end telomeres were minute. The heterochromatin bands with high AT content varied irregularly, and the course of aberrant spermatogenesis was evident.

Conclusions

Diversification of homologues is a unique phenomenon very likely caused by the long-term absence of a recombination and consequential accumulation of chromosome rearrangements in the genome of A. huronensis during species evolution. Unalterable asexual reproduction of the tapeworm, along with international trade in its host (carp), is facilitating its ongoing spread.

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.