The genus Pedicularis L. (Orobanchaceae) comprises approximately 600 species of root-hemiparasitic plants that are distributed mostly in high-latitude or alpine habitats of the Northern Hemisphere (Yang et al., 1998). Within the traditional family Scrophulariaceae, Pedicularis was originally placed in tribe Rhinantheae Benth. Since then, this genus, along with other hemiparasitic rhinanthoids, has been transferred to Orobanchaceae based on molecular evidence and pollen morphology (Minkin and Eshbaugh, 1989; dePamphilis et al., 1997; Young et al., 1999; Olmstead et al., 2001). It is characterized by its diversification of corolla morphology, including variations in the galea (beaked, curved, toothed, or crested) and the length of the corolla tube, as a result of adaptive radiation (Li, 1951). Pedicularis ishidoyana Koidz. & Ohwi is a Korean endemic that is distinguished from its congeners by its long pedicels and undeveloped stems. This species is listed as Vulnerable (VU) in the Korea Red Data Book (Ministry of the Environment of Korea, 2012). As of 2012, it is also protected under the Endangered Species Act within Korean law. Populations are restricted to fewer than 10 locations in the lowlands of cool valleys, and are threatened by anthropogenic disturbances such as waterfront development. Pedicularis ishidoyana is considered a potentially important natural resource for medicinal products. The genus Pedicularis is known as “pseudo-ginseng” and is used for traditional medicines in East Asia. New pharmaceutical iridoids have also been discovered in its congener P. artselaeri Maxim. (Su et al., 1998). As part of the effort to preserve these threatened plants, we developed microsatellites using next-generation sequencing technology so that they can serve as valuable molecular tools for understanding population dynamics based on genetic diversity.
METHODS AND RESULTS
We collected 26 individuals of P. ishidoyana from a natural population at Mt. Geomma, Gyeongbuk, Korea, and deposited voucher specimens in the Inha University herbarium (IUI), Incheon, Korea (voucher no. Cho. 105024; 36°43′27″N, 128°14′52″E). Whole genomic DNA was extracted from silica gel-dried leaf tissues by a protocol that used the DNeasy Plant Mini Kit (QIA-GEN, Seoul, Korea). Measurements were made with a NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, Delaware, USA). High-quality DNA (concentration 186 ng·µL−1; A260/280 = 1.87; A260/230 = 2.3) was sequenced using the Illumina MiSeq platform (BML Co., Daejeon, Korea).
A total of 17,537,200 reads (300 × 300) were produced by Illumina paired-end sequencing and then trimmed and read by Trimmomatic 0.32 (Bolger et al., 2014). To identify the microsatellites from those reads, we screened them using SSR_pipeline version 0.951 (Miller et al., 2013). The parameters were set for detection of di-, tri-, tetra-, and pentanucleotide motifs with flanking regions larger than 40 bp and a minimum of 10, seven, five, and four repeats, respectively. In all, we found 63,531 microsatellite loci meeting the above criteria. Primer pairs were designed with Primer3 in Geneious R7 (Biomatters, available from http://www.geneious.com) and labeled via the M13 sequence tag method. To design primers efficiently, we attempted reference mapping of total reads to each microsatellite-containing singleton using Geneious R7. After discarding putative multicopy loci with exceptionally high coverage, we selected fragments with unique patterns that had two separate alleles, few variations at the site to which a primer was attached, and no additional single nucleotide polymorphisms (SNPs) in the flanking region. For the 26 tested individuals of P. ishidoyana, we designed 74 primer pairs and successfully amplified 32 of them.
We conducted PCR with 10-µL reaction volumes containing 5 µL of 2× PCR Plus Mix (400 µM dNTP, 4 mM MgCl2, 0.4 units of Taq DNA polymerase), 10 ng of DNA, 0.01 µM forward M13(−21)-tagged primer, 0.1 µM reverse primer, and 0.1 µM M13(−21)-labeled fluorescent marker (NED, PET, VIC, 6-FAM). Reactions were performed in a GeneAmp PCR System 2700 Thermal Cycler (Applied Biosystems, Foster City, California, USA) under the following conditions: denaturation at 94°C for 2 min; then 35 cycles at 94°C/30 s, 52°C to 60°C/1 min, and 72°C/1 min; and a final extension at 72°C for 7 min. Fluorescently labeled PCR products were resolved to genotype on an ABI 3730XL sequencer with GeneScan 500 LIZ Size Standard (Applied Biosystems). Sizes were determined with GeneMapper 3.7 (Applied Biosystems) and Geneious R7. Null allele frequencies were calculated by MICRO-CHECKER version 2.2 (van Oosterhout et al., 2004). The number of alleles plus values for expected heterozygosity (He) and observed heterozygosity (Ho) were determined in GenAlEx 6 (Peakall and Smouse, 2006). Deviations from Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were estimated with GENEPOP 4.0 (Rousset, 2008).
Characteristics of 18 microsatellite markers developed in Pedicularis ishidoyana.
Finally, 18 primer pairs proved polymorphic while the remaining 14 either were monomorphic or produced inconsequential peaks (Table 1). Alleles per locus numbered two to six (average 3.72), while values for He and Ho ranged from 0.142 to 0.703 and from 0.077 to 0.615, respectively (Table 2). Due to heterozygote deficiencies, two loci (Pi040 and Pi073) significantly deviated from HWE values after Bonferroni correction (P < 0.05). No significant LD was detected among locus pairs at the population level (P > 0.05).
In conclusion, we developed 18 microsatellite markers for the endangered P. ishidoyana. These markers will be informative tools for investigating genetic structure and diversity among populations of this species, and will help facilitate effective strategies for its conservation. They will also be useful in future studies to increase understanding of the phylogeographic history of the species based on gene flow and spatial genetic patterns.
Genetic properties of 18 newly developed, polymorphic microsatellites for Pedicularis ishidoyana.a