Abstract
We investigated the mitochondrial genetic diversity of Echinococcus multilocularis in red foxes from northern and northeastern Poland, and tested the hypothesis that the presence of Asian haplotypes in Poland stems from their introduction from Northeast Asia to Eastern Europe. We targeted the cytochrome c oxidase subunit 1 gene (cox1) of 263 tapeworms from 59 red foxes hunted in selected districts of Pomorskie, Warmińsko-Mazurskie, and Podlaskie Voivodships. We calculated several measures of genetic variation, and constructed median-joining haplotype networks and a RelTime timetree. In addition, a novel alignment-free computational method was applied: the 4D-Dynamic Representation of DNA/RNA Sequences combined with K-means clustering. PCR and sequencing were successful for 252 tapeworms. We identified six haplotypes, three of which were novel. The mismatch distribution was bimodal. The fixation index (FST) between the European and Asian haplogroups was 0.97. Our findings suggest secondary contact between two previously isolated metapopulations with distinct genetic lineages. Clustering patterns appear to support the hypothesis that Asian variants of E. multilocularis spread into northeastern Poland as a consequence of multiple introductions of raccoon dogs from the Soviet Far East to Eastern Europe.
Data availability
The 252 complete coding sequences of the cytochrome c oxidase subunit 1 gene (1,608 bp) obtained in this study were submitted to GenBank and received accessions PV055184–PV055435.
References
Sulima, M. et al. Liver transplantation for incurable alveolar echinococcosis: an analysis of patients hospitalized in department of tropical and parasitic diseases in gdynia. Transplant. Proc. 48, 1708–1712 (2016).
Wen, H. et al. Echinococcosis: advances in the 21st Century. Clin. Microbiol. Rev. 32, 10 (2019).
Knapp, J. et al. Genetic diversity of the Cestode Echinococcus multilocularis in red foxes at a continental scale in Europe. PLoS Negl. Trop. Dis. 3, e452 (2009).
Laurimaa, L. et al. First report of the zoonotic tapeworm Echinococcus multilocularis in raccoon dogs in Estonia, and comparisons with other countries in Europe. Vet. Parasitol. 212, 200–205 (2015).
Santoro, A. et al. Mitochondrial genetic diversity and phylogenetic relationships of Echinococcus multilocularis in Europe. Int. J. Parasitol. 54, 233–245 (2024).
McDevitt, A. D. et al. Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes). Europe. Mol. Ecol. 31, 993–1006 (2022).
Dalén, L. et al. Ancient DNA reveals lack of postglacial habitat tracking in the arctic fox. Proc. Natl. Acad. Sci. 104, 6726–6729 (2007).
Nakao, M., Sako, Y. & Ito, A. Isolation of polymorphic microsatellite loci from the tapeworm Echinococcus multilocularis. Infect. Genet. Evol. 3, 159–163 (2003).
Haag, K., Araújo, A., Gottstein, B. & Zaha, A. Selection, Recombination and History in a Parasitic Flatworm (Echinococcus) Inferred from Nucleotide Sequences. Mem. Inst. Oswaldo Cruz 93, 695–702 (1998).
Bartel, M., Seesee, F. & Worley, D. Comparison of Montana and Alaska isolates of Echinococcus multilocularis in gerbils with observations on the cyst growth, hook characteristics, and host response. J. Parasitol. 78, 529–532 (1992).
Konyaev, S. V. et al. Genetic diversity of Echinococcus spp.in Russia. Parasitol. 140, 1637–1647 (2013).
Umhang, G. et al. Asian admixture in European Echinococcus multilocularis populations: new data from poland comparing EmsB microsatellite analyses and mitochondrial sequencing. Front. Vet. Sci. 7, 620722 (2021).
Karamon, J. et al. Genetic diversity of Echinococcus multilocularis in red foxes in Poland: The first report of a haplotype of probable Asian origin. Folia Parasitol (Praha) 64, 10 (2017).
Karamon, J. et al. Echinococcus multilocularis genetic diversity based on isolates from pigs confirmed the characteristic haplotype distribution and the presence of the Asian-like haplotype in Central Europe. J. Vet. Res. 67, 567–574 (2023).
Gładysz, P. & Lass, A. Detection of Asian genetic components in autochthonous human Echinococcus multilocularis infections from endemic Warmia-Masuria (north-eastern Poland). One Health 17, 100623 (2023).
Gładysz, P. et al. Genetic diversity of Echinococcus multilocularis from red foxes and humans in northern and northeastern Poland investigated using the microsatellite EmsB. Infect. Genet. Evol. 135, 105843 (2025).
Hofer, S. et al. High prevalence of Echinococcus multilocularis in urban red foxes (Vulpes vulpes) and voles (Arvicola terrestris) in the city of Zürich Switzerland. Parasitology 120, 135–142 (2000).
Gładysz, P., Świątalska, A., Samorek-Pieróg, M., Karamon, J. & Lass, A. An update on the prevalence of Echinococcus multilocularis in red foxes (Vulpes vulpes) in a Central European focus: Northern and northeastern Poland (2022–2024). One Health 21, 101151 (2025).
Nakao, M. et al. Geographic pattern of genetic variation in the fox tapeworm Echinococcus multilocularis. Parasitol Int. 58, 384 (2009).
Nei, M. Molecular Evolutionary Genetics (Columbia University Press, 1987). https://doi.org/10.7312/nei-92038.
Nei, M. & Li, W. H. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. 76, 5269–5273 (1979).
Rozas, J. et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302 (2017).
Jombart, T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405 (2008).
Jombart, T. & Ahmed, I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics 27, 3070–3071 (2011).
Paradis, E. pegas: an R package for population genetics with an integrated–modular approach. Bioinformatics 26, 419–420 (2010).
R Core Team. R: A Language and environment for statistical computing. Preprint at https://www.R-project.org/ (2023).
Bandelt, H. J., Forster, P. & Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).
Kumar, S. et al. MEGA12: molecular evolutionary genetic analysis version 12 for adaptive and green computing. Mol. Biol. Evol. 41, msae263 (2024).
Nakao, M., McManus, D. P., Schantz, P. M., Craig, P. S. & Ito, A. A molecular phylogeny of the genus Echinococcus inferred from complete mitochondrial genomes. Parasitology 134, 713–722 (2006).
Knapp, J. et al. Phylogenetic relationships within Echinococcus and Taenia tapeworms (Cestoda: Taeniidae): An inference from nuclear protein-coding genes. Mol. Phylogenet. Evol. 61, 628–638 (2011).
Marshall, L. G., Webb, S. D., Sepkoski, J. J. & Raup, D. M. Mammalian Evolution and the Great American Interchange. Science 1979(215), 1351–1357 (1982).
Bielińska-Wąż, D. & Wąż, P. Non-standard bioinformatics characterization of SARS-CoV-2. Comput. Biol. Med. 131, 104247 (2021).
Bielińska-Wąż, D., Wąż, P., Lass, A. & Karamon, J. 4D-dynamic representation of DNA/RNA sequences: studies on genetic diversity of echinococcus multilocularis in red foxes in Poland. Life 12, 877 (2022).
Bielińska-Wąż, D., Wąż, P. & Błaczkowska, A. 20D-dynamic representation of protein sequences combined with K-means clustering. Comb. Chem. High Throughput Screen https://doi.org/10.2174/0113862073359729250220131623 (2025).
Umhang, G. et al. A step forward in the understanding of the presence and expansion of Echinococcus multilocularis in Eastern Europe using microsatellite EmsB genotyping in Poland. Infection Genet. Evol. 54, 176–182 (2017).
Avise, J. C. Phylogeography: the history and formation of species (Harvard University Press, UK, 2000).
Guo, B. et al. The genetic variation of mitochondrial sequences and pathological differences of Echinococcus multilocularis strains from different continents. Spectr. Microbiol. https://doi.org/10.1128/spectrum.01318-24 (2025).
Hayashi, N. et al. Mitogenomic exploration supports the historical hypothesis of anthropogenic diffusion of a zoonotic parasite Echinococcus multilocularis. iScience 26, 107741 (2023).
Kauhala, K. & Kowalczyk, R. Invasion of the raccoon dog Nyctereutes procyonoides in Europe: History of colonization, features behind its success, and threats to native fauna. Curr. Zool. 57, 584–598 (2011).
Nowak, E. & Pielowski, Z. Die Verbreitung des Marderhundes in Polen im Zusammenhang mit seiner Einbürgerung und Ausbreitung in Europa. Acta Theriol. (Warsz) 9, 81–110 (1964).
Gesy, K., Hill, J. E., Schwantje, H., Liccioli, S. & Jenkins, E. J. Establishment of a European-type strain of Echinococcus multilocularis in Canadian wildlife. Parasitology 140, 1133–1137 (2013).
Massolo, A. et al. European Echinococcus multilocularis Identified in Patients in Canada. N. Engl. J. Med. 381, 384–385 (2019).
Santa, M. A., Musiani, M., Ruckstuhl, K. E. & Massolo, A. A review on invasions by parasites with complex life cycles: the European strain of Echinococcus multilocularis in North America as a model. Parasitology 148, 1532–1544 (2021).
Conlon, C. L., Schuler, K. L., Lejeune, M. & Whipps, C. M. Novel Report of the European Variant of Echinococcus multilocularis in Coyotes (Canis latrans) in New York State. Journal of Parasitology 109, 357 (2023).
Conlon, C. L. et al. Echinococcus multilocularis in New York Wildlife: distribution and genetic diversity of an emerging pathogen. J. Parasitol. 110, 697 (2024).
Rogers, A. R. & Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552–569 (1992).
Tamura, K. & Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 10, 512–526 (1993).
Felsenstein, J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution (N Y). 39, 783–791 (1985).
Saitou, N. & Nei, M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Evol. Biol. https://doi.org/10.1093/oxfordjournals.molbev.a040454 (1987).
Tamura, K. et al. Estimating divergence times in large molecular phylogenies. Proc. Natl. Acad. Sci. 109, 19333–19338 (2012).
Tamura, K., Tao, Q. & Kumar, S. Theoretical foundation of the reltime method for estimating divergence times from variable evolutionary rates. Mol. Biol. Evol. 35, 1770–1782 (2018).
Nei, M. & Kumar, S. Molecular evolution and phylogenetics (Oxford University Press, 2000).
Tao, Q., Tamura, K., Mello, B. & Kumar, S. Reliable confidence intervals for reltime estimates of evolutionary divergence times. Mol. Biol. Evol. 37, 280–290 (2020).
R Core Team. R: A Language and Environment for Statistical Computing. Preprint at https://www.R-project.org/ (2024).
Acknowledgements
The authors would like to express their utmost gratitude to Adrianna Kilikowska, WSB Merito University in Gdańsk, for insightful comments and suggestions that enabled us to improve the original draft. Many thanks to Ioannis Kontogeorgos for creating Figure 1.
Funding
This work was supported by the National Science Centre grant no. 2020/37/B/NZ7/03934.
Author information
Authors and Affiliations
Contributions
P. Gładysz : Conceptualisation, Data Curation, Formal Analysis, Investigation, Methodology, Validation, Visualisation, Writing – Original Draft, Writing – Review & Editing. D. Bielińska-Wąż : Formal Analysis, Methodology, Software, Visualisation, Writing – Review & Editing. P. Wąż : Formal Analysis, Methodology, Software, Visualisation, Writing – Review & Editing. M. Samorek-Pieróg : Resources. J. Karamon : Resources, Writing – Review & Editing. A. Lass : Conceptualisation, Funding Acquisition, Project Administration, Resources, Supervision, Writing – Review & Editing.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Supplementary Information 1.
Supplementary Information 2.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissions
About this article
Cite this article
Gładysz, P., Bielińska-Wąż, D., Wąż, P. et al. Euro-Asian hybrids of Echinococcus multilocularis from red foxes in northern and northeastern Poland result from secondary contact between long-isolated populations.
Sci Rep (2026). https://doi.org/10.1038/s41598-026-40313-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-026-40313-z
Keywords
- Alignment-free method
- Alveolar echinococcosis
- Asian origin
- Phylogeography
Source: Ecology - nature.com

