NORTH-WESTERN JOURNAL OF ZOOLOGY 13 (1): 186-191 ©NwjZ, Oradea, Romania, 2017 Article No.: e179501 http://biozoojournals.ro/nwjz/index.html

In the Memory of Andrey STOYANOV & Nikolay TZANKOV:

Individual variation of pileus scalation characteristics in bosniensis Boettger, 1889 (Reptilia: Sqaumata: )

Andrey STOYANOV and Nikolay TZANKOV

National Museum of Natural History, Bulgarian Academy of Sciences, 1 Tsar Osvoboditel Blvd., 1000 Sofia, Bulgaria.

Received: 31. August 2015 / Accepted: 30. March 2017 / Available online: 30. March 2017 / Printed: June 2017

Editorial justification. This paper was submitted to evaluation on 31st August 2015. After the start of the evaluation process, we found out from an that both of the authors died in a car accident, during execution of herpetological studies in northern Greece (Acta Zoologica Bulgarica 68: 453-453). The editors of North-Western Journal of Zoology and the two anonymous referees agreed to publish the manuscript in its last stage in the memory of the authors.

Abstract. In the present study of mountainous populations of Vipera berus bosniensis a total of 122 individuals were photo-documented for future identification. When examining the material a set of anomalous manifestations and deviations from the standard shield arrangement were observed. Additionally 37 newborn juveniles were examined and compared with adults. The prevalence of different types of shield deviations as well as percentage of specimens with different types of arrangements were analyzed. Samples were subjected to analyses with the aim of elucidating any presumptive differences in scalation between age and sex groups. Results indicated a similar pattern of group frequencies among sexes as well among age groups. In both adults and juvenile , individuals with deviations from the accepted normal configuration of shields were prevalent.

Key words: ontogenetic change, Balkan adder, scutellation, variability, Vitosha Mt.

Introduction fiable scales. Several factors hinder a definite con- clusion, taking into account that in some Shield arrangements of the dorsal aspect of the groups (such as Boidae, Pythonidae and Viperi- head, i.e. the pileus region in , have been dae) exhibit a large variability of pileus scalation. accepted as a group, respectively species and in- In a comprehensive study. Werner (1898) regarded dividually specific diagnostic traits (Werner 1899, taxa with large shields as evolutionarily younger Schreiber 1912, Ursel 1978, Benson 1999, Krupitz and subsequently, more specialized forms evolved 1999, Schwarz 1997). In a set of publications de- to form increasingly fragmented small scales. In scribing the adder Vipera berus (Linnaeus, 1758), contrast, the converse hypothesis, with large specific schemes that illustrated the commonly ac- shields stated to be the evolutionarily more recent cepted standard: number, size and arrangement of and progressive type, was purported by Zacharias pileus shields have been mentioned (Schreiber (1898) concerning members of Boidae and Py- 1912, Schiemenz 1987, Benson 1999, Toth et al. thonidae. According to Ursel (1978), in whose 2005). Examples of individual deviations from the opinion head scale evolution conformed to the lat- typical arrangement have been presented, albeit ter model in viperids and crotalids, the fragmenta- without establishing causal inferences (Schreiber tion is a secondary phenomenon, the reasons for 1912, Saint Girons 1978, Benson 1999, Biella 1988, which remain unknown. The proposition was Biella 1990). Saint Girons (1978) and Schiemenz about the range of attachment, facilitated by the (1987) presented differences in pholidosis among head skin extensibility in this type of scalation. the subspecies of V. berus. It is known that the formation of the individ- To a large extent still, the issue about the pri- ual pattern, number and arrangement of head mary type of head scalation is open to debate. shields and scales takes place during late embryo- Both types include the oligomerization, i.e. the re- genesis (Lorioux et al. 2013). Of particular interest duction of scales to a small set of large shields, or are the postnatal changes in pileus scalation, such the fragmentation into a large number of unidenti- as has been rarely observed in different snake Individual variation in Balkan adder pileus scutellation 187 families (Werner 1899, Tomović et al. 2008; frontal (FFr); fragmentation of parietals (FPa); fragmenta- Hodges & Seabrook 2014). tion of interoculars (FInoc); reduction of intercanthals The relationship between changes in a com- (RInca); reduction of interoculars (RInoc); reduction of canthals (RCa); reduction of parietals (RPa). The speci- plete body scalation of snakes has been found to mens with mixed characters, when at the same time both be dependent on thermoregulation of individuals, types of deviation were found, were also mentioned evaporation of body humidity and climatic condi- (MCh). Additionally, all shields with incomplete frag- tions and habitat characteristics in species of the mentation (starting sulci) in different degrees of fusion genus Echis (Cherlin 1983) and in Cerastes vipera were noted, where there were incompletely closed sulci. (Jooris & Fourmy 1996). This would certainly have Eleven types of such changes were defined. Partly sepa- a physiological effect when it comes to a large rated frontal (PSFr); partly separated parietals (PSPa); partly separated interoculars (PSInoc); partly separated body surface, i.e. concerning changes in whole- intercanth (PSInca); partly separated supraoculars body scalation in either direction. Regarding only (PSSuoc); partly separated canthals (PSCa); fusion of in- the insignificant total area covered by the pileus, tercanthals (FuInca); fusion of canthals (FuCa); fusion of the hypothesis about the increased head kinesis, as interoculars (FuInoc); fusion of frontal and interoculars proposed by Ursel (1978) seems to be more plausi- (FuFr+Inoc). ble. All elements situated posteriorly of an imaginary line The present study aims to characterize the pil- denoting the anterior edges of the supraoculars were ac- cepted to be part of the frontal formed due to a fragmen- eus shield arrangement in an open population of tation (Fig. 1). Elements closed between this line and can- V. berus bosniensis Boettger, 1889 with respect to thals, respectively apicals were defined to be intercan- the standardized model as defined for example by thals. Commonly occurring fragmentation at the end of Schreiber (1912) and Toth et al. (2005) of this re- the acute finishing[?] frontal in a caudal direction, in all gion as well as in regards to the individual cases was treated as a part of this element. As separate uniqueness of number, dimension, form and ar- from the frontal we also accepted those elements that rangement of shields. Additionally, the question were forming a second row of interoculars. about the possible effects of age and sex were tested.

