Primula halleri J. F. Gmel. (Primulaceae) is a herbaceous perennial species with a distributional range limited to the Alps, Carpathians, and the Rila and Pirin mountains of the Balkan region. In these Central European mountains, P. halleri populations are relatively rare and restricted to high alpine habitats between 1500 and 2700 m (Richards, 2003; de Vos et al., 2012). Due to insufficient information for much of the geographic range, P. halleri is not formally listed as vulnerable or threatened in the IUCN Red List (Walter and Gillett, 1998), although it is defined as vulnerable or regionally extinct in Switzerland (Moser et al., 2002). Together, P. halleri, P. farinosa L., P. frondosa Janka, and P. scotica Hook. form a European clade within Primula sect. Aleuritia (Guggisberg et al., 2006). Primula farinosa has a wide but patchy distribution across Europe, P. frondosa is endemic to northeastern Bulgaria, and P. scotica is endemic to northern Scotland and the Orkney Islands (Richards, 2003). Primula sect. Aleuritia exhibits considerable variation in ploidy, but the species sampled here do not show evidence of intraspecific ploidy variation (Guggisberg et al., 2006). Primula farinosa and P. frondosa are both diploid (2n = 2x = 18), P. scotica is hexaploid (2n = 6x = 54), and P. halleri is tetraploid (2n = 4x = 36). These four species also exhibit different breeding systems; P. farinosa and P. frondosa are distylous, while P. halleri and P. scotica are homostylous (de Vos et al., 2012). Populations of distylous species have individuals of two distinct floral morphs, which differ in the reciprocal positioning of sexual organs and pollen mating type. Homostylous species are composed of individuals carrying only one floral morph and pollen mating type (Richards, 2003).
Here we report the development of 12 polymorphic novel microsatellite loci for P. halleri, which will be used in future studies to characterize the geographic structure of genetic diversity and mating system of this species. Additionally, we evaluate the amplification and polymorphism of these loci in the related species of the European clade of Primula sect. Aleuritia.
METHODS AND RESULTS
Genomic DNA was isolated from one P. halleri individual from the Binnergale population in Switzerland (voucher Z-88501; 46°22′53.1″N, 08°13′40.8″E) using the DNeasy Plant Mini Kit (QIAGEN, Hombrechtikon, Switzerland) according to the manufacturer's standard protocol. This DNA sample was sent to the Cornell Evolutionary Genetics Core Facility ( https://cores.lifesciences.cornell.edu/brcinfo/) for the development of microsatellite loci. Briefly, genomic DNA was digested with two restriction enzymes, adapters were ligated, and the library was cloned into a pUC19 plasmid. The library was enriched for genomic fragments containing microsatellite loci using a standard pull-down method with streptavidin-coated magnetic beads and 3′-biotinylated oligonucleotides carrying repeat motifs (TA)17, (GA)17, (TC)17, (TTC)7, (TACA)7, and (TCAA)7. The enriched library was amplified and cloned, and 32 positive library colonies were selected for sequencing. From the sequence of these 32 colonies, 28 microsatellite loci were identified and primers were designed to amplify products ranging from 100 to 500 bp using Primer3 version 3.0 software (Rozen and Skaletsky, 2000).
Each primer pair was tested for amplification using a test panel of four P. halleri individuals from each of three populations in Switzerland (see Table 1). Genomic DNA was isolated as described above, and PCR was performed in 25-µL reactions containing: 10–30 ng of DNA, 2.5 µL of l0× reaction buffer, 2 µL of MgCl2 (50 mM), 1.25 µL of each dNTP (10 mM), 0.5 µL of both forward and reverse primer (10 µM), and 0.25 µL of Taq DNA polymerase (5 U/µL; Bioline GmbH, Luckenwalde, Germany). PCR was performed in a Bio-Rad C1000 Thermal Cycler (Bio-Rad Laboratories, Hercules, California, USA) under the following conditions: initial denaturation at 95°C for 3 min followed by 30 cycles of 95°C for 1 min, 55–57°C for 1 min (see Table 1), and 72°C for 1 min, with a final extension at 72°C for 5 min. Successful amplification and quantification of fragment length was determined on a QIAxcel Advanced electrophoresis system using the DNA High Resolution kit and ScreenGel software (QIAGEN).
TABLE 1.
Primer sequences and characteristics of microsatellite loci developed in Primula halleri.
Thirteen of the 28 primer pairs tested consistently amplified a single locus and were thus selected to evaluate polymorphism in the 12 individuals of the test panel used above. PCR was performed in 10-µL reactions containing 10–30 ng of genomic DNA, 2.0 µL of 5× GoTaq buffer, 0.5 µL of each dNTP (10 mM), 0.8 µL of MgCl2 (50 mM), 0.2 µL of the forward primer with an M13(−21) tail at the 5′ end (10 µM), 0.5 µL of the reverse primer (10 µM), 0.5 µL of fluorescently labeled M13(−21) primer (FAM, NED, PET; 10 µM), and 0.1 µL of GoTaq DNA polymerase. Amplifications were performed in a BioRad C1000 Thermal Cycler (Bio-Rad Laboratories) following the two-step method described by Schuelke (2000) with cycling conditions as described above and an additional 10 cycles of 95°C for 30 s, 53°C for 30 s, and 72°C for 1 min prior to the final extension. Amplified products from as many as three loci were multiplexed for capillary electrophoresis on an ABI 3130xl Genetic Analyzer (Applied Biosystems, Foster City, California, USA) using GeneScan 500 LIZ as internal standard. Fragment lengths were determined using GeneMapper software version 4.1 (Applied Biosystems).
TABLE 2.
Summary statistics of microsatellite polymorphism determined by screening 45 Primula halleri samples from three populations in Switzerland.a
TABLE 3.
Cross-amplification of 11 novel microsatellite loci in Primula farinosa, P. frondosa, and P. scotica.
Amplicons from 11 of the 13 microsatellite loci showed clear evidence of polymorphism, while two loci were monomorphic in the 12 individuals tested. Only one (H4) of the 11 polymorphic loci was found to carry more than two alleles per individual. In all of the samples genotyped at this locus, no more than two alleles were observed within two distinct allelic size ranges (110–126 bp and 130–142 bp). Thus, following Saltonstall (2003), we treated this locus as if it were two distinct diploid loci (H4-1 and H4-2; see Table 1). The remaining 10 loci appear to be amplifying only one of the two diploid copies of the tetraploid genome, suggesting that homeologous loci in this species may be quite divergent at the sequence level (Soltis et al., 2010).
The 12 polymorphic loci were used to estimate genetic diversity parameters and deviation from Hardy–Weinberg equilibrium in 45 individuals chosen randomly from the three Swiss populations described in Table 1. Genetic diversity statistics (e.g., observed [Ho] and expected heterozygosity [He], respectively) were calculated in GENEPOP version 4.0 (Rousset, 2008; Table 2). The total number of observed alleles per locus in all populations ranged from two to eight, per locus Ho ranged from 0 to 0.733, and He ranged from 0 to 0.867.
Cross-amplification was tested in closely related species using 16 individuals of P. farinosa (voucher no. Z-88504), eight individuals of P. frondosa (voucher no. Z-88505), and 12 individuals of P. scotica (voucher no. Z-88506). Because locus H4 does not show more than two alleles per individual in these three species (unlike the results for P. halleri), it is treated as one locus in the cross-amplification results. The results presented in Table 3 show that of the 11 loci tested, only eight consistently amplify across species boundaries.
CONCLUSIONS
We describe here the development of 12 polymorphic microsatellite loci for the high alpine endemic species P. halleri. These loci make possible future studies focused on estimating parameters of the mating system and genetic structure of P. halleri populations. Successful cross-amplification suggests that some of these loci could be applied to estimate the genetic diversity in other species of Primula sect. Aleuritia.