BACKGROUND: The application of molecular-phylogenetic approaches to taxonomy has had a dramatic effect on our understanding of the diversity of reptiles. These approaches have allowed researchers to reveal previously hidden lineages as well as taxonomic overestimation in morphologically plastic taxa. Slow worms, legless lizards of the genus Anguis (Squamata: Anguidae), were previously considered to comprise either one or two species, and morphology-based intraspecific taxonomy of Anguis fragilis remained controversial throughout the 20th century. After the discovery of deep genetic divergences within the genus, its taxonomy was reconsidered, and as a result, five extant species have been recognized. In order to better understand the patterns of their interspecific differentiation, here we studied phenotypic differences between the two most widespread of them-A. fragilis and A. colchica, and their putative hybrids across the contact zone of both species in Central Europe. METHODS: We used multivariate and univariate statistics and analyzed ten metric, eleven meristic, and six categorical phenotypic variables in material comprising a total of 326 individuals. We also genotyped individuals from the contact zone for one mitochondrial and two nuclear DNA fragments in order to delineate the distribution of individuals of hybrid and non-hybrid origin. The clines in morphological traits were studied using HZAR. RESULTS: We show that the two species are morphologically differentiated. Anguis fragilis has a less robust head, fewer scales covering the body, lower frequency of the external ear opening presence, lower frequency of separated prefrontal scales, higher frequency of prefrontal scales in contact with each other, and body coloration more similar to the juvenile coloration than A. colchica. Slow worms from the contact/hybrid zone are characterized by an intermediate morphology, with more similarities to A. fragilis than to A. colchica. DISCUSSION: None of the analyzed characters alone proved to be fully diagnostic, although more than 90% of all individuals could be successfully assigned to one or another species based on numbers of scales around the body. Our results indicate concordant, coincident, and steep clines in character states change. We present several hypotheses on the origin and evolutionary maintenance of the morphological divergence between both species and suggest that different evolutionary histories of the taxa rather than recently acting selection explain the observed morphological variation.

Department of Zoology Comenius University in Bratislava Bratislava Slovak Republic

Department of Zoology National Museum Prague Czech Republic

Institute of Vertebrate Biology of the Czech Academy of Sciences Brno Czech Republic

See more in PubMed

Adnan N, Ahmad MH, Adnan R. A comparative study on some methods for handling multicollinearity problems. Matematika. 2006a;22:109–119.

Adnan N, Ahmad MH, Adnan R. A comparative study on some methods for handling multicollinearity problems. Matematika. 2006b;22:109–119.

Aguilera A, Escabias M. Functional and operatorial statistics. Physica-Verlag; Heidelberg: 2008. Solving multicollinearity in functional multinomial logit models for nominal and ordinal responses; pp. 7–13.

Arnold ML. Natural hybridization and evolution. Oxford University Press; New York: 1997.

Baker RJ, Bradley RD. Speciation in mammals and the genetic species concept. Journal of Mammalogy. 2006;87:643–662. doi: 10.1644/06-MAMM-F-038R2.1. PubMed DOI PMC

Baran I. Türkiye’de Anguidae familyasi türlerinin taksonomisi. Ege Üniversitesi Fen Fakültesi Dergisi Series B, C.I. S.2. 1977;1977:145–153.

Barton NH, Hewitt GM. Analysis of hybrid zones. Annual Review of Ecology and Systematcs. 1985;16:113–148. doi: 10.1146/annurev.es.16.110185.000553. DOI

Bergmann PJ, Morinaga G. The convergent evolution of snake-like forms by divergent evolutionary pathways in squamate reptiles. Evolution. 2019;73:481–496. doi: 10.1111/evo.13651. PubMed DOI

Beschkov V. Untersuchungen über Systematik und Verbreitung der Blindschleiche (Anguis fragilis L.) in Bulgarien. Bulletin de L’Institut De Zoologie Et Musée (Sofia) 1966;21:185–201. [in Bulgarian with German summary]

Blosat B. Morphometrische und ökologische Feldstudien an Reptilien im Bergischen Land (Nordrhein-Westfalen) –I. Blindschleiche (Anguis f. fragilisLinnaeus, 1758) Salamandra. 1997;33:161–174.

Böhme W. Notes on the combat behaviour of male slow worms, Anguis fragilis. Zeitschrift Für Feldherpetologie. 2006;13:240–244.

