in

Euro-Asian hybrids of Echinococcus multilocularis from red foxes in northern and northeastern Poland result from secondary contact between long-isolated populations


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

  1. 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).

    Google Scholar 

  2. Wen, H. et al. Echinococcosis: advances in the 21st Century. Clin. Microbiol. Rev. 32, 10 (2019).

    Google Scholar 

  3. 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).

    Google Scholar 

  4. 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).

    Google Scholar 

  5. Santoro, A. et al. Mitochondrial genetic diversity and phylogenetic relationships of Echinococcus multilocularis in Europe. Int. J. Parasitol. 54, 233–245 (2024).

    Google Scholar 

  6. 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).

    Google Scholar 

  7. 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).

    Google Scholar 

  8. Nakao, M., Sako, Y. & Ito, A. Isolation of polymorphic microsatellite loci from the tapeworm Echinococcus multilocularis. Infect. Genet. Evol. 3, 159–163 (2003).

    Google Scholar 

  9. 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).

    Google Scholar 

  10. 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).

    Google Scholar 

  11. Konyaev, S. V. et al. Genetic diversity of Echinococcus spp.in Russia. Parasitol. 140, 1637–1647 (2013).

    Google Scholar 

  12. 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).

    Google Scholar 

  13. 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).

    Google Scholar 

  14. 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).

    Google Scholar 

  15. 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).

    Google Scholar 

  16. 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).

    Google Scholar 

  17. 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).

    Google Scholar 

  18. 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).

    Google Scholar 

  19. Nakao, M. et al. Geographic pattern of genetic variation in the fox tapeworm Echinococcus multilocularis. Parasitol Int. 58, 384 (2009).

    Google Scholar 

  20. Nei, M. Molecular Evolutionary Genetics (Columbia University Press, 1987). https://doi.org/10.7312/nei-92038.

    Google Scholar 

  21. Nei, M. & Li, W. H. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. 76, 5269–5273 (1979).

    Google Scholar 

  22. Rozas, J. et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302 (2017).

    Google Scholar 

  23. Jombart, T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405 (2008).

    Google Scholar 

  24. Jombart, T. & Ahmed, I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics 27, 3070–3071 (2011).

    Google Scholar 

  25. Paradis, E. pegas: an R package for population genetics with an integrated–modular approach. Bioinformatics 26, 419–420 (2010).

    Google Scholar 

  26. R Core Team. R: A Language and environment for statistical computing. Preprint at https://www.R-project.org/ (2023).

  27. Bandelt, H. J., Forster, P. & Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).

    Google Scholar 

  28. Kumar, S. et al. MEGA12: molecular evolutionary genetic analysis version 12 for adaptive and green computing. Mol. Biol. Evol. 41, msae263 (2024).

    Google Scholar 

  29. 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).

    Google Scholar 

  30. 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).

    Google Scholar 

  31. 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).

    Google Scholar 

  32. Bielińska-Wąż, D. & Wąż, P. Non-standard bioinformatics characterization of SARS-CoV-2. Comput. Biol. Med. 131, 104247 (2021).

    Google Scholar 

  33. 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).

    Google Scholar 

  34. 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).

    Google Scholar 

  35. 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).

    Google Scholar 

  36. Avise, J. C. Phylogeography: the history and formation of species (Harvard University Press, UK, 2000).

    Google Scholar 

  37. 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).

    Google Scholar 

  38. Hayashi, N. et al. Mitogenomic exploration supports the historical hypothesis of anthropogenic diffusion of a zoonotic parasite Echinococcus multilocularis. iScience 26, 107741 (2023).

    Google Scholar 

  39. 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).

    Google Scholar 

  40. 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).

    Google Scholar 

  41. 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).

    Google Scholar 

  42. Massolo, A. et al. European Echinococcus multilocularis Identified in Patients in Canada. N. Engl. J. Med. 381, 384–385 (2019).

    Google Scholar 

  43. 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).

    Google Scholar 

  44. 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).

    Google Scholar 

  45. Conlon, C. L. et al. Echinococcus multilocularis in New York Wildlife: distribution and genetic diversity of an emerging pathogen. J. Parasitol. 110, 697 (2024).

    Google Scholar 

  46. Rogers, A. R. & Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552–569 (1992).

    Google Scholar 

  47. 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).

    Google Scholar 

  48. Felsenstein, J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution (N Y). 39, 783–791 (1985).

    Google Scholar 

  49. 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).

    Google Scholar 

  50. Tamura, K. et al. Estimating divergence times in large molecular phylogenies. Proc. Natl. Acad. Sci. 109, 19333–19338 (2012).

    Google Scholar 

  51. 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).

    Google Scholar 

  52. Nei, M. & Kumar, S. Molecular evolution and phylogenetics (Oxford University Press, 2000).

    Google Scholar 

  53. 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).

    Google Scholar 

  54. R Core Team. R: A Language and Environment for Statistical Computing. Preprint at https://www.R-project.org/ (2024).

Download references

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

Authors

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

Correspondence to
Paweł Gładysz or Anna Lass.

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

Download citation

  • 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

Trend analysis of dam inflow data using the trend accuracy index and the potential-evapotranspiration correction factor

Mapping the marine distribution of eulachon (Thaleichthys pacificus) in the Northeast Pacific using environmental DNA

Back to Top