Skip to main content
Have a personal or library account? Click to login
Molecular detection and COX1-based phylogenetic placement of Dipylidium caninum caninum in stray dogs from Iğdır Province, eastern Türkiye Cover

Molecular detection and COX1-based phylogenetic placement of Dipylidium caninum caninum in stray dogs from Iğdır Province, eastern Türkiye

Open Access
|Sep 2026

Full Article

Introduction

Dipylidium caninum (Cestoda: Dipylidiidae) is a globally distributed tapeworm that inhabits the small intestine of domestic dogs and cats and may occasionally infect humans, particularly children (Bowman et al., 2002; Jiang et al., 2017). Given the documented genetic heterogeneity and unresolved taxonomy within the Dipylidium caninum group, the lineage detected in the present study is hereafter referred to as Dipylidium caninum within the D. caninum species complex, unless otherwise specified. Transmission occurs through an indirect life cycle in which cysticercoids develop in flea intermediate hosts, mainly Ctenocephalides felis and Ctenocephalides canis, and definitive hosts become infected by ingesting infected fleas (Bowman et al., 2002).

Although D. caninum is among the most frequently reported cestodes in companion animals, human dipylidiasis has traditionally been considered rare. However, confirmed human cases reported in Asia and Eastern Europe indicate that zoonotic transmission continues in regions with high flea exposure (Ramana et al., 2011; Jiang et al., 2017; Bronstein et al., 2020). A recent global synthesis further supports the ongoing public health relevance of this parasite (Rousseau et al., 2022).

Recent molecular studies have challenged the long-standing view of D. caninum as a genetically uniform species. Substantial mitochondrial divergence has been documented, and at least two host-associated genotypes, adapted to dogs and cats, respectively, have been proposed (Labuschagne et al., 2018; Beugnet et al., 2018). In addition, morphologically similar taxa such as D. carracidoi may have been historically misidentified as D. caninum, highlighting unresolved taxonomic complexity within Dipylidiidae (Millán & Casanova, 2009; Enitez-Bolivar et al., 2022).

In Türkiye, numerous coprological surveys have reported D. caninum in both owned and shepherd dog populations across multiple provinces, with prevalence varying by geography, diagnostic approach and dog population structure (Yaman et al., 2006; Orhun & Ayaz, 2006; Kozan et al., 2007; Yıldırım et al., 2007; Öter et al., 2011; Işık et al., 2014; Nas & Biçek, 2018; Karakuş & Denizhan, 2019; Ceylan et al., 2024). These findings are consistent with earlier studies dating back to the late 1970s, indicating long-term endemic transmission in both rural and urban settings (Güralp et al., 1977; Taşan, 1983; Tınar et al., 1989; Doğanay & Öge, 1993; Şahin et al., 1993; Aydenizöz, 1997; Ayçiçek et al., 1998).

Despite this extensive parasitological background, molecular data on Dipylidium caninum in dogs from Türkiye remain scarce. Although recent work has demonstrated the presence of D. caninum DNA in flea vectors (Özdemir et al., 2025), molecular prevalence and sequence-based characterization studies in canine hosts are still limited, particularly at the regional level. In eastern Türkiye, previous studies have mainly focused on the molecular epidemiology of zoonotic and parasitic agents in animals and environmental samples, highlighting their prevalence and potential public health significance (Ertaş et al., 2022; Ertaş & Ayan, 2022; Kılınç et al., 2023; Koca et al., 2023; Akkuş & Ertaş Oğuz, 2024; Ertaş Oğuz et al., 2025). However, molecular information on Dipylidiidae cestodes in shepherd dog populations from this region remains limited.

Therefore, the present study aimed to investigate the molecular presence of Dipylidiidae cestodes in shepherd dogs from Iğdır Province, eastern Türkiye, using a COX1-based PCR approach and to characterize the detected lineage through partial COX1 gene sequencing. The resulting data provide new molecular evidence from an under-studied region and contribute to a better understanding of Dipylidium diversity and transmission dynamics in shepherd dog populations.

Materials and Methods

Sample Collection

In 2023, a total of 100 fecal samples were collected from shepherd dogs in Iğdır Province, eastern Türkiye. Before sampling, dog owners were informed about the study, and voluntary consent was obtained. During field visits, dogs were evaluated individually with assistance from their owners. Fresh fecal samples were collected non-invasively either directly from the ground when they could be attributed to the relevant dog or immediately after observed defecation. To prevent sample mixing, each sample was collected using disposable gloves, placed into a separate leak-proof container, and labeled with a unique sample code. Age, sex, and flock/farm information were also recorded. The samples were first transported to the laboratory of Iğdır University Animal Hospital and then transferred under cold-chain conditions to the Department of Genetics, Faculty of Veterinary Medicine, Van Yüzüncü Yıl University, for molecular analyses.

