Skip to main content

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.)



Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish.


In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions.


These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.


The advent of next-generation sequencing (NGS) has greatly advanced our opportunities to obtain abundant sequences data by resequencing plant genome and transcriptome in greater depth [1, 2]. RNA-Seq technology had been used to rapidly isolate and identify alternative splicing (AS) [3, 4], novel transcripts [5] and structural variations [6]. Surprisingly, the transcriptome data offers rich resources to identify and develop large number of single nucleotide polymorphism (SNP) and insertion/deletion (InDel) markers based on comparing transcriptome sequences to the reference genome sequences [7, 8]. This approach offers a direct, reliable and high-efficient strategy for identification and development of SNPs and InDels using data mining by bioinformatic methods. Recently, using comparative transcriptomic analysis, a large number of SNP and InDel markers have been identified in many plant species including Arabidopsis [7], Brassica napus [9], rice [10] and tomato [11]. However, a limited number of SNP and InDel markers have been developed in transcriptome data of radish [12].

Molecular markers related to phenotypes provide us with valuable tools for elucidation of genetic variations. In the past two decades, several molecular marker systems, such as restriction fragment length polymorphism (RFLP), random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP) and simple sequence repeats (SSR) markers had been widely employed in genetic diversity analysis, and genetic mapping studies in many plant species [10, 13,14,15]. Recently, SNP and InDel markers has increasingly being received attention owing to their high frequency, wide distribution, co-dominant inheritance and abundant DNA polymorphisms [16, 17]. Previous studies indicated that the SNPs and InDels within or near the coding sequences could be influenced by certain phenotype, transcription factor binding and alternative splicing [18, 19]. SNP and InDel markers have been extensively utilized in several crop improvement programs including genetic diversity analysis, quantitative trait locus (QTL) mapping and genome-wide association analyses (GWAS) [20, 21].

Alternative splicing, generating multiple transcript variants from a single pre-mRNA and increasing the transcriptome complexity, is one of the most important regulatory mechanisms for gene expression and functional diversity of proteome [22, 23]. In plants, it has been estimated that 60% of Arabidopsis genes [24], 33% of rice genes [23], 40% of soybean genes [25] and 59.3% of tomato genes [26] could be alternatively spliced. Exon-skipping (ES) events were the predominant AS type in animals, while intron retention (IR) was more frequent in plants [27]. Moreover, AS played crucial roles in regulation of biological functions, particularly in developmental processes and stress responses [24]. Genome-wide analysis of AS events had been surveyed in tomato, which showed that the splice variants of multi-exon genes in the seedlings and flowers were lower than those in early growth fruits [26]. In six teosinte and 10 maize transcriptomes, a large number of genes were induced in AS complexity during domestication from teosinte to maize [3]. These results provided useful information for identification and characterization of AS in other plant species.

Radish (Raphanus sativus L., 2n = 2× = 18) is an economically annual or biennial root vegetable crop of the Brassicaceae family. Although a set of SNP and/or InDel markers had been developed from the public radish EST database or whole genome resequencing data [28, 29], the number of SNP and InDel markers are far from saturated for linkage mapping and association study in radish. Recently, the release of radish reference genome [30] and availability of de novo transcriptome datasets [31,32,33] generated a great number of transcriptome sequences, which provided useful sequence information for identification and development of SNP, InDel and SSR markers in radish.

To better systematically characterize the overall transcriptional landscape in radish, three available radish transcriptomes data were collected in this study. The alignment of the short reads to the radish genome dataset were conducted for identification of AS events, novel mRNA transcripts as well as SNP, InDel and SSR markers among these three radish transcriptomes. Additionally, the frequency, distribution and polymorphic of SNPs and InDels were estimated and characterized in three radish genotypes. The newly developed InDel markers were applied to investigate the genetic diversity among 32 accessions. The outcomes of this study could enhance our understanding of radish gene structure and facilitate marker-assisted selection and genetic improvement of some important horticulture traits in radish breeding programs.


Overview of the RNA-Seq data

To identify alternative splicing events and molecular marker in radish, three transcriptome datasets of the radish genotypes ‘NAU-RG’, ‘NAU-LB’ and ‘NAU-YH’ were collected and used for further analysis. After filtering the low-quality reads and adapter sequences, 64,000,056, 49,631,196 and 50,386,238 clean reads were obtained in ‘NAU-RG’, ‘NAU-LB’ and ‘NAU-YH’ libraries, respectively, from which 50,047,249 (78.2%), 34,770,143 (70.06%) and 37,512,238 (74.45%) clean reads in ‘NAU-RG’, ‘NAU-LB’ and ‘NAU-YH’ libraries, respectively, were mapped onto the reference genome sequence of radish (Table 1). Totally, 49,481,475 (77.31%), 34,413,111 (69.34%) and 37,111,575 (73.65%) clean reads in the ‘NAU-RG’, ‘NAU-LB’ and ‘NAU-YH’ libraries, respectively, were mapped to unique genome locations. The uniquely mapped reads were aligned to the exon, intron and intergenic regions of reference genome (Fig. 1).

Table 1 Summary of clean reads and mapped to the radish reference genome
Fig. 1
figure 1

Overview of RNA-Seq reads mapping to radish reference genome. The percentage of clean reads mapped to reference genome regions in the NAU-RG, NAU-LB and NAU-YH libraries

Identification of AS events and discovery of novel mRNA transcripts

In all, 170,086 unique splice junctions sites were identified using TopHat program [34] (Additional file 1: Table S2); among which 131,106 (77.1%) were annotated in the radish genome and 38,980 (22.9%) were novel splice junctions. As shown in Additional file 1: Table S2, the ‘NAU-LB’ had the largest number of junctions, followed by ‘NAU-YH’ and ‘NAU-RG’. Among the dinucleotides at the intron border, a total of 218,574 GT-AG splice sites (95.9%), 163 GC-AG splice sites (0.1%), and 9256 splice sites of other types (4.0%) were identified (Fig. 2a).

Fig. 2
figure 2

Overall characteristics of AS events. a The proportions for the three types of dinucleotides at the splicing border. b Number and frequency of the different AS types. (IR) Intron retention, (A3SS) alternative 3′ splice site, (A5SS) alternative 5′ splice site, (ES) exon skipping. c Four genes representing four types of AS events were validated by RT-PCR. The lanes from left to right are DNA ladder, leaf and root of ‘NAU-RG’, ‘NAU-LB’, ‘NAU-YH’, and genomic DNA of ‘NAU-YH’. The RNA-Seq read density of the four genes was viewed by IGV browser (green peaks). The red arrows indicate alternative splice sites. The blue arcs indicate splice junction reads that support the junctions

