Isolation and characterization of fourteen polymorphic microsatellite markers in the viperine snake Natrix maura

Abstract Nineteen polymorphic microsatellite loci were identified and developed for Natrix maura. Polymorphism was assessed for 120 individuals sampled across four sampling sites from the French Pyrenees Mountains. The number of alleles per locus ranged from 3 to 15, and expected heterozygosity per locus ranged from 0.227 to 0.863. We tested for deviation from Hardy–Weinberg equilibrium and linkage disequilibrium and assessed the presence of null alleles for all loci, resulting in a selection of 14 high‐quality polymorphic markers. These markers will be extremely useful in identifying fine‐scale genetic structures and providing insight into conservation management plans of this species.


| INTRODUC TI ON
Anthropogenic activities have already led to massive species extinction, and this loss of biodiversity is expected to continue at an unprecedented pace (Ceballos, Ehrlich, & Dirzo, 2017). Global warming is likely the most preoccupying threat given the potential synergy with many other environmental changes (Cahill et al., 2013;Thomas et al., 2004), impacting organisms at both the individual and population levels and resulting in local increase in extinction risks, species redistribution and community reshuffling (Aubret & Shine, 2010;Pauls, Nowak, Bálint, & Pfenninger, 2013;Walther et al., 2002).
Ectotherms represent more than 98% of animal species and are the more likely to be affected because of direct physiological sensitivity to climate conditions (Deutsch et al., 2008;Dupoué et al., 2017;Sinervo et al., 2010). When the conservation status of a given population is uncertain, genetic studies constitute an indirect and valuable approach to assess the impacts of these environmental threats on levels of population genetic diversity and structure, effective dispersal, demographic status and possible past and future responses to global change (Segelbacher et al., 2010).
The viperine snake (Natrix maura) is a common Mediterranean snake inhabiting natural and artificial aquatic environments in Northwestern Africa, Iberian Peninsula, Southern France and Northern Italy. Although some localities may exhibit high snake densities, populations are generally considered as declining (Santos & Llorente, 2009). The viperine snake is threatened by multiple factors, such as aquatic pollution, habitat loss and fragmentation, direct destruction by humans because of confusion with venomous vipers and climate change (Santos & Fernández Cardenete, 2015;Santos & Llorente, 2009). All of these environmental threats are likely to interact significantly, impacting viperine snake populations (Gangloff, Sorlin, Cordero, Souchet, & Aubret, 2019;Muthoni, 2010). In this context, the development of polymorphic genetic markers is critical for this species in order to study patterns of genetic diversity and understand population structure and functioning. Here, we isolated and characterized 19 new polymorphic microsatellite markers for N. maura using Illumina high-throughput sequencing.

| MATERIAL & ME THODS
We sampled DNA from 120 viperine snakes from four populations in the southwestern France (Ariège , Table 1), using buccal swabs as a noninvasive sampling method (Beebee, 2008 A total of 108 primer pairs, each targeting a different locus, were identified and organized into 31 multiplexes using Multiplex Manager (Holleley & Geerts, 2009). The computer-designed multiplexes were validated and checked for polymorphism using DNA samples from an additional set of seven individuals. The polymerase chain reactions (PCRs) were carried out following Schuelke (2000).
As oligonucleotide tails, we used the universal sequences M13 (GGA AAC AGC TAT GAC CAT), CAG (CAG TCG GGC GTC ATC), and T3 (AATTAA CCC TCA CTA AAGGG) labeled with the HEX dye, the FAM dye, and the TAMRA dye, respectively. PCRs were performed in a final reaction volume of 12.5 µl, containing around 10 ng of DNA, Type-it Multiplex PCR Master Mix (Qiagen), and Primer Mix 1× (0.2 µM forward primers and labeled tails, and 0.02 µM reverse primers). The optimal PCR protocol consisted in an initial denaturation step at 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 57°C for 90 s, 72°C for 30 s; 8 cycles of 95°C for 30 s, 53°C for 90 s, 72°C for 30 s; and a final extension step at 68°C for 30 min.
All PCR rounds included a negative control to check for potential cross-contamination. PCR products were subsequently subjected to fragment analysis. Allele calling was performed using Geneious 11.1.2 (Biomatters). Finally, the 19 primer pairs with the highest polymorphism were organized into seven multiplexes according to dye colors and expected amplicon sizes (Table 2). We finally applied this genotyping protocol to the 120 viperine snake samples to assess markers' quality.
To avoid any bias in further analyses, we first identified populations showing HWE across the maximum number of markers. We used the test_HW function from the R-distribution of the Genepop software (Rousset, 2008) to assess HWE for each locus and each population, and only retained loci showing HWE.
Considering each populations independently, we then estimated for each locus the number of alleles (n a ), observed and expected heterozygosity (H O and H E, respectively) using FSTAT 2.9.3.2 (Goudet, 1995). We also computed null allele frequency along with 95% confidence intervals using the null.all function from the R-package PopGenReport (Adamack & Gruber, 2014). Loci showing a lower bound exceeding a null allele frequency of 5% where discarded. We finally assessed linkage disequilibrium across populations using the test_LD function (Rousset, 2008). Tests for HWE and linkage disequilibrium were all conducted using false discovery rate FDR-correction to account for multiple-related tests (Benjamini & Hochberg, 1995).

| RE SULTS AND D ISCUSS I ON
All loci amplified well (from 0% to 5.8% of missing values). Five loci (NM_013, NM_170, NM_346, NM_384, and NM_462) did not conform to HWE. Using a subset of the 14 remaining markers, we found that all populations conformed to HWE.
In each considered population, the 19 loci were found to be polymorphic with n a ranging from 3 to 15 and H E ranging from 0.227 to 0.863 (Table 2). Yet, the presence of null alleles was detected in five loci (NM_013, NM_170, NM_245, NM_346, and NM_384). These markers should therefore be used with caution as they can significantly affect the results of genetic analyses (

CO N FLI C T O F I NTE R E S T
None declared.

AUTH O R CO NTR I B UTI O N S
All authors contribute significantly to the present study and to the

DATA AVA I L A B I L I T Y S TAT E M E N T
The microsatellite data are available on Dryad: https ://doi. org/10.5061/dryad.0vd1fj3