Materials and methods

Field work was carried out from 2010 to 2014 in the higher northern part of Vitosha Mt. near Sofia city. The study polygon is situated between 1750 and 1850 m a.s.l., covering approximately 1.7 km2. A total of 122 specimens (56 males and 66 females) was captured and photo- documented individually. A smaller set of 8 specimens comes from another open subpopulation situated in the southeastern part of Vitosha Mt. located between 1450- 1500 m a.s.l. Both source territories belong to one meta- population. Pregnant females (6) when caught in late Au- gust and early September were kept in captivity until they gave birth, and then were released at their respective point of capture. In this manner, data from 37 newborn (17 male and 20 females) were obtained. In order to establish individual characteristics of the pileus as well as the overall deviation from the standard scutellation, we accepted one particular type as number, Figure 1. Example of standard shields arrangement of the size and arrangement of shields as the standard to which pileus: small open circles – apicals; small open squares – the corresponding changes were registered. Based on canthals; large open pentagons – intercanthals; small head form and the symmetrical arrangement of scales, we filed pentagons – interoculars; large open squares – su- followed the model proposed by Toth et al. (2005). Similar praoculars; large filled circles – parietals; star – frontal. schemes were given by Schreiber (1912), where the in- Line denotes the anterior edges of supraoculars. terocular scales were with one more. For every caught specimen a set of 9 types of changes (fragmentation, re- spectively reduction and fusion of the main elements) of Four groups were designated containing: adult fe- head scutellation in a so-called pileus region were traced. males (AdF), adult males (AdM), juvenile females (JuvF) Fragmentation of intercanthals (FInca); fragmentation of and juvenile males (JuvM), respectively. The presence of 188 A. Stoyanov & N. Tzankov the above mentioned changes was analyzed with respect Results to individual characteristics in the number of the com- pleted fragmentation, and reductions respectively, while All of defined 19 types of changes in scutellation partial changes in the shields-started fusions or fragmen- were found to be presented both individually and tations were similarly analyzed. Differences were tested with Welch Anova (data were with heteroscedastic distri- most frequently in various combinations among bution) set at significance p<0.05 level. For post-hoc different individuals. Some of the changes are pre- analyses Tukey’s honest significant difference test (THSD) sented in Fig 2. The predominant types of changes for unequal sample size was performed. When comparing in all studied groups concerned frontal and parie- the qualitative states between groups a Chi square test tal shields (Table 1). was performed. Statistical procedures were performed Concerning the fragmentation respectively, with Statistica (StatSoft, Inc. 2011) software.