To analyze the AS events in radish, the ASTALAVISTA software was employed to extract the AS events and classify different types of alternative splicing [35]. In total, 56,530 AS events were identified in three transcriptome datasets. These AS events were distributed in 36,546 genes, which accounted for 59.4% of the total radish genes, with an average of 1.68 AS events per gene. In this study, 25.0% of the radish AS-related genes undergo multiple AS events resulting in a variety of transcripts from a single gene (Additional file 1: Figure S1). The splice junctions of these genes were used to identify and categorize the AS events into five main types of AS events as intron retention (IR), exon skipping (ES), alternative 3′ splice site donor (A3SS), alternative 5′ splice site acceptor (A5SS) and other events. Among the different AS types, IR events were the predominant AS type (27,254, 48.2%), followed by other (19,674, 34.8%), A3SS (4228, 7.5%), and A5SS (3995, 7.1%) events (Fig. 2b, Additional file 2). To verify the accuracy of AS events, a reverse transcription-PCR (RT-PCR) analysis was performed on the leaf and root of ‘NAU-RG’, ‘NAU-LB’ and ‘NAU-YH’, respectively. Nine AS events from nine genes were exhibited to produce splice variants, which were consistent with the RNA-Seq data results. The RT-PCR results indicated that the expression of AS events was tissue-specific, for example, Rsa1.0_01846.1_g00001.1 was expressed in the root of ‘NAU-LB’, while not expressed in the leaf (Fig. 2c, Additional file 1: Table S3).

To detect the differential splicing events, the AS events consisted of A3SS, A5SS, IR, ES and mutually exclusive exons (MXE) were analyzed using the Multivariate Analysis of Transcript Splicing (MATS) program [36]. The maximum number of differential AS (1360) was identified in ‘NAU-LB’ vs ‘NAU-RG’, of which RI (459 events in 422 genes) possessed the predominant events (Fig. 3a, Additional file 3). Gene ontology (GO) enrichment analysis showed that a total of 336 (79.6%) differentially spliced genes were assigned into 43 GO terms in the comparison between ‘NAU-LB’ and ‘NAU-RG’ (Fig. 3b). In terms of three main categories, ‘cellular process’, ‘cell’ and ‘binding’ represented the most abundant GO terms in biological process, cellular component and molecular function, respectively. In ‘NAU-LB’ and ‘NAU-RG’ comparison, the Rsa1.0_01282.1_g00001.1 with IR event was selected and validated (Fig. 3c). A novel isoform was found in ‘NAU- RG’, but not in ‘NAU- LB’ by RT-PCR (Fig. 3d). In this study, the novel and unknown genes were defined as novel transcripts, a total of 10,657 novel transcripts were detected based on the transcriptome data (Additional file 4). The newly identified transcripts would be useful for updating the radish genome annotation.

Fig. 3
figure 3

Characterization of differential splicing events in three genotypes. a Distribution of differential splicing events in three pairwise comparisons. b Gene ontology classification of genes with differential IR events in ‘NAU-LB’ vs ‘NAU-RG’. c Reads distribution of Rsa1.0_01282.1_g00001.1 derived from RNA-Seq data of ‘NAU-LB’ and ‘NAU-RG’. The horizontal and vertical axis represents per-base expression and genomic coordinates, respectively. mRNA isoform quantified is given at the bottom, black box and line with arrow heads represents exon and intron, respectively. d Validation of Rsa1.0_01282.1_g00001.1 with differential IR event in the ‘NAU-LB’ and ‘NAU-RG’ comparison by RT-PCR

Identification and characterization of SNPs, InDels and SSRs in radish

To systematically identify genetic variations, three radish transcriptomes were utilized for SNPs and InDels discovery using Genome Analysis Toolkit [37]. Overall, 22,412 putative SNPs were detected among these three genotypes with an average frequency of 1 SNP/17.9 kb in comparison with the reference sequence (Additional file 5). Of these, 19,013, 20,547 and 20,551 SNPs were identified in ‘NAU-LB’, ‘NAU-RG’ and ‘NAU-YH’, respectively. A total of 16,997 SNPs were shared by all three genotypes, while 439, 571 and 385 SNPs were specific for ‘NAU-YH’, ‘NAU-LB’ and ‘NAU-RG’, respectively (Fig. 4a). The identified SNPs were distributed on the 8248 scaffolds, and the majority of scaffolds only had one SNP (Fig. 4b).

Fig. 4
figure 4

Summary of SNPs and InDels between three genotypes. Number of SNPs (a) and InDels (c) presented in three genotypes; Number of SNPs (b) and InDels (d) on the scaffold. The intersecting portions represent the common loci between three genotypes

SNPs were further classified based on their zygosity and nucleotide substitution. Between these three pairwise comparisons, 14,583 SNPs (13,241 homozygous and 1342 biallelic SNPs) detected between ‘NAU-RG’ and ‘NAU-YH’ were more than those in the other two transcriptome comparisons (Table 2). A total of 3838 (17.1%) and 3209 (14.3%) homozygous SNPs were presented in ‘NAU-LB’ vs ‘NAU-RG’ comparison and ‘NAU-LB’ vs ‘NAU-YH’ comparison, respectively. There were 1342 homozygous SNPs between ‘NAU-RG’ vs ‘NAU-YH’ comparison, which comprised only 6.0% of the total SNPs, indicating that the self-incompatible ‘NAU-YH’ presented a larger degree of heterozygosity than other two accessions (Table 2). Moreover, the frequency of transitions (58.5%) was higher than transversions (41.5%) with transition to transversion ratio (Ts/Tv) of 1.41:1 (Table 2). Across all pairwise comparisons, the Ts/Tv of bi-allelic SNPs was slightly more than the homozygous SNPs, the overall Ts/Tv being relatively stable between genotypes with a range from 1.29 to 1.4.

Table 2 SNPs zygosity and corresponding transition to transversion ratio (Ts/Tv) in pairwise transcriptome comparisons

Compared to the reference sequence, a total of 9436 InDels were identified with an average density of 1 InDel/42.5 kb (Fig. 4c, Additional file 6). These InDels were distributed on 3661 scaffolds with the length ranged from 1 to 13 bp, and for 1873 InDels (51.2%), each was located in one scaffold (Fig. 4d). The single-nucleotide InDels (82.7%) was the dominant type, followed by bi- (14.7%) and tri- (1.3%) nucleotide InDels (Additional file 1: Table S4). A comparative analysis showed that the maximum (8508) and minimum (7986) number of InDels were presented in ‘NAU-RG’ vs ‘NAU-YH’ and ‘NAU-LB’ vs ‘NAU-YH’ comparison, respectively, of which 2066 (21.9%) and 2835 (30.0%) InDels were homozygous in ‘NAU-RG’ vs ‘NAU-YH’ and ‘NAU-LB’ vs ‘NAU-YH’ comparison, respectively.

