First establishment of microsatellite markers in clausiliid snails (Mollusca: Gastropoda: Clausiliidae)

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.


Findings
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 [1]. Their eponymous closing apparatus shows several morphological structures (clausiliumplate, 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 [4] and even population genetic questions [5], 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 [6], 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).
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 SSRenriched 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 [7]. 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 [8] 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 (H o ) and the expected (H e ) heterozygosity as well as linkage disequilibria and the deviation from Hardy-Weinberg-equilibrium (HWE) was calculated per population in Genepop 4.2 [9] ( 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 [10]. Thus, the potential presence of null alleles was tested using MICRO-CHECKER 2.2.3 [11]. 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.