DNA extraction

Genomic DNA was extracted from approximately 200 mg of each fecal sample using the GeneMATRIX Stool DNA Purification Kit (Eurx, Poland), following the manufacturer’s protocol. DNA concentration and purity were assessed spectrophotometrically, then all extracts were stored at −20 °C until PCR analysis.

PCR amplification

A 446-bp fragment of the mitochondrial cytochrome c oxidase sub-unit I (COX1) gene was amplified using the primer pair JB3 and JB4.5 described by Bowles et al. (1992). The primer sequences were as follows:

  • JB3 (forward): 5′-TTTTTTGGGCATCCTGAGGTTTAT-3′

  • JB4.5 (reverse): 5′-TAAAGAAAGAACATAATGAAAATG-3′

PCR reactions were performed in a final volume of 25 μL containing 200 μM of each dNTP, 1.5 mM MgCl2, 10 pmol of each primer, 1.25 U Taq DNA polymerase, 1× PCR buffer, nuclease-free water, and 2 μL of template DNA. Amplification was carried out in a Veriti™ Thermal Cycler (Applied Biosystems) with the following cycling conditions: initial denaturation at 95 °C for 15 min, followed by 35 cycles of denaturation at 95 °C for 1 min, annealing at 45 °C for 1 min, and extension at 72 °C for 1 min, with a final extension step at 72 °C for 10 min. PCR products were resolved by electrophoresis on 1.5 % agarose gels stained with Safe-T-Stain, run at 90 V for 60 min, and visualized under UV light using a Syngene Bio-Imaging System.

Sequencing and sequence quality control

PCR amplicons that produced clear bands of the expected size were purified and bidirectionally sequenced by BM Labosis Laboratory (Ankara, Türkiye). Forward and reverse chromatograms (.ab1 files) were examined in MEGA11 (Tamura et al., 2021), and base calls were manually inspected. Only sequences with clear chromatogram peaks and agreement between forward and reverse reads at all variable positions were used for further analyses. A 393-bp consensus COX1 sequence was generated from the forward and reverse reads.

BLAST analysis

The consensus COX1 sequence was searched against the NCBI nucleotide database using BLASTn to assess molecular similarity with reference sequences. BLASTn hits were evaluated according to query coverage, percentage identity, E-value, accession number, host, country of origin, submission date, and reference/source information. The relevant BLASTn matches to Dipylidium caninum and related cestode sequences are presented in Supplementary Table 1.

Supplementary Table S1.

BLASTn comparison and phylogenetic dataset of Dipylidium caninum COX1 sequences

No.AccessionIsolate / descriptionHostCountryIdentity (%)Query coverSequence length (bp)
1PZ299077.1D. caninum isolate Igdir_Dog1Shepherd dog / Canis lupus familiarisTürkiye100100 %396
2PX765914.1D. caninum isolate Hp03DogGhana92.1293 %373
3OR251824.1D. caninum isolate Sv_PS_Amz_22Speothos venaticus / bush dogColombia91.92100 %1602
4OR251823.1D. caninum isolate Sv_PS_Amz_21Speothos venaticus / bush dogColombia91.92100 %1602
5OK523384.1D. caninum isolate Dcan_canine_FL1, complete mitogenomeCanis lupus familiarisUSA91.92100 %14296
6PZ165824.1D. caninum isolate DICAACanis lupus familiarisCroatia91.92100 %396
7PX765913.1D. caninum isolate Hp02DogGhana91.8593 %373
8PX765912.1D. caninum isolate Hp01DogGhana91.8593 %373
10MN099047.1D. caninum complete mitogenomeStray dogChina91.67100 %14226
11PP842203.1D. caninum isolate DogSRTDc2DogTürkiye91.2795 %383
12PP842202.1D. caninum isolate DogSRTDc1DogTürkiye91.2795 %383
13MG587892.1D. caninum isolate R166, complete mitogenomeFelineSouth Africa89.899 %13598
14OQ281679.1D. caninum isolate VCG6Felis catusMexico89.899 %416
15OK523385.1D. caninum isolate Dcan_feline_KS1, complete mitogenomeFelis catusUSA89.899 %13598
16MT806359.1D. caninum voucher 4245P18Vulpes vulpes / red foxItaly89.899 %411
17PX765915.1D. caninum isolate Hp04DogGhana89.6793 %373
18OR511473.1D. caninum isolate Alex-2Stray catEgypt89.5392 %363
19PZ165825.1D. caninum isolate DICAB.1Felis catusCroatia89.3990 %358
20ON506044.1D. caninum isolate Isb1CatPakistan89.3492 %384
21PP054312.1D. caninum haplotype AFelis catus familiarisEthiopia89.2796 %963
22OR511472.1D. caninum isolate Alex-1Stray catEgypt88.7192 %363
23ON954760.1D. caninum isolate dog1DogIraq88.0192 %384
24ON954628.1D. caninum isolate cat1CatIraq88.0192 %384
25ON514129.1D. caninum isolate Isb2CatPakistan88.0192 %384
26ON954763.1D. caninum isolate cat1CatIraq87.9192 %383
27ON954761.1D. caninum isolate dog2DogIraq86.5892 %384
28PP054311.1Dipylidium sp. 01-SA haplotype AProteles cristata / aardwolfSouth Africa85.86100 %1647
29PZ165826.1D. caninum isolate DICAB.2Canis lupus familiarisCroatia358
30ON954762.1D. caninum isolate dog3DogIraq384
31ON954764.1D. caninum isolate cat2CatIraq384
32ON954765.1D. caninum isolate cat3CatIraq384