Figure 2. Examples of different type of changes in scutellation of pileus in V. b. bosniensis. (A) fragmentation of frontal with the formation of double row of intraoculars; (B) fragmenta- tion of frontal (anteriorly) and fusion of parietals to one shield; (C) mixed characters: reduction of apicals to one shield, fragmentation of parietals; (D) partial reduction of frontal and full re- duction of parietals; (E) reduction of intercanthals and interoculars; (F) fragmentation of frontal and parietals but fusion of interoculars; (G) fragmentation in most of elements of the pileus; (H) reduction in interparetals and partial fragmentation of parietals; (I) nearly compete frag- mentation of pileus elements. Individual variation in Balkan adder pileus scutellation 189

Table 1. Frequency of occurrence [%] of different types of dency for adults was opposite and males had changes, fragmentation respectively reduction and fu- higher scores (Fig 4a,b). sion of the main elements, incomplete fragmentation (starting sulci) as well as terminated fusion among the four designated groups. (n.p.) not presented, for other abbreviations see material and methods. Sum per group and per type of changes gives more them 100% because of the occurrence of more than one type per specimen.

JuvM JuvF AdM AdF (n=17) (n=20) (n=56) (n=66) FInca 18 45 14 39 FFr 65 90 75 70 FPa 29 30 23 65 FInoc 18 60 5 11 RInca 47 20 25 21 RInoc 24 5 29 26 RCa 24 15 13 9 RPa 6 5 n.p. 2 MCh 53 40 38 39 PSFr 24 25 30 41 PSPa 68 50 89 79 PSInoc n.p. 9 18 14 PSSuoc n.p. n.p. 2 3 PSCa n.p. n.p. 9 2

PSInca n.p. n.p. 11 3 FuInca 4 n.p. 2 n.p. Figure 3. Frequency of occurrence of the number FuCa 4 9 9 5 of fragmentations respectively reductions and FuInoc n.p. n.p. 7 n.p. fusions per specimens in adults (A) and juve- FuFe+Inoc n.p. 6 2 2 niles (B). Light gray bars – males, black – females, white – juvenile males, dark gray – juvenile females. reduction and fusion of the main elements, all four designated groups exhibited similar tendencies of distribution, but for FPa the frequency of occur- rence was nearly three times higher in adult fe- males than the adult males. Juveniles females had similarly higher scores of FInca. Statistically sig- nificant differences were found among adults (p=0.027) and among female age classes (p=0.002). When stating the incomplete fragmentations, the observed overall distribution among groups was similar but a set of traits were not presented in juveniles (PSSuoc, PSCa, PSInca) or it was with lower scores (PSFr, PSPa) (Table 1). These results could be interpreted as an indirect evidence for the presence of individual changes that appeared during their growth. Differences in numbers of fusions per speci- men when comparing both sexes were statistically significant (F Welch(3,155)=12.95, p=0.000003) for both adults (pTHSD=0.0002) and juveniles (pTHSD=0.005) (Fig. 3a,b). In adults, females were more fre- Figure 4. Frequency of occurrence of the number quently associated with higher scores. No signifi- of incomplete fragmentations per specimens in cant differences were found between age groups adults (A) and juvenils (B). Light gray bars – and when concerning the number of incomplete males, black – females, white – juvenile males, fragmentation (results not shown), but the ten- dark gray – juvenile females. 190 A. Stoyanov & N. Tzankov