To identify the potential SSRs, a total of 104,801 assembled unigenes were obtained from three transcriptome datasets. In total, 43,680 potential SSRs were identified in 31,604 unigenes with a density of 1 SSR/2.5 kb, among which 8009 (25.0%) possessed more than one SSR (Additional file 7). Among the identified SSRs, the mono-nucleotide repeats were the most abundant (21,599, 49.4%), followed by tri- (10,529, 25.6%), di- (11,164, 24.1%), and tetra- (287, 0.7%) nucleotide repeats, respectively (Fig. 5). For the mono-nucleotide motifs, A/T was the primary motif. AG/CT and GA/TC were the most common motifs among the di-nucleotide motifs (Additional file 1: Table S5). The SSRs identified in this study would facilitate genetic linkage map construction and marker-assisted selection (MAS) in radish breeding programs.

Fig. 5
figure 5

Distribution of microsatellites with different motif lengths in radish

Analysis of genetic variation between three genotypes

To investigate the distributions of SNPs and InDels in genic and intergenic regions, the location of genetic variation was compared with the annotation of the reference gene models. Among the 22,412 SNPs, 11,236 (50.1%) SNPs were located in intergenic regions, 5483 (24.5%) in intron, and 5693 (25.4%) in coding sequence (CDS); and the SNP frequencies were 1 SNP/35.8, 73.3 and 70.6 kb, respectively (Table 3). For SNPs in the CDS regions, the maximum number of 4895 was detected in ‘NAU-RG’ and ‘NAU-YH’ comparison, with a frequency of 1 SNP/82.1 kb; the lowest number of 4536 SNPs was observed in ‘NAU-LB’ vs ‘NAU-RG’ comparison with a frequency of 1 SNP/88.6 kb. Moreover, 1413 (15%) InDels were located in intergenic regions, 1135 (12%) in intron, 6888 (73%) in CDS, with a frequency of 1 InDel/58.4, 354.2 and 284.5 kb, respectively. The variations of InDel in the CDS regions, the maximum (6301, 1 InDel/0.06 Mb) and minimum (6015, 1 InDel/0.07 Mb) numbers of InDels also existed in ‘NAU-RG’ vs ‘NAU-YH’ comparison and ‘NAU-LB’ vs ‘NAU-RG’ comparison, respectively.

Table 3 Distribution and frequency of SNPs and InDels in functional regions of the radish genome

In addition, all the detected SNPs in the coding region were categorized as synonymous or non-synonymous. A set of 5318 coding SNPs were obtained, of which 2592 (48.74%) were synonymous SNPs, whereas 2726 (51.26%) were nonsynonymous SNPs and caused amino acid changes in protein sequences. The ratio of nonsynonymous vs. synonymous SNPs was 1.05:1. The majority of synonymous SNPs were the leucine (L) (13.35%), while methionine (M) and tryptophan (W) were absent in synonymous changes. Serine (S) (9.21%) and arginine (R) (8%) were the top two types of the nonsynonymous SNPs change. For the nonsynonymous mutation, the alanine (A) (7.3%) was mutated easier than the original amino acid loci (6.09%) (Fig. 6). Moreover, the predominant proteins correlation to chromatin binding, signal transducer activity and motor activity had more variations as their encoding genes had more nonsynonymous SNPs. Among these genes containing synonymous SNPs, the predominant molecular functions of proteins were involved in structural molecular activity, DNA binding, kinase, and transporter activity (Additional file 1: Figure S2).

Fig. 6
figure 6

Amino acid mutations of SNP in coding regions. The vertical axis represents the frequencies of amini acid locus. The horizontal axis show the class of amnio acid identified in this study

Validation of markers and linkage group (LG) locations in radish

To evaluate the accuracy and efficiency of putative polymorphisms, a total of 17 DNA sequences spanning 35 SNPs in ‘NAU-LB’, ‘NAU-YH’ and ‘NAU-RG’ were randomly selected for PCR amplification and Sanger sequencing (Additional file 8). In total, 15 primer pairs could amplify expected sizes, containing a total of 31 SNP loci, of which 26 (83.9%) SNPs revealed identifiable polymorphisms among three genotypes. A total of 180 candidate InDels with lengths greater than or equal to 3 were designed for PCR validation on polyacrylamide gels electrophoresis (PAGE) gels, of which 131 primer pairs showed distinguishable polymorphisms among three genotypes (Fig. 7a-b). However, the other 25 primer pairs showed monomorphism, 14 primer pairs exhibited unexpected size and 10 primer pairs could not successfully be amplified in three genotypes. As for the InDels with ≤2 bp could not be unambiguously discriminated on PAGE gels, a total of 20 InDels with lengths of less than 3 bp were randomly selected for Sanger sequencing (Additional file 1: Table S6), from which 9 loci were polymorphic among three genotypes. Accordingly, 140 of the 200 primer pairs (70%) revealed identifiable polymorphisms among three genotypes.

Fig. 7
figure 7

Genetic diversity analysis of 32 radish accessions with InDel markers. Amplification of 32 genotypes with InDel primers NAU InDel15 (a) and NAU InDel27 (b) by polyacrylamide gel electrophoresis, the UPGMA dendrogram of 32 radish genotypes based on 40 new InDel markers (c). M: 500 bp DNA ladder, the genotype name of numbers (1–32) were listed in Additional file 1: Table S1

To investigate the linkage group (LG) locations of SNP and InDel markers developed in this study, the scaffolds contained SNPs and InDels were separately anchored onto the approximate location of linkage groups. For the validated SNPs and InDels, 15 SNPs and 125 InDels were aligned on radish linkage groups. A total of 140 markers were found to be unevenly distributed on 9 linkage groups (R1-R9) (Additional file 1: Figure S3). Among the 9 linkage groups, R2 contained the largest number of markers (26 markers), while only 5 markers were mapped to R3. Moreover, there were 15 and 9 markers clustered in the bottom of chromosome R2 and R4, respectively.

Application of InDels in genetic diversity analysis of radish germplasm

