Abstract:Background Single nucleotide polymorphisms (SNP) have been applied as important molecular markers in genetics and breeding studies. The rapid advance of next generation sequencing (NGS) provides a high-throughput means of SNP discovery. However, SNP development is limited by the availability of reliable SNP discovery methods. Especially, the optimum assembler and SNP caller for accurate SNP prediction from next generation sequencing data are not known. Results Herein we performed SNP prediction based on RNA-seq data of peach and mandarin peel tissue under a comprehensive comparison of two paired-end read lengths (125 bp and 150 bp), five assemblers (Trinity, IDBA, oases, SOAPdenovo, Trans-abyss) and two SNP callers (GATK and GBS). The predicted SNPs were compared with the authentic SNPs identified via PCR amplification followed by gene cloning and sequencing procedures. A total of 40 and 240 authentic SNPs were presented in five anthocyanin biosynthesis related genes in peach and in nine carotenogenic genes in mandarin. Putative SNPs predicted from the same RNA-seq data with different strategies led to quite divergent results. The rate of false positive SNPs was significantly lower when the paired-end read length was 150 bp compared with 125 bp. Trinity was superior to the other four assemblers and GATK was substantially superior to GBS due to a low rate of missing authentic SNPs. The combination of assembler Trinity, SNP caller GATK, and the paired-end read length 150 bp had the best performance in SNP discovery with 100% accuracy both in peach and in mandarin cases. This strategy was applied to the characterization of SNPs in peach and mandarin transcriptomes. Conclusions Through comparison of authentic SNPs obtained by PCR cloning strategy and putative SNPs predicted from different combinations of five assemblers, two SNP callers, and two paired-end read lengths, we provided a reliable and efficient strategy, Trinity-GATK with 150 bp paired-end read length, for SNP discovery from RNA-seq data. This strategy discovered SNP at 100% accuracy in peach and mandarin cases and might be applicable to a wide range of plants and other organisms.

Improved Annotation of the Peach (Prunus persica) Genome and Identification of Tissue- or Development Stage-Specific Alternative Splicing through the Integration of Iso-Seq and RNA-Seq Data

Deep RNA-Seq Uncovers the Peach Transcriptome Landscape.

Global Identification of Alternative Splicing Via Comparative Analysis of Smrt- and Illumina-Based Rna-Seq in Strawberry

A high-throughput SNP discovery strategy for RNA-seq data

Peach genomics and genome-wide association study: a review]

Transcriptomic Complexity in Strawberry Fruit Development and Maturation Revealed by Nanopore Sequencing

[Peach Genomics and Genome-Wide Association Study: a Review].

Differential transcript abundance and genotypic variation of four putative allergen-encoding gene families in melting peach

Alternative Splicing During Fruit Development among Fleshy Fruits

Genome-wide identification and characterization of long noncoding RNAs during peach (Prunus persica) fruit development and ripening

Multiple-statistical Genome-Wide Association Analysis and Genomic Prediction of Fruit Aroma and Agronomic Traits in Peaches.

An integrated peach genome structural variation map uncovers genes associated with fruit traits.

Whole-Genome Analysis of Diversity and SNP-Major Gene Association in Peach Germplasm.

Genome-Wide Analysis of the Gene Structure, Expression and Protein Interactions of the Peach (<i>Prunus persica</i>) <i>TIFY</i> Gene Family

Genome-Wide Analysis of the Gene Structure, Expression and Protein Interactions of the Peach (Prunus persica) TIFY Gene Family

De Novo Assembly, Functional Annotation, and Marker Development of Asian Pear (pyrus Pyrifolia) Fruit Transcriptome Through Massively Parallel Sequencing.

Development of a 775 SNP Array for Peach Based on Whole-Genome Resequencing Data, and Assessment of the Potential of Its Application

Transcriptome Analysis of Differentially Expressed Genes Relevant to Variegation in Peach Flowers

Genome re-sequencing reveals the evolutionary history of peach fruit edibility

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Detection and Application of Genome-Wide Variations in Peach for Association and Genetic Relationship Analysis