Cabela A, Grillitsch H. Zum systematischen Status der Blindschleiche (Anguis fragilis Linnaeus, 1758) von Nordgriechenland und Albanien (Squamata: Anguidae) Herpetozoa. 1989;2:51–69.

Calsbeek R, Knouft JH, Smith TB. Variation in scale numbers is consistent with ecologically based natural selection acting within and between lizard species. Evolutionary Ecology. 2006;20:377–394. doi: 10.1007/s10682-006-0007-y. DOI

Capula M, Anibaldi C, Filippi E, Luiselli L. Sexual combats, matings, and reproductive phenology in an Alpine population of the slow worm, Anguis fragilis. Herpetological Natural History. 1998;6:33–39.

Capula M, Luiselli L, Capanna E. The blue-spotted morph of the slow worm, Anguis fragilis: colour polymorphism and predation risks. Italian Journal of Zoology. 1997;64:147–153. doi: 10.1080/11250009709356188. DOI

Carranza S, Arnold EN. A review of the geckos of the genus Hemidactylus (Squamata: Gekkonidae) from Oman based on morphology, mitochondrial and nuclear data, with descriptions of eight new species. Zootaxa. 2012;3378:1–95. doi: 10.11646/zootaxa.3378.1.1. DOI

Costello AB, Osborne JW. Best practices in exploratory factor analysis: four recommendations for getting the most from your analysis. Practical Assessment, Research & Evaluation. 2005;10:1–9.

D’Ambra A, Sarnacchiaro P. Some data reduction methods to analyze the dependence with highly collinear variables: a simulation study. Asian Journal of Mathematics and Statistics. 2010;3:69–81. doi: 10.3923/ajms.2010.69.81. DOI

Dely OG. Über die Unterarten der Blindschleiche. Anguis fragilis L. Vertebrata Hungarica. 1974;15:11–37.

Dely OG. Anguis fragilis Linnaeus 1758 –Blindschleiche. In: Böhme W, editor. Handbuch der Reptilien und Amphibien Europas. Band 1. Echsen (Sauria) 1. AULA-Verlag; Wiesbaden: 1981. pp. 241–258.

Derryberry EP, Derryberry GE, Maley JM, Brumfield RT. hzar: hybrid zone analysis using an R software package. Molecular Ecology Resources. 2014;14:652–663. doi: 10.1111/1755-0998.12209. PubMed DOI

Dixon JR. A systematic review of the teiid lizards, genus Bachia, with remarks on Heterodactylus and Anotosaura. University of Kansas publications, Museum of Natural History. 1973;57:1–47.

Džukić G. Taxonomic and biogeographic characteristics of the slow-worm (Anguis fragilis Linnaeus 1758) in Yugoslavia and on the Balcan Peninsula. Scopolia. 1987;12:1–47.

Everit BS, Dunn G. Applied multivariate data analysis (2nd edition) John Wiley & Sons Ltd; West Sussex: 2001.

Galis F, Arntzen JW, Lande R. Dollo’s law and the irreversibility of digit loss in Bachia. Evolution. 2010;64:2466–2476. PubMed

Games PA. Curvilinear transformations of the dependent variable. Psychological Bulletin. 1983;93:382–387. doi: 10.1037/0033-2909.93.2.382. DOI

Graham JM, Guthrie AC, Thompson B. Consequences of not interpreting structure coefficients in published CFA research: a reminder. Structural Equation Modeling. 2003;10:142–153. doi: 10.1207/S15328007SEM1001_7. DOI

Greer AE. The Loss of the external ear opening in scincid lizards. Journal of Herpetology. 2002;36:544–555. doi: 10.1670/0022-1511(2002)036[0544:TLOTEE]2.0.CO;2. DOI

Grillitsch H, Cabela A. Zum systematischen Status der Blindschleichen (Squamata: Anguidae) der Peloponnes und der südlichen Ionischen Inseln (Griechenland) Herpetozoa. 1990;2:131–153.

Groth JG, Barrowclough GF. Basal divergences in birds and the phylogenetic utility of the nuclear RAG-1 gene. Molecular Phylogenetics and Evolution. 1999;12:115–123. doi: 10.1006/mpev.1998.0603. PubMed DOI

Günther A. Description of a new Saurian (Hyalosaurus) allied to Pseudopus. Annals and Magazine of Natural History (4) 1873;11:351.

Gvoždík L, Van Damme R. Evolutionary maintenance of sexual dimorphism in head size in the lizard Zootoca vivipara: a test of two hypotheses. Journal of Zoology. 2003;259:7–13. doi: 10.1017/S0952836902003308. DOI