To assess the value of the InDel markers, 40 InDel marker primers were further used for genetic diversity analysis among 32 worldwide radish accessions. The amplification of 40 polymorphic InDels markers generated 86 alleles across 32 germplasm accessions. The total number of alleles (Na) ranged from 2 to 4 with an average of 2.15. Expected homozygosty (HO) per locus ranged from 0 to 0.53, with an average of 0.24. Expected heterozygosity (HE) per locus ranged from 0.32 to 0.75 with an average of 0.44. The polymorphic information content (PIC) values varied from 0.09 to 0.69 and the mean PIC value was 0.40 (Table 4). The dendrogram showed that the 32 radish accessions could be divided into four major clusters with similarity coefficient varying from 0.55 to 0.81 (Fig. 7c). The cluster I comprised four Chinese accessions, most of which had late maturity characters, and the root shapes were long cylindrical. The cluster II included 22 accessions, and most of the accessions had different geographical origins, taproot colors and maturities. In this cluster, LLYB showed a high similarity with XBY, both had white phloem and xylem with medium maturity. PI358483 originated from Macedonia and PI263262 collected from Japan were divided into a subgroup, and both of them had similar taproot color and maturity. The cluster III contained four wild accessions, the RKZ was collected from China, while the PI271451, PI436536 and X50 with different taproot colors were originated from India, Guatemala and Germany, respectively. The cluster IV consisted of one Russian and Turkey accession with medium maturity.

Table 4 Genetic diversity analysis with InDel markers developed in radish


Alternative splicing emerged as a key scheme in multicellular eukaryotes to enhance the proteome diversity [38]. SNP and InDel markers played important roles in molecular breeding program, which had been widely utilized in high-density genetic map construction, genome-wide association mapping and marker-assisted selection [20, 21]. Although several molecular markers had been successfully developed in radish [28, 29], the number of SNP and InDel markers were still limited for genetic diversity analysis and genetic mapping. To better systematically characterize the overall transcriptional landscape in radish, the identification and development of large-scale AS, SNP and InDel markers is needed for future genetic studies. In this study, a total of 177,540 AS events, 22,412 SNPs, 9436 InDels and 43,680 SSRs were successfully detected from three transcriptomes sequencing data of radish. Moreover, 40 InDel markers were further applied for genetic diversity analysis among 32 different radish accessions. To the best of our knowledge, this is the first report on identification of AS events, SNPs and InDels using comparative transcriptomic analysis in radish.

Alternative splicing plays a crucial role in transcriptional regulation

RNA sequencing had been widely used for identification of alternative splicing levels in plants [23, 39]. In this study, a genome-wide identification and characterization of alternative splicing in radish was conducted by three transcriptome datasets. It was found that 59.4% of the total radish genes were alternatively spliced, and the frequency was corroborated with previous findings in Arabidopsis [24] and tomato [26]. It was commonly accepted that IR is the major type of AS in plants [23, 25]. IR is the most frequent AS type in radish, while ES is the least type (2.4% of the total events). In Arabidopsis, the fewer number of ES was caused by multiple exons skipped together and exon skipping utilized alternative 5′ and/or 3′ splice sites [4]. The UDP-galactose 4-epimerase as a member of the short chain dehydrogenase superfamily were underwent splicing variants in Aspergillus flavus [40]. In this study, Rsa1.0_16819.1_g00001.1, encoding the UDP-D-galactose 4-epimerase 1 also showed alternative splicing. Rsa1.0_01846.1_g00001.1, encoding the 2-oxoglutarate (2OG) and Fe (II)-dependent oxygenase superfamily protein showed tissue-specific expression, which was coincided with previously studies that a few genes exhibited tissue-specific splicing variants in plants [23, 26]. Recent studies revealed that different transcripts produced from individual multi-exon genes could accelerate genome evolution through generating new functions [3]. Three BrFLC1 alleles with alternative splicing patterns greatly contributed to flowering-time variations in B. rapa [41], indicating that the AS has important biological consequences in plants. Together, the identification of AS patterns in genes would facilitate understanding the transcriptional complexity in root vegetable crops.

Features of SNPs and InDels by comparative transcriptomic analysis

It has increasingly been proven that SNPs and InDels can potentially be linked to functional genes involved in important agronomic traits [8, 17]. Comparative transcriptomic analysis has been widely applied for identification of SNPs and InDels in many plants, such as B. napus [9], tomato [11] and Petunia [42]. In this study, the SNP and InDel density was 1 SNP/17.9 kb and 1 InDel/42.6 kb, respectively, which was lower than previously described average occurrence of 3.9/kb by resequencing 93 radish F2 individuals [29]. Generally, genomic DNA sequences had higher SNP frequency than those observed in the transcribed regions. Additionally, the SNPs and InDels within the coding regions are more likely to contribute to agronomic phenotypes. Among pairwise comparisons, the largest numbers of SNPs and InDels were present in the ‘NAU-RG’ vs ‘NAU-YH’ comparison, which suggested that ‘NAU-RG’ vs ‘NAU-YH’ had a far relationship than the other two comparisons. Moreover, the transition-transversion ratio (1.41:1) presented in this study was similar to the previous findings in melon [43] and Petunia [42]. The Ts/Tv bias could be explained by the fact that cytosine-guanine (CpG) dinucleotides exhibited the high transition frequencies after methylation [44]. Notably, a total of 13 InDel types were found in these three genotypes, and the two most common types were single- and bi-nucleotide InDels, which have also been reported in several plant species including B. rapa, cotton and sesame [45,46,47].

In the current study, the identified SNPs/InDels heterozygotes (81.5%/80.5%) were predominate as compared with homozygotes (18.5%/19.5%), confirming that the high heterozygosity rate was found in three accessions. Meanwhile, the clean reads uniquely mapped to the radish reference genome were less than 80% in each transcriptome, which was similar to previously reported results that heterozygous polymorphisms had been supported by a low proportion of uniquely mapped reads [48]. Previous studies have identified a number of heterozygous single-nucleotide variants (SNVs) in Puccinia striiformis f. sp. Tritici (Pst)-CY32 genome, and some of the predicted genes were located in the heterozygous regions [49]. Wang et al. [50] revealed that heterozygous alleles played significant roles in restoring male fertility of cytoplasmic male-sterile in radish. Therefore, the results would provide useful information for the research of genome heterozygosity and its molecular function in radish.

Based on the genome position of SNPs and InDels, the sequence region types were classified. The SNPs frequency was enriched almost equally between genic (49.9%) and intergenic (50.1%) regions. Approximately 85% of the total InDels was identified in genic regions, which was similar to previously reported in rice and sesame [47, 51]. It has been extensively shown that SNPs/InDels in the coding region were most influential variations, and mainly lead to non-function of protein [48, 52]. Notably, 73% of all InDels were identified in the CDS regions, which could cause frameshift mutation that results in the changes of gene function. Numerous studies demonstrated that non-synonymous played significant roles in amino acid variation in the protein product of genes [53]. The total numbers of nonsynonymous SNPs and synonymous SNPs were 2592 and 2726, the nonsynonymous to synonymous ratio of 1.05 exhibited a medium value compared with Arabidopsis (0.83) [54] and rice (1.29) [55]. Together, these results illustrated the identified SNPs and InDels would be valuable molecular tool for identification of causal gene mutations in radish.

