Characterization of the complete mitochondrial genome of the cloacal tapeworm Cloacotaenia megalops (Cestoda: Hymenolepididae)

First molecular insights into gastrointestinal helminths of domestic birds in the Caspian Sea Littoral of Iran with an emphasis on the One Health concern

The current investigation was carried out during the period from July 2022 to March 2023, aiming to investigate the prevalence of gastrointestinal helminths in domestic birds collected from traditional markets in Guilan province. One hundred forty-eight domestic birds, including chickens (Gallus gallus domesticus), domestic ducks (Anas platyrhynchos domesticus), greylag geese (Anser anser), and domestic turkeys (Meleagris gallopavo domesticus) were examined. Totally, 42.56% of the investigated birds were positive for helminthic parasites. Morphological analysis revealed varying infection rates among birds: Echinostoma revolutum (5.40%), Hypoderaeum conoideum (2.02%), Cloacotaenia megalops (0.67%), Hymenolepididae family (4.05%), Ascaridia galli (16.89%), and Heterakis gallinarum (4.72%). The investigation involved molecular analysis of the 18S and ITS1 + 5.8S + ITS2 rRNA gene regions. The findings indicated that the 18S region of nematode isolates exhibited a similarity of 92 to 100% with sequences in the GenBank, whereas trematode and cestode isolates showed a gene similarity ranging from 88 to 99%. The ITS regions of nematode, trematode, and cestode isolates exhibited genetic similarities ranging from 87 to 100%, 73-99%, and 75-99%, respectively. Furthermore, phylogenetic analysis confirmed the categorization of the identified species within the Ascaridiidae, Heterakidae, Hymenolepididae, and Echinostomatidae families, indicating their close affinity with previously documented species. Implementing precise control measures such as consistent monitoring, adequate sanitation protocols, and administering anthelmintic treatments is crucial for effectively managing parasitic infections in free-range and backyard poultry farms. Additionally, conducting further surveys is advisable to assess the impact of these parasites on the health and productivity of poultry in the investigated area.

Genomic differences and species delimitation: a case for two species in the zoonotic cestode Dipylidium caninum

Dipylidium caninum (Linnaeus, 1758) is a common zoonotic cestode of dogs and cats worldwide. Previous studies have demonstrated the existence of largely host associated canine and feline genotypes based on infection studies, genetic differences at the nuclear 28S rDNA gene and complete mitochondrial genomes. There have been no comparative studies at a genome-wide scale. Here, we sequenced the genomes of a dog and cat isolate of Dipylidium caninum from the United States using the Illumina platform and conducted comparative analyses with the reference draft genome. Complete mitochondrial genomes were used to confirm the genotypes of the isolates. D. caninum canine and feline genomes generated in this study had mean coverage depths of 45x and 26x and an average identity of 98% and 89% respectively when compared to the reference genome. SNPs were 20 times higher in the feline isolate. Comparison and species delimitation using universally conserved orthologs and protein coding mitochondrial genes revealed that the canine and feline isolates are different species. Data from this study builds a base for future integrative taxonomy. Further genomic studies from geographically diverse populations are necessary to understand implications for taxonomy, epidemiology, veterinary clinical medicine, and anthelmintic resistance.

Heteroplasmic mitochondrial genomes of a Raillietina tapeworm in wild Pangolin

Background

Raillietina species belong to the family Davaineidae, which parasitizes in a wide variety of mammals and birds, causing stunted growth, lethargy, emaciation, and digestive tract obstruction. However, only a limited number of Raillietina species have been identified in wild animals.

Methods

We analyzed and annotated the complete mitochondrial (mt) genome of a worm from the intestine of a wild pangolin using Illumina sequencing of whole genomic DNA.

Results

These findings showed the presence of two mtDNA sequences in Raillietina sp., designated as mt1 and mt2, with the lengths of 14,331 bp and 14,341 bp, respectively. Both the mts genomes of Raillietina sp. comprised 36 genes, containing 12 protein-coding genes (PCGs), 2 ribosomal RNAs, and 22 transfer RNAs. Gene arrangements of both mt genomes of Raillietina sp. were similar to those of most flatworms, except for taeniids, which shift positions between tRNAL1 and tRNAS2 genes. Twenty of 22 tRNA secondary structures of Raillietina sp. had a typical cloverleaf structure similar to Raillietina tetragona. Sequence differences between the mt1 and mt2 genomes were 4.4%, and this difference arises from the mtDNA heteroplasmic mutations. Moreover, heteroplasmic mtDNA mutations were detected in PCGs, tRNAs, rRNAs, NCRs, and intergenes, but the highest proportion of heteroplasmy of 79.0% was detected in PCGs, indicating the occurrence of mtDNA heteroplasmy in Raillietina sp. To our knowledge, this is the first report of mtDNA heteroplasmy in tapeworm parasites. Phylogenetic analyses of 18S rRNA, ITS2, and 12 PCG sequences demonstrated that the worm was clustered with other Raillietina species in the Davaneidae family.

Conclusions

We found a novel Raillietina species in wild pangolin with the existence of mitochondrial DNA heteroplasmy. Thus, these findings provide insights into the heterogeneity of the mt genome in parasitic cestodes, and mt genome data contributes to the understanding of pangolin-parasitic cestodes in terms of their molecular biology, epidemiology, diagnosis, and taxonomy.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05301-y.