Gvoždík V, Benkovský N, Crottini A, Bellati A, Moravec J, Romano A, Sacchi R, Jandzik D. An ancient lineage of slow worms, genus Anguis (Squamata: Anguidae), survived in the Italian Peninsula. Molecular Phylogenetics and Evolution. 2013;69:1077–1092. doi: 10.1016/j.ympev.2013.05.004. PubMed DOI

Gvoždík V, Harca Z, Hánová A, Jablonski D, Pupins M, Čeirãns A, Paasikunnas T. Two species of slow worm (Anguis fragilis, A. colchica) present in the Baltic region. Amphibia-Reptilia. 2021;42:383–389. doi: 10.1163/15685381-bja10055. DOI

Gvoždík V, Jandzik D, Lymberakis P, Jablonski D, Moravec J. Slow worm, Anguis fragilis (Reptilia: Anguidae) as a species complex: Genetic structure reveals deep divergences. Molecular Phylogenetics and Evolution. 2010;55:460–472. doi: 10.1016/j.ympev.2010.01.007. PubMed DOI

Gvoždík V, Moravec J. Anguis fragilis Linnaeus, 1758–slepýš křehký [Anguis fragilis Linnaeus, 1758–Common European slow worm] In: Moravec J, editor. Fauna ČR, Plazi–Reptilia [Fauna of the Czech Republic, Reptiles –Reptilia] Praha: Academia; 2015. pp. 237–261. [in Czech with an English summary]

Gvoždík V, Moravec J, Zavadil V, Jeřábková L. Slepýš křehký a slepýš východní –výskyt v České republice [Common European slow worm and Eastern slow worm –distribution in the Czech Republic] In: Moravec J, editor. Fauna ČR, Plazi –Reptilia [Fauna of the Czech Republic, Reptiles –Reptilia] Praha: Academia; 2015. pp. 275–278. [in Czech with an English summary]

Hawkins DM. The problem of overfitting. Journal of Chemical Information and Computer Sciences. 2004;44:1–12. doi: 10.1021/ci0342472. PubMed DOI

Henson RK, Roberts JK. Use of exploratory facor analysis in published research: common errors and some comment on improved practice. Educational and Psychological Measurement. 2006;66:393–416. doi: 10.1177/0013164405282485. DOI

Herrel AK, Huyghe B, Vanhooydonck T, Backeljau K, Breugelmans I, Grbac R, Van Damme, et al. Rapid large-scale evolutionary divergence in morphology and performance associated with exploitation of a different dietary resource. Proceedings of the National Academy of Sciences of the United States of America. 2008;105(12):4792–4795. doi: 10.1073/pnas.0711998105. PubMed DOI PMC

Jablonski D, Jandzik D, Mikulíček P, Džukić G, Ljubisavljević K, Tzankov N, Jelić D, Thanou E, Moravec J, Gvoždík V. Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula. BMC Evolutionary Biology. 2016;16:99. doi: 10.1186/s12862-016-0669-1. PubMed DOI PMC

Jablonski D, Najbar B, Grochowalska R, Gvoždík V, Strzała T. Phylogeography and postglacial colonization of Central Europe by Anguis fragilis and Anguis colchica. Amphibia-Reptilia. 2017;38:562–569. doi: 10.1163/15685381-00003133. DOI

Jablonski D, Sillero N, Oskyrko O, Bellati A, Čeirãns A, Cheylan K, Cogălniceanu D, Crnobrnja-Isailović J, Crochet P-A, Crottini A, Doronin I, Džukić G, Geniez P, Ç Ilgaz, Iosif R, Jandzik D, Jelić D, Litvinchuk S, Ljubisavljević K, Lymberakis P, Mikulíček P, Mizsei E, Moravec J, Najbar B, Pabijan M, Pupins M, Sourrouille P, Strachinis I, Szabolcs M, Thanou E, Tzoras E, Vergilov V, Vörös J, Gvoždík V. The distribution and biogeography of slow worms (Anguis, Squamata) across the Western Palearctic, with an emphasis on secondary contact zones. Amphibia-Reptilia. 2021;42:519–530. doi: 10.1163/15685381-bja10069. DOI

Jennrich RI. A simple general method for oblique rotation. Psychometrika. 2002;67:7–20. doi: 10.1007/BF02294706. DOI

