Development and characterization of EST-SSR markers in Ostryopsis (Betulaceae)1

• Premise of the study: A set of expressed sequence tag (EST) microsatellite markers were developed and characterized using next-generation sequencing technology for the Chinese genus Ostryopsis (Betulaceae). • Methods and Results: A total of 38 high-quality simple sequence repeat (SSR) primers were identified, of which 15 could be successfully amplified. Subsequently, we selected 80 individuals to represent the three species of the genus to evaluate the efficacy of these markers for examining genetic diversity of each species in the future. We found that the number of alleles per locus ranged from one to nine, with an average of 3.8. The expected heterozygosity and observed heterozygosity per locus varied from 0 to 0.829 and from 0 to 1, respectively, with their respective mean values as 0.483 and 0.416. • Conclusions: These EST-SSR markers will be useful for evaluating the range-wide genetic diversity of each species and examining genetic divergence and gene flow between the three species.

. Ostryopsis davidiana is mainly distributed in northern China, while O. nobilis and O. intermedia are limited to southwestern China. The northern vs. southwestern distributions of the three species in this genus suggest that Ostryopsis is a good model system to explore species divergence of plants in response to both habitat and temperature change across China. Until now, a total of 10 simple sequence repeat (SSR) markers have been developed for O. davidiana ( Qiu et al., 2009 ). However, some of these markers could not be successfully amplifi ed with the other two species. Furthermore, the development of expressed sequence tag (EST)-SSR markers is particularly attractive because they represent coding regions of the genome. EST-derived SSRs developed from one species can be easily amplifi ed and used in closely related species ( Wünsch, 2009 ). Here, we report 15 new polymorphic EST-SSR loci for Ostryopsis , which will facilitate the characterization of genetic diversity of each species in the genus and examination of gene fl ow and genetic divergence between the three species.