Discussion hypothesis is the partial instability in the pileus region in V. berus was expressed in fragmentation Some differences between the results presented started with grooving to complete splitting of one here and some other similar studies (Schreiber of the elements of the interocular shields in a (1912, Toth et al. 2005, Korosov 2008) were found, traced individual that happened within one sea- but the reason for this is mainly due to method- son observed by Hodges & Seabrook (2014). A set ologies, and their respective interpretations. In this of such individual postnatal changes in the pileus study a more stringently defined model of the region of V. ursinii macrops were reported by To- number and arrangement of shields in the pileus mović et al. (2008). Our current knowledge on ad- region served as a basis of comparison against der populations in Bulgaria showed that this type which changes/deviations were recorded. In gen- of postnatal change should be considered to be an eral, the observed changes among different shield extremely rare abnormality. For example, it could groups were expressed to various extent. The most be present in a recaptured female specimen after common occurrence was the fragmentation of the the third year from the initial capture when a pit largest elements of the pileus, such as frontal and tag was implanted. In such cases no trace of parietal shields. The observed statistically signifi- changes was observed (Fig. 5). It must be pointed cant differences between sexes for some elements that in other similar researches in recaptured indi- deserve further studies. Similar results and ten- viduals no traces of changes have been observed dencies have been observed for the juveniles but despite the relatively large sample size, as exam- the increased values for certain parameters must ple in Korosov (2008) (n=3570), Hodges & Sea- be considered with caution because of the smaller brook (n=594) for V. berus and Üveges et al. (2012) sample size. Generally the pileus elements show for V. ursinii rakosiensis (n=430). much less pronounced tendency toward fragmen- tation. Elements with smaller degrees of reduction were the canthals and the parietals and with re- gard to fragmentation - the canthals, respectively apical shield. However, in most cases (38,52% for adults and 45,95% for the juveniles), specimens with both types of changes have been registered. Although sampling originated from a large free-ranging population, the results show an ex- tremely high percentage of individualism in re- spect to model pileus arrangement, as in adults it was close to 96% and for the juveniles it was 97.3%. Similar high values of individual variation in pileus scalation of adders were already men- tioned from different regions of the species range and proposed as a useful tool for individual rec- ognition (Biella 1988, Krupitz 2009, Sheldon & Bradley 1989).

When tracing the incomplete fragmentations (starting sulci) as well as terminated fusions where Figure 5. Comparison of head scutellation in female there were still incompletely closed sulci, the fol- specimen first caught in 2011 (A) and then recaptured in 2013 (B). lowing peculiarities have been observed: all of the designated changes were presented in adults (fe- males being exception only for some of the rarest – The phylogenetic importance of this trait was FuInca and FuInoc); in juveniles some of the suggested by Marx et al. (1988) while directional changes were completely missing (PSInca, PSSuoc, line as proposed by them was from large scales in PSCa and RInoc). The observed differences be- ancestral state (presented in genus Causus for ex- tween the age groups could be discussed in terms ample) thought partial to complete fragmentation of smaller sample size in juveniles but also the of pileus region. It could be also well traced later ontogenetic initiation of the mentioned among the main clades within group changes cannot be excluded. In support of this based on group phylogeny as presented by Pyron Individual variation in Balkan adder pileus scutellation 191 et al. (2013). This character seems to be also well Hodges, R., Seabrook K. (2014): Head-scale instability and the apparent heritability of a head scale anomaly in the northern related to climatic zonation, and viperid taxa dis- viper (Vipera berus). The Herpetological Bulletin 130: 16-17. tributed in tropical and subtropical regions in Jooris, R., Fourmy R. (1996): An analysis of scutellation in general have more fragmented head scales. How- populations of Cerastes vipera (Linnaeus, 1758): Scale characters co-vary with environmental temperature. African Journal of ever this directional trend was observed also Herpetology 45: 59-67. within closely related tropical genera, as example Korosov, A.V. (2008): The dynamics of an island Vipera berus in the circumrostral region scalation reduction in population. Zhoologicheskii Journal 87: 1235-1249. (in Russian) Krupitz, W. (2009): Raumnutzung männlicher Kreuzottern (Vipera African forest-dwelling arboreal viperids referred berus L.) während der Paarungszeit. Diplomarbeit. 73 p. to tribe Atherini (Broadley 1996). Within the berus Lorioux, S., Vaugoyeau, M., DeNardo, D.F., Clobert, J., Guillon, M., species group there is also well-expressed ten- Lourdais O. (2013): Stage Dependence of Phenotypical and Phenological Maternal Effects: Insight into Squamate dency of fragmentation of head scalation toward Reproductive Strategies. The American Naturalist 182: 223-233. the southern edge of areal distribution and in Marx, H., Ashe, J.S., Watrou, L.E. (1988): Phylogeny of viperine lower altitude, as for example in taxa V. barani, V. snakes (Viperinae): Part I. Character analysis. Fieldiana Zoology (N.S.) 51: 1-16. b. bosniensis, V. nikolskii and in an extinct popula- Milto, K.D., Zinenko, O.I. (2005): Distribution and morphological tion from Po River valley in northern Apennine variability of Vipera berus in Eastern Europe. pp. 64-73. In: Peninsula (Scali & Gentilli 1999, Franzen & Heckes Ananyeva, N., Tsinenko O. (eds.), Herpetologia Petropolitana: Proceedings of 12th Ord. Gen. Meet. S.E.H., St. Petersburg– 2000, Milto & Zinenko 2005). Moscow. Some of the here-presented observed changes, Pyron, R.A., Burbrink F.T., Wiens J.J. (2013): A phylogeny and especially the sexual and age differences in pileus revised classification of , including 4161 species of lizards and snakes. BMC Evolutionary Biology 13: 93. scalation, definitively require future investigation Saint Girons, H. (1978): Morphologie externe comparėe et in order to reveal the reasons and driving forces of systematique des Viperes d’Europe (Reptilia, Viperidae). Revue these processes. For such study a larger set of in- suisse de Zoologie 85: 565-595. Scali, S., Gentilli A. (1999): Morphometric analysis and sexual dividually traced specimens is required. Addi- dimorphism of extinct adders (Vipera berus) of the Po plane tionally similar researches on other viperid taxa (Northern Italy). pp. 391-396. In: Miaud C., Guyetant R. (eds.), will provide the opportunity to designate the phy- Current Studies in Herpetology. Proc. of the 9th Ordinary General Meeting of the Soc. Eur. Herpetol., Soc. Eur. Herpetol., logenetic value of the observed processes and their Le Bourget du Lac. evolutionary direction. Schiemenz, H. (1987): Die Kreuzotter: Vipera berus. Die Neue Brehm-Bücherei 332. A. Ziemsen Verlag, Wittenberg Lutherstadt. Schreiber, E. (1912): Herpetologia europaea. Jena Verlag von Gustav Fischer. Acknowledgments (by the authors). We thank Schwarz, A. (1997): Möglichkeiten der Ermittlung von Alexander Westerström for the English improvement and Raumnutzung und Populationsdichte bei der Kreuzotter (Vipera valuable suggestions. b. berus L.). Mertensiella 7: 247-260. Acknowledgments (by the editors). We want to thank Sheldon, S., Bradley C. (1989): Identification of individual adders the two anonymous referees for the collaboration and (Vipera berus) by their head markings. Herpetological Journal 1: 392-396. Steve Anderson for the spelling editing of the manuscript. StatSoft, Inc. (2011): Statistica (data analysis software system), version 10, Tomović, L., Carretero, M.A., Ajtic, R., Crnobrnja-Isailovic, J. (2008): Evidence for post-natal instability of head scalation in the References meadow viper (Vipera ursinii) – patterns and taxonomic Benson, P. (1999): Identifying individual adders, Vipera berus, implications. Amphibia-Reptilia 29: 61-70. within an isolated colony in east Yorkshire. British Tóth, T., Covács, T., Herczeg, G. (2005): Morphological anomalies in Herpetological Society Bulletin 67: 21-27. a finnish subarctic population of the adder (Vipera berus berus). Biella, H.J. (1988): Untersuchung zur Varistion der Litteratura serpentium 25: 206-215. Kopfbeschilderung von Kreuzottern aus dem Süden der DDR. Ursel, F. (1978): Der Pileus der Squamata. Stuttgarter Beiträge zur Abhandlungen und Berichte des Naturkundemuseums Görlitz Naturkunde, Serie A 62: 29-38. Üveges, B., Halpern, B., Péchy, T., Posta, J., Komlósi, I. (2012): Biella, H.J. (1990): Zur Variation eidonomischer Merkmale bei Characteristics and heritability analysis of head scales of the Kreuzottern (Vipera b. berus) aus dem Gesamtgbiet der DDR. Hungarian meadow viper (Vipera ursinii rakosiensis Méhely, Zoologische Abhandlungen Staatliches Museum für Tierkunde 1893). Amphibia-Reptilia 33: 393-400. Dresden 45: 127-135. Werner, F. (1898): Phylogenetische Studien über die Homologien Broadley, D.G. (1996): A review of the tribe Atherini (Serpentes: und Veränderungen der Kopfschilder bei den Schlangen. Viperidae), with the descriptions of two new genera. African Arbeiten aus den Zoologischen Instituten der Universität Wien Journal of Herpetology 45(2): 40-48. und der Zoologischen Station in Triest 11: 117-162. Cherlin, V.А. (1983): Dependence of scutellation of snakes of genus Zacharias, H.C.E. (1898): Die Phylogenie der Kopfschilder bei den Echis from the climatic factors. Zoologicheskii Zhurnal 62: 252– Boiden. Zoologische Jahrbücher 10: 57-90. 257. [in Russian] Franzen, M., Heckes, U. (2000): Vipera barani Böhme & Joger, 1983 aus dem östlichen Pontus-Gebirge, Türkei: Differentialmerkmale, Verbreitung, Habitate. Spixiana 23: 61-70.