Kaliantzopoulou A, Pinho C, Harris JD, Carretero MA. When cryptic diversity blurs the picture: a cautionary tale from Iberian and North African Podarcis wall lizards. Biological Journal of the Linnean Society. 2011;103:779–800. doi: 10.1111/j.1095-8312.2011.01703.x. DOI

Keating JP, Hensley OL. A note on the critical values of the f-max test. Communications in Statistics –Simulation and Computation. 1983;12:257–263. doi: 10.1080/03610918308812317. DOI

Lác J. K systematike slepúcha lámavého (Anguis fragilis L.) a jeho rozšírenie na Slovensku. Biológia. 1967;22:908–921.

Lance RF, Kennedy ML, Leberg PL. Classification bias in discriminant function analysis used to evaluate putatively different taxa. Journal of Mammalogy. 2000;81:245–249. doi: 10.1644/1545-1542(2000)081<0245:CBIDFA>2.0.CO;2. DOI

Lavin BR, Girman DJ. Phylogenetic relationships and divergence dating in the Glass Lizards (Anguinae) Molecular Phylogenetics and Evolution. 2019;133:128–140. doi: 10.1016/j.ympev.2018.12.022. PubMed DOI

Macey RJ, Schulte IIJA, Larson A, Tuniyev BS, Orlov N, Papenfuss TJ. Molecular phylogenetics, tRNA evolution, and historical biogeography in anguid lizards and related taxonomic families. Molecular Phylogenetics and Evolution. 1999;12:250–272. doi: 10.1006/mpev.1999.0615. PubMed DOI

Mayer W, Grillitsch H, Cabela A. Proteinelektrophoretische Untersuchungen zur Systematik der südgriechischen Blindschleiche (Squamata: Anguidae) Herpetozoa. 1991;4:157–165.

McLachlan GJ. Discriminant analysis and statistical pattern recognition. John Wiley & Sons, Inc; Hoboken: 2004.

Mikulíček P, Jandzik D, Fritz U, Schneider C, Široký P. AFLP analysis shows high incongruence between genetic differentiation and morphology-based taxonomy in a widely distributed tortoise. Biological Journal of the Linnean Society. 2013;108:151–160. doi: 10.1111/j.1095-8312.2012.01999.x. DOI

Moravec J, Gvoždík V. Anguis colchica (Nordmann, 1840) –slepýš východní [Anguis colchica (Nordmann, 1840) –Eastern slow worm] In: Moravec J, editor. Fauna ČR, Plazi –Reptilia [Fauna of the Czech Republic, Reptiles –Reptilia] Praha: Academia; 2015. pp. 262–275. [in Czech with an English summary]

Musters CJM, In den Bosch HAJ. Einige Bemerkungen zu den Unterarten von Anguis fragilis L. mit berucksichtigung niederlandisher exemplare (Reptilia: Sauria: Anguidae) Salamandra. 1982;18:196–204.

Oufiero CE, Gartner GEA, Adolph SC, Garland Jr T. Latitudinal and climatic variation in body size and dorsal scale counts in Sceloporus lizards: a phylogenetic perspective. Evolution. 2011;65:3590–3607. doi: 10.1111/j.1558-5646.2011.01405.x. PubMed DOI

Parizanganeh A, Lakhan VC, Yazdani M, Ahmad SR. Modelling categorical data to identify factors influencing concern for the natural environment in Iran. Journal of Environmental Management. 2011;92:2836–2843. PubMed

Petzold H-G. Blindschleiche und Scheltopusik. Die Familie Anguidae. Die Neue Brehm-Bücherei 448. Wittenberg-Lutherstadt: A. Ziemsen; 1971.

Pyron RA, Burbrink FT, Wiens JJ. A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. BMC Evolutionary Biology. 2013;13:1–53. doi: 10.1186/1471-2148-13-1. PubMed DOI PMC

R-Core-Team A language and environment for statistical computing. http://www.R-project.org/ 2020

R Studio Team . RStudio, Inc. Boston, MA; 2015.

Rafter JA, Abell ML, Braselton JP. Multiple comparison methods for means. SIAM Review. 2002;44:259–278. doi: 10.1137/S0036144501357233. DOI

Ščerbak NN, Ščerbaň MI. Zemnovodnye i presmykayushchiesya Ukrainskikh Karpat. Naukova Dumka; Kiev: 1980.

Ščerbaň MI. K vnitrovidovoj sistematike veretenicy lomkoj (Reptilia, Sauria, Anguis fragilis L.) Sbornik Trudov Zoologicheskogo Muzeia (Kiev) 1976;36:81–83.