Multiple sequence alignment

Selected reference COX1 sequences of D. caninum and related dipylidiid taxa were retrieved from GenBank. All sequences were combined into a single FASTA dataset. The sequences were first checked, and multiple sequence alignment was then performed using MAFFT with the “adjust direction” option. Poorly aligned, ambiguous, and non-overlapping regions were removed to improve alignment quality, resulting in a final alignment of 358 bp. The accession numbers and metadata of these sequences are listed in Supplementary Table 1. ON981095 (Joyeuxiella pasqualei) was used as the outgroup.

Phylogenetic analysis

Phylogenetic analysis was performed in IQ-TREE v3 using the Maximum Likelihood method. The best-fit nucleotide substitution model was selected with ModelFinder Plus, and the GTR+I+G model was used for tree construction. Branch support was assessed with 1,000 ultrafast bootstrap replicates and 1,000 SH-aL-RT replicates. The resulting tree was rooted using Joyeuxiella pasqualei ON981095 as the outgroup and visualized in FigTree v1.4.5. Branch values are shown as SH-aLRT support (%) / ultra-fast bootstrap support (%).

Literature data extraction

To place the findings from Iğdır in a broader epidemiological context, published prevalence data for Dipylidium caninum infections in dogs were compiled from the literature. For each country, mean prevalence values were calculated across available studies irrespective of diagnostic method, and these country-level estimates were subsequently grouped to generate continent-level summaries. Although this approach does not constitute a formal meta-analysis, it provides a comparative framework for interpreting regional variation and situating the Iğdır data within global transmission patterns. A descriptive overview of the global prevalence of D. caninum in dogs is presented in Supplementary Figure S1. Data handling and visualization were performed using R statistical software (R Core Team, 2025).

Supplementary Fig. S1.

Global prevalence of Dipylidium caninum in dogs.

Reported canine D. caninum prevalence data from previously published studies were descriptively summarized by continent and country to provide broad epidemiological context. (A) Prevalence values across continents. Box plots show the median, interquartile range, and individual study points. (B) Country-level prevalence grouped by continent. Colored points represent individual studies according to approximate mid-study year. No molecular sequence data were included in this overview. Because the included studies differed in sampling design, dog population, geographic setting, and diagnostic method, this figure should be interpreted only as descriptive context and not as a formal quantitative comparison among regions.

Ethical Approval and/or Informed Consent

Permission was obtained from the animal owners for this study. Additionally, fecal samples were collected from the ground without coming into contact with the animals. According to Turkish HADYEK (Local Ethics Committee for Animal Research) regulation, fecal sample collection is not subject to ethics committee approval when no animal handling or intervention is performed. Therefore, an ethics committee approval number was not required for this study.

Results

Out of 100 screened fecal samples, eight (8 %) yielded a PCR product of the expected size for the COX1 gene fragment. All positive samples produced identical sequences across the analyzed region, indicating the presence of a single mitochondrial haplotype among the infected dogs.