Characterization and polymorphism of SNP and InDel markers

In this study, a total of 35 SNPs and 200 InDels identified from three transcriptome sequences of the ‘NAU-RG’, ‘NAU-LB’ and ‘NAU-YH’ were employed for designing PCR primer. The polymorphism rate of SNP/InDel (83.9%/70%) among these three radish genotypes was in accordance with previously reported results [42, 45]. Obviously, the polymorphism ratios for transcriptome derived SNPs was 84.3% by selecting a 51 SNPs set in Petunia spp. [42]. A large-scale ESTs sequence data from different cotton species were used to identify putative InDels, of which 71.1% InDels were showed loci polymorphisms [46]. Several identified markers amplified with unexpected bands and failed amplification, which might be related to the presence of large introns in genomic sequence [56]. Most of the InDel (87.5%, 35/40) markers produced two allele products, which is similar to the previous results in tomato [8]. Moreover, the average PIC value of these 40 InDels markers was 0.40, which is relatively higher than the PIC value calculated in sesame [47] but lower than that in trifoliate orange [48], indicating that the moderate level of polymorphism of the developed markers would greatly facilitate genetic diversity analysis in radish.

Availability of InDels in radish genetic diversity analysis

InDels are important tools for a wide range of genetics and genomic studies. In this study, phylogenetic analysis revealed that 32 radish accessions were grouped into four distinct clusters, thereby demonstrating the availability of the InDel markers to differentiate the cultivars. Besides, the distribution of radish cultivars was not strictly base on the geographical origins, which is consistent with previous studies in genus Oryza [57]. It is generally accepted that InDel markers display higher interspecies differentiation than SSR markers [58]. Currently, the efficient differentiation of the cultivars and wild radish by InDel markers in third and fourth main cluster was higher than previous reported SSR marker in 32 radish accessions [59]. However, most of Chinese cultivars were clustered in the same subgroup in the dendrogram, indicating that the genetic relationship of Chinese radish cultivars is still narrow. Therefore, it is important to introduce more exotic genotype for expanding the genetic basis of radish.


In summary, this is the first report on transcriptome-based characterization of AS events, SNP and InDel markers in radish. Totally, 56,530 AS events were accounted for 59.4% of the total radish genes among these three transcriptomes. Moreover, 22,412 SNPs, 9436 InDels and 43,680 SSRs were found in all the three genotypes. Among the total SNPs located in coding regions, the percentage of nonsynonymous and synonymous SNPs were 48.74% and 51.26%, respectively. The polymorphism rate of SNPs and InDels was verified among three radish accessions, and 83.9% of SNPs and 70% of InDels exhibited polymorphism. Totally, 15 SNPs and 125 InDels were found to be unevenly distributed on radish nine linkage groups. Germplasms were successfully differentiated by several newly developed InDel markers. Taken together, our results could provide a better understanding of the radish transcriptome complexity, and these newly developed molecular markers would be excellent tool for construction of linkage map, identification of candidate genes for interest traits, and marker-assisted selection in radish breeding programs.


Plant materials and DNA extraction

Seeds of radish advanced inbred line, ‘NAU-RG’, ‘NAU-LB’, ‘NAU-YH’ and other 32 radish accessions with different root colors and origins (Additional file 1: Table S1) were soaked, surface-sterilized and incubated for 3 days. The germinated seeds were grown in plastic pots and cultured in greenhouse under 25 °C/14 h light and 18 °C/10 h dark. Young leaves of 30-day-old seedlings were collected for DNA extraction. Genomic DNA of all accessions were extracted from young leaves using a modified CTAB procedure [60]. Then DNA samples were diluted to a final concentration of 10 ng μl−1 with 1 × TE buffer for further use.

Transcriptome sequences and mapping reads to reference genome

