Open Access
Translator Disclaimer
5 June 2015 Isolation and Characterization of Microsatellite Markers for an Endemic Tree in East Asia, Quercus variabilis (Fagaceae)
Xian Wang, Jing Li, Yong Li
Author Affiliations +

Quercus variabilis Blume (Fagaceae) is a tree native to East Asia found mainly in warm-temperate and subtropical regions (Chen et al., 2012). Previous studies on this species have focused on its physiology and ecology (Lei et al., 2013; Du et al., 2014). To date, only a few studies have reported on the genetic diversity and population structure of this species using isozyme, chloroplast DNA, and microsatellite markers (Zhou et al., 2003; Xu et al., 2004; Chen et al., 2012). However, the polymorphic isozyme loci are limited in number and the chloroplast DNA is maternally inherited, which limits the comprehensive understanding of its genetic diversity and population structure. Xu et al. (2004) developed 16 microsatellite loci for Q. variabilis using primer pairs from three other Quercus L. species; however, these loci do not meet the requirement of high-density markers in landscape genetics (Hall and Beissinger, 2014). Thus, developing a greater number of microsatellite loci for this species is necessary for population genetic and landscape genetic studies of Q. variabilis to progress. Microsatellite markers, or simple sequence repeats (SSRs), are widely used in population genetic and landscape genetic studies, providing useful information to plant natural resources management and promoting future long-term strategies (Ohtsuki et al., 2014; Lei et al., 2015). Next-generation sequencing allows for rapid development of a large number of SSR markers (Huang et al., 2014; Kumar et al., 2014). Here, we report the development of microsatellite markers for Q. variabilis using Roche 454 pyrosequencing combined with a magnetic bead enrichment protocol.


Fresh leaves from 24 individuals of Q. variabilis were collected in two wild populations in Henan and Shandong provinces (Song Mountain: 34°28′08″N, 113°05′07″E; Meng Mountain: 35°29′34″N, 117°47′58″E) during May 2012. The samples were stored in silica gel until genomic DNA extraction and genotyping. Vouchers were deposited at the herbarium of Henan Agricultural University (Herbarium accession no. HHAU-2101-2112 [Song Mountain population] and HHAU-2113-2124 [Meng Mountain population]). Genomic DNA was extracted from approximately 30 mg of dried leaf tissue using a Tiangen Biotech (Beijing, China) DNA extraction kit (Plant no. DP305) following the manufacturer's protocol. A shotgun library was prepared by shearing 1 µg of genomic DNA using the DNA Library Preparation Kit (Roche Applied Science, Indianapolis, Indiana, USA) following the GS FLX+ library preparation protocol. The shotgun library was further enriched by eight 5′-biotinylated probes, namely, (AG)10, (AC)10, (AAC)8, (ACG)8, (AAG)8, (AGG)8, (ACAT)8, and (ATCT)6, and eventually sequenced on a Roche 454 GS FLX platform. A 1/16th run was performed using this pyrosequencing system. MISA software (Thiel et al., 2003) was used to search for microsatellites, while Primer3 software (Rozen and Skaletsky, 1999) was used to design primer pairs. Primer design parameters were set as follows: primer size = 18−27 nucleotides with 20 as optimum; primer melting temperature (Tm) = 55−65°C with 60°C as optimum; primer GC% = 40−60%; product size = 100−500 bp with 200 bp as optimum; and CG clamp = 3. PCRs were performed with a single microsatellite primer pair in a 30-µL reaction mixture containing 30 ng of genomic DNA, 0.2 mmol/L of each dNTP, 0.3 µmol/L of each primer, 3 µL 10× polymerization buffer, and 1 unit of Taq polymerase (Tiangen Biotech). PCR reactions were performed in a PTC-200 Thermal Cycler (MJ Research, Watertown, Massachusetts, USA) as follows: an initial denaturation at 95°C for 5 min; followed by 35 cycles at 94°C for 40 s, locus-specific annealing temperature (Table 1) for 45 s, and 72°C for 50 s; and a final extension at 72°C for 5 min. The PCR products were resolved on 8% native polyacrylamide gel and visualized by silver staining. The band size was determined in comparison with a 50-bp DNA ladder (TaKaRa Biotechnology Co., Dalian, Liaoning, China). Analyses of the number of alleles per locus (A), observed heterozygosity (Ho), expected heterozygosity (He), Hardy-Weinberg equilibrium (HWE), and linkage disequilibrium (LD) were performed using GENEPOP version 4.2 (Rousset, 2008). Analysis of the inbreeding coefficients (FIS) and its significance testing were performed using FSTAT version (Goudet, 1995).

Table 1.

Primer sequences and characterization for 34 microsatellite loci isolated from Quercus variabilis.


Table 2.

Genetic diversity parameters for two populations of Quercus variabilis.a


Our results suggest that next-generation sequencing technology combined with the magnetic bead enrichment protocol is a rapid and cost-effective means to identify large numbers of polymorphic microsatellite loci in Q. variabilis. By taking advantage of the Roche 454 sequencing platform, 63,851 individual sequence reads with an average length of 327 bp were obtained. A total of 2121 microsatellites were identified that contained enough flanking sequences for primer design. These sequences were deposited in GenBank (KP854041–KP856161). To test the amplification efficiency of these microsatellite loci, we randomly selected 100 loci to design primer pairs and perform PCRs. Sixty-six primer pairs were discarded for unsuccessful or nontarget amplification, whereas the remaining 34 primer pairs were tested for polymorphisms (Table 1). The 34 microsatellite sequences were compared with others in GenBank through BLAST searches. Qv040 and Qv089 matched significantly with the microsatellite CsCAT35 sequence (GenBank accession no. FJ868967) in Castanea sativa Mill. (Fagaceae) and REO_390 sequence (GenBank accession no. KJ207438) in Q. rubra L. (Fagaceae), respectively. Qv078 matched significantly with the mRNA sequence (GenBank accession no. XM_011077301) in Sesamum indicum L. (Pedaliaceae).