The obtained 396 bp partial mitochondrial COX1 sequence of Dipylidium caninum isolate Igdir_Dog1 from a shepherd dog in Iğdır Province, Türkiye, was deposited in GenBank under accession number PZ299077.1.

BLASTn analysis showed that the sequence matched D. caninum entries deposited in GenBank, including partial COX1 sequences and complete mitochondrial genome records. The highest identity was observed with a dog-derived sequence from Ghana (PX765914.1; 92.12 % identity, 93 % query cover), followed by bush dog-derived Colombian sequences (OR251824.1 and OR251823.1), a canine complete mitogenome from the USA (OK523384.1), and a dog-derived Croatian sequence (PZ165824.1), each showing 91.92 % identity and 100 % query cover. Turkish dog-derived sequences from Elazığ (PP842203.1 and PP842202.1) showed 91.27 % identity with 95 % query cover. Cat-associated sequences were also included in the comparison and showed lower identity values, ranging approximately from 88.71 % to 89.80 %. The full BLASTn comparison is provided in Supplementary Table 1.

Discussion

The present study provides molecular evidence of Dipylidium caninum in dogs from Iğdır Province, eastern Türkiye, and represents one of the first COX1-based molecular characterizations of Dipylidiidae in this region. Although D. caninum has long been reported in Türkiye using conventional parasitological methods, molecular data remain limited. Two unpublished GenBank records from Elazığ (PP842202.1 and PP842203.1; Aslan Celik et al., GenBank direct submissions) and the molecular detection of D. caninum DNA in fleas from Tokat indicate that molecular evidence from Türkiye is gradually increasing (Özdemir et al., 2025). In this context, the Iğdır sequence reported here adds new COX1 data from eastern Türkiye and helps expand the limited molecular data-set available for Türkiye.

In the phylogenetic tree, the analyzed D. caninum sequences were arranged into three broad groups, but the separation between dog- and cat-associated lineages was unclear (Fig. 1). The first group mainly contained cat-associated and mixed-host sequences from South Africa, Mexico, the USA, Italy, Egypt, Croatia, Pakistan and Ethiopia (Labuschagne et al., 2018; Jesudoss Chelladurai et al., 2023; Camacho-Giles et al., 2024; Citterio et al., 2021; Satour et al., 2026; Beck & Šikić, 2026; Alvi et al., 2022; Dumendiak et al., 2024). This group also included some non-cat hosts, suggesting that host origin alone may not fully explain the observed phylogenetic pattern. Dogs dominated the second group and other canine-derived sequences from Ghana, Colombia, the USA, Croatia, China and Türkiye (Xie et al., 2019; Uribe et al., 2023; Jesudoss Chelladurai et al., 2023; Addy et al., 2026; Beck & Šikić, 2026; Aslan Çelik et al., 2024). The Iğdır isolate PZ299077.1 was positioned within this broad canine-derived group, together with the previously submitted Turkish dog sequences from Elazığ and several dog or wild-canid sequences from other regions. This placement indicates that the Iğdır sequence is not isolated from the global D. caninum dataset, but falls within a group largely represented by canine hosts. The third group included mainly Iraqi and Pakistani sequences obtained from both dog and cat hosts, forming a mixed regional cluster (Al-Ardi, 2022; Alvi et al., 2022). This pattern may reflect regional sequence similarity rather than a strict host-associated structure. Overall, the phylogenetic pattern supports the identification of the Iğdır isolate D. caninum and shows that it contributes new eastern Turkish COX1 data to the limited molecular records available from Türkiye. However, because several branches showed low to moderate support values and the analyzed COX1 fragment was relatively short, the observed clustering should be interpreted cautiously rather than as definitive evidence of host-specific lineages.

Fig. 1.

Maximum Likelihood phylogenetic tree based on a 363-bp COX1 alignment of Dipylidium caninum and related Dipylidiidae sequences.

The tree was reconstructed in IQ-TREE v3 using the GTR+I+G substitution model and rooted with Joyeuxiella pasqualei ON981095 as the outgroup. Branch values are shown as SH-aLRT support (%) / ultrafast bootstrap support (%) based on 1,000 replicates. The Turkish isolate generated in this study, PZ299077.1, is highlighted in red.

