- Technical Note
- Open Access
PSR: polymorphic SSR retrieval
BMC Research Notesvolume 8, Article number: 525 (2015)
With the advent of high-throughput sequencing technologies large-scale identification of microsatellites became affordable and was especially directed to non-model species. By contrast, few efforts have been published toward the automatic identification of polymorphic microsatellites by exploiting sequence redundancy. Few tools for genotyping microsatellite repeats have been implemented so far that are able to manage huge amount of sequence data and handle the SAM/BAM file format. Most of them have been developed for and tested on human or model organisms with high quality reference genomes.
In this note we describe polymorphic SSR retrieval (PSR), a read counter and simple sequence repeat (SSR) length polymorphism detection tool. It is written in Perl and was developed to identify length polymorphisms in perfect microsatellites exploiting next generation sequencing (NGS) data. PSR has been developed bearing in mind plant non-model species for which de novo transcriptome assembly is generally the first sequence resource available to be used for SSR-mining. PSR is divided into two modules: the read-counting module (PSR_read_retrieval) identifies all the reads that cover the full-length of perfect microsatellites; the comparative module (PSR_poly_finder) detects both heterozygous and homozygous alleles at each microsatellite locus across all genotypes under investigation. Two threshold values to call a length polymorphism and reduce the number of false positives can be defined by the user: the minimum number of reads overlapping the repetitive stretch and the minimum read depth. The first parameter determines if the microsatellite-containing sequence must be processed or not, while the second one is decisive for the identification of minor alleles. PSR was tested on two different case studies. The first study aims at the identification of polymorphic SSRs in a set of de novo assembled transcripts defined by RNA-sequencing of two different plant genotypes. The second research activity aims to investigate sequence variations within a collection of newly sequenced chloroplast genomes. In both the cases PSR results are in agreement with those obtained by capillary gel separation.
PSR has been specifically developed from the need to automate the gene-based and genome-wide identification of polymorphic microsatellites from NGS data. It overcomes the limits related to the existing and time-consuming efforts based on tools developed in the pre-NGS era.
Simple sequence repeat (SSR), also known as short tandem repeats (STRs), variable number tandem repeats (VNTRs) or microsatellite markers, are repetitive stretches of 1–6 nucleotide units randomly spread in eukaryotic genomes. Microsatellites typically vary in length between 5 and 40 repeat units and can be classified as perfect or imperfect: in the first case, the repeat unit is arranged in head to tail manner without any interruption; in the second case the repeat unit includes insertions, deletions or substitutions of bases. In addition, SSRs can be of two types: simple or compound (i.e. two adjacent distinct SSRs are separated by a certain number of nucleotides). A further categorization allows SSRs to be recognized as genomic SSRs (gSSRs) and expressed sequence tag (EST)-SSRs; the latter result from transcribed regions of a genome, are associated with functionally annotated genes and are supposed to be more conserved among related species/genera. SSRs are very polymorphic due to the high mutation rate affecting the number of repeat units and have several advantages over other molecular markers: they allow the identification of multiple alleles at single locus, are equally distributed all over the genome, show a co-dominant inheritance pattern and could be mapped in different populations becoming the “backbone” of high resolution linkage maps. The development of EST-SSR markers was historically based on collection of sequence data of complementary DNAs (cDNAs) generated using the Sanger sequencing technology. Some useful tools have been made available for the detection and localization of both perfect and imperfect microsatellites along genome and mRNA sequences [1–6] (http://www.bioinformatics.org/ftp/pub/msatfinder/; http://pgrc.ipk-gatersleben.de/misa/).
Generally, these tools are equipped with additional programs able to design primers for microsatellite flanking regions. One of the most common approaches to the detection of length polymorphisms relies on the random selection of SSRs and on the analysis of PCR products on high-resolution gels . Since the application of this method is labour-intensive and time-consuming, it is more advisable to adopt an in silico approach to identify polymorphic SSRs and subsequently validate them through capillary gel separation. Relatively, little work has been done till now toward the automatic identification of polymorphic SSRs by exploiting EST sequence redundancy . With the advent of next-generation sequencing (NGS) methods, the large-scale identification of microsatellites became rapid and cost-effective  and was especially directed to non-model species with un-sequenced genomes [9–14]. But, high-throughput NGS technologies posed new challenges to existing bioinformatics applications that need to face a huge amount of sequence data and to handle new file formats. Until a few years ago, no adequate tools have been developed for the identification of polymorphic SSRs from NGS data with the exception of custom scripts still based on the detection of pad characters in the output of the final CAP3 assembly  or of the use of graphical viewer for NGS assemblies as starting point to manually explore and estimate microsatellite variability at each SSR locus [12, 15, 16]. It is clear, therefore, that it is no longer possible to work with human-readable text files including multiple sequence alignments of reads to the reference; rather it is essential to consider compressed text files in SAM/BAM format  to analyse sequence alignments and extract useful patterns. As far as we know, lobSTR , RepeatSeq  and the recently developed STR-FM  are the mainly publically available tools for profiling microsatellites from SAM/BAM data. All of them were conceived to identify gSSR alleles at each locus bearing in mind the short size of NGS reads and the difficulties this type of sequence data arises in the correct identification of polymorphic SSRs. An additional effort has been freshly published , but the method described is not publically available. However, these tools, even though capable of genotyping also EST-SSRs, were developed for and tested on human or model organisms for which the genome is available and represents a good reference for SSR-mining. Indeed, it seems that their adaptability is not likely to allow the easy use of custom reference transcriptomes and/or automatic comparative analysis at microsatellite loci. In this paper, we present polymorphic SSR retrieval (PSR) a Perl package developed to identify polymorphic SSRs from NGS data and provide quantitative information to each call. Unlike the aforementioned tools, PSR has been developed in a context where de novo transcriptome assembly is generally the first sequence resource for non-model plant species to be used for SSR-mining and downstream investigations.
PSR is a Perl package, conceived as a modular and flexible tool, developed to identify polymorphic SSRs from NGS data. The PSR package and the user manual are available at the following URL: http://sourceforge.net/projects/polyssr/.
The user guide provides general information on software dependencies and installation procedure as well as detailed instructions for running the application. PSR workflow is shown in Fig. 1 and the text that follows describes the key points. Psr_read_retrieval aims at the identification of all the reads that align to the reference sequences covering the full-length of perfect microsatellites. We decided to focus on perfect microsatellite only, because the total length of the reads obtained by the Illumina sequencing instruments rarely exceeds the 100 high quality nucleotides. These constraints appear to limit polymorphism discovery only to short polymorphic SSRs, but as it is evident from data in Table 1, the maximum length of the microsatellite is always slightly lower than the one of the sequenced read. In addition, the number of reads covering each microsatellite does not reflect the total number of reads at each locus. Indeed, in case SSRs are located at the ends of the read, the microsatellite-containing sequences are discarded since they can strongly affect the call of polymorphic sites (Fig. 2).
In addition, PSR filters out all the reads that match twice or more on the reference sequence as well as non-overlapping paired-end reads that aligned on the same microsatellite. This filtering step allows the user to exclude misalignments that can greatly influence downstream analysis.
Sometimes it could happen that the most frequent microsatellite identified by PSR presents a number of repeat units that differs from the one detected on the microsatellite-containing reference sequence. This is not surprising since the reference sequence was assembled considering all the reads; by contrast PSR takes into account only those reads that entirely overlap the repeating pattern.
Psr_read_retrieval takes two files as input: a 7-column tab-delimited file and a BAM file. The first one is a MISA-like output file (http://pgrc.ipk-gatersleben.de/misa/misa.html) that includes the type and the size of identified microsatellites and their localization within the reference sequences. The second input file stores read alignments against reference sequences for each genotype under investigation. The latter can be genome sequences or de novo assembled transcripts depending on the type of microsatellites under study.
Based on our experience and on results by Highnam et al. , we suggest bowtie2  as reference sequence alignment mapping tool. However, we are well-aware that read mapping tools are constantly being developed; for that reason we believe more convenient to stay independent from alignment algorithm.
Psr_read_retrieval should be run as many times as the number of genotypes under investigation (N in Fig. 1). It generates two different outputs: a collection of FASTA files including all the reads that successfully aligned along each SSR and different 3-column tab-delimited files, one for each genotype under investigation, where for each microsatellite the number of repeat units and the number of supporting reads are reported. In case of genomic DNA, results from this step of the analysis allowed ambiguities associated with repetitive stretches to be fixed as well as sequence assembly quality to be assessed.
Then, all the table-wise text files are fed into the psr_poly_finder module that, exploiting MySQL DBMS, creates a query table joining all the info recorded in text files in order to detect both heterozygous and homozygous alleles at each microsatellite locus. Users can define two different parameters, namely the minimum number of supporting reads (−t) and the minimum read depth (−p), for length polymorphism to be called. By default the values of these parameters are fixed to 10 and 30, respectively.
To give a clear example, a microsatellite is analysed in case at least ten reads overlap the repetitive stretch and it is called as heterozygous if 30 % or more of the reads supports the identification of the minor allele. In case of polyploid genotypes, the user can fine-tune the value of the “p” parameter in order to support the detection of all the putative alleles at each locus.
Filtering by a minimum number of reads supports variant calls and by minimum read depth within the microsatellite interval helps to reduce false positive calls.
Psr_poly_finder produces as final output a table-wise text file where the first two columns report the identifiers of microsatellite-containing sequences and the type of the identified microsatellite; the remaining columns include one or more (in case of heterozygosity) SSR allele sizes per genotype.
In order to evaluate PSR reliability and effectiveness, the package has been tested in the identification of polymorphic SSRs in (1) a collection of de novo assembled transcripts defined by RNA-sequencing of two different plant genotypes; (2) a set of newly sequenced chloroplast genomes.
In the first case, ten randomly selected polymorphic microsatellites, as they were identified by PSR, were assessed by capillary electrophoresis and eight of them showed the expected length polymorphism (Table 2). Only one SSR locus is heterozygous in case of genotype 1 and results from in silico and wet-lab experiments are in agreement. As for genotype 2, seven out of ten microsatellite loci resulted heterozygous even if PSR identifies both the alleles (major and minor) only in five cases.
In the latter case PSR combined with capillary gel separation was used to assess whether in silico identified length polymorphisms represented sequencing/assembling errors or, alternatively, real genotype-specific sequence variations. In seven cases out of eight, the size of the microsatellite supported by the highest number of reads was confirmed by capillary electrophoresis (Table 3). The use of PSR proved to be particularly useful to fix sequencing errors in monomeric microsatellites.
It is very likely that differences from in silico and wet-lab experiments depend on the amount of reads PSR uses to determine the number of repeated unit and microsatellite length. Indeed, for each SSR region PSR discards all the reads (on average the 60 % of the total) that do not completely span the microsatellite. Furthermore, in case of EST-SSRs, the quality of the reference transcriptome may negatively affect read alignment as a first step and, subsequently, read count. Finally, differences can be the result of PCR slippage products of repetitive stretches or can be ascribed to errors generated by PCR-based NGS technologies, especially in case of monomeric microsatellites .
The ability of next generation sequencing technologies to produce large amounts of sequence data accelerated microsatellite identification and facilitated the discovery of polymorphic SSRs. In addition, NGS resulted in the development of new specialised tools and file formats that support the management of huge amount of sequence data. However, as far as we know, few automatic and efficient systems have been developed so far for detecting both gene-based and genome-wide polymorphic microsatellites from NGS data. In this note we proposed PSR, a Perl Package developed for the detection of length polymorphisms and for the automatic comparative analysis at microsatellite loci. Indeed, the availability of polymorphic SSR markers is an essential requirement for plant community since it serve as basic resources for QTL mapping, molecular breeding, genetic diversity analysis and SSR fingerprinting for food traceability. All these and further applications, motivated us to develop and distribute PSR within plant community.
Availability and requirements
Project name: PSR—polymorphic SSR retrieval
Project home page: http://sourceforge.net/projects/polyssr/
Operating system(s): tested on Linux and Mac OS X
Programming language: Perl
Other requirements: BioPerl, MySQL
License: GNU GPL 3
Restrictions to use by non-academics: None.
polymorphic SSR retrieval
simple sequence repeat
short tandem repeat
variable number tandem repeats
expressed sequence tag-SSR
next generation sequencing
quantitative trait loci
database management system
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.
Kolpakov R, Bana G, Kucherov G. mreps: efficient and flexible detection of tandem repeats in DNA. Nucleic Acids Res. 2003;31(13):3672–8.
da Maia LC, Palmieri DA, de Souza VQ, Kopp MM, de Carvalho FI, Costa de Oliveira A. SSR locator: tool for simple sequence repeat discovery integrated with primer design and PCR simulation. Int J Plant Genom. 2008;2008:412696. doi:10.1155/2008/412696.
Tang J, Baldwin SJ, Jacobs JM, Linden CG, Voorrips RE, Leunissen JA, et al. Large-scale identification of polymorphic microsatellites using an in silico approach. BMC Bioinform. 2008;9:374. doi:10.1186/1471-2105-9-374.
Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001;11(8):1441–52. doi:10.1101/gr.184001.
Churbanov A, Ryan R, Hasan N, Bailey D, Chen H, Milligan B, et al. HighSSR: high-throughput SSR characterization and locus development from next-gen sequencing data. Bioinformatics. 2012;28(21):2797–803. doi:10.1093/bioinformatics/bts524.
Chandra A, Grisham MP, Pan YB. Allelic divergence and cultivar-specific SSR alleles revealed by capillary electrophoresis using fluorescence-labeled SSR markers in sugarcane. Genome/Natl Res Counc Can = Genome/Conseil national de recherches Can. 2014;57(6):363–72. doi:10.1139/gen-2014-0072.
Zalapa JE, Cuevas H, Zhu H, Steffan S, Senalik D, Zeldin E, et al. Using next-generation sequencing approaches to isolate simple sequence repeat (SSR) loci in the plant sciences. Am J Bot. 2012;99(2):193–208. doi:10.3732/ajb.1100394.
Iorizzo M, Senalik DA, Grzebelus D, Bowman M, Cavagnaro PF, Matvienko M, et al. De novo assembly and characterization of the carrot transcriptome reveals novel genes, new markers, and genetic diversity. BMC Genom. 2011;12:389. doi:10.1186/1471-2164-12-389.
Zhang H, Wei L, Miao H, Zhang T, Wang C. Development and validation of genic-SSR markers in sesame by RNA-seq. BMC Genom. 2012;13:316. doi:10.1186/1471-2164-13-316.
Shirasawa K, Koilkonda P, Aoki K, Hirakawa H, Tabata S, Watanabe M, et al. In silico polymorphism analysis for the development of simple sequence repeat and transposon markers and construction of linkage map in cultivated peanut. BMC Plant Biol. 2012;12:80. doi:10.1186/1471-2229-12-80.
D’Agostino N, Golas T, van de Geest H, Bombarely A, Dawood T, Zethof J, et al. Genomic analysis of the native European Solanum species, S. dulcamara. BMC Genom. 2013;14:356. doi:10.1186/1471-2164-14-356.
Liu Z, Chen T, Ma L, Zhao Z, Zhao PX, Nan Z, et al. Global transcriptome sequencing using the Illumina platform and the development of EST-SSR markers in autotetraploid alfalfa. PLoS One. 2013;8(12):e83549. doi:10.1371/journal.pone.0083549.
Xiao Y, Zhou L, Xia W, Mason AS, Yang Y, Ma Z, et al. Exploiting transcriptome data for the development and characterization of gene-based SSR markers related to cold tolerance in oil palm (Elaeis guineensis). BMC Plant Biol. 2014;14:384. doi:10.1186/s12870-014-0384-2.
Sio C-P Lu Y-L, Chen C-M, Pai T-W, Chang H-T, editors. Mining polymorphic SSRs from individual genome sequences. In: The seventh international conference on complex, intelligent, and software intensive systems (CISIS), Taichung; 2013.
Hoffman JI, Nichols HJ. A novel approach for mining polymorphic microsatellite markers in silico. PLoS One. 2011;6(8):e23283. doi:10.1371/journal.pone.0023283.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. doi:10.1093/bioinformatics/btp352.
Gymrek M, Golan D, Rosset S, Erlich Y. lobSTR: a short tandem repeat profiler for personal genomes. Genome Res. 2012;22(6):1154–62. doi:10.1101/gr.135780.111.
Highnam G, Franck C, Martin A, Stephens C, Puthige A, Mittelman D. Accurate human microsatellite genotypes from high-throughput resequencing data using informed error profiles. Nucleic Acids Res. 2013;41(1):e32. doi:10.1093/nar/gks981.
Fungtammasan A, Ananda G, Hile SE, Su MS, Sun C, Harris R, et al. Accurate typing of short tandem repeats from genome-wide sequencing data and its applications. Genome Res. 2015;25(5):736–49. doi:10.1101/gr.185892.114.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9. doi:10.1038/nmeth.1923.
CC and NDA designed the core of the Perl package. CC developed the package and contributed to the writing of the manuscript. NDA supervised PSR implementation and wrote the manuscript. All authors read and approved the final manuscript.
We sincerely thank Nunzia Scotti (CNR-IBBR), for sharing plastid genome-scale data and Andrea Allavena (CREA-FSO) for sharing RNA-sequencing data. We also thank Teodoro Cardi (CREA-ORT) for critical reading of the manuscript. This work was carried out in the frame of the “GenHORT—adding value to elite Campania horticultural crops by advanced genomic technologies” Project (PON02_00395_3215002) and is supported by the PON R&C 2007–2013 grant funded by the Italian Ministry of Education, University and Research in cooperation with the European Funds for the Regional Development (FESR).
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.