Fourteen of the 34 tested loci in 24 individuals from two Q. variabilis wild populations were polymorphic, with A ranging from two to six. He and Ho ranged from 0.042 to 0.721 and 0.042 to 0.750, respectively (Table 2). Locus Qv026 deviated significantly from HWE (P < 0.05) because of excessive homozygosity. FIS varied from −0.571 to 0.759 (Table 2), and no significant negative FIS values were observed. Three pairs of loci (Qv034 and Qv060, Qv034 and Qv041, Qv028 and Qv049) showed significant LD, indicating significant allelic associations between them. These new microsatellite loci will provide useful tools for further population genetic and landscape genetic studies on Q. variabilis.


In this study, we isolated 2121 microsatellite loci from Q. variabilis. One hundred microsatellite loci were randomly selected to design primer pairs and perform PCRs in two wild populations. Among these tested loci, 14 microsatellite markers were polymorphic. These new microsatellite loci increase the density of the microsatellite markers for this species, which will be helpful for further population genetic and landscape genetic studies on Q. variabilis.



D. Chen , X. Zhang , H. Kang , X. Sun , S. Yin , H. Du , N. Yamanaka , et al. 2012. Phylogeography of Quercus variabilis based on chloroplast DNA sequence in East Asia: Multiple glacial refugia and mainland-migrated island populations. PLoS ONE 7: e47268. Google Scholar


Q. R. Du , J. P. Lei , J. F. Liu , P. C. Wang , W. F. Xiao , and L. Pan . 2014. Eco-physiological response of Quercus variabilis seedlings to increased atmospheric CO2 and N supply. Chinese Journal of Applied Ecology 25: 24–30. Google Scholar


J. Goudet 1995. FSTAT (version 1.2): A computer program to calculate F-statistics. Journal of Heredity 86: 485–486. Google Scholar


L. A. Hall , and S. R. Beissinger . 2014. A practical toolbox for design and analysis of landscape genetics studies. Landscape Ecology 29: 1487–1504. Google Scholar


D. N. Huang , Y. Q. Zhang , M. D. Jin , H. K. Li , Z. P. Song , Y. G. Wang , and J. K. Chen . 2014. Characterization and high cross-species transferability of microsatellite markers from the floral transcriptome of Aspidistra saxicola (Asparagaceae). Molecular Ecology Resources 14: 569–577. Google Scholar


S. Kumar , N. Shah , V. Garg , and S. Bhatia . 2014. Large scale in-silico identification and characterization of simple sequence repeats (SSRs) from de novo assembled transcriptome of Catharanthus roseus (L.) G. Don. Plant Cell Reports 33: 905–918. Google Scholar


J. P. Lei , W. Xiao , J. F. Liu , D. Xiong , P. Wang , L. Pan , Y. Jiang , and M. H. Li . 2013. Responses of nutrients and mobile carbohydrates in Quercus variabilis seedlings to environmental variations using in situ and ex situ experiments. PLoS ONE 8: e61192. Google Scholar


Y. K. Lei , W. Wang , Y. P. Liu , D. He , and Y. Li . 2015. Adaptive genetic variation in the smoke tree (Cotinus coggygria Scop.) is driven by precipitation. Biochemical Systematics and Ecology 59: 63–69. Google Scholar


T. Ohtsuki , T. Shoda , Y. Kaneko , and H. Setoguchi . 2014. Development of microsatellite markers for Vitex rotundifolia (Verbenaceae). an endangered coastal plant in Lake Biwa, Japan. Applications in Plant Sciences 2: 1300100. Google Scholar


F. Rousset 2008. GENEPOP'007: A complete reimplementation of the GENEPOP software for Windows and Linux. Molecular Ecology Resources 8: 103–106. Google Scholar


S. Rozen , 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


T. Thiel , W. Michalek , R. K. Varshney , and A. Graner . 2003. Exploiting EST database for the development and characterization of gene-derived SSR markers in barley (Hordeum vulgare L.). Theoretical and Applied Genetics 106: 411–422. Google Scholar


X. L. Xu , L. A. Xu , M. R. Huan , and Z. R. Wang . 2004. Genetic diversity of microsatellites (SSRs) of natural populations of Quercus variabilis. Hereditas 26: 683–688. Google Scholar


J. Y. Zhou , J. Z. Guo , Z. S. Yang , and W. H. Zhang . 2003. Variation of peroxidase isozyme on natural populations of Quercus variabilis. Journal of Northwest Forestry University 18: 33–36. Google Scholar


[1] Funding for this project was provided by the National Natural Science Foundation of China (31100272) and Henan Education Department Nature Science Research (2011B220004).

Xian Wang, Jing Li, and Yong Li "Isolation and Characterization of Microsatellite Markers for an Endemic Tree in East Asia, Quercus variabilis (Fagaceae)," Applications in Plant Sciences 3(6), (5 June 2015).
Received: 24 March 2015; Accepted: 1 April 2015; Published: 5 June 2015

454 pyrosequencing
magnetic bead enrichment
microsatellite markers
Quercus variabilis
Get copyright permission
Back to Top