Three transcriptome data, sequencing from roots of genotypes ‘NAU-YH’ and ‘NAU-RG’ (NCBI accession ID No. SRX707630 and SRX316199) [31, 32], leaves of ‘NAU-LB’ (NCBI accession ID No. SRX1671013) [33], were collected from our previously reports. The radish reference genome sequence and gene model annotation files were downloaded from the available genome website ( Index of the reference genome was built using Bowtie v2.2.3, and paired-end clean reads mapped to the radish reference genome using TopHat v2.0.12 [34]. TopHat can generate a database of splice junctions based on the gene model annotation file.

Identification and validation of AS events

Splice junctions were identified using TopHat program, and then filtered the junctions with reads number less than three. The known and novel transcripts were identified by Cufflinks v2.1.1 Reference Annotation Based Transcript (RABT) assembly method [61]. The final transcript annotations in a GTF format file were used to identify AS events using ASTALAVISTA program with default settings [35]. Five types including IR, A3SS, A5SS, ES and others that contain more than one of the four basic types were analyzed as previously described [4]. The mapping results of AS events were visualized using the Integrative Genomics Viewer (IGV) [62]. Alternative splicing isoforms were quantified between the three genotypes using the MATS program [36]. The differentially spliced events were defined with at least 10% change in exon inclusion level and false discovery rate (FDR) of less than 0.05. The splicing events were represented with the sashimi plots [63].

Total RNA was isolated from leaves and roots of ‘NAU-RG’, ‘NAU-LB’, and ‘NAU-YH’ with Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The RNAs were treated with RNase-free DNase I (Takara, Japan), and reverse-transcribed into cDNAs using SuperScript II reverse transcriptase (Invitrogen). The primer information is listed in Additional file 1: Table S3. RsActin was used as the internal control, and PCR products were visualized in agarose gel stained by ethidium bromide (EB).

Identification of SNPs, InDels and SSRs

The Picard-tools v1.96 and samtools v0.1.18 programs were used to sort and mark duplicated reads as well as reorder the alignment results of each sample [64]. SNP and InDel calling were carried out using the Genome Analysis Toolkit (GATK, version v3.2) [37], and the variant were filtered as follows: low quality scores (QUAL) is less than 30 and quality by depth (QD) is less than 5. Only the putative SNPs or InDels with 10 or more reads coverage were recorded to reduce the numbers of false positive SNPs. The localization of SNPs and InDels was performed based on the annotation of gene models of the radish genome. Synonymous and non-synonymous substitutions of the SNPs were determined using the SNPEff program [65]. Blast2GO program was applied to achieve Gene Ontology analysis on the genes contained SNPs and InDels.

To generate the reference transcriptome sequences, a de novo assembly of three transcriptome data was performed using Trinity program [66]. The potential SSRs were identified from the assembled unigenes using the MISA (microsatellite identification tool) program [67]. The number of repetitive units with a minimum of 10, 6, 5, 5, 5 and 5 were defined for the mono- to hexa-nucleotide, respectively.

Primer design and validation of DNA polymorphism

Primer pairs flanking the SNPs or InDels were designed using Primer Premier 5.0 program. The length of primer pairs ranged from 18 to 24 bp, while the predicted product sizes were varied from 80 bp to 400 bp. The Tm value was restricted to between 50 and 60 °C. PCR was carried out in a 15 μl reaction volume containing 10 ng of template DNA, 2.0 mM of MgCl2, 0.2 mM dNTPs, 0.1 μM of each primer and 0.5 U Taq DNA polymerase (TaKaRa Bio Inc., Dalian, China). The PCR procedure comprised an initial denaturation at 94 °C for 3 min, followed by 35 cycles of 94 °C for 50 s, 56 °C for 50 s and 72 °C for 1 min, with a final extension at 72 °C for 10 min. The SNP amplified PCR products were sequenced by Sanger method and analyzed by BioEdit v ( bioedit.html). The InDel amplified products were separated on 8% PAGE and visualized with silver staining as described previously [68]. The validated SNPs and InDels were mapped on the prior dense genetic maps of radish according to their approximate positions on the linkage group locations using MapInspect software ( [30].

Data analysis

To estimate the allelic variation of InDel markers in 32 radish accessions, POPgene program (v1.32) was employed to estimate the total number of alleles, expected homozygosty and expected heterozygosity [48]. The polymorphic information content value for the InDel marker was calculated using Power Marker v3.0 [69]. The 0 – 1 data matrix was further used to calculate the coefficients of genetic similarity among all the accessions using the SIMQUAL program of NTSYS-pc software. The dendrogram was constructed basing on the unweighted pair-group method with arithmetic averages (UPGMA) in the SAHN module of NTSYS-pc software [59].



Alternative 3′ splice site donor


Alternative 5′ splice site donor


Alternative splicing


Coding sequence


Exon skipping


Expected heterozygosity


Expected homozygosity




Intron retention


Linkage group


Number of alleles


National center for biotechnology information


Polyacrylamide gel electrophoresis


Polymorphic information content


Single nucleotide polymorphism


Simple sequence repeat






Unweighted pair-group method with arithmetic averages


  1. Lu T, Lu G, Fan D, Zhu C, Li W, Zhao Q, et al. Function annotation of the rice transcriptome at single-nucleotide resolution by RNA-seq. Genome Res. 2010;20(9):1238–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. van Dijk EL, Auger H, Jaszczyszyn Y, Thermes C. Ten years of next-generation sequencing technology. Trends Genet. 2014;30(9):418–26.

    Article  PubMed  Google Scholar 

  3. Huang J, Gao YJ, Jia HT, Liu L, Zhang D, Zhang ZX. Comparative transcriptomics uncovers alternative splicing changes and signatures of selection from maize improvement. BMC Genomics. 2015;16(1):1.

    Article  Google Scholar 

  4. Marquez Y, Brown JW, Simpson C, Barta A, Kalyna M. Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res. 2012;22(6):1184–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011;27(17):2325–9.

    Article  CAS  PubMed  Google Scholar 

  6. Maher CA, Kumar-Sinha C, Cao X, Kalyana-Sundaram S, Han B, Jing X, et al. Transcriptome sequencing to detect gene fusions in cancer. Nature. 2009;458(7234):97–101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Barbazuk WB, Emrich SJ, Chen HD, Li L, Schnable PS. SNP discovery via 454 transcriptome sequencing. Plant J. 2007;51(5):910–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Yang J, Wang Y, Shen H, Yang W. In Silico identification and experimental validation of Insertion–Deletion polymorphisms in tomato genome. DNA Res. 2014;21(4):429–38.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Trick M, Long Y, Meng J, Bancroft I. Single nucleotide polymorphism (SNP) discovery in the polyploid Brassica napus using Solexa transcriptome sequencing. Plant Biotechnol J. 2009;7(4):334–46.

    Article  CAS  PubMed  Google Scholar 

  10. Varshney RK, Graner A, Sorrells ME. Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005;23(1):48–55.

    Article  CAS  PubMed  Google Scholar 

  11. Koenig D, Jiménez-Gómez JM, Kimura S, Fulop D, Chitwood DH, Headland LR, et al. Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato. Proc Natl Acad Sci U S A. 2013;110(28):E2655–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Wu G, Zhang LB, Yin YT, Wu JS, Yu LJ, Zhou YH, et al. Sequencing, de novo assembly and comparative analysis of Raphanus sativus transcriptome. Front Plant Sci. 2015;6:198.

    PubMed  PubMed Central  Google Scholar 

  13. Kalia RK, Rai MK, Kalia S, Singh R, Dhawan AK. Microsatellite markers: an overview of the recent progress in plants. Euphytica. 2011;177(3):309–34.

    Article  CAS  Google Scholar 

  14. Xu L, Wang LJ, Gong YQ, Dai WH, Wang Y, Zhu XW, et al. Genetic linkage map construction and QTL mapping of cadmium accumulation in radish (Raphanus sativus L.). Theor Appl Genet. 2012;125(4):659–70.

    Article  CAS  PubMed  Google Scholar 

  15. Yang H, Li C, Lam HM, Clements J, Yan G, Zhao S. Sequencing consolidates molecular markers with plant breeding practice. Theor Appl Genet. 2015;128(5):779–95.

    Article  CAS  PubMed  Google Scholar 

  16. Zou X, Shi C, Austin RS, Merico D, Munholland S, Marsolais F, et al. Genome-wide single nucleotide polymorphism and Insertion-Deletion discovery through next-generation sequencing of reduced representation libraries in common bean. Mol Breed. 2014;33(4):769–7.

    Article  CAS  Google Scholar 

  17. Thakur S, Singh PK, Rathour R, Variar M, Prashanthi SK, Gopalakrishnan S, et al. Genotyping and development of single-nucleotide polymorphism (SNP) markers associated with blast resistance genes in rice using GoldenGate assay. Mol Breed. 2014;34(3):1449–63.

    Article  CAS  Google Scholar 

  18. Salem M, Vallejo RL, Leeds TD, Palti Y, Liu S, Sabbagh A, et al. RNA-Seq identifies SNP markers for growth traits in rainbow trout. PLoS One. 2012;7(5):e36264.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. de Souza LM, Toledo-Silva G, Cardoso-Silva CB, Da Silva CC, de Araujo Andreotti IA, Conson ARO, et al. Development of single nucleotide polymorphism markers in the large and complex rubber tree genome using next-generation sequence data. Mol Breed. 2016;36(8):115.

    Article  Google Scholar 

  20. Pfender W, Saha M, Johnson E, Slabaugh M. Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. Theor Appl Genet. 2011;122(8):1467–80.

    Article  CAS  PubMed  Google Scholar 

  21. Chen W, Gao Y, Xie W, Gong L, Lu K, Wang W, et al. Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism. Nat Genet. 2014;46(7):714–21.

    Article  CAS  PubMed  Google Scholar 

  22. Wang BB, Brendel V. Genome wide comparative analysis of alternative splicing in plants. Proc Natl Acad Sci U S A. 2006;103(18):7175–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zhang G, Guo G, Hu X, Zhang Y, Li Q, Li R, et al. Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome. Genome Res. 2010;20(5):646–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Filichkin SA, Priest HD, Givan SA, Shen R, Bryant DW, Fox SE, et al. Genome-wide mapping of alternative splicing in Arabidopsis thaliana. Genome Res. 2010;20(1):45–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Shen Y, Zhou Z, Wang Z, Li W, Fang C, Wu M, et al. Global dissection of alternative splicing in paleopolyploid soybean. Plant Cell. 2014;26(3):996–1008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Sun Y, Xiao H. Identification of alternative splicing events by RNA sequencing in early growth tomato fruits. BMC Genomics. 2015;16(1):1.

    Article  Google Scholar 

  27. Mandadi KK, Scholthof K-BG. Genome-wide analysis of alternative splicing landscapes modulated during plant-virus interactions in Brachypodium distachyon. Plant Cell. 2015;27(1):71–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Shen D, Sun H, Huang M, Zheng Y, Li X, Fei Z. RadishBase: a database for genomics and genetics of radish. Plant Cell Physiol. 2013;54(2):e3.

    Article  CAS  PubMed  Google Scholar 

  29. Mun JH, Chung H, Chung WH, Oh M, Jeong YM, Kim N, et al. Construction of a reference genetic map of Raphanus sativus based on genotyping by whole-genome resequencing. Theor Appl Genet. 2015;128(2):259–72.

    Article  CAS  PubMed  Google Scholar 

  30. Kitashiba H, Li F, Hirakawa H, Kawanabe T, Zou Z, Hasegawa Y, et al. Draft sequences of the radish (Raphanus sativus L.) genome. DNA Res. 2014;21(5):481–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wang Y, Pan Y, Liu Z, Zhu XW, Zhai LL, Xu L, et al. De novo transcriptome sequencing of radish (Raphanus sativus L.) and analysis of major genes involved in glucosinolate metabolism. BMC Genomics. 2013;14(1):1.

    Article  Google Scholar 

  32. Yu RG, Xu L, Zhang W, Wang Y, Luo XB, Wang RH, et al. De novo taproot transcriptome sequencing and analysis of major genes involved in sucrose metabolism in radish (Raphanus sativus L.). Front Plant Sci. 2016;7:585.

    PubMed  PubMed Central  Google Scholar 

  33. Nie SS, Li C, Xu L, Wang Y, Huang DQ, Muleke EM, et al. De novo transcriptome analysis in radish (Raphanus sativus L.) and identification of critical genes involved in bolting and flowering. BMC Genomics. 2016;17(1):1.

    Google Scholar 

  34. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Foissac S, Sammeth M. ASTALAVISTA: Dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007;35:W297–9.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Shen SH, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, et al. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A. 2014;111(51):E5593–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Vitulo N, Forcato C, Carpinelli EC, Telatin A, Campagna D, D'Angelo M, et al. A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype. BMC Plant Biol. 2014;14(1):1.

    Article  Google Scholar 

  39. Bao H, Li E, Mansfield SD, Cronk QC, El-Kassaby YA, Douglas CJ. The developing xylem transcriptome and genome-wide analysis of alternative splicing in Populus trichocarpa (black cottonwood) populations. BMC Genomics. 2013;14(1):1.

    Article  Google Scholar 

  40. Chang KY, Muddiman DC. Identification of alternative splice variants in Aspergillus flavus through comparison of multiple tandem MS search algorithms. BMC Genomics. 2011;12(1):358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Yuan YX, Wu J, Sun RF, Zhang XW, Xu DH, Bonnema G, et al. A naturally occurring splicing site mutation in the Brassica rapa FLC1 gene is associated with variation in flowering time. J Exp Bot. 2009;60(4):1299–308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Guo Y, Wiegert-Rininger KE, Vallejo VA, Barry CS, Warner RM. Transcriptome-enabled marker discovery and mapping of plastochron-related genes in Petunia spp. BMC Genomics. 2015;16(1):1.

    Article  Google Scholar 

  43. Blanca JM, Cañizares J, Ziarsolo P, Esteras C, Mir G, Nuez F, et al. Melon transcriptome characterization: Simple sequence repeats and single nucleotide polymorphisms discovery for high throughput genotyping across the species. Plant Genome. 2011;4(2):118–31.

    Article  CAS  Google Scholar 

  44. Keller I, Bensasson D, Nichols RA. Transition-transversion bias is not universal: a counter example from grasshopper pseudogenes. PLoS Genet. 2007;3(2):e22.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Liu B, Wang Y, Zhai W, Deng J, Wang H, Cui Y, et al. Development of InDel markers for Brassica rapa based on whole-genome re-sequencing. Theor Appl Genet. 2013;126(1):231–9.

    Article  CAS  PubMed  Google Scholar 

  46. Lv YD, Zhao L, Xu XY, Wang L, Wang C, Zhang TZ, et al. Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species. BMC Genomics. 2013;14(1):1.

    Article  Google Scholar 

  47. Wei L, Miao H, Li C, Duan Y, Niu J, Zhang T, et al. Development of SNP and InDel markers via de novo transcriptome assembly in Sesamum indicum L. Mol Breed. 2014;34(4):2205–17.

    Article  CAS  Google Scholar 

  48. Zhang JZ, Liu SR, Hu CG. Identifying the genome-wide genetic variation between precocious trifoliate orange and its wild type and developing new markers for genetics research. DNA Res. 2016;23(4):403–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Zheng W, Huang L, Huang J, Wang X, Chen X, Zhao J, et al. High genome heterozygosity and endemic genetic recombination in the wheat stripe rust fungus. Nat Commun. 2013;4:2676.

    Google Scholar 

  50. Wang ZW, Wang C, Gao L, Mei SY, Zhou Y, Xiang CP, et al. Heterozygous alleles restore male fertility to cytoplasmic male-sterile radish (Raphanus sativus L.): a case of overdominance. J Exp Bot. 2013;64(7):2041–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Arai-Kichise Y, Shiwa Y, Nagasaki H, Ebana K, Yoshikawa H, Yano M, et al. Discovery of genome-wide DNA polymorphisms in a landrace cultivar of japonica rice by whole-genome sequencing. Plant Cell Physiol. 2011;52(2):274–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Sathya B, Dharshini AP, Kumar GR. NGS meta data analysis for identification of SNP and INDEL patterns in human airway transcriptome: A preliminary indicator for lung cancer. Appl Transl Genom. 2015;4:4–9.

    Article  Google Scholar 

  53. Ramensky V, Bork P, Sunyaev S. Human non-synonymous SNPs: server and survey. Nucleic Acids Res. 2002;30(17):3894–900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Clark RM, Schweikert G, Toomajian C, Ossowski S, Zeller G, Shinn P, et al. Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana. Science. 2007;317(5836):338–42.

    Article  CAS  PubMed  Google Scholar 

  55. Xu X, Liu X, Ge S, Jensen JD, Hu F, Li X, et al. Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes. Nat Biotechnol. 2012;30(1):105–11.

    Article  CAS  Google Scholar 

  56. Zeng SH, Xiao G, Guo J, Fei ZJ, Xu YQ, Roe BA, et al. Development of a EST dataset and characterization of EST-SSRs in a traditional Chinese medicinal plant, Epimedium sagittatum (Sieb. Et Zucc.) Maxim. BMC Genomics. 2010;11(1):1.

    Article  Google Scholar 

  57. Yamaki S, Ohyanagi H, Yamasaki M, Eiguchi M, Miyabayashi T, Kubo T, et al. Development of INDEL markers to discriminate all genome types rapidly in the genus Oryza. Breed Sci. 2013;63(3):246–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. García-Lor A, Luro F, Navarro L, Ollitrault P. Comparative use of InDel and SSR markers in deciphering the interspecific structure of cultivated citrus genetic diversity: a perspective for genetic association studies. Mol Gen Genomics. 2012;287(1):77–94.

    Article  Google Scholar 

  59. Zhai LL, Xu L, Wang Y, Cheng H, Chen YL, Gong YQ, et al. Novel and useful genic-SSR markers from de novo transcriptome sequencing of radish (Raphanus sativus L.). Mol Breed. 2014;33(3):611–24.

    Article  CAS  Google Scholar 

  60. Liu L, Guo W, Zhu X, Zhang T. Inheritance and fine mapping of fertility restoration for cytoplasmic male sterility in Gossypium hirsutum L. Theor Appl Genet. 2003;106(3):461–9.

    Article  CAS  PubMed  Google Scholar 

  61. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Katz Y, Wang ET, Airoldi EM, Burge CB. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010;7:1009–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Cingolani P, Platts A, Wan LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.

    Article  CAS  Google Scholar 

  66. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.

    Article  CAS  PubMed  Google Scholar 

  67. Thiel T, Michalek W, Varshney RK, Graner A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003;106:411–22.

    Article  CAS  PubMed  Google Scholar 

  68. Liu LW, Zhao LP, Gong YQ, Wang MX, Chen LM, Yang JL, et al. DNA fingerprinting and genetic diversity analysis of late-bolting radish cultivars with RAPD. ISSR and SRAP markers Sci Hort. 2008;116(3):240–7.

    Article  CAS  Google Scholar 

  69. Liu K, Muse SV. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005;21(9):2128–9.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


This work was in part supported by grants from the Natural Science Foundation of China (31,372,064, 31,501,759), National Key Technology R&D Program of China (2016YFD0100204; 2017YFD0101803), Jiangsu Agricultural Science and Technology Innovation Fund (JASTIF, CX(16)1012] and Key Technology R&D Program of Jiangsu Province (BE2016379, BE2013429). These funding sources had no role in study design and data collection, analysis, interpretation of data and in writing the manuscript.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

Author information

Authors and Affiliations



XL performed the data analysis and wrote the manuscript. MT, DL and HJ conducted validation of experiments. YW, DL and WZ contributed powerful analytical tools. XZ and YZ helped with the revise of manuscript. LL and LX conceived and designed the research. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Liwang Liu.

Ethics declarations

Ethics approval and consent to participate

In this study, all the radish lines were collected from germplasm bank of National Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University. Collection of plant materials complied with the institutional, national and international guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests for this research.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Radish genotypes used for genetic diversity analysis in this study. Table S2. Numbers of identified splice junction reads in three transcriptomes. Table S3. RT-PCR experimental validation of AS events. Table S4. Distribution and frequency of InDels in radish. Table S5. Distribution of mono- to tetranucleotide repeats in radish ranscriptome. Table S6. The InDels validated by Sanger sequencing among three genotypes. Figure S1. Distribution of AS events in radish genome. Figure S2. Synonymous and non-synonymous SNP distribution in 16 major protein molecular functions. Figure S3. Linkage group (LG) localization of SNP and InDel markers in the radish genome. The bar on the left shows the marker positions [cM], marker names are shown on the left of each linkage group. (DOCX 1013 kb)

Additional file 2:

The information of detected alternative splicing (AS) events in radish. (XLSX 5752 kb)

Additional file 3:

The differential alternative splicing events between three genotypes (XLSX 740 kb)

Additional file 4:

List of novel transcripts detected in the genotypes ‘NAU-YH’, ‘NAU-LB’ and ‘NAU-RG’. (XLSX 523 kb)

Additional file 5:

Putative single nucleotide polymorphisms (SNPs) in the genotypes ‘NAU-YH’, ‘NAU-LB’ and ‘NAU-RG’. (XLSX 1184 kb)

Additional file 6:

Putative insertion/deletion (InDel) in the genotypes ‘NAU-YH’, ‘NAU-LB’ and ‘NAU-RG’. (XLSX 517 kb)

Additional file 7:

The detailed information of putative SSR markers identified in radish. (XLSX 2112 kb)

Additional file 8:

Primer sequences used for SNPs validation between ‘NAU-YH’, ‘NAU-LB’ and ‘NAU-RG’ by Sanger sequencing. (XLSX 54 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Luo, X., Xu, L., Liang, D. et al. Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.). BMC Genomics 18, 505 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: