Zoosyst. Evol. 99 (2) 2023, 535-543 | DOI 10.3897/zse.99.110610 pg MuseuM TOR BERLIN The Caucasus is neither a cradle nor a museum of diversity of the land snail genus Helix (Gastropoda, Stylommatophora, Helicidae), while Crimea is home to an ancient lineage Ondéej Korabek!*, Igor Balashov°, Marco T. Neiber*, Frank Walther*, Bernhard Hausdorf* Leibniz Institute for the Analysis of Biodiversity Change, Zoological Museum, Martin-Luther-King-Platz 3, 20146 Hamburg, Germany Department of Zoology, Faculty of Science, Charles University, Vinicéna 7, 12844 Praha 2, Czech Republic Schmalhausen Institute of Zoology, National Academy of Sciences of Ukraine, B. Khmelnitsky str. 15, Kyiv 01030, Ukraine Universitat Hamburg, Martin-Luther-King-Platz 3, 20146 Hamburg, Germany FB WM fF https://zoobank. org/7F4E305A-A732-470C-AEB5-S8FSES8605BF Corresponding author: Ondiej Korabek (ondrej.korabekK@gmail.com) Academic editor: Frank Kohler # Received 6 August 2023 # Accepted 10 October 2023 @ Published 15 November 2023 Abstract The Caucasus and the adjacent Pontic Mountains in north-eastern Anatolia are home to numerous endemic land snail genera and species. The diversity of the region is the result of both intra-regional speciation and the persistence of relict lineages. The same seemed to be true for the genus Helix, which has been present in the Greater Caucasus since the Miocene. In the Caucasus region, there are three Helix species. Helix buchii (Pontic Mountains and Georgia) and Helix albescens (southern Ukraine to northern Lesser Caucasus) are both separated by deep splits from the major Helix clades in the mitochondrial phylogeny. In contrast, Helix luco- rum belongs to the Anatolian radiation of Helix. At least part of its intraspecific diversification may have occurred in north-eastern Anatolia and the adjacent parts of the Caucasus. Here, we report new evidence suggesting that the Caucasus and the Pontus regions were less important as a refugium of ancient He/ix lineages or as a diversification centre than previously hypothesised. Helix luco- rum probably diversified more westwards, while H. buchii is a less ancient lineage than previously thought. Helix albescens had its long-term refugium on the Crimean Peninsula in southern Ukraine, not in the Caucasus. The Caucasus is close to the eastern limit of the distribution range of the genus and, although the fossil record shows that He/ix was present there as early as the Miocene, the current diversity of the genus there is the result of much later colonisation. Key Words Anatolia, fossil, land snail, phylogeography, refugia, Ukraine Introduction The Caucasus ecoregion (Zazanashvili et al. 2020) is one of the global biodiversity hotspots, areas of high endem- ic diversity threatened by ongoing habitat loss (Myers et al. 2000; Marchese 2015). Mountainous topography and climatic stability can explain accumulation of diver- sity (Harrison and Noss 2017; Rahbek et al. 2019) and these two factors combine in the region. There are high mountains and there are climatic refugia as well. For ex- ample, some representatives of the late Neogene woody flora (Browicz 1989; Denk et al. 2001; Nakhutsrishvili et al. 2015; Magri et al. 2017; Maharramova et al. 2018), as well as some endemic land snails (Pokryszko et al. 2011; Walther et al. 2014; Neiber and Hausdorf 2017; Neiber et al. 2017; Mumladze et al. 2023) and other taxa (Tarkhnishvili et al. 2000), survived in the Colchis re- fugium in the eastern part of the Euxine floristic prov- ince of humid temperate forests, which extends along the southern Black Sea coast to the western Greater Caucasus (Browicz 1989). Quaternary glacial forest re- fugia extended from the western Caucasus along north- ern Anatolia approximately to Ordu in northern Turkey (Tarkhnishvili et al. 2012). Copyright Korabek, O. et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 536 Phylogenetic relationships between some Pontic land snail taxa on the one hand and those from southern Ana- tolia on the other hand suggest that the Pontic region con- tains relicts of lineages that were eliminated in central Anatolia by the uplift and aridification of the Anatolian Plateau and its surrounding mountain ranges (Neiber and Hausdorf 2017; Neiber et al. 2021). The uniqueness of the land snail fauna of the Caucasus and the adjacent Pon- tus Mountains (north-eastern Turkey) is demonstrated by the many species and genera that occur only in the east- ernmost part of the Pontic Mountains, especially on the northern seaward side around Artvin (Cook 1997; Schutt 2005). Other unique lineages live in the Greater Caucasus (e.g. Walther et al. (2018)). In light of this, the presence of two species of the genus He/ix without close relatives in the Caucasus and Pontus regions seemed to add to the list of ancient lineages preserved in this long-term refugium (Korabek et al. 2015), the more that there are Micoene and Pliocene Helix fossils in the Greater Caucasus (Stek- lov 1966). The first, Helix buchii Dubois de Montpéreux, 1840, is mainly distributed in the Pontus Mountains and in Georgia. The available data suggest that the diversity of mitochondrial lineages is greatest in the western part of its range (Korabek et al. 2022). The other species, Helix albescens Rossmassler, 1839 is widely distributed in the Greater Caucasus and southern Ukraine, as far west as Odesa Oblast. Until recently, the species was thought to live as far west as Bulgaria (e.g. Irikov and Eréss (2008)), but this was based on misidentifications (Korabek et al. 2022) and the confusion with Helix philibinensis Ross- massler, 1839 (as in Schileyko 1978). The distribution of its intraspecific genetic diversity is unknown. In addition to the two possibly relictual species, there is a third Helix species in the Caucasus, Helix lucorum Linnaeus, 1758. It is distributed all over the Caucasus region, Anatolia and the southern Balkan Peninsula. The occurrence of most of the main intraspecific mitochondri- al clades in the Caucasus region or adjacent north-eastern Anatolia, with endemic clades in the south-eastern Cau- casus (Korabek et al. 2022) may indicate an origin of the Species in north-eastern Anatolia or the Caucasus. It was hypothesised that the Pontus region could be a source re- gion of its Holocene expansion to the western Caucasus as well as to the Balkans (Korabek et al. 2018). Based on our expanded sampling, we tested the hypoth- esis that the Caucasus region and the Pontic Mountains were the area where the intraspecific diversity of H. /uco- rum originated and that the other two species are old, phy- logenetically isolated lineages that persisted in the Cauca- sus (H. albescens) or the Pontic Mountains (H. buchii). Methods Partial sequences of the mitochondrial genes for cyto- chrome c oxidase subunit I (cox/), 16S rRNA (rrnL) and 12S rRNA (rrnS) were analysed. Laboratory methods were as described in Korabek et al. (2022) with some zse.pensoft.net Korabek, O. et al.: Helix in the Caucasus modifications. DNA was extracted as described by Scheel and Hausdorf (2012). The primers 16Sar, 16Sbr (Palum- bi et al. 1991), 16S_MN3R (Neiber et al. 2017), 16Scs1 and 16Scs2 (Chiba 1999) were used for amplification and sequencing of the rrnL gene, a modified pair LCO1490 and HCO2198 (Hausdorf et al. 2003) for cox/ and 12SGast_fwd2 and 12SGast_rev3 (Cadahia et al. 2014) for the rrnS. Internal primers were designed for PCR and sequencing of cox/ in a number of H. buchii samples: COL buchii_F: TATTTGGCGTGTGATGTGGC, COL_ buchit_ RR: TAATAGCTCCGGCCAAAACA. The new data were combined with previously pub- lished sequences (Korabek et al. 2022 and references therein). Only seven samples of H. albescens, one of them from outside its autochthonous distribution range (Bulgaria), have been analysed previously (Korabek et al. 2022): we added 40 more from 15 additional lo- calities within the presumed native range and two from introduced populations in Ukraine (Kyiv, Poltava). For H. buchii, we analysed 57 (27 new) individuals from 32 (13 new) localities (Fig. 1B). Finally, 91 newly-added individuals from 88 sites were analysed for H. /ucorum, mainly from south-eastern Balkans and the Caucasus re- gion. All the sequences used are listed, along with locality information, in Suppl. material 1: table S1. Alignments were done with MAFFT 7.487 (E-INS-1 algorithm; Katoh and Standley (2013)). Substitution model selection was performed with IQ-TREE 1.6.12 (Nguyen et al. 2015; Kalyaanamoorthy et al. 2017; par- titioning into 3 cox/ codon positions and the two rRNA genes). Concatenated alignments of cox], rrnL and, in H.. buchii, rrnS sequences were analysed with BEAST 1.10.4 (Suchard et al. 2018) for each species separate- ly. Helix nicaeensis Férussac, 1821, the sister species of H. lucorum, was used as an outgroup for H. /ucorum; the other two species lack close relatives. BEAST was used to produce ultrametric trees where comparably old clades may be defined. Tree prior was set to Bayesian Skyline (a flexible prior not imposing strong constraints on the distribution of branch lengths), using a lognormal relaxed clock and no calibration. However, we specified the same mean rate for all species for approximate comparabili- ty. We ran the analyses in two replicates and collected a total of 10,000 post-burn-in trees for each of the three analysed species, making sure that effective sample sizes were > 800 for all parameters, except the skyline model. PCR and sequencing with the primer pair LCO1490 and HCO219 targeting the “barcoding” fragment of cox/ produced for two H. buchii individuals an incomplete pseudogene sequence containing an in-frame stop codon and frame-shift indels. Two phylogenetically distant mi- tochondrial lineages are known from H. buchii, one of which was probably acquired from an unknown, possi- bly extinct species (Korabek et al. 2015). A phylogenet- ic analysis of the pseudogene could reveal which of the mitochondrial lineages is more likely to be the original one for H. buchii. We used a representative set of Helix and Maltzanella Hesse, 1917 cox! and rrnL sequences Zoosyst. Evol. 99 (2) 2023, 535-543 537 ES H. albescense ' é H. lugorum a - BH. buchii a ——— Md ala —— EH. albescens Ee @ posterior probability >0.98 ed | F H. buchii majority lineage Mi ES ee minority lineage ; —a Figure 1. Intraspecific mitochondrial lineages of the Caucasian Helix species and their distribution. A—C. Distribution of major clades (i.e. the branches labelled with correspondingly coloured circles in the trees). Non-native or potentially non-native occur- rences are not distinguished. A. H. albescens, B. H. buchii, C. H. lucorum; D. Classification of the discussed regions: southern Ukraine (mainland, orange), Crimea (yellow), the Caucasus region (red), eastern and western part of Anatolia (dark and light blue) and the Balkans (green). The east-west division of Anatolia is for plotting purposes only and coincides with the western distribution limit of H. buchii, which also marks the limit of the wettest area of the Pontic Mountains; E-G. Bayesian mitochondrial trees based on concatenated alignments of partial rrnL, rrnS (only in H. buchii) and cox] sequences. Supported nodes are indicated by small circles at nodes (omitted from the shallowest nodes), large circles on branches denote the clades used for plotting the maps along with the corresponding numbers of the clades as given in Korabek et al. (2022). Tips are labelled with colour indicating the area of origin as in D (samples from outside the covered area are coloured white). The species are illustrated by shells from Georgia: Patara Dmanisi (H. buchii, ZMH 100320), Armenia: Gegha Valley by Kavchut (1. /ucorum, ZMH 101308) and Russia: Kamennomostkty in Adygea (H. albescens, ZMH 86704). zse.pensoft.net 538 (from a dataset used by Kordabek et al. (2022): fig 3) to analyse the phylogenetic position of the cox/ nuclear pseudogene sequences obtained from H. buchii (merged for the analysis into a single sequence at their 404 bp overlap). The sequences used are listed in Suppl. material 1: table S2. The merged pseudogene sequence was added manually to the alignment from Korabek et al. (2022) and the data were analysed with Maximum Likelihood phylo- genetic analysis in IQ-TREE following partitioning and automated model selection as above. The support was as- sessed with 500 bootstraps. The sequence data underlying this article are available from the NCBI Nucleotide database (available at: www. ncbi.nlm.nih.gov/Genbank). Accession numbers are pro- vided in the online supplementary material along with sampling localities and voucher information. Results Helix albescens Rossmiassler, 1839 We found two major clades, each of them with only shal- low internal divergences (Fig. 1E). Only one was found in the Caucasus, but both are present in southern Ukraine. In the Crimean Mountains, the two clades are parapat- ric, One occurring in the west and the other in the east. The eastern clade has its basal-most haplotype there. The western clade is surprisingly the one found also in the Caucasus. We observed three groups of haplotypes in this clade, one (without strong support) in southern Ukraine inclusive of Crimea and the other two on the southern (Georgia, Armenia) and northern (Russia) side of the Greater Caucasus, respectively. Helix buchii Dubois de Montpéreux, 1840 Two divergent lineages have been found in H. buchii. All but one of the analysed individuals yielded haplo- types belonging to the ‘majority lineage’, which is deep- ly separated from the major He/ix radiations (Figs 1F, 2). The data show generally low intraspecific diversity in Georgia, whereas multiple lineages including the ‘mi- nority lineage’, which 1s likely the sister lineage of Helix pathetica Mousson, 1854 within the Anatolian radiation of Helix (Fig. 2), were found in north-eastern Turkey. The oldest split of the majority lineage separated a clade (47 as numbered by Korabek et al. (2022); Fig. 1F) from the westernmost sampling site from all other haplotypes. Almost all Georgian samples are closely related to each other and the adjacent eastern Turkish samples. An ex- ception to these are samples from the Goderdzi Pass in south-western Georgia (Mumladze et al. 2013), which were described as Helix goderdziana Mumladze, Tarkh- nishvili & Pokryszko, 2008, sharing a mitochondrial lineage (45) with some individuals from the west of the species range. zse.pensoft.net Korabek, O. et al.: Helix in the Caucasus Partial, but overlapping cox] pseudogene sequences, were obtained from two individuals (GenBank acces- sions: isolate SP2 O0Q148366, isolate 6916 OQ148367). These samples otherwise yielded haplotypes of the clades 44 and 45 of the majority lineage. Phylogenetic analysis placed the pseudogene sequences into the Anatolian clade of Helix (Korabek et al. 2015), along with the sole sam- ple of the minority mitochondrial lineage collected near Stmela Monastery in Turkey (Fig. 2). Helix lucorum Linnaeus, 1758 In Armenia and Azerbaijan, all but one sample (from Tatev Monastery, Armenia) yielded haplotypes from two co-occurring clades (Fig. 1G), which as yet have not been found elsewhere, while almost all samples from Georgia and Russia belong to an unrelated clade that is otherwise distributed in northern Anatolia and accounts also for the majority of the occurrences in the Balkans, as well as introduced populations in various parts of Europe. This widespread intraspecific clade (number 89 in Korabek et al. (2022), blue in the map in Fig. 1C), is associated with conical to flattened conical shells in contrast to more globular ones typical for more southerly populations in Anatolia and the Lesser Caucasus (Korabek et al. 2018). The haplotypes of the clade 89 form four well-support- ed groups. One of them was found only in a non-native sample from Spain, but one was recorded so far only in the eastern Balkans and none of the remaining two is unique to the Caucasus and Pontus. Furthermore, there are nearly identical haplotypes (1 bp difference over 1469 bp of the rrnL+cox/ alignment) in the Caucasus on the one hand and the eastern Balkans on the other hand. The sister lineage of the whole clade (red in Fig. 1C) was found exclusively in the eastern Balkans. Populations with globular shells, which share mito- chondrial clade 93 (green in Fig. 1C), are widespread in northern Anatolia, but are also found in Crimea and north-eastern Bulgaria. The geographically isolated popu- lations of this lineage in Bulgaria differ in shell shape and colouration from the majority of Bulgarian H. /ucorum. Discussion Helix albescens Rossmiassler, 1839 Both major mitochondrial clades within H. albescens occur in the Crimean Mountains and the lowlands of southern Ukraine, while only one of them occurs in the Caucasus (Fig. 1A). It is possible that H. albescens was present in the Caucasus already before the Last Glacial, as we found some differentiation of Caucasian popula- tions on both sides of the Greater Caucasus ridge. How- ever, it is unlikely that H. a/lbescens has persisted in the Caucasus for an extended period of time. Instead, the Crimean Mountains, where its two main clades appear Zoosyst. Evol. 99 (2) 2023, 535-543 539 H. pathetica i> H. pomacella D H, buchil cox1 pseudogene 1S H. lucorum, H. nicaeensis 7 H. asemnis 1A H. "escherichi" |@ j oO H. asemnis 1 H. albescens Mediterranean clade European clade H. buchii majority lineage Helix (Pelasga) H. godetiana H, philibinensis bootstrap support 0.2 @ 97-100% Figure 2. Maximum Likelihood tree based on rrnL and cox] sequences showing the relationships of H. buchii. Both mitochondri- al lineages and the cox/ pseudogene (presumed NUMT) sequence obtained from two of its individuals with the majority lineage mitochondrion are included. The samples listed in Suppl. material 1: table S2 were used; the tree is rooted with Maltzanella (not shown). Support values were obtained by standard bootstrapping. to be parapatric (in contrast to the rest of Ukraine) are the most probable long-term refugium. The Crimean Mountains show some biogeographic affinities with the western Caucasus (Cameron et al. 2013; Turbanov and Balashov 2015; Neiber and Hausdorf 2017), but were a glacial refugium in their own right (e.g. Kukushkin et al. (2020)) and host a number of endemic land snail spe- cies and genera (Hausdorf 1994; Balashov and Baidash- nikov 2013; Cameron et al. 2013; Balashov 2014, 2016; Neiber et al. 2019). The long-term stronghold of H. albescens in southern Crimea supported by our data, therefore, corresponds to known patterns. Furthermore, within Ukraine, the variability in shell shape was found to be the greatest in Crimea (Kramarenko 2016) and the diversity of banding patterns is greater in the Crimean Mountains than in the Crimean lowlands (Kramarenko and Leonov 2011). The extent of the native distribution of H. albescens on the East European Plain is unclear. It is relatively continuously distributed across the steppe ecoregion of southern Ukraine and in adjacent regions of Russia south to the Caucasus, both in natural and transformed habi- tats. In Ukraine, there are no Quaternary fossil records of this species outside Crimea (Balashov 2016: 85). Since at least 2006, H. albescens has been present in Kyiv, cen- tral Ukraine, marking the first case of an evident invasion outside its native range far from its continuous distribu- tion. Since then, it has also reached several other settle- ments in central and eastern Ukraine. Thus, H. albescens is NOW spreading northwards and it 1s possible that much of its distribution in the East European Plain is the re- sult of recent expansion, possibly largely human-assisted. Similarly, several other stylommatophoran gastropods that originated or were isolated for a long time in the Crimean Mountains have relatively recently spread to the southern part of the East European Plain. To what extent this range extension is anthropogenic remains unclear. In particular, this applies to Brephulopsis cylindrica (Men- ke, 1828) (Enidae), Helicopsis filimargo (Krynicki, 1833) (Geomitridae) and Monacha fruticola (Krynicki, 1833) (Hygromiidae) (Balashov 2016; Neiber and Hausdorf 2017; Balashov et al. 2021). All these species, along with H. albescens, form a group of inhabitants of dry open and semi-open habitats that spread to the south of the East European Plain postglacially from Crimea and anthropo- genic factors allow for continuation of this expansion. Helix buchii Dubois de Montpéreux, 1840 Korabek et al. (2015) found two divergent mitochondri- al lineages within H. buchii. One sample from Stmela (Trabzon Province of Turkey) yielded a lineage belong- ing to a group of largely Anatolian species (“Anatolian clade” sensu Korabek et al. (2015)). However, all other samples yielded a different, highly divergent lineage with no relatives and unaffiliated with any major clade within the genus. The majority lineage indicates that H. buchii could be an old species persisting in a long term refu- gium from the Neogene. One or the other lineage might be the result of an introgression from another species, but the source of the introgression remained unclear be- cause neither lineage is closely related to other known extant species. It is unlikely that one of them 1s a nucle- ar pseudogene (NUMT), because the results for the rrnL and cox/ genes were replicated here also with the rrnS gene, which is in Helicidae located in a different part of the mitogenome and coded on the other strand (Terrett et al. 1996; Groenenberg et al. 2012; Gaitan-Espitia et al. 2013; Korabek et al. 2019). zse.pensoft.net 540 Based on the distribution of individual clades from the majority lineage, Helix buchii did not originate in the Caucasus hotspot or the centre of land snail diversity in the mountains around Artvin in Turkey. The data show generally low intraspecific diversity in Georgia, whereas multiple lineages (including the minority lineage) were found in north-eastern Turkey. Our results obtained with improved sampling thus confirm earlier observations (Mumladze et al. 2013; Korabek et al. 2022). The ori- gins of H. buchii can be traced to the western part of the Eastern Pontic Mountains and the additional data confirm that its current presence in the Caucasus is the result of a recent range expansion from the west. The phylogenetic placement of the cox/ pseudogene (presumed NUMT) sequenced here in the Anatolian clade indicates that the minority mitochondrial lineage is the original mitochondrial lineage of H. buchii. This means that H. buchii is younger and less isolated than suggested by the sequence data representing the majority mitochon- drial lineage. The origin of the majority lineage remains unclear. Perhaps the inferred early branching is an arte- fact caused by its long branch. In any case, the closest relatives of H. buchii appear to be species from the north- west of Anatolia. Helix lucorum Linnaeus, 1758 With two sympatric lineages restricted to the south-east- ern Caucasus region (Azerbaijan, Armenia, adjacent Iran), the species is not a complete newcomer to the area, although it is still only a very young group in the context of the Helix phylogeny (Korabek et al. 2022). However, the phylogenetic relationships of these two lineages sug- gest that they originated in Anatolia. The sister lineages of both are found there, one of them being distributed all over the north of Anatolia. As the sister group of H. /u- corum, Helix nicaeensis, lives in north-western Anatolia, the likelihood that the species diversified near the Cauca- sus region seems low and the two endemic lineages are likely the result of an earlier, pre-Holocene colonisation. The geographic origin of the widespread clade 89 (see Fig. 1C, G), which with one exception accounts for the distribution of H. /ucorum in Russia and Georgia, is un- clear, as sampling of Anatolia remains inadequate and the distribution ranges of the H. /ucorum lineages are likely to have been affected by anthropogenic translocations. Mumladze (2015) raised reasonable doubts as to whether H. lucorum is native to Georgia and the new data add to these. The fact that the sister clade of clade 89 has, so far, been found only in Bulgaria suggests that the origin of this clade is located more westerly than the Pontus re- gion, which was previously suggested as its area of origin (Korabek et al. 2018). Korabek et al. (2018) found one haplotype group in this clade to be unique to the Pontus Mountains and Georgia (apart from an introduced popu- lation in Moscow), but it is now also known from locali- ties in the Provinces of Tokat and Sakarya, 1.e. westwards zse.pensoft.net Korabek, O. et al.: Helix in the Caucasus up to the westernmost Anatolia. Very similar haplotypes at the western and eastern extremes of the range of clade 89 (1.e. eastern Balkans versus the western Caucasus), given the distance, suggest human-mediated dispersal in at least one direction, although natural dispersal cannot be excluded at present. The current distribution of different H. /ucorum \in- eages within the species’ range is influenced by anthro- pogenic translocations (Korabek et al. 2018). The oc- currence near the Tatev Monastery in Armenia could be such a case, as the lineage present there differs from the regional background. Similarly, globular-shelled pop- ulations of clade 93 (green in Fig. 1C, G) were found in north-eastern Bulgaria in an area disconnected from the main range of the clade in Anatolia. The presence of clade 93 in Crimea (a globular form previously separat- ed as H. lucorum taurica Krynicki, 1833; see Korabek et al. (2018)) is also most easily explained by introduction (Korabek et al. 2018). The most conchologically similar individuals have so far been found in north-western Ana- tolia in the provinces of Bolu and Kastamonu. It is likely that translocations have also occurred within Anatolia. The co-occurrence of several lineages in a region must, therefore, be interpreted with caution, as they may have been brought together by humans, especially if samples come from synanthropic habitats. Conclusion It is likely that Helix buchii is not an ancient and isolat- ed lineage of He/ix in the Caucasus region, but part of a group of species that diversified in Anatolia (the ‘Anato- lian clade’ of Korabek et al. (2015)). Our new data sug- gest that the species originated in the Pontus region in northern Anatolia and colonised the Caucasus region rel- atively recently. Helix albescens has also not been present in the Caucasus region for a long time, as its crown group probably originated in Crimea. The third Caucasian spe- cies, H. Iucorum, is another member of the Anatolian clade of Helix. Populations of H. /ucorum from Armenia and Azerbaijan have mitochondrial lineages that may be endemic to the region, suggesting that this species may have been present there for a longer period of time than the other two species. However, the more northerly popu- lations in Georgia and Russia belong to a recently-spread clade, for which the new data suggest an origin elsewhere, possibly as far west as the Balkans, and which probably colonised the region only with human help. The hypothesis of the Caucasus and the adjacent Pon- tus region of Turkey as a place where ancient lineages of Helix (not belonging to the European, Mediterranean or Anatolian clades of Helix) persisted for a long period of time (Korabek et al. 2015) can thus be rejected. Within the wider region, deeper intraspecific diversity was preserved or endemic species exist in the western Pontus, in the Less- er Caucasus and in Crimea, but the Greater Caucasus and the territory of Georgia were colonised relatively recently Zoosyst. Evol. 99 (2) 2023, 535-543 (comparison with the central-European species Helix po- matia suggests as late as Late Pleistocene or Holocene; Korabek et al. (2023)). There appears to be no direct link between the Helix species that lived in the northern foot- hills of the Greater Caucasus (Chechnya, Dagestan) from the Middle Miocene and Late Pliocene to the Early Pleis- tocene (Steklov 1966; Beluzhenko 2014, 2015) and the Helix species currently living in the region. This is bad news for the potential use of the remarkable fossil record of this region for evolutionary analyses of Helix. Wheth- er there is a link between Miocene Helix fossils from Crimea (Mayer 1856; Andrusov 1896; Wenz 1923) and H. albescens remains to be established, but it would be consistent with a previous estimate of the divergence time of H. albescens and the European clade of Helix at around 10 million years (Neiber and Hausdorf 2015). Acknowledgements We thank colleagues who provided us samples for anal- ysis or helped with collecting, especially S. Kramarenko (Mykolayiv National Agrarian University, Ukraine), M. Son (Institute of Marine Biology NAS Ukraine, Odesa, Ukraine), E. Rybalchenko (Poltava, Ukraine) and P. Ro- manov (Kryvyi Rih, Ukraine). The work was funded by the Alexander von Humboldt Foundation. The publica- tion of this article was funded by the Open Access Fund of the Leibniz Association. References Andrusov N (1896) Die stidrussischen Neogenablagerungen (Eine kurze Uebersicht). 1-ter Theil. Alteres Miocan. Verhandlungen der Russisch-Kaiserlichen mineralogischen Gesellschaft zu St. Petersburg. Zweite Serie 34: 195-242. https://books.google.cz/ books? id=J2tjAAAAIAAJ&redir_esc=y Balashov JIA (2014) Taurinellushka babugana gen. nov., sp. nov. (Sty- lommatophora: Pristilomatinae) from the Crimean Mountains (Ukraine) and revision of Crimean Mediterranea (Oxychilinae). Journal of Conchology 41: 575-584. Balashov IA (2016) Fauna Ukrainy. Tom 29, Mollyuski. Vypusk 5, Stebel’ chatoglazye (Stylommatophora). Naukova Dumka, Kyiv. Balashov IA, Baidashnikov AA (2013) Nazemnye mollyuski red- kolesii mozhzhevel’nika vysokogo v Krymskykh gorakh. Zoolog- icheskiy jurnal [3oonormyeckuii xx ypHan| 92: 257-263. https://doi. org/10.7868/S00445 13413030033 Balashov IA, Neiber MT, Hausdorf B (2021) Phylogeny, species delimita- tion and population structure of the steppe-inhabiting land snail genus Helicopsis in Eastern Europe. Zoological Journal of the Linnean So- ciety 193(3): 1108-1125. https://doi.org/10.1093/zoolinnean/zlaa156 Beluzhenko EV (2014) Neogene-Eopleistocene shelly limestones in the North Caucasus and Ciscaucasia regions. Stratigraphy and Geological Correlation 22(4): 406-425. https://doi.org/10.1134/ S08695938 14040029 Beluzhenko EV (2015) Stratigraficheskoe polozhenie i lateral’noe raspredelenie mshakovo-vodoroslevykh izvestnyakov miotsena 541 severnogo Kavkaza. Byulleten’ Moskovskogo obschestva ispYtatelei prirody. Otdel geologicheskii 90: 50-61. Browicz K (1989) Chorology of the Euxinian and Hyrcanian element in the woody flora of Asia. Plant Systematics and Evolution 162(1—4): 305-314. https://doi.org/10.1007/BF00936923 Cadahia L, Harl J, Duda M, Sattmann H, Kruckenhauser L, Fehér Z, Zopp L, Haring E (2014) New data on the phylogeny of Ariantinae (Pulmonata, Helicidae) and the systematic position of Cylindrus ob- tusus based on nuclear and mitochondrial DNA marker sequences. Journal of Zoological Systematics and Evolutionary Research 52(2): 163-169. https://doi.org/10.1111/jzs.12044 Cameron RAD, Pokryszko BM, Horsak M (2013) Forest snail faunas from Crimea (Ukraine), an isolated and incomplete Pleistocene refugium. Biological Journal of the Linnean Society, Linnean Society of London 109(2): 424-433. https://doi.org/10.1111/ bij. 12040 Chiba S (1999) Accelerated evolution of land snails Mandarina in the oceanic Bonin Islands: Evidence from mitochondrial DNA sequenc- es. Evolution; International Journal of Organic Evolution 53(2): 460-471. https://doi.org/10.2307/2640782 Cook LM (1997) Geographic and ecological patterns in Turkish land snails. Journal of Biogeography 24(4): 409-418. https://doi. org/10.1111/j.1365-2699.1997.00139.x Denk T, Frotzler N, Davitashvili N (2001) Vegetational patterns and distribution of relict taxa in humid temperate forests and wetlands of Georgia (Transcaucasia). Biological Journal of the Linnean Society. Linnean Society of London 72(2): 287-332. https://doi. org/10.1111/j.1095-8312.2001.tb01318.x Gaitan-Espitia JD, Nespolo RF, Opazo JC (2013) The complete mi- tochondrial genome of the land snail Cornu aspersum (Helicidae: Mollusca): intra-specific divergence of protein-coding genes and phylogenetic considerations within Euthyneura. PLoS ONE 8(6): e67299. https://doi.org/10.1371/journal.pone.0067299 Groenenberg DS, Pirovano W, Gittenberger E, Schilthuizen M (2012) The complete mitogenome of Cylindrus obtusus (Helicidae, Arian- tinae) using Illumina next generation sequencing. BMC Genomics 13(1): 114. https://doi.org/10.1186/1471-2164-13-114 Harrison S, Noss R (2017) Endemism hotspots are linked to stable climatic refugia. Annals of Botany 119(2): 207-214. https://doi. org/10.1093/aob/mcw248 Hausdorf B (1994) Additive typogenesis in Thoanteus (Gastropoda: Bulminidae). Zoological Journal of the Linnean Society 112(3): 353-361. https://doi.org/10.1111/).1096-3642.1994 tb00325.x Hausdorf B, Ropstorf P, Riedel F (2003) Relationships and origin of endemic Lake Baikal gastropods (Caenogastropoda: Rissooidea) based on mitochondrial DNA sequences. Molecular Phylogenet- ics and Evolution 26(3): 435-443. https://doi.org/10.1016/S1055- 7903(02)00365-2 Irikov A, Eréss Z (2008) An updated and annotated checklist of Bulgarian terrestrial gastropods (Mollusca: Gastropoda). Fo- lia Malacologica 16(4): 197-205. https://doi.org/10.12657/fol- mal.016.015 Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nature Methods 14(6): 587-589. https://doi.org/10.1038/ nmeth.4285 Katoh K, Standley DM (2013) MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. zse.pensoft.net 542 Molecular Biology and Evolution 30(4): 772-780. https://doi. org/10.1093/molbev/mst010 Korabek O, Petrusek A, Neubert E, Jurfi¢kova L (2015) Molecular phy- logeny of the genus Helix (Pulmonata: Helicidae). Zoologica Scripta 44(3): 263-280. https://doi.org/10.1111/zsc.12101 Korabek O, Juri¢kova L, Balashov I, Petrusek A (2018) The contribu- tion of ancient and modern anthropogenic introductions to the colo- nization of Europe by the land snail He/ix /ucorum Linnaeus, 1758 (Helicidae). Contributions to Zoology (Amsterdam, Netherlands) 87(2): 61-74. https://doi.org/10.1163/18759866-08702001 Korabek O, Petrusek A, Rovatsos M (2019) Mitogenome of Helix pomatia and the basal phylogeny of Helicinae (Gastropoda, Sty- lommatophora, Helicidae). ZooKeys 827: 19-30. https://doi. org/10.3897/zookeys.827.33057 Korabek O, Juri¢kova L, Petrusek A (2022) Diversity of land snail tribe Helicini (Gastropoda: Stylommatophora: Helicidae): where do we stand after 20 years of sequencing mitochondrial markers? Diversity (Basel) 14(1): 24. https://doi.org/10.3390/d14010024 Korabek O, Adamcova T, Prockow M, Petrusek A, Hausdorf B, Juriékova L (2023) In both directions: Expansions of European land snails to the north and south from glacial refugia. Journal of Bioge- ography 40(4): 654-668. https://doi.org/10.1111/jbi.14531 Kramarenko SS (2016) Patterns of spatio-temporal variation in land snails: A multi-scale approach. Folia Malacologica 24(3): 111-177. https://doi.org/10.12657/folmal.024.008 Kramarenko SS, Leonov SV (2011) Phenetic population structure of the land snail Helix albescens (Gastropoda, Pulmonata, Helicidae) in the Crimea. Russian Journal of Ecology 42(2): 170-177. https://doi. org/10.1134/S1067413611020068 Kukushkin OV, Ermakov OA, Ivanov AY, Doronin IV, Sviridenko EY Sviridenko, Simonov EP, Gorelov RA, Khramova MA, Blokhin IG (2020) Cytochrome b mitochondrial gene analysis-based phyloge- ography of a Sand lizard in the Crimea: ancient refugium at the pen- insula, late expansion from the North, and first evidence of Lacerta agilis tauridica and L. a. exigua (Lacertidae: Sauria) hybridization. Proceedings of the Zoological Institute RAS 324: 56-99. https://doi. org/10.31610/trudyzin/2020.324.1.56 Magri D, di Rita F, Aranbarri J, Fletcher W, Gonzalez-Sampériz P (2017) Quaternary disappearance of tree taxa from Southern Eu- rope: Timing and trends. Quaternary Science Reviews 163: 23-55. https://doi.org/10.1016/j.quascirev.2017.02.014 Maharramova E, Huseynova I, Kolabaia S, Gruenstaeud! M, Borsch T, Muller LAH (2018) Phylogeography and population genetics of the riparian relict tree Pterocarya fraxinifolia (Juglandaceae) in the South Caucasus. Systematics and Biodiversity 16(1): 14-27. https:// doi.org/10.1080/14772000.2017.1333540 Marchese C (2015) Biodiversity hotspots: A shortcut for a more com- plicated concept. Global Ecology and Conservation 3: 297-309. https://doi.org/10.1016/j.gecco.2014.12.008 Mayer MC (1856) Description de coquilles fossiles des terrains tertiaires du midi de la Russie. Journal de Conchyliologie 5: 96-113. [pl. 4] https://www. biodiversitylibrary.org/item/5387 1#page/102/mode/lup Mumladze L (2015) Species of the genus Helix (Mollusca, Gastropoda) in Georgia. Proceedings of the Institute of Zoology (Ilia State Uni- versity, Tbilisi) 24: 148-157. http://eprints. iliauni.edu. ge/4295/ Mumladze L, Tarkhnishvili D, Pokryszko BM (2008) A new species of the genus Helix from the Lesser Caucasus (SW Georgia). Journal of Conchology 39: 483-485. zse.pensoft.net Korabek, O. et al.: Helix in the Caucasus Mumladze L, Tarkhnishvili D, Murtskhvaladze M (2013) Sys- tematics and evolutionary history of large endemic snails from the Caucasus (Helix buchii and H. goderdziana) (Helicidae). American Malacological Bulletin 31(2): 225-234. https://doi. org/10.4003/006.03 1.0202 Mumladze L, Grego J, Szekeres M (2023) The land snail family Clausi- liidae (Gastropoda, Pulmonata, Stylommatophora) in Georgia: Overview, novel records and a new species. Caucasiana 2: 29-61. https://doi.org/10.3897/caucasiana.2.e101013 Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J (2000) Biodiversity hotspots for conservation priorities. Nature 403(6772): 853-858. https://doi.org/10.1038/35002501 Nakhutsrishvili G, Zazanashvili N, Batsatsashvili K, Montalvo Man- cheno CS (2015) Colchic and Hyrcanian forests of the Caucasus: Similarities, differences and conservation status. Flora Mediterranea 25: 185-192. https://doi.org/10.7320/FlMedit25SI.185 Neiber MT, Hausdorf B (2015) Molecular phylogeny reveals the poly- phyly of the snail genus Cepaea (Gastropoda: Helicidae). Molecular Phylogenetics and Evolution 93: 143-149. https://doi.org/10.1016/). ympev.2015.07.022 Neiber MT, Hausdorf B (2017) Molecular phylogeny and biogeography of the land snail genus Monacha (Gastropoda, Hygromiidae). Zoo- logica Scripta 46(3): 308-321. https://doi.org/10.1111/zsc.12218 Neiber MT, Razkin O, Hausdorf B (2017) Molecular phylogeny and biogeography of the land snail family Hygromiidae (Gastropoda: Helicoidea). Molecular Phylogenetics and Evolution 111: 169-184. https://doi.org/10.1016/j.ympev.2017.04.002 Neiber MT, Helfenrath K, Walther F, Hausdorf B (2019) Ecological spe- cialization resulting in restricted gene flow promotes differentiation in door snails. Molecular Phylogenetics and Evolution 141: 106608. https://doi.org/10.1016/j.ympev.2019.106608 Neiber MT, Walther F, Kijashko PV, Mumladze L, Hausdorf B (2021) The role of Anatolia in the origin of the Caucasus biodiversity hotspot illustrated by land snails in the genus Oxychilus. Cladistics 38(1): 83-102. https://doi.org/10.1111/cla.12479 Neuyen L-T, Schmidt HA, von Haeseler A, Minh BQ (2015) IQ-TREE: A fast and effective stochastic algorithm for estimating maxi- mum-likelihood phylogenies. Molecular Biology and Evolution 32(1): 268-274. https://doi.org/10.1093/molbev/msu300 Palumbi S, Martin A, Romanco S, McMillan WO, Stice L, Grabowski G (1991) The simple fool’s guide to PCR. Department of Zoology, University of Hawaii, Honolulu. Pokryszko BM, Cameron RAD, Mumladze L, Tarkhnishvili D (2011) Forest snail faunas from Georgian Transcaucasia: Patterns of diver- sity in a Pleistocene refugium. Biological Journal of the Linnean Society, Linnean Society of London 102(2): 239-250. https://doi. org/10.1111/j.1095-8312.2010.01575.x Rahbek C, Borregaard MK, Colwell RK, Dalsgaard B, Holt BG, Morueta-Holme N, Nogues-Bravo D, Whittaker RJ, Fjeldsa J (2019) Humboldt’s enigma: What causes global patterns of mountain bio- diversity? Science 365(6458): 1108-1113. https://doi.org/10.1126/ science.aax0149 Scheel BM, Hausdorf B (2012) Survival and differentiation of subspe- cies of the land snail Charpentieria itala in mountain refuges in the Southern Alps. Molecular Ecology 21(15): 3794-3808. https://doi. org/10.1111/j.1365-294X.2012.05649.x Schileyko AA (1978) Nazemnye mollyuski nadsemeistva Helicoidea. In: Skarlato OA (Ed.) Fauna SSSR, Tom 3 vyp. 6. Nauka, Leningrad. Zoosyst. Evol. 99 (2) 2023, 535-543 Schutt H (2005) Turkish land snails. 4", revised and enlarged edition. Verlag Natur & Wissenschaft, Solingen. Steklov AA (1966) Nazemnye mollyuski neogena Predkavkaz’ya 1 ikh stratigraficheskoe znachenie. Nauka, Moskva. Suchard MA, Lemey P, Baele G, Ayres DL, Drummond AJ, Rambaut A (2018) Bayesian phylogenetic and phylodynamic data integra- tion using BEAST 1.10. Virus Evolution 4(1): vey016. https://doi. org/10.1093/ve/vey016 Tarkhnishvili DN, Thorpe RS, Arntzen JW (2000) Pre-Pleistocene re- fugia and differentiation between populations of the Caucasian sal- amander (Mertensiella caucasica). Molecular Phylogenetics and Evolution 14(3): 414422. https://doi.org/10.1006/mpev. 1999.0718 Tarkhnishvili DN, Gavashelishvili A, Mumladze L (2012) Palaeocli- matic models help to understand current distribution of Caucasian forest species. Biological Journal of the Linnean Society, Linnean Society of London 105(1): 231-248. https://doi.org/10.1111/j.1095- 8312.2011.01788.x Terrett JA, Miles S, Thomas RH (1996) Complete DNA sequence of the mitochondrial genome of Cepaea nemoralis (Gastropoda: Pulmo- nata). Journal of Molecular Evolution 42(2): 160-168. https://doi. org/10.1007/BF02198842 Turbanov I, Balashov IA (2015) A second record of Selenochlamys (Stylommatophora: Trigonochlamydidae) from Crimea. Malacolog- ica Bohemoslovaca 14: 1-4. https://do1.org/10.5817/MaB2015-14-1 Walther F, Kiyashko P, Harutyunova L, Mumladze L, Neiber MT, Hausdorf B (2014) Biogeography of the land snails of the Caucasus region. Tenta- cle 22: 3—5. http://www.hawaii.edu/cowielab/Tentacle/Tentacle_22.pdf Walther F, Neiber MT, Hausdorf B (2018) Systematic revision of the Caucasigenini (Gastropoda: Hygromiidae) from the Caucasus re- gion. Archiv fur Molluskenkunde 147(1): 129-169. https://doi. org/10.1127/arch.moll/147/129-169 543 Wenz W (1923) Fossilium Catalogus I: Animalia. Pars 18: Gastropoda extramarina tertiaria. II. W. Junk, Berlin. https://www.biodiversityli- brary.org/item/125702#page/1/mode/1up Zazanashvili N, Manvelyan K, Askerov E, Mousavi M, Krever V, Kalem S, Garforth M (2020) The boundaries and bio-physical fea- tures of the Caucasus Ecoregion. In: Zazanashvili N, Garforth M, Bitsadze M (Eds) Ecoregional Conservation Plan for the Cauca- sus, 2020 Edition: Supplementary Reports. WWF, KfW, Tbilisi, 9-20. _ https://wwfint.awsassets.panda.org/downloads/ecp_2020_ part_2.pdf Supplementary material | List of sequences Authors: Ondfej Korabek, Igor Balashov, Marco T. Neiber, Frank Walther, Bernhard Hausdorf Data type: xlsx Explanation note: table S1. List of sequences used for intraspecific phylogenetic analyses presented in Fig. 1. table S2. List of sequences used for inferring the phy- logenetic position of a cox/ pseudogene in Fig. 2. Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons. org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow us- ers to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://do1.org/10.3897/zse.99.110610.suppl1 zse.pensoft.net