Previous molecular and genomic studies have shown that D. caninum contains largely host-associated canine and feline lineages, with marked genetic differentiation between these groups (Labuschagne et al., 2018; Rousseau et al., 2022; Jesudoss Chelladurai et al., 2023). However, host association is not strict, and dogs may carry parasites related to different genotype groups. Recent data from Ghana also reported both canine and feline–canine-associated genotypes, although most dog-derived isolates were assigned to the canine genotype (Addy et al., 2026). Therefore, the placement of PZ299077.1 near Ghanaian dog-derived sequences and other canine-host records supports its affinity with the canine-derived group. Considering the known mitochondrial diversity within the D. caninum complex, the relatively low COX1 identity values observed here do not preclude this placement. However, because the analysis was based on a partial COX1 fragment, the result should be interpreted as molecular placement within the canine-derived group rather than as a final genotype-level conclusion.

Earlier parasitological surveys conducted in different provinces of Türkiye have repeatedly reported D. caninum in dogs using copromicroscopic examination or necropsy-based approaches, with prevalence values varying widely depending on geographic region, dog population structure, and diagnostic methodology (Güralp et al., 1977; Taşan, 1983; Tınar et al., 1989; Doğanay & Öge, 1993; Aydenizöz, 1997; Ayçiçek et al., 1998; Yaman et al., 2006; Kozan et al., 2007; Yıldırım et al., 2007; Öter et al., 2011; Işık et al., 2014; Nas & Biçek, 2018; Karakuş & Denizhan, 2019; Ceylan et al., 2024). However, most of these investigations lacked molecular confirmation, leaving uncertainties regarding species identity and the extent of genetic diversity within Dipylidiidae. The national overview presented in Figure 2 supports this pattern. It shows that D. caninum has been reported from both the western and eastern regions of Türkiye, using both dog-based surveys and flea-based molecular detection. The 5 % molecular prevalence detected in Iğdır therefore adds PCR-based evidence from eastern Türkiye and supports the need for further molecular surveillance. Additionally, the Tokat flea-based PCR record was considered as supporting evidence for transmission rather than as directly comparable dog prevalence data (Özdemir et al., 2025).

Fig. 2.

Geographic and temporal distribution of Dipylidium caninum in Türkiye.

The map and accompanying bar plot summarize published prevalence data from 2006 to 2025 across different provinces of Türkiye. Coprological and necropsy-based surveys are shown in light shading, whereas PCR-based studies are shown in dark shading. The figure illustrates the long-term and widespread circulation of D. caninum in both western and eastern regions of the country.

In conclusion, this study provides a new COX1-based molecular record of D. caninum in dogs from Iğdır Province, eastern Türkiye. It adds new mitochondrial data to the limited molecular dataset currently available from Türkiye. The analyzed sequence showed affinity with canine-derived D. caninum sequences in the phylogenetic analysis and contributes additional information to the growing body of evidence for genetic diversity within the D. caninum complex. In line with recent studies, the present findings also highlight the need for broader sampling and the inclusion of additional mitochondrial and nuclear markers from underrepresented regions to understand better the evolutionary structure and global diversity of this parasite. From a One Health perspective, continued molecular surveillance of dogs, together with monitoring of flea intermediate hosts, may contribute to an improved understanding of the transmission dynamics and zoonotic risk associated with canine dipylidiasis.

Acknowledgements

The authors thank the Department of Genetics, Faculty of Veterinary Medicine, Van Yüzüncü Yıl University, for providing technical support and access to laboratory facilities.

Notes

[1] Contributed by Author Contributions

Fatma Ertaş Oğuz: Study design, field sampling, laboratory work, and final review of the manuscript. Hasret Öztürk: GenBank data handling, data analysis, figure preparation, and manuscript writing and revision. Adnan Ayan: Laboratory work, data analysis, and final review of the manuscript.

[2] Data Availability

The COX1 sequence generated in this study has been deposited in GenBank under accession number PZ299077.1.

[3] Conflicts of interest Conflict of Interest

The authors declare no conflicts of interest.

[4] Financial disclosure Funding

This study was conducted using the authors’ own resources and did not receive external funding.

DOI: https://doi.org/10.2478/helm-2026-0011 | Journal eISSN: 1336-9083 | Journal ISSN: 0440-6605
Language: English
Page range: 99 - 108
Submitted on: Jan 3, 2026
Accepted on: Jun 15, 2026
Published on: Sep 7, 2026
Published by: Slovak Academy of Sciences, Institute of Parasitology
In partnership with: Paradigm Publishing Services
Publication frequency: Volume open

© 2026 F. E. Oğuz, A. Ayan, H. Öztürk, A. H. Kirmizigül, published by Slovak Academy of Sciences, Institute of Parasitology
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.