Primula L. section Petiolares Pax (Primulaceae) is mainly distributed in the Hengduanshan—Himalaya Mountains with only a few members occurring in Kashmir, central China, and other regions (Hu and Kelso, 1996). Primula ovalifolia Franch. and P. tardiflora (C. M. Hu) C. M. Hu are two closely related species in the section. Primula ovalifolia is widely distributed in southwestern and adjacent central China, around the Sichuan Basin, mainly growing in shaded habitats in broad-leaved forests and ravines, with altitudes ranging from 600 to 2500 m. Primula tardiflora is morphologically similar to P. ovalifolia but with a more limited distribution, known only from a single locality in E'mei Mountain. It was first considered as a subspecies of P. ovalifolia, and later was regarded as distinct from P. ovalifolia by its higher altitude habitat, later floral phenology, and some neglected vegetative traits (Hu and Kelso, 1996). A phylogeographic study based on chloroplast DNA data suggested that P. tardiflora is genetically close to the E'mei Mountain population of P. ovalifolia (Xie et al., 2012). Development of highly polymorphic nuclear markers (e.g., simple sequence repeats [SSRs]) will help to delineate between the two species. In Primula, only a few genomic SSR markers have been developed for several species thus far, such as P. vulgaris Huds. (Van et al., 2006), P. obconica Hance (Yan et al., 2010), P. sieboldii E. Morren (Ueno et al., 2011), P. veris L. (Bickler et al., 2013), P. poissonii Franch., and P. wilsonii Dunn (Zhang et al., 2013). Considering that these markers demonstrate low polymorphism and limited transferability, and that species of section Petiolares are phylogenetically distant from the above-named species, expressed sequence tag (EST)-SSRs specific to section Petiolares could provide useful tools for evolutionary and ecological studies in this section. In this study, we first obtained transcriptome data for P. ovalifolia using the Illumina platform and then designed marker pairs based on SSR loci. A subset of the markers was selected to investigate their polymorphism and transferability in congeneric species of Primula.
METHODS AND RESULTS
Transcriptome sequencing —Plants of P. ovalifolia were collected from Mount E'mei, Sichuan Province, China. Leaves from one individual were sampled, frozen immediately in liquid nitrogen, and stored at -80°C for RNA extraction and transcriptome sequencing. RNA was extracted using TRIzol Reagent (QIAGEN, Dusseldorf, Germany) and delivered to Genepioneer Technologies Corporation (Nanjing, China) for construction of cDNA libraries and sequencing. The cDNA libraries were sequenced on the Illumina HiSeq 2500 platform (Illumina, San Diego, California, USA) according to the manufacturer's recommendations. A total of 49,094,910 raw reads were obtained and deposited in the National Center for Biotechnology Information (NCBI) Short Read Archive under BioProject ID PRJNA379052 (accession no. SRP102475). Raw reads were first cleaned by trimming adapters and removing ambiguous reads (‘N’ > 10%) and low-quality reads (Phred score < 30). Clean reads were assembled into 142,468 transcripts using Trinity tools with default parameters (Haas et al., 2013) and were then clustered into 67,577 unigenes with TGICL version 2.1 (Pertea et al., 2003).
Development of EST-SSR markers—The EST-SSR loci were identified from unigenes longer than 1 kb using MIcroSAtellite identification tool (MISA) based on the Perl language (Thiel et al., 2003). We searched for SSRs with motifs ranging from mono- to hexanucleotides in size, and 4753 primer pairs were designated from 5139 putative loci using Primer3 web version 0.4.0 (Rozen and Skaletsky, 1999). A total of 220 markers comprising two nucleotides with at least eight contiguous repeat units were chosen for screening, among which 102 primers produced clear bands with suitable fragment lengths (<500 bp) during the preliminary test with four individuals of P. ovalifolia.
Results of initial primer screening in populations of Primula species.a
These 102 loci were further tested with eight individuals of P. ovalifolia. PCR reactions were performed with three primers: a sequence-specific forward primer with an M13(−21) tail at its 5′ end, a sequence-specific reverse primer, and the universal fluorescent-labeled M13(−21) primer (FAM, ROX, HEX, or TAMRA; Invitrogen, Guangzhou, Guangdong, China) (Schuelke, 2000). Genomic DNA was extracted using a modified cetyltrimethylammonium bromide (CTAB) protocol (Doyle, 1991). The amplified 10-µL mixture for SSRs included 5 µL of Master Mix (Generay Biotech, Guangzhou, China), 0.4 mM of each primer pair, 3.2 µL of deionized water, and 30–50 ng of genomic DNA. PCRs were run following a touchdown procedure with initial denaturation for 4 min at 94°C; followed by 10 cycles of 94°C for 35 s, 35 s at 60°C with an increment of −1°C per cycle, 45 s at 72°C; followed by 28 cycles of 94°C for 35 s, 35 s at 50°C, 45 s at 72°C; ending with an extra extension of 10 min at 72°C. PCR products were scanned by an ABI PRISM 3100 Genetic Analyzer using GeneScan 500 LIZ internal size standard (Invitrogen). Allele binning and calling were done using GeneMarker version 2.4.0 (SoftGenetics, State College, Pennsylvania, USA), and 38 primer pairs were selected for further polymorphism and transferability detection (Table 1). All of these SSR sequences have been deposited in GenBank (Table 1).
Polymorphism and transferability assessment—To assess the polymorphism level of these 38 loci, we genotyped 20–24 individuals in each of five populations from three species (Appendix 1). DNA extraction, PCR amplification, and length assessment of PCR products were performed following the procedures described above. Linkage disequilibrium among loci per population and deviation from Hardy–Weinberg equilibrium were tested using FSTAT version 2.9.3 (Goudet, 2001). We used GenAlEx 6.5 (Peakall and Smouse, 2012) to calculate the number of observed alleles per locus (A), expected heterozygosity (He), and observed heterozygosity (Ho).
No significant linkage disequilibrium was detected among loci after Bonferroni correction at α = 0.05 confidence level, and some loci showed significant deviations from Hardy—Weinberg equilibrium (Table 2). The 38 EST-SSRs displayed varied genetic diversity in three populations of P. ovalifolia (Table 2). A, Ho, and He for each locus ranged from one to 19, 0 to 0.938, and 0 to 0.915, respectively (Table 2). Excluding monomorphic loci, the polymorphic EST-SSR markers showed an average A of 5.2, 6.1, and 4.0; He of 0.587, 0.582, and 0.51; and Ho of 0.389, 0.421, and 0.40, in each population, respectively (Table 2). Out of the 38 SSR markers, 36 loci were successfully amplified in P. tardiflora and 31 loci showed polymorphism, with A ranging from two to six (Table 2). Similarly, 34 loci were successfully amplified in P. epilosa Craib, among which 23 loci showed polymorphism, with A ranging from two to six (Table 2). Overall, most of the EST-SSR markers developed for P. ovalifolia could be successfully cross-amplified, leading to a high transferability in the two congeneric species.
We developed and characterized 38 EST-SSR markers based on transcriptome sequencing of P. ovalifolia, a widely distributed species of Primula section Petiolares. These markers demonstrated high polymorphism in P. ovalifolia, with A ranging from one to 19, Ho from 0.000 to 0.938, and He from 0.000 to 0.915. Most of the markers could be successfully cross-amplified in congeneric species. These SSR makers are found to be useful tools for investigation of genetic structure and interspecific gene flow in this section.
The authors are indebted to Mingzhi Li and Jing Deng for their kind help in the laboratory. This work was supported by the National Natural Science Foundation of China (grant no. 31400209, U1603231) and the Ministry of Science and Technology of China (grant no. 2013FY111200).
- Bickler, C., S. A'Hara, J. Cottrell, L. Rogers, and J. Bridle. 2013. Characterisation of thirteen polymorphic microsatellite markers for cowslip (Primula veris L.) developed using a 454 sequencing approach. Conservation Genetics Resources 5: 1185–1187. Google Scholar
- Doyle, J. 1991. DNA protocols for plants: CTAB total DNA isolation. In G. M. Hewitt and A. Johnston [eds.], Molecular techniques in taxonomy, 283–293. Springer-Verlag, Berlin, Germany. Google Scholar
- Haas, B. J., A. Papanicolaou, M. Yassour, M. Grabherr, P. D. Blood, J. Bowden, M. B. Couger, et al. 2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nature Protocols 8: 1494–1512. Google Scholar
- Hu, C. M., and S. Kelso. 1996. Primulaceae. In Z. Y. Wu and P. H. Raven [eds.], Flora of China, vol. 15. Science Press, Beijing, China, and Missouri Botanical Garden Press, St. Louis, Missouri, USA. Google Scholar
- Peakall, R., and P. Smouse. 2012. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics (Oxford, England) 28: 2537–2539. Google Scholar
- Pertea, G., X. Huang, F. Liang, V. Antonescu, R. Sultana, S. Karamycheva, Y. Lee, et al. 2003. TIGR Gene Indices clustering tools (TGICL): A software system for fast clustering of large EST datasets. Bioinformatics (Oxford, England) 19: 651–652. Google Scholar
- Rozen, S., and H. Skaletsky. 1999. Primer3 on the WWW for general users and for biologist programmers. In S. Misener and S. A. Krawetz [eds.], Methods in molecular biology, vol. 132: Bioinformatics: Methods and protocols, 365–386. Humana Press, Totowa, New Jersey, USA. Google Scholar
- Schuelke, M. 2000. An economic method for the fluorescent labeling of PCR fragments. Nature Biotechnology 18: 233–234. Google Scholar
- Thiel, T., W. Michalek, R. Varshney, and A. Graner. 2003. Exploiting EST databases for the development and characterization of genederived SSR-markers in barley (Hordeum vulgare L.). Theoretical and Applied Genetics 106: 411–422. Google Scholar
- Ueno, S., Y. Yoshida, Y. Taguchi, T. Ujino-Ihara, N. Kitamoto, M. Honjo, R. Ohsawa, et al. 2011. Generation of expressed sequence tags, development of microsatellite and single nucleotide polymorphism markers in Primula sieboldii E. Morren (Primulaceae) for analysis of genetic diversity in natural and horticultural populations. Breeding Science 61: 234–243. Google Scholar
- Van, G. A., R. F. Van, I. Stiers, T. Sierens, J. Barker, and L. Triest. 2006. Isolation and characterization of microsatellite loci in primrose (Primula vulgaris). Belgian Journal of Botany 139: 261–264. Google Scholar
- Xie, X. F., H. F. Yan, F. Y. Wang, X. J. Ge, C. M. Hu. and G. Hao. 2012. Chloroplast DNA phylogeography of Primula ovalifolia in central and adjacent southwestern China: Past gradual expansion and geographical isolation. Journal of Systematics and Evolution 50: 284–294. Google Scholar
- Yan, H., X. Ge, C. M. Hu. and G. Hao. 2010. Isolation and characterization of microsatellite loci for the ornamental plant Primula obconica Hance (Primulaceae). HortScience 45: 314–315. Google Scholar
- Zhang, L., H. F. Yan, W. Wu, H. Yu, and X. J. Ge. 2013. Comparative transcriptome analysis and marker development of two closely related Primrose species (Primula poissonii and Primula wilsonii). BMC Genomics 14: 329. Google Scholar