METHODS AND RESULTS
In this study, 20 or 30 individuals per species were collected for O. davidiana (Yijun, Shaanxi Province; Diebu, Gansu Province), O. nobilis (Tangdui, Yunnan Province; Jirenhe, Yunnan Province), and O. intermedia (Judian, Yunnan Province; Deqin, Yunnan Province) from their native distributions in China, and the voucher specimens were deposited in the herbarium of Lanzhou University (LZU), Lanzhou, China (Appendix 1). Genomic DNA was extracted using a cetyltrimethylammonium bromide (CTAB) procedure from leaves of each individual per species ( Ghangal et al., 2009 ). The total RNA of O. davidiana and O. nobilis was also extracted using a CTAB procedure ( Ghangal et al., 2009 ) and then the complementary DNA (cDNA) libraries were constructed and sequenced, respectively.
Approximately 4 μ g of RNAs were purifi ed using poly (dT)-conjugated beads (Life Technologies, Carlsbad, California, USA) to clear poly(A)-tagged mRNA. These RNAs were then broken into ~200-bp fragments under divalent cations at 75 ° C. We synthesized the fi rst strand of cDNA by the reverse transcriptase with random hexmer primers, and the second strand of cDNA by RNase H (Invitrogen, Ghent, Belgium) and DNA polymerase I (New England BioLabs, Ipswich, Massachusetts, USA). We sequenced the transcriptome on an Illumina (Solexa) Genome Analyzer II (Illumina Inc., San Diego, California, USA). After removing adapter sequences, we fi ltered and assembled two data sets of raw reads from O. davidiana and O. nobilis as described in Qiu et al. (2011) . To identify orthologous genes between two species, their reads were mapped back against the assembled unigenes using Bowtie 2 ( Langmead and Salzberg, 2012 ). We recalled single-nucleotide polymorphisms (SNPs) and indel calling with SAMtools ( Li et al., 2009 ). We identifi ed EST-SSRs using MISA (http://pgrc.ipk-gatersleben.de/misa) based on the orthologous unigene sequences. Because the de novo assembly introduced multiple indels, we therefore removed the shared ones by two orthologous unigenes. We used only the indels with end-to-end alignment as candidate SSR regions for further primer design. The fi nal analyses resulted in a total of 72 and 69 SSR indels according to the above reference unigenes from two species, respectively. We then designed the paired primers with Primer3 software ( Rozen and Skaletsky, 2000 ). Primers were not retained when they targeted to the SNP region of the unigenes. In this way, we obtained a total of 38 primer pairs from both species. • Premise of the study: A set of expressed sequence tag (EST) microsatellite markers were developed and characterized using next-generation sequencing technology for the Chinese genus Ostryopsis (Betulaceae). • Methods and Results: A total of 38 high-quality simple sequence repeat (SSR) primers were identifi ed, of which 15 could be successfully amplifi ed. Subsequently, we selected 80 individuals to represent the three species of the genus to evaluate the effi cacy of these markers for examining genetic diversity of each species in the future. We found that the number of alleles per locus ranged from one to nine, with an average of 3.8. The expected heterozygosity and observed heterozygosity per locus varied from 0 to 0.829 and from 0 to 1, respectively, with their respective mean values as 0.483 and 0.416. • Conclusions: These EST-SSR markers will be useful for evaluating the range-wide genetic diversity of each species and examining genetic divergence and gene fl ow between the three species.
temperature was decreased by 0.7 ° C for every cycle to a fi nal temperature of 53 ° C. This annealing temperature was employed for the last 30 cycles of the amplifi cation, followed by one cycle of 72 ° C for 10 min. PCR products were assayed on 1.5% (w/v) agarose gels to test the utility of the primers. Finally, a total of 15 primer pairs were successfully amplifi ed across all three species and all of them displayed clear polymorphisms ( Table 1 ) . Subsequently, we selected a total of 80 individuals from the three species to evaluate the potential value of these markers for estimating genetic diversity of each species. Fluorescence-based SSR genotyping was performed using a To evaluate polymorphisms of these primer pairs and their possible amplifications in the other species, we selected fi ve individuals from each of the three species. PCR reactions were carried out in a 20-μ L solution containing 20 ng of DNA template, 5 pmol of each primer, 100 μ M each of dNTPs, 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl 2 , and 0.2 U of Taq polymerase (Intron Biotechnology, Seongnam, South Korea). After a denaturing step of 5 min at 95 ° C, a touchdown amplifi cation program was performed. This profi le included a denaturing step of 45 s at 95 ° C and an extension step of 30 s at 72 ° C. The initial annealing step was 45 s at 60 ° C for one cycle, and subsequently the T ABLE 1. Characteristics of 15 polymorphic EST-SSRs developed in Ostryopsis .

Locus
Primer sequences (5 ′ -3 ′ ) Repeat motif Allele size (bp) T a ( ° C) GenBank accession no.  Hayden et al. (2008) . Briefl y, the forward primers of the 15 EST-SSRs were labeled with 6-FAM, VIC, or NED fl uorescent tags (Applied Biosystems, Foster City, California, USA). The PCR reactions were carried out separately for each microsatellite as described above. Amplifi cation products for which size and color did not overlap were pooled together for simultaneous detection of the amplifi ed alleles.
To characterize each EST-SSR marker, we calculated three genetic diversity statistics using POPGENE version 1.31 ( Yeh et al., 1999 ): number of alleles per locus, observed heterozygosity, and expected heterozygosity ( Table 2 ) . We found that the number of alleles ranged from one to nine, with an average of 3.8 alleles per locus. The expected heterozygosity and observed heterozygosity ranged from 0 to 0.829 and 0 to 1, respectively, with their respective mean values as 0.483 and 0.416 ( Table 2 ).

CONCLUSIONS
We developed 15 polymorphic EST-SSR markers for Ostryopsis from two cDNA libraries. The polymorphisms of these markers were further evaluated with 80 individuals representing the three species. These newly developed EST-SSRs have a high degree of universality between species and will be useful for studying the genetic diversity of each species and genetic divergence between the three species.