Characterization of 12 polymorphic SSR markers in Veronica subsect. Pentasepalae (Plantaginaceae) and cross-amplification in 10 other subgenera1

Premise of the study: Microsatellite primers were developed in the perennial herbs of the diploid-polyploid complex Veronica subsect. Pentasepalae (Plantaginaceae) to investigate the role that hybridization has played in the evolution of the group, which includes several endangered species. Methods and Results: Twelve pairs of primers leading to polymorphic and readable markers were identified and optimized from V. jacquinii and V. orbiculata using a microsatellite-enriched library method and 454 GS-FLX technique. The set of primers amplified dinucleotide to pentanucleotide repeats, and the number of alleles per locus ranged from one to six, one to 11, and one to nine for V. orsiniana, V. javalambrensis, and V. rosea, respectively. Transferability analyses were performed in 20 species corresponding to 10 different subgenera. Conclusions: These results indicate the utility of the newly developed microsatellites across Veronica subsect. Pentasepalae, which will help in the study of gene flow patterns and genetic structure.

The genus Veronica L. (Plantaginaceae) comprises ca. 450 species, which are grouped into 12 subgenera with between two and 180 species each ( Albach et al., 2004 ;Garnock-Jones et al., 2007 ). It includes some perennials of relative economic importance in ornamental horticulture and others that are well-known widespread weeds. Additionally, several species of Veronica are registered on the International Union for Conservation of Nature Red List ( http://www.iucnredlist.org/ ) and other regional catalogs of endangered plants (e.g., Peñas de Giles et al., 2004 ), or are threatened plants with narrow distribution areas (e.g., Petrova and Vladimirov, 2009 ).
Veronica subsect. Pentasepalae Benth. is a monophyletic diploid-polyploid complex and one of the four subsections currently recognized within the also monophyletic Veronica subgen. Pentasepalae M. M. Mart. Ort., Albach & M. A. Fischer ( Albach et al., 2008 ). This subsection comprises ca. 20 perennial taxa and is represented in the temperate regions of Eurasia with one species in North Africa. The complex seems to be of recent origin and divergence, as many diploid representatives are still extant and short branches are found in the phylogenetic analyses based on ITS and plastid DNA sequence data ( Rojas-Andrés et al., 2015 ). Although the diploid species are characterized by subtle morphological differences, each has been recovered as monophyletic in previous studies. Hybridization and polyploidization are widespread in the group, and several authors ( Lehmann, 1937 ;Scheerer, 1949 ;Rojas-Andrés et al., 2015 ) have concluded that gene fl ow and complex relationships among polyploids and their diploid relatives might exist. Interestingly, some of the diploid and polyploid species belonging to Veronica subsect. Pentasepalae are Mediterranean orophytes that face a high risk of extinction with climate warming and/ or grow in Important Plant Areas (IPAs; IPA online database: http://www.plantlifeipa.org/reports.asp ), regions that display exceptionally rich fl oras of biogeographic interest ( Rojas-Andrés et al., 2015 ). Given that current gene fl ow and introgression may have blurred species limits, particularly in hybrid zones, accurate investigations of gene fl ow patterns within and among Veronica subsect. Pentasepalae populations are necessary for conservation and species delimitation purposes.

METHODS AND RESULTS
Microsatellite development -For the microsatellite library, silica gel-dried leaves of 12 diploid individuals of V. jacquinii Baumg. and V. orbiculata A. Kern. were selected from eight different populations (Appendix 1). Ploidy level was checked using fl ow cytometry. A microsatellite library was prepared by Genoscreen (Lille, France) using a 454 GS-FLX (Roche Diagnostics, Meylan, France) high-throughput DNA sequencer ( Malausa et al., 2011 ). Genomic DNA was extracted using the cetyltrimethylammonium bromide method described in Doyle and Doyle (1987 a See Appendix 1 for locality and voucher information for each population . b Deviations from HWE were not statistically signifi cant (ns) and statistically signifi cant at * P < 0.05, ** P < 0.01, and *** P ≤ 0.001.
AAC, AAG, AGG, ACG, ACAT, and ACTC motifs. A total of 32,052 high-quality sequences were obtained. Analyses of these sequences with QDD software ( Meglécz et al., 2010 ) revealed 3010 sequences with microsatellite motifs, for which 195 pairs of primers were obtained. Given that it is too time consuming and not affordable to check all of the primer pairs obtained, 54 of them with low primer pair penalty and different lengths and repeat motifs were selected. These primers were ordered (Eurofi ns, Ebersberg, Germany) to evaluate polymorphic loci on 12 individuals from the complex V. jacquinii-V. orbiculata . PCRs were performed in a total volume of 15 μ L, which contained 1 × PCR Green GoTaq Buffer (Promega Corporation, Madison, Wisconsin, USA), 0.25 mM of each dNTP (Life Technologies, Carlsbad, California, USA), 0.33 mM of each primer, 0.5 units GoTaq DNA Polymerase (Promega Corporation), and 18.2 ng of DNA template. PCRs used the following conditions: an initial step at 94 ° C for 2 min; followed by 35 cycles of 1 min at 94 ° C, 1 min at 50-58 ° C, and 50 s at 72 ° C; and a fi nal extension of 15 min at 72 ° C. All the reactions were conducted on a Mastercycler pro S thermocycler (Eppendorf, Hamburg, Germany). The PCR products were separated by electrophoresis on a 2.5% agarose gel and sent to Macrogen Europe sequencing service (Amsterdam, The Netherlands). In a second step, those primers that were polymorphic in the V. jacquinii-V. orbiculata complex were tested in two individuals from three species, each from a different clade ( V. orsiniana Ten. [core clade], V. javalambrensis Pau [Iberian clade], and V. rosea Desf. [North African clade]), using the same PCR conditions. Twelve polymorphic primer pairs were selected (see Appendix 2 for additional primers). Following the procedure developed by Schuelke (2000) , the sequence-specifi c forward primers were marked at the 5 ′ end with an M13 tail (5 ′ -TGTAAAACGACGGCCAGT-3 ′ ) (Eurofi ns), which was then labeled with 5-FAM, VIC, NED, or PET fl uorescent dyes ( Table 1 )  and 50 ng of DNA template in a total volume of 15 μ L. Conditions of the PCR amplifi cation were as described above, adding 10 cycles of 1 min at 94 ° C, 1 min at 53 ° C, and 50 s at 72 ° C before the fi nal extension. PCR products were analyzed with GeneMarker AFLP/Genotyping Software version 1.8 (SoftGenetics, State College, Pennsylvania, USA).

Population genetics parameters in three further species from
Veronica subsect. Pentasepalae -The fi rst comprehensive phylogenetic analysis of Veronica subsect. Pentasepalae based on DNA sequence data revealed four main clades each corresponding to a broad geographic area ( Rojas-Andrés et al., 2015 ). Thus, for the characterization of the microsatellite markers, diploid populations corresponding to species from different clades were selected (Appendix 1): V. orsiniana (core clade), V. javalambrensis (Iberian clade), and V. rosea (North African clade). The Central Asian clade was not considered because no material was available. The mean number of alleles per locus, observed and expected heterozygosities, possible deviations from Hardy-Weinberg equilibrium (HWE; Table 2 ) , and tests for linkage disequilibrium between markers in each population were estimated using Arlequin version 3.5.1.2 ( Excoffi er and Lischer, 2010 ).
The number of alleles per locus ranged from one to six, one to 11, and one to nine in the V. orsiniana , V. javalambrensis , and V. rosea populations, respectively. Loci 26, 49, and 52 were monomorphic in V. orsiniana , loci 10 and 52 were monomorphic in V. javalambrensis , and in V. rosea , loci 8 and 13 were monomorphic and locus 49 did not amplify. The observed and expected heterozygosities for all populations are shown in Table 2 . Signifi cant deviation from HWE ( P < 0.05) was seen for loci 8, 10, 13, and 54 in V. orsiniana , for locus 50 in V. javalambrensis , and for loci 10 and 50 in V. rosea . Linkage disequilibrium showed signifi cance levels below 0.05 after false discovery rate (FDR) correction in two pairwise comparisons (pair 20-52 in V. rosea and pair 27-54 in V. orsiniana ).

Cross-amplifi cation in other species from
Veronica subsect. Pentasepalae and 10 subgenera of Veronica -Cross-amplifi cation performed for these 12 polymorphic loci showed successful results within the expected allele size in two additional species from Veronica subsect. Pentasepalae : V. austriaca L. and V. dentata F. W. Schmidt. Tests were also performed for 20 additional species from 10 different subgenera within the large genus Veronica ( Table 3 ) . The tests were carried out with the original PCR protocol. The 12 loci tested in agarose gel showed successful amplifi cation of at least several bands. Six of these (8, 10, 13, 19, 26, and 35) showed good amplifi cation results in most samples.

CONCLUSIONS
A set of polymorphic microsatellite markers for Veronica subsect. Pentasepalae is reported. Amplifi cation success for these markers in the cross-transferability tests extends their potential usefulness to other subgenera. These markers will be useful for investigating genetic parameters, which may provide essential information for the conservation of threatened species, as well as data on the role of interspecifi c hybridization in the evolution of the genus.  c Populations used to generate the data included in Appendix 2.