- Short Report
- Open Access
First establishment of microsatellite markers in clausiliid snails (Mollusca: Gastropoda: Clausiliidae)
BMC Research Notesvolume 10, Article number: 137 (2017)
Clausiliidae (door snails) are gastropods with a very high diversity concerning shell morphology, especially of their complex closing apparatus, which provides the most important diagnostic traits for classification of taxa. Due to the high variability, a high number of taxa has been described, though their systematics and taxonomy is partially controversially discussed. Montenegrina is the second most speciose door snail genus in Europe. It is an obligate rock-dwelling land snail and has, compared to its complex systematics, a rather small distribution range in the western parts of the Balkan Peninsula. The different taxa themselves show a very narrow and patchy distribution range. As Montenegrina is comprehensively sampled over the whole distribution range, it is a perfect study system for general questions on speciation and morphological differentiation in land snails. To study the amount of gene flow between geographically close or co-occurring populations, highly polymorphic markers are needed.
Thirteen microsatellite loci with a tetranucleotid repeat were isolated and tested in three geographically close Montenegrina populations (two populations of M. dofleini prespaensis from the Prespa Lake, n = 35 and one population from M. stankovici from the Ohrid Lake, n = 20). The number of alleles per locus ranged from 2 to 27. No significant linkage disequilibria between the same two loci were found in all three tested populations. The deviation from Hardy–Weinberg equilibrium reveal only for two loci a significant deviation from HWE in more than one population (Mont_5483 and Mont_4477).
The 13 newly established genetic markers will help to gain better insights to the population genetic structure of Montenegrina and might reveal new results about speciation processes in co-occurring taxa. Furthermore, these microsatellite loci could also be tested in other clausiliid species.
The gastropod family Clausiliidae (door snails) is one of the most speciose among land snails and is distributed in four continents with diversity hot-spots in southern Europe, southeastern Asia and the northwestern part of South America . Their eponymous closing apparatus shows several morphological structures (clausilium-plate, lunella-complex, lamellae and plicae), which are very important for discriminating taxa within this family. Due to the high morphological variability and the high number of described taxa, systematics and taxonomy is quite complicated and controversial. Although several clausiliid taxa have been intensively studied concerning phylogenetic [2, 3], phylogeographic  and even population genetic questions , no microsatellite markers have been applied so far in any of the door snail genera.
Montenegrina is an obligate rock-dwelling door-snail genus comprising 106 taxa , which are morphologically very diverse. It is the second speciose door snail genus in Europe and has, compared to its complex systematics, a rather small distribution range in the western parts of the Balkan Peninsula. The different taxa show very patchy distribution ranges due to low dispersal ability and the insular occurrence of their calcareous rock habitats. Due to this narrow distribution range and well-known occurrences of many populations, a quite comprehensive sampling is possible, which makes Montenegrina a perfect study system for general questions on speciation and morphological diversification of rock-dwelling gastropods. Due to the high morphological variability species delimitation is sometimes challenging in Montenegrina, especially in the case of co-occurring morphotypes. Therefore, estimating the amount of gene flow is crucial.
To test for possible gene flow between populations or morphotypes and to get insights into population structure in general, highly polymorphic nuclear genetic markers are needed, which can be applied to a high number of individuals at moderate cost. Therefore, we established 13 microsatellite loci for Montenegrina on the basis of next generation sequencing data and tested them in two species, which occur very close together (about 25 km apart). (1) Montenegrina dofleini prespaensis Nordsieck, 1988 occurs along the western and southern shores of the Prespa Lake (Albania, Greece). For this study we selected specimen from two populations that live in close vicinity: MprP from Psarades (15 specimens) and MprPE from Panagia Eleousa Cave (20 specimens). (2) Montenegrina stankovici (Urbański, 1960) lives along the eastern shore of Lake Ohrid (Macedonia). The studied population MstSV was collected in Sveti Naum (20 specimens).
For DNA extraction we used the Qiagen spin column kit (Qiagen DNeasy Blood & Tissue Kit, Germany) and followed the manufacturer’s protocol. Microsatellite markers were established by ecogenics GmbH (Schlieren, Switzerland) based on seven individuals. Size-selected fragments from genomic DNA were enriched for SSR content by using magnetic streptavidin beads and biotin-labeled GTAT, GATA, AAAC and AAAG repeat oligonucleotides. The SSR-enriched library was analysed on an Illumina MiSeq platform using the Nano 2 × 250 v2 format. After assembly, 13,656 contigs or singlets contained a microsatellite insert with a tetra- or a trinucleotide of at least six repeat units or a dinucleotide of at least ten repeat units. Suitable primer design was possible in 2905 microsatellite candidates, of which 36 were positively tested for functionality and polymorphism. This was done with fluorescent labelled M13 tails that were added to the forward primers . The primers for these 36 microsatellite loci were tested for consistency in PCR amplification and finally the 13 best loci were chosen for further testing in a larger set of individuals. Using the software Multiplex Manager 1.2  these primer sets were then combined in three multiplex PCR reactions (Set1/2/3). For conducting multiplex analyses the forward primers were marked with fluorescent dyes. PCR was performed in 10 µl reaction volume using the Qiagen Multiplex-Kit (1 µl DNA, 2.5 µl Qiagen Multiplex PCR Mastermix, 1 µl primer mix, [1 µl Q-Solution in Set2], 5.5 µl AD [4.5 µl AD in Set2]). The two-step PCR profile started with an initial denaturation at 95 °C for 15 min, two cycles denaturation at 94 °C for 30 s, annealing at 61 °C (Set1) and 58 °C (Set2 and 3) for 90 s, extension at 72 °C for 60 s. For the following 35 cycles the same profile was used, except an altered annealing temperature of 58 °C (Set1) and 51 °C (Set2 and 3), respectively, and a final extension at 60 °C for 30 min. Positive PCR products were diluted 1:10 with AD. 1 µl of the PCR dilution was added to 9 µl of a Hi-Di formamide-Size Standard mixture (0.25 µl, Gene Scan 500 LIZ, Applied Biosystems, USA) and denatured at 94 °C for 4 min. The products were then analysed on a 3130xl Sequence Analyzer (Applied Biosystems, USA). The alleles obtained from the electropherograms were identified and binned using the software GENEMAPPER 5.0 (Applied Biosystems, USA) and checked manually. The multiplex PCR conditions for the 13 loci, the repeat motif and the primer sequences as well as the respective GenBank Accession numbers are given in Table 1.
In total 55 individuals—35 of M. d. prespaensis and 20 of M. stankovici—were genotyped. For each locus the observed (Ho) and the expected (He) heterozygosity as well as linkage disequilibria and the deviation from Hardy–Weinberg-equilibrium (HWE) was calculated per population in Genepop 4.2  (Table 2). The number of alleles per locus ranged from 2 to 27. No significant linkage disequilibria between the same two loci were found in all three tested populations, two loci pairs show linkage disequilibrium but only in one population (Mont_2349/Mont_4477 in MprPE; Mont_17419/Mont_5483 in MprP). The deviation from the Hardy–Weinberg equilibrium is significant for four loci in MprP, four in MprPE and three in MstSV. However, checking the results per locus in all populations, reveal only in two cases a significant deviation from HWE over two (Mont_5483) or even all three populations (Mont_4477). The high deficiency of heterozygotes over several loci has already been reported in other gastropod species and could be due to self-fertilization or the patchy distribution leading to drift effects due to restricted gene flow, single loci could be affected by allele scoring artefacts or the presence of null alleles . Thus, the potential presence of null alleles was tested using MICRO-CHECKER 2.2.3 . The results revealed three potential null alleles in the MprP population five in the MprPE population and two in the MstSV population. The loci Mont_4477 and Mont_5483 have potential null alleles in the MprP and in the MprPE population. However, in MstSV these loci give no hint for null alleles and were therefore not excluded from the analysis. At the locus Mont_5717 only one allele per population was found, yet in preliminary tests with other Montenegrina species, further alleles are found at this locus. Consequently, the 13 microsatellite markers are suitable for population genetic studies in the rock-dwelling land snail genus Montenegrina.
The newly established genetic markers will be useful in many further investigations within Montenegrina. We expect to test gene flow between syntopically occurring congeners, as well as among parapatric and allopatric populations of conspecific taxa. We hope that this will promote the better understanding of speciation patterns and speciation processes in obligate rock-dwelling gastropods. Furthermore, they could also be tested in other clausiliid species.
studied population of Montenegrina dofleini prespaensis Nordsieck, 1988, collected in Psarades
studied population of Montenegrina dofleini prespaensis Nordsieck, 1988, collected in the Panagia Eleousa Cave
studied population of Montenegrina stankovici (Urbański, 1960), collected in Sveti Naum
- Ho :
- He :
Nordsieck H. Worldwide door snails. Hackenheim: ConchBooks; 2007.
Hausdorf B, Naegele K-L. Systematics of Strobeliella from the Southern Alps and its relationships within Clausilia (Gastropoda: Clausiliidae). J Molluscan Stud. 2016;82:31–6.
Uit de Weerd DR, Gittenberger E. Phylogeny of the land snail family Clausiliidae (Gastropoda: Pulmonata). Mol Phylogenet Evol. 2013;67:201–16.
Schilthuizen M, Gutteling E, Van Moorsel CHM, Welter-Schultes FW, Haase M, Gittenberger E. Phylogeography of the land snail Albinaria hippolyti (Pulmonata: Clausiliidae) from Crete, inferred from ITS-1 sequences. Biol J Linn Soc. 2004;83:317–26.
Scheel BM, Hausdorf B. Survival and differentiation of subspecies of the land snail Charpentieria itala in mountain refuges in the Southern Alps. Mol Ecol. 2012;21:3794–808.
Fehér Z, Szekeres M. Taxonomic revision of the rock-dwelling door snail genus Montenegrina Boettger, 1877. Zookeys. 2016;599:1–137.
Schuelke M. An economic method for the fluorescent labeling of PCR fragments—A poor man’s approach to genotyping for research and high-throughput diagnostics. Nat Biotechnol. 2000;18:233–4.
Holleley CE, Geerts PG. Multiplex Manager 1.0: a cross-platform computer program that plans and optimizes multiplex PCR. Biotechniques. 2009;46:511–7.
Rousset F. GENEPOP’ 007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.
Armbruster GFJ, Alvarez C, Pesaro M, Baur B. Polymorphic microsatellite DNA markers in the endangered land snail. Trichia caelata (Gastropoda. Stylommatophora). Mol Ecol Notes. 2007;7:1123–4.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.
KJ, ZF and LK conceived the study. ZF provided the samples. KJ and LK analysed the data. KJ carried out the laboratory work and drafted the manuscript. EH contributed with technical suggestions and improving the draft. All authors read and approved the final manuscript.
We are grateful to Helmut Sattmann for providing material and for fruitful discussions. We thank Sandra Kirchner for critical comments on the manuscript and Sonja Bamberger for help with the laboratory work.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
This project was funded by the Austrian Science Foundation (FWF. Project-number P26581-B25).