Open Access
Translator Disclaimer
19 July 2017 Development of SSR Markers from Transcriptomes for Orinus (Poaceae), an Endemic of the Qinghai–Tibetan Plateau
Yu-Ping Liu, Xu Su, Wen-Chun Luo, Ting Lv, Ke-Long Chen, AJ Harris, Sayed Afzal Shah
Author Affiliations +

Simple sequence repeats (SSRs) have been widely used for DNA fingerprinting, molecular-assisted breeding, detecting gene locations, genetic diversity analyses, and evolutionary studies because they are codominant, are highly polymorphic, can be amplified repeatably, and provide many informative sites distributed throughout the genome (Agarwal et al., 2008; Izzah et al., 2014). SSRs that are developed from transcribed RNA sequences, known as expressed sequence tag SSRs (EST-SSRs), can be developed cheaply and efficiently using next-generation sequencing technology (Simon et al., 2009). Previously, microsatellite markers developed from transcriptomes have primarily been for woody and medicinal plants (e.g., Liu et al., 2014; Kim et al., 2017). Few microsatellite markers have been developed to study endemic alpine grasses.

Here, we show the utility of SSR markers derived from transcriptome assemblies to detect genomic variation in the three species of Orinus Hitchc. (Poaceae), which are endemic grasses that occur at high elevations in xeric, alpine areas of the Qinghai–Tibetan Plateau (QTP) (Su et al., 2015). Orinus comprises O. thoroldii (Stapf ex Hemsl.) Bor, which is mainly distributed in the western QTP, O. kokonoricus (K. S. Hao) Tzvelev, which occurs in the eastern QTP, and O. intermedius X. Su & J. Quan Liu, which is a hybrid of O. thoroldii and O. kokonoricus and has a range in the southeastern QTP (Su et al., 2015, 2017). Orinus has ecological and conservation value for its role in soil stabilization in the QTP, especially due to its expansive root system (Su et al., 2013), and this genus may represent a good system for elucidating the timing and mechanisms of desertification in the QTP (Su et al., 2015). Therefore, developing SSR markers for Orinus is a first step toward utilizing the population history of its species to better understand the origins of the genus and to improve conservation efforts for desert habitats in the QTP. Here, we develop microsatellite markers in Orinus using transcriptomes obtained via Illumina paired-end sequencing.

Table 1.

Characteristics of 16 polymorphic SSR markers developed in three species of Orinus.



Twenty-five to 30 individuals were collected from three populations each of the three Orinus species (248 total collections) (Appendix 1). In addition, representative individuals of O. thoroldii and O. kokonoricus were collected (Appendix 1), from which fresh leaves were obtained; these were immediately frozen in liquid nitrogen in the field and later stored at −80°C prior to RNA extraction. For all collections, voucher specimens were deposited at the Herbarium of the Northwest Plateau Institute of Biology (HNWP) (Appendix 1).

Total RNA and genomic DNA were extracted using a cetyltrimethylammonium bromide (CTAB) procedure (Ghangal et al., 2009). The RNA was quantified and its quality assessed using NanoDrop 2000 (Thermo Fisher Scientific, Waltham, Massachusetts, USA) and the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California, USA). The total RNA samples were purified to remove poly(A) tags using approximately 5 µg of RNAs and (dT)-conjugated beads (Life Technologies, Carlsbad, California, USA), and purified RNAs were divided into 200-bp fragments using divalent cations at 75°C. The first strand of cDNA was synthesized with reverse transcriptase and random hexamer primers, and the second strand was synthesized by RNase H (Invitrogen, Ghent, Belgium) and Taq DNA polymerase I (New England BioLabs, Ipswich, Massachusetts, USA). Finally, the cDNAs representing the transcriptome were sequenced on an Illumina (Solexa) Genome Analyzer II (Illumina, San Diego, California, USA). All sequence information has been deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (Bioproject no. PRJNA385721).