Phylogeny of hymenolepidids (Cestoda: Cyclophyllidea) from mammals: sequences of 18S rRNA and COI genes confirm major clades revealed by the 28S rRNA analyses

The aim of the study is to test a hypothesis for the phylogenetic relationships among mammalian hymenolepidid tapeworms, based on partial (D1-D3) nuclear 28S ribosomal RNA (rRNA) genes, by estimating new molecular phylogenies for the group based on partial mitochondrial cytochrome c oxidase I (COI) and nuclear 18S rRNA genes, as well as a combined analysis using all three genes. New sequences of COI and 18S rRNA genes were obtained for Coronacanthus integrus, C. magnihamatus, C. omissus, C. vassilevi, Ditestolepis diaphana, Lineolepis scutigera, Spasskylepis ovaluteri, Staphylocystis tiara, S. furcata, S. uncinata, Vaucherilepis trichophorus and Neoskrjabinolepis sp. The phylogenetic analyses confirmed the major clades identified by Haukisalmi et al. (Zoologica Scripta 39: 631-641, 2010): Ditestolepis clade, Hymenolepis clade, Rodentolepis clade and Arostrilepis clade. While the Ditestolepis clade is associated with soricids, the structure of the other three clades suggests multiple evolutionary events of host switching between shrews and rodents. Two of the present analyses (18S rRNA and COI genes) show that the basal relationships of the four mammalian clades are branching at the same polytomy with several hymenolepidids from birds (both terrestrial and aquatic). This may indicate a rapid radiation of the group, with multiple events of colonizations of mammalian hosts by avian parasites.

Dipylidium caninum draft genome - a new resource for comparative genomic and genetic explorations of flatworms

Here, we present a draft genome of the tapeworm Dipylidium caninum (family Dipylidiidae) and compare it with other cestode genomes. This draft genome of D. caninum is 110 Mb in size, has a repeat content of ~13.4% and is predicted to encode ~10,000 protein-coding genes. We inferred excretory/secretory molecules (representing the secretome), other key groups of proteins (including peptidases, kinases, phosphatases, GTPases, receptors, transporters and ion-channels) and predicted potential intervention targets for future evaluation. Using 144 shared single-copy orthologous sequences, we investigated the genetic relationships of cestodes for which nuclear genomes are available. This study provides first insights into the molecular biology of D. caninum and a new resource for comparative genomic and genetic explorations of this and other flatworms.

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.

Mitochondrial genome of Paruterina candelabraria (Cestoda: Paruterinidae), with implications for the relationships between the genera Cladotaenia and Paruterina

The taxonomic concept for the family Paruterinidae is controversial, especially concerning the position of the genus Cladotaenia, since the latter genus has been placed sometimes in other families, i.e. in the Taeniidae or in the distinct family Cladotaeniidae; finding a solution based on morphological data is difficult and molecular data on paruterinids and related groups are scarce. In this study, the complete mitochondrial (mt) genome sequence of the type-species of the type-genus of the Paruterinidae, Paruterina candelabraria, was determined and annotated. Gene arrangements are identical with those of Cladotaenia vulturi but differing from those of species of the family Taeniidae by the order change between tRNA-SerUCN and tRNA-LeuCUN. Phylogenetic tree was constructed by Bayesian Inference (BI) analysis using the concatenated amino acid sequences of 12 protein-coding genes. The analysis clearly shows that the Paruterinidae and Taeniidae are sister-groups, and Cladotaenia is a sister taxon of Paruterina. This supports the position of the genus Cladotaenia in the family Paruterinidae and reveals the necessity for sequencing additional taxa of the Paruterinidae for better understanding of phylogenetic relationships within the group.

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.

Moniezia benedeni and Moniezia expansa are distinct cestode species based on complete mitochondrial genomes

Moniezia spp. parasitize the intestines of ruminants, causing monieziasis. In this study, the complete mitochondrial (mt) genomes of M. benedeni and M. expansa have been determined, characterized and employed to test the hypothesis that M. benedeni and M. expansa are distinct species by phylogenetic analysis based on the concatenated amino acid sequences derived from 12 protein-coding genes, inferred with Bayesian and Maximum-likelihood methods. The complete mt genomes of M. benedeni and M. expansa were 13,958bp and 13,934bp in size, respectively. Nucleotide sequence identity between the two mt genomes was 83.4%. Each of the two circular mt genomes encodes 36 genes including two ribosomal RNA genes, 22 transfer RNA genes and 12 protein-coding genes, which are transcribed from the same direction. The gene orders of the two mt genomes are identical to those of Anoplocephala spp. (Anoplocephalidae), Hymenolepis spp. (Hymenolepididae) and Dipylidium caninum (Dipylidiidae), but distinct from the species of the Taeniidae family. Phylogenetic analysis confirmed that M. benedeni and M. expansa are taxonomically valid species and have a sister relationship, regardless of the analytical method employed. Furthermore, comparing the cox1 gene sequences of Moniezia spp. from the NCBI deposited sequences and the ones obtained in the present study revealed that the nucleotide sequence differences were 12.5% for M. benedeni and 6.2% for M. expansa, respectively, suggesting the existence of cryptic species in these parasites. The complete mt genome sequences reported herein will be valuable in further studies of diagnoses, molecular ecology and population genetics of Moniezia spp. of socio-economic importance.