Shine R. Ecological causes for sexual dimorphism: a review of the evidence. The Quarterly Review of Biology. 1989;64:419–461. doi: 10.1086/416458. PubMed DOI

Siler CD, Brown RM. Evidence for repeated acquisition and loss of complex body-form characters in an insular clade of Southeast Asian semi-fossorial skinks. Evolution. 2011;65:2641–2663. doi: 10.1111/j.1558-5646.2011.01315.x. PubMed DOI

Simpson GG. The species concept. Evolution. 1951;5:285–298. doi: 10.1111/j.1558-5646.1951.tb02788.x. DOI

Sos T. Evaluating the accuracy of morphological traits used in Anguis (sub)species differentiation. Herpetologica Romanica. 2010;4:29–44.

Sos T. Spot polymorphism in Anguis colchica Nordmann, 1840 (Reptilia: Anguidae): inter-size class variation. North-Western Journal of Zoology. 2011;7:171–175.

Sos T, Herczeg G. Sexual size dimorphism in eastern slow-worm (Anguis fragilis colchica, Reptilia: Anguidae) Russian Journal of Herpetology. 2009;16:304–310.

Soulé M, Kerfoot CW. On the climatic determination of scale size in a lizard. Systematic Zoology. 1972;21:97–105. doi: 10.2307/2412262. DOI

Speybroeck J, Beukema W, Dufresnes C, Fritz U, Jablonski D, Lymberakis P, Martínez-Solano I, Razzetti E, Vamberger M, Vences M, Vörös J, Crochet P-A. Species list of the European herpetofauna –2020 update by the Taxonomic Committee of the Societas Europaea Herpetologica. Amphibia-Reptilia. 2020;41:139–189. doi: 10.1163/15685381-bja10010. DOI

Stugren B, Fuhn J, Popovici N. Untersuchungen über die Systematik des Blindschleiches (Anguis fragilis L.) in Rumänien. Zoologischer Anzeiger. 1962;169:460–466.

Szabó K, Vörös J. Distribution and hybridization of Anguis fragilis and A. colchica in Hungary. Amphibia-Reptilia. 2014;35:135–140. doi: 10.1163/15685381-00002927. DOI

Tabachnick BG, Fidell LS. Using multivariate statistics (5th Edition) Allyn & Bacon; Boston, Pearson: 2007.

Thorpe RS, Baez M. Geographic variation within an island: univariate and multivariate contouring of scalation, size, and shape of the lizard Gallotia galloti. Evolution. 1987;41:256–268. doi: 10.1111/j.1558-5646.1987.tb05795.x. PubMed DOI

Valakos ED, Pafilis P, Sotiropoulos K, Lymberakis P, Marangou P, Foufopoulos J. Chimaira Publications; Frankfurt am Mainz: Reptiles and Amphibians of Greece; p. 463.

Voipio P. Multiple phaneromorphism in the European slow-worm (Anguis fragilis) and the distributional and evolutionary history of the species. Annales Zoologici Societatis Zoologicae-Botanicae Fennicae Vanamo. 1962;23:1–20.

Völkl W, Alfermann D. Zeitschrift für Feldherpetologie, Beiheft 11. Bielefeld: Laurenti-Verlag; 2007. Die Blindschleiche–die vergessene Echse.

Wermuth H. Variationsstatistische Untersuchung der Rassen- und Geschlechtsmerkmale bei der Blindschleiche (Anguis fragilis Linné) Deutsche Zoologische Zeitschrift. 1950;1:81–121.

Wiens JJ, Slingluff JL. How lizards turn into snakes: aphylogenetic analysis of body-form evolution in anguid lizards. Evolution. 2001;55:2303–2318. doi: 10.1111/j.0014-3820.2001.tb00744.x. PubMed DOI

Wiley EO. Evolutionary species concept reconsidered. Systematic Zoology. 1978;27:17–26. doi: 10.2307/2412809. DOI

Winne CT, Willson JD, Andrews KM, Reed RN. Efficacy of marking snakes with disposable medical cautery units. Herpetological Review. 2006;37:52–54.

Zadravec M, Golub A. Longer than the longest –two record-breaking specimens of Anguis fragilis from north-western Croatia. Zeitschrift Für Feldherpetologie. 2018;25:105–107.

Contact zone of slow worms Anguis fragilis Linnaeus, 1758 and Anguis colchica (Nordmann, 1840) in Poland