SSR primers were designed using the transcriptomes of O. thoroldii and O. kokonoricus. To accomplish this, 79,111,257 and 67,617,602 raw reads were first filtered from O. thoroldii and O. kokonoricus. As a result, 60,124,556 and 50,037,026 high-quality reads, respectively, were obtained for de novo assembly. De novo assembly was performed in Trinity 2.2 with default parameters (Grabherr et al., 2011). From Trinity, 23,029 high-quality contigs were obtained from O. thoroldii and 24,086 from O. kokonoricus, representing total lengths of 20,808,832 and 22,281,570 bp, respectively, with an average size of 903 bp (N50 = 1188 bp) and 925 bp (N50 = 1203 bp). The assembly of each species was used to map the reads of the other using Bowtie 2 (Langmead and Salzberg, 2012) to identify an orthologous set of genes. Within the orthologous set, we searched for candidate SSRs using MISA 4.0 ( and also identified single-nucleotide polymorphisms (SNPs) and insertion/deletion polymorphisms (indels) using SAMtools 1.4 (Li et al., 2009). Indels were selected as SSR candidates only if the indels appeared informative between the two species and could be confidently aligned at their 3′ and 5′ ends. In total, 58 polymorphic candidate loci were recovered from O. thoroldii and 52 from O. kokonoricus. For each locus, primers were designed in Primer3 (Rozen and Skaletsky, 2000), and primers with binding sites containing SNPs were rejected. Thus, a total of 50 primers suitable for both species of Orinus were found. Ten individuals were selected from each of three species of Orinus to test the 50 primers. Amplification was performed using a standard 25-µL PCR reaction containing 1.00 µL of template DNA, 0.20 mM MgCl2, 0.25 mM dNTPs, 2.00 µM of each primer, 2.50 µL of 10× PCR buffer, 0.25 µL of Taq DNA polymerase (5 U·µL-1; TaKaRa Biotechnology Co., Dalian, Liaoning, China), and distilled water up to the final volume. The amplification was carried out under the following thermocycling protocol: enzyme initiation at 94°C for 5 min; followed by 36 cycles of denaturation at 94°C for 50 s, annealing at 49–58°C for 50 s (Table 1), and 72°C for 1 min; and a final extension at 72°C for 10 min. Amplification success was determined by viewing the results in a 1.5% (w/v) agarose gel. Sixteen primers showed successful amplification in all three Orinus species, and all of them possessed clear polymorphisms (Table 1).

Subsequently, the 16 primers were used to carry out fluorescence-based genotyping for all 248 sampled individuals of Orinus. Fluorescence-based genotyping was performed using a modification of the method presented in Hayden et al. (2008). In brief, the forward primers for the 16 SSRs were labeled with 6-FAM fluorescent tags (Applied Biosystems, Foster City, California, USA) and PCR reactions were performed as described above. The labeled products were detected on an ABI 3730XL sequencer with the GeneScan 500 LIZ Size Standard (Applied Biosystems). The profiles of the amplified loci were examined using GeneMapper 3.7 (Applied Biosystems), and peaks were scored manually by visual inspection. Each SSR marker was characterized by calculating three measures of genetic diversity in GENEPOP 4.2 (Rousset, 2008): number of alleles per locus, observed heterozygosity, and expected heterozygosity (Table 2). By these measures, the markers were highly polymorphic, with the number of alleles ranging from one to seven within this genus, with an average of 2.6 alleles per locus. The expected and observed heterozygosity ranged from 0.00 to 0.83 and 0.00 to 1.00, respectively, with respective mean values of 0.32 and 0.34 (Table 2).

Table 2.

Genetic diversity statistics for each sampled population of the three Orinus species based on 16 pairs of SSR primers.a



We have developed 16 polymorphic SSR markers from two cDNA libraries for investigating population structure in Orinus. These markers amplified easily and showed considerable polymorphisms for 248 individuals from three populations each of the three species of Orinus. These markers represent valuable new tools that will facilitate the development of Orinus as a model for understanding the origins and phylogeographic processes of the alpine desert of the QTP.


The authors thank Prof. Jianquan Liu for providing laboratory space and funding and thank the reviewers for their comments on the manuscript. This study was supported by the National Natural Science Foundation of China (no. 31260052), the Natural Science Foundation of Qinghai Province (no. 2016-ZJ-937Q, 2014-ZJ-947Q), the State Scholarship Fund of China Scholarship Council (no. 201508635003), and the Young and Middle-Aged Research Foundation of Qinghai Normal University (no. 2015-19).



Agarwal, M., N. Shrivastava, and H. Padh. 2008. Advances in molecular marker techniques and their applications in plant sciences. Plant Cell Reports 27: 617–631. Google Scholar


Ghangal, R., S. Raghuvanshi, and P. C. Sharma. 2009. Isolation of good quality RNA from a medicinal plant sea buckthorn, rich in secondary metabolites. Plant Physiology and Biochemistry 47: 1113–1115. Google Scholar


Grabherr, M. G., B. J. Haas, M. Yassour, J. Z. Levin, D. A. Thompson, I. Amit, X. Adiconis, et al. 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology 29: 644–652. Google Scholar


Hayden, M. J., T. M. Nguyen, A. Waterman, G. L. McMichael, and K. J. Chalmers. 2008. Application of multiplex-ready PCR for fluorescencebased SSR genotyping in barley and wheat. Molecular Breeding 21: 271–281. Google Scholar


Izzah, N. K., J. Lee, M. Jayakodi, S. Perumal, M. Jin, B. S. Park, K. Ahn, et al. 2014. Transcriptome sequencing of two parental lines of cabbage (Brassica oleracea L. var. capitata L.) and construction of an EST-based genetic map. BMC Genomics 15: 149. Google Scholar


Kim, B. Y., H. S. Park, S. Kim, and Y. D. Kim. 2017. Development of microsatellite markers for Viscum coloratum (Santalaceae) and their application to wild populations. Applications in Plant Sciences 5: 1600102. Google Scholar


Langmead, B., and S. Salzberg. 2012. Fast gapped-read alignment with Bowtie 2. Nature Methods 9: 357–359. Google Scholar


Li, H., B. Handsaker, A. Wysoker, T. Fennell, J. Ruan, N. Homer, G. Marth, et al. 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25: 2078–2079. Google Scholar


Liu, B. B., T. Bin, H. Ma, Z. Q. Lu, Q. Qiu, K. S. Mao, and J. Q. Liu. 2014. Development and characterization of EST-SSR markers in Ostryopsis (Betulaceae). Applications in Plant Sciences 2: 1300062. Google Scholar


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


Rozen, S., and H. Skaletsky. 2000. 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


Simon, S. A., J. Zhai, R. S. Nandety, K. P. McCormick, J. Zeng, D. Mejia, and B. C. Meyers. 2009. Short-read sequencing technologies for transcriptional analyses. Annual Review of Plant Biology 60: 305–333. Google Scholar


Su, X., W. Yue, and J. Q. Liu. 2013. Germplasm collection and preservation of Orinus (Poaceae) in the Qinghai-Tibet Plateau. Plant Diversity and Resources 35: 343–347. Google Scholar


Su, X., G. L. Wu, L. L. Li, and J. Q. Liu. 2015. Species delimitation in plants using the Qinghai-Tibetan Plateau endemic Orinus (Poaceae: Tridentinae) as an example. Annals of Botany 116: 35–48. Google Scholar


Su, X., Y. P. Liu, G. L. Wu, W. C. Luo, and J. Q. Liu. 2017. A taxonomic revision of Orinus (Poaceae) with a new species, O. intermedius, from the Qinghai-Tibet Plateau. Novon 25: 206–213. Google Scholar


Appendix 1.

Locality information for the populations of Orinus used in this study.

Yu-Ping Liu, Xu Su, Wen-Chun Luo, Ting Lv, Ke-Long Chen, AJ Harris, and Sayed Afzal Shah "Development of SSR Markers from Transcriptomes for Orinus (Poaceae), an Endemic of the Qinghai–Tibetan Plateau," Applications in Plant Sciences 5(7), (19 July 2017).
Received: 28 March 2017; Accepted: 1 May 2017; Published: 19 July 2017

next-generation sequencing
orthologous gene
simple sequence repeat (SSR) marker
Get copyright permission
Back to Top