Transcriptome sequencing of two parental lines of cabbage (Brassica oleracea L. var. capitata L.) and construction of an EST-based genetic map
- Nur Kholilatul Izzah†1, 4,
- Jonghoon Lee†1,
- Murukarthick Jayakodi1,
- Sampath Perumal1,
- Mina Jin2,
- Beom-Seok Park2,
- Kyounggu Ahn3 and
- Tae-Jin Yang1Email author
© Izzah et al.; licensee BioMed Central Ltd. 2014
Received: 28 October 2013
Accepted: 17 February 2014
Published: 22 February 2014
Expressed sequence tag (EST)-based markers are preferred because they reflect transcribed portions of the genome. We report the development of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers derived from transcriptome sequences in cabbage, and their utility for map construction.
Transcriptome sequences were obtained from two cabbage parental lines, C1184 and C1234, which are susceptible and resistant to black rot disease, respectively, using the 454 platform. A total of 92,255 and 127,522 reads were generated and clustered into 34,688 and 40,947 unigenes, respectively. We identified 2,405 SSR motifs from the unigenes of the black rot-resistant parent C1234. Trinucleotide motifs were the most abundant (66.15%) among the repeat motifs. In addition, 1,167 SNPs were detected between the two parental lines. A total of 937 EST-based SSR and 97 SNP-based dCAPS markers were designed and used for detection of polymorphism between parents. Using an F2 population, we built a genetic map comprising 265 loci, and consisting of 98 EST-based SSRs, 21 SNP-based dCAPS, 55 IBP markers derived from B. rapa genome sequence and 91 public SSRs, distributed on nine linkage groups spanning a total of 1,331.88 cM with an average distance of 5.03 cM between adjacent loci. The parental lines used in this study are elite breeding lines with little genetic diversity; therefore, the markers that mapped in our genetic map will have broad spectrum utility.
This genetic map provides additional genetic information to the existing B. oleracea map. Moreover, the new set of EST-based SSR and dCAPS markers developed herein is a valuable resource for genetic studies and will facilitate cabbage breeding. Additionally, this study demonstrates the usefulness of NGS transcriptomes for the development of genetic maps even with little genetic diversity in the mapping population.
The genus Brassica includes some of the most economically important crops with wide-ranging adaptability for cultivation under various agro-climatic conditions. Among these are six widely cultivated species, including three monogenomic diploids: B. rapa (AA, 2n = 20), B. nigra (BB, 2n = 16), B. oleracea (CC, 2n = 18); and three amphidiploids (allotetraploids): B. juncea (AABB, 2n = 36), B. napus (AACC, 2n = 38) and B. carinata (BBCC, 2 n = 34), which evolved through hybridization between different monogenomic diploids as described by U’s triangle . Among those cultivated species, B. oleracea exhibits the largest genetic and morphological diversity , encompassing many common vegetables such as cabbage, broccoli, cauliflower, Brussels sprout, kale, kohlrabi, and kai-lan.
In recent years, development of massive sequencing technology has led to new possibilities for high-throughput genome analysis . Transcriptome sequencing has become an attractive approach compared to whole-genome sequencing because it allows efforts to be targeted to genic regions . The advent of transcriptome sequencing has yielded a huge amount of transcribed sequence data, such as expressed sequence tags (ESTs) that can be exploited for gene expression profiling, genome annotation, comparative genomics and physical mapping. ESTs offer a simple strategy to study the transcribed portions of genomes and provide a robust sequence resource from which to develop functional markers [5, 6].
Among next-generation sequencing (NGS) technologies, the 454 sequencing platform has been successfully employed for de novo transcriptome sequencing of many plant species, including American ginseng (Panax quinquefolius) , A. thaliana, maize (Zea mays)  and olive (Olea europaea) . A single-plate run on the 454 GS-FLX titanium platform typically produces around million reads with an average length of 400 bp, and is faster and cheaper than traditional Sanger sequencing methods .
Molecular DNA markers have greatly contributed to the development of plant genetics and breeding studies. The use of DNA markers has become essential for crop improvement programs, such as for cultivar identification, genetic diversity, linkage map construction and identification of quantitative trait loci (QTL) . Among the many types of DNA markers, simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) are the preferred marker types for many genetic applications. SSRs are efficient co-dominant anchor markers with high levels of polymorphism and can easily be amplified by polymerase chain reaction (PCR) using primers designed from flanking sequences of the SSR motifs. Meanwhile, SNPs are abundant in virtually all populations, with the majority being biallelic, and can be tightly linked to or are the actual cause of allelic (phenotypic) differences in traits .
A total of 19 different genetic linkage maps of B. oleracea have been published so far . However, most of those maps were developed using RFLP markers, which could not be easily transferred to other genetic maps. In addition, only a few functional markers have been mapped in B. oleracea. The latest cabbage genetic map constructed by Wang et al.  contained 80 EST-SSR markers. However, publicly available EST-based SSR and EST-based dCAPS markers are still limited in this plant species. Hence, development of more functional markers is desirable in order to facilitate the mapping, tagging and identification of important trait loci.
To generate large-scale EST data and to develop functional markers using ESTs, we performed large-scale transcriptome sequencing of two cabbage parental lines, C1184 and C1234, which have been used as elite breeding lines for development of commercial F1 varieties by the Joeun Seed company in Korea, using the 454 sequencing platform, Roche GS FLX Titanium series. These newly developed EST-based SSR and dCAPS markers were mapped along with SSR markers and intron-based polymorphism (IBP) markers that were previously reported for Brassica species. This genetic map will promote QTL mapping and breeding.
Sequence assembly and functional annotation
Summary of 454 transcriptome sequencing and assembly data
Response to black rot disease
No. of HQ reads
Total number of assembled contigs
No. of singletons
No. of unigenes
Largest contig length (bp)
Average contig length (bp)
N50 length (bp)
Prior to functional annotation, the singletons from both lines that were less than 200 bp in length were removed. The remaining 33,244 and 38,088 assembled unique transcripts found in cabbage lines C1184 and C1234, respectively, were compared against the NCBI non-redundant (nr) protein database using the BLASTX algorithm. Of these, 27,740 (79%) of the C1184 unigenes and 31,458 (76%) of the C1234 unigenes had significant hits. Among the transcripts with hits, more than 76% matched known functional genes in both lines.
NBS-encoding genes in the black rot-resistant line
Most disease resistance R genes encode putative proteins containing nucleotide binding sites (NBS) and leucine-rich repeat (LRR) domains. NBS-LRR genes can be divided into the Toll-interleukin-1 receptor (TIR-NBS-LRR) and coiled-coil (CC-NBS-LRR) sub-families. Plant NBS-LRR-encoding genes play an important role in the responses of plants to various pathogens, including viruses, bacteria and fungi. Hence, the black rot-resistant line C1234 was searched for NBS-encoding genes. A total of 29 ESTs related to the NBS-LRR domain were identified in the C1234 line (Additional file 1). Among them, 22 were classified as TIR-NBS-LRRs and 7 were categorized as CC-NBS-LLRs. Thus, we conclude that many of the NBS genes in B. oleracea belong to the TIR-NBS-LRR type. Further, we compared the abundance/expression profile of NBS genes between C1234 and C1184 lines based on the mapping read count. As expected, we found a higher number of mapping counts for NBS genes in the resistant line (C1234) when compared to the susceptible (C1184) line (Additional file 1).
Frequency and distribution of SSRs in cabbage ESTs
Summary of EST-based SSRs identified
Total EST sequences examined
- Dinucleotide motif
- Trinucleotide motif
- Tetranucleotide motif
- Pentanucleotide motif
- Hexanucleotide motif
ESTs containing SSRs (SSR-ESTs)
- ESTs containing 1 SSR
- ESTs containing more than 2 SSRs
- ESTs with SSRs in compound formation
SSR primer design
- ESTs used for SSR primer design
- ESTs excluded for primer design
· exact match between both parents
· short flanking sequence insufficient to design primers
Characteristics of EST-SSRs and efficiency of marker development in cabbage
No. of EST-SSRs (%)
No. of designed primer pairs
No. of primer pairs amplifying product (%)a
No. of polymorphic primers (%)b
Development of EST-based SSR markers
From the 2,214 ESTs containing SSRs, we designed a total of 937 EST-based SSR markers. The remaining 1,677 ESTs were excluded for SSR marker development due to short flanking sequence insufficient for primer design (624 ESTs) or because they had identical sequences in both parents (740 ESTs). The 937 EST-based SSR markers were used for a parental polymorphism survey between lines C1184 and C1234. Successful amplification was obtained from 829 primer sets (88.47%). Polymorphism was identified from 116 (13.99%) primer sets (Additional file 2), and we used 99 SSR markers for further mapping after excluding those giving rise to unclear band patterns and dominant marker types (Table 3). The majority of the primer pairs amplified a single polymorphic locus, except for BoESSR045, which had two polymorphic loci.
Development of EST-based dCAPS markers
Summary of SNPs between homologous EST pairs from two cabbage lines
Type of SNP
A < − > G
C < − > T
A < − > T
G < − > T
G < − > C
A < − > C
Construction of genetic linkage map
Distribution of molecular markers on the cabbage genetic map
No. of mapped loci
Average distance between adjacent loci (cM)
Number of gaps (> 20 cM)
We identified some large gaps throughout the LGs. Twelve gaps with >20 cM between adjacent markers were identified in eight LGs (Table 5). C05 and C09 were each found to have three gaps in their LGs. The largest gaps were detected in C03, with 30.6 cM between BodCAPS22 and CB10267. This indicates that the marker loci were unevenly distributed in the nine LGs of the cabbage genetic map.
Segregation distortion of polymorphic markers
Features of the molecular markers used for mapping
No. of primers used for polymorphic survey
No. of polymorphic primers
No. of mapped markers
No. of unlinked markers
No. of distorted markers
Distribution of molecular markers in the segregation distortion regions of the cabbage linkage map
No. of distorted markers
No. of distorted markers located in SDRsa
Length of SDRs (cM)
Transcriptome sequencing, assembly and gene annotation
Transcriptome sequencing has proven to be an important tool for gene discovery, allele mining and marker development. In this study, the 454 GS-FLX platform was utilized due to its longer read length, which enables high-quality de novo assembly of the transcriptome without a characterized reference genome . Additionally, Newbler v.2.3 software, which is currently the most robust software for 454 transcriptome assembly , was chosen for assembling the sequence reads. Consequently, a vast quantity of potential unique ESTs were generated, representing a large fraction of the cabbage transcriptome, and were further used for development of SSR and SNP markers. The quality of the sequence information obtained here was confirmed by the high percentage of unigenes matching to known proteins by BLASTX and the high rate of successful PCR amplifications.
In order to examine the potential functions represented in the cabbage transcriptome, BLASTX analyses were performed against NCBI non-redundant protein database. At first, we removed singletons shorter than 200 bp because the percentage of BLAST hits usually shows positive correlation with sequence length, as has been confirmed previously in sweet potato  and celery . The BLASTX results revealed that more than 76% of the cabbage transcripts had similarity to known unique proteins. In addition, a large proportion of transcripts (± 89%) were assigned to a remarkable range of GO categories (Figure 2), indicative of the diversity of genes represented by the cabbage transcriptome. Of all assignments made in three categories, most mapped into the GO terms of binding activity (49.5%) and catalytic activity (35.1%) under the molecular function category. Our results for GO term distribution are in agreement with previous studies in rubber tree and pea transcriptome sequencing [43, 44]. Further, we identified and analyzed NBS-encoding disease resistance genes. Although 454 technologies have low sequencing depth, we estimated the NBS-LRR gene expression profiles based on read count. In general, the resistant line showed higher expression of NBS-LRR genes than did the susceptible line. These findings contribute to understanding the evolution of NBS-encoding genes in Brassica species.
General features of EST-SSRs in the cabbage genome
The large collection of EST sequences generated in the present study facilitates the identification of SSRs by in silico mining, which requires relatively little time and has been applied in a variety of plant species . Approximately 5.41% of the 40,947 cabbage unigenes possessed at least one SSR, which is in accord with values reported for other species, ranging from ~2 to ~16% [46–48]. The EST-SSR frequency observed here was higher than previous reports for A. thaliana, maize, tomato, cotton, poplar, and flax [49, 50]. However, it is important to note that values for SSR abundance and frequency among different plant ESTs significantly depend on the parameters used to detect SSRs, the size of the dataset, the database mining tools, and the EST sequence redundancy [45, 51].
Our investigation revealed that trinucleotide repeats are the most common repeat motif attributed to the fact that they can generate non-frame-shift mutations in the coding region , and perhaps result in variation of amino acid residue number at the protein level . Earlier studies demonstrated that AG/CT and AAG/CTT were the predominant di- and trinucleotide SSR motifs, respectively, in plant dicot ESTs [2, 14, 50, 51, 53, 54], which is similar to our observation. These findings suggest that AG and AAG motifs can be considered common features of EST-SSRs in dicot plants.
Marker development and polymorphism level of EST-based SSR markers
EST-SSRs are known to have high level of transferability across taxa and could be useful as anchor markers for comparative mapping and evolutionary studies . In present study, a new set of 937 EST-based SSR markers was developed and 88.47% of them yielded amplification products. The amplification rate observed here is slightly higher than that reported for sweet potato (84.6%)  and tomato (83%) , but lower than that for B. rapa (97.74%) . The success rate for SSR amplification generally ranges between 60-90%, as previously reported for several crop plants .
The EST-based markers designed here showed low polymorphism, which reflects the fact that EST-SSR markers have less polymorphism than genomic SSRs due to highly conserved DNA sequences in genic regions [45, 56]. Another factor likely contributing to the low level of polymorphism is that the parental lines used in this study are elite breeding lines that have a close genetic relationship. Even though EST-based SSR markers exhibit relatively low polymorphism, they may be linked to candidate genes or a trait of interest , and as such can be more valuable than anonymous markers. Moreover, markers designed from compound formations exhibited the highest level of polymorphism. This could be due to compound formations containing more than one SSR motif, which could increase the probability of polymorphism. In addition, previous studies reported that the level of polymorphism of SSR markers is usually correlated with SSR length, as observed in pepper and rice [53, 58]. Likewise, in this study, higher polymorphism was also observed when EST-SSR markers included more than 20 bp of SSR length.
Validation and polymorphism analysis of EST-based dCAPS markers
A considerable number of SNPs were successfully identified and used to design dCAPS markers. The dCAPS method creates polymorphism from SNPs by restriction endonuclease digestion of the PCR products . Approximately 54.44% of the amplified primers showed polymorphism between the two parental lines. The remaining 46% did not show any polymorphism that might be derived from sequencing errors or mis-alignment between paralogous genes of the triplicated Brassica genome [60, 61]. The polymorphism analysis obtained in this study was in accordance with a previous study in cabbage by Wang et al. , in which it was also observed that SNP markers had higher polymorphism than SSR markers. The combination of SSR and SNP markers designed here allowed a larger number of EST-based markers to be mapped onto the cabbage genetic linkage map.
Linkage map construction for cabbage
We effectively constructed a genetic linkage map for cabbage spanning a total 1,331.88 cM, which is slightly larger than the earlier cabbage genetic map (1197.9 cM) generated by Wang et al. . The difference in length between these two maps could be contributed by the difference in the chromosomal recombination frequency caused by environmental factors and the genetic distance between mapping parents, and also the different size of the mapping population as well as the number and types of markers used for map construction .
Some of the mapped EST-based marker loci were found to be clustered in narrow regions, e.g. C03 (6 loci within 6.3 cM), CO7 (3 loci within 2.2 cM) and C08 (3 loci within 1.1 cM). This clustering might correspond to the gene-rich regions of cabbage. Clustered markers in genetic maps were also reported in soybean , pepper , and B. rapa, as well as a previous cabbage map . In addition, we observed 12 gaps in total along this map that varied in size (>20 cM). These gaps were detected in all LGs except LG 2, suggesting that such gaps are not restricted to a particular region of the chromosomes. The presence of these gaps may have negative effects on the application of mapped DNA markers. As mentioned by Cregan et al. , genomic regions that lack DNA markers will make detection of quantitative trait loci (QTL) difficult. Therefore, we plan to develop more markers in the near future to fill in gaps between markers and achieve a high-density genetic linkage map. We also observed that six markers that we used as anchor markers were mapped to the different LGs than previously reported. This could be due to some of those markers producing multiple bands, which could lead to a band derived from a paralogous locus being mapped in our population. Another possible explanation is due to genomic rearrangement which may be observed in genetically unstable population such as in F2 population that we used in this study. As also reported by Wang et al.  that F2 population are temporary and difficult to maintain for long term period.
Segregation distortion phenomenon in the cabbage genetic map
Segregation distortion is a common fact in segregated populations generated from crosses between diverse genotypes [64, 65]. In plants, segregation distortion was first reported in maize , and subsequently in many species including rice , wheat , Arabidopsis , and cabbage . Segregation distortion can have important implications for the construction of a genetic map and QTL mapping, but if addressed properly, distorted markers can also be helpful for QTL mapping . Many factors such as mapping population type, marker type, and genetic relationships of the parents are closely related to the extent of segregation distortion .
A recent study identified 26 SDRs on seven LGs of a cabbage genetic map . In the present study, we detected six SDRs located on six LGs, specifically C01, C03, C05, C07, C08, and C09. The difference in SDR numbers may be related to the different types of mapping population used for constructing the maps, since for the earlier map they used a double haploid (DH) population, whereas in this study an F2 population between elite breeding lines was used. Thus, this finding is in strong agreement with those of Zhang et al.  that segregation distortion is more frequent in DHs and RILs than F2 populations. Moreover, we found that the longest distorted regions were on C05. This result reflects that the distorted markers were non-randomly distributed throughout the genome . Additionally, the fact that markers with segregation distortion are clustered in particular regions indicates that segregation distortion in the F2 population is most likely caused by genetic factors and unlikely to be due to statistical bias, genotyping or scoring errors . The existence of SDRs suggests that there has been a selective process in gametophytes or sporophytes . Also, based on studies in other crops, SDR loci may be linked to sterility genes and pollen-suppressed genes that can affect the selection of partial gametophytes or sporophytes . Overall, the results represent an initial finding of segregation distortion in cabbage; therefore, further investigation is needed to understand better the mechanism underlying the segregation distortion phenomenon in the cabbage genetic map.
The 454 GS-FLX platform has been established to be a powerful tool for de novo transcriptome sequencing due to its long read length. A large number of cabbage EST sequences were generated and used as a reliable source for marker development and discovery of a new candidate disease resistance gene. Subsequently, a novel set of 937 SSR and 97 dCAPS markers were successfully developed and validated using two parental lines of cabbage. Of these, 99 SSR and 21 dCAPS markers revealed clear polymorphism between the two cabbage parental lines, and together with previously developed markers were used to construct a genetic linkage map for cabbage. The map generated herein will facilitate the identification of candidate QTL for economically important traits. In addition, these newly developed markers increase the publicly available EST-based markers in cabbage, which readily can be utilized for other Brassica species. This result demonstrates that transcriptome sequencing using the 454 GS-FLX Titanium sequencer can be a fast and efficient approach for gene discovery and marker development, especially for species without reference genome sequence.
Plant materials and genomic DNA extraction
Young leaf samples of two cabbages parental lines, C1184 and C1234, were collected for RNA extraction. The samples were immediately frozen in liquid nitrogen and stored at −70°C until use. For construction of a genetic linkage map, 97 F2 plants were developed from a cross between C1184 as the female parent and C1234 as the male parent. These two cabbage inbred lines were selected because they are relatively diverse among 16 inbred lines bred for F1 cultivar development in the Joeun Seed company in Korea after a study on their genetic distance based on SSR markers used in a previous report . Additionally, they show different responses to black rot disease: C1184 is susceptible, while C1234 is resistant. All plant materials used in this study were kindly provided by Joeun Seeds, Chungcheongbuk-Do, Korea.
The total genomic DNA was extracted from the leaves of each F2 plant according to the modified cetyltrimethylammonium bromide (CTAB) method . The quality and quantity of the extracted DNA were estimated with a NanoDrop ND-1000 (NanoDrop Technologies, Inc., Wilmington, DE, USA). The final concentration of each DNA sample was adjusted to 10 ng/μL for PCR analysis.
454 transcriptome sequencing and assembly
Total RNA was extracted from approximately 5 g leaf tissue of cabbage C1184 and C1234 using the SV Total RNA Isolation Kit (Promega, Madison, WI) according to the manufacturer’s instructions. cDNA synthesis and library construction from 5 μg extracted mRNAs was then performed as described in the cDNA Rapid Library Preparation Method Manual provided with the Roche GS FLX Titanium Series. Total RNAs were fragmented using a 96 ring Magnetic Particle Concentrator (MPC), and double-stranded cDNA was then synthesized with the cDNA Synthesis System Kit (Roche, IN, USA). Constructed libraries were amplified using emPCR kits (Roche, IN, USA), and sequencing was then performed by 1/8 lane of the 454 GS FLX Titanium Sequencer at the National Instrumentation Center for Environmental Management (NICEM, Seoul National University). The sequence data generated in this study have been deposited at NCBI in the Short Read Archive database under the accession number SRA098802 (experiment accession number SRX338064). The data sets supporting the results of this article can be downloading at http://www.ncbi.nlm.nih.gov/sra/?term=SRA098802. The raw sequence reads generated were assembled by Newbler2.3 software (Roche) with 98% sequence similarity threshold.
To assess the quality of the de novo assembly, a similarity search against the NCBI nr protein database (ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz), was conducted using the BLASTx algorithm with an E value threshold of 10-5. Further, all unigenes were searched against the NCBI non-redundant (nr) protein database (http://www.ncbi.nlm.nih.gov) for functional annotation using BLASTx with an e-value cutoff of 1e-5. The resulting BLAST hits were analyzed for the mapping step in order to retrieve Gene Ontology (GO) terms associated with the hits from the BLAST results. Subsequently, a GO annotation step to select GO terms from the GO pool obtained from the mapping step was performed by the Blast2GO program .
Selection of EST sequences containing SSRs and primer design
The MIcroSAtellite identification tool (MISA) at http://pgrc.ipk-gatersleben.de/misa/misa.html was used for detection of simple sequence repeats (SSRs). The criteria used for detection of EST sequences containing SSRs was a minimum of six repeats for dinucleotide motifs, five repeats for trinucleotide motifs and four repeats for tetra-, penta- and hexa- nucleotide motifs. EST sequences containing SSRs (SSR-ESTs) of cabbage C1234 was BLAST searched against EST sequences of cabbage C1184 (susceptible to black rot disease) using our local database (http://im-crop.snu.ac.kr/). After comparison, only C1234-unique SSR-ESTs, found in C1234 but not in C1184, were used for primer design. Primer pairs were designed for all selected SSR-ESTs from the flanking sequences of SSR motif using the Primer3 program (http://primer3.wi.mit.edu/). The parameters used for primer design were: 55–65°C melting temperature (Tm) with an optimum Tm of 60°C, primer length ranging from 18–24 nt with an optimum size of 20, GC content between 40% and 70% with an optimum set to 50% and product size estimated from 100 to 350 bp. The newly developed EST-SSR markers were designated with the BoESSR (Brassica oleracea EST-based SSR) prefix (e.g., BoESSR001, BoESSR002, BoESSR003, etc.).
SNP discovery and primer design
SNP identification was accomplished by CLC mapping of two cabbage parental lines, C1184 and C1234. Raw reads of C1234 were mapped onto C1184 contigs that were used as reference. In order to improve the accuracy of SNPs, the detected SNPs were then filtered based on the criteria of a minimum 70% of read depth. The selected SNPs were used to develop dCAPS markers using the dCAPS Finder 2.0 program (http://helix.wustl.edu/dcaps) for generation of nearly matched primers including SNP positions . After designing mismatched primers for each SNP, the opposite primers were designed using the Primer3 program (http://primer3.wi.mit.edu/). All of the primers were synthesized by Macrogen (Seoul, Korea).
Molecular marker analysis
A total of 3,570 markers were screened for detection of polymorphisms between the parental lines C1184 and C1234. Of these, 1,034 were EST-based markers comprising 937 EST-based SSR and 97 EST-based dCAPS markers that were developed in this study. Also included were 1,841 intron-based polymorphism (IBP) markers that were developed from B. rapa genome sequences . Furthermore, 695 publically reported SSR markers were used to integrate the reference genetic map: 264 primers derived from the public domain [16, 17], 94 primers from Wang et al. , 71 primers designed from publicly available B. napus genome survey sequences (GSSs) , 45 primers isolated from B. napus[23, 29], 41 primers from Agriculture and Agri-Food Canada , 35 primers obtained from Burgess et al. , 27 primers designed from a microsatellite-enriched genomic library of B. rapa, 24 BAC-derived SSR primers from Kim et al. , 21 primers developed from B. rapa by Choi et al. , 18 primers from EST sequences of B. napus, 14 primers from Long et al. , 12 primers developed by Iniguez-Luy et al. , 11 primers originally isolated from B. napus[33, 34], 9 database sequence-derived primers from Louarn et al. , 4 primers isolated from an A. thaliana library , and one primer each from Lagercrantz et al. , Sebastian et al. , Smith and King , HRI (unpublished data) and Saito et al. .
SSR and IBP analysis
PCR amplifications were performed in a total volume of 10 μL containing 10 ng DNA template, 1X PCR reaction buffer (Inclone Biotech), 0.2 mM each dNTP (Inclone Biotech), 0.2 μM each primer and 1 unit Taq DNA polymerase (Inclone Biotech). The PCR profile was as follows: initial denaturation at 94°C for 4 min, and then 35 cycles of 30 s denaturation at 94°C, 30 s annealing at 55°C - 60°C, 30 s extension at 72°C, and 10 min at 72°C for final extension. The PCR-amplified products were separated by 6% non-denaturing polyacrylamide gel electrophoresis using 1X TBE buffer. The gels were stained with ethidium bromide for 20 min and DNA bands were visualized under UV light using a gel documentation system. The PCR products of some markers were genotyped using Fragment Analyzer, an automated capillary electrophoresis system (Advanced Analytical Technologies Inc., USA), in order to obtain clear separation. The genotyping results were analyzed using PROSizeTM 2.0 analytical software, which can easily screen electropherograms or digital images using the flagging feature (Advanced Analytical Technologies Inc., USA).
PCR reactions were performed in a total volume of 25 μL containing 20 ng DNA template, 0.2 μM each primer set, 1 × PCR buffer, 0.2 mM each dNTP, 1 unit Taq DNA polymerase (VIVAGEN, Korea). Amplification was conducted as described above. The amplified PCR products were digested with appropriate restriction enzymes (3 units) in the presence of the appropriate 1 × buffer, 1 × BSA (if necessary), and distilled water, with incubation at 37°C for more than 3 hours. The products were analyzed using 9% non-denaturing polyacrylamide gel electrophoresis and visualized on a UV trans-illuminator after ethidium bromide staining.
Linkage analysis and map construction
Reproducible polymorphic markers were scored in the F2 population. Linkage analysis and map construction were performed using JoinMap version 3.0 . Linked loci were grouped in the LOD (logarithm of odds) with minimum scores of 2.0, and linkage groups were assigned as C01 to C09, corresponding to the formerly reported map of this species. Locus order within the LOD grouping was generated for each linkage group using a recombination frequency below 0.4 and an LOD score above 0.5 for all marker pairs within each linkage group. The Kosambi function was used to convert recombinant values to genetic distances between the markers . Seventy-five SSR markers derived from the reference map were used as anchor markers in order to assign the newly designed markers in this study to specific linkage groups (LGs). The final genetic linkage map was drawn using MapChart .
This work was supported by the Technology Development Program (No. 309008–05) for Agriculture and Forestry, Ministry of Food, Agriculture, Forestry and Fisheries, Republic of Korea. Nur Kholilatul Izzah is supported by the Islamic Development Bank (IDB), Saudi Arabia.
- UN: Genome-analysis in Brassica with special reference to the experimental formation of Brassica napus and its peculiar mode of fertilization. Jpn J Bot. 1935, 7: 389-452.Google Scholar
- Ramchiary N, Nguyen VD, Li X, Hong CP, Dhandapani V, Choi SR, Yu G, Piao ZY, Lim YP: Genic microsatellite markers in brassica Rapa: development, characterization, mapping, and their utility in other cultivated and wild Brassica relatives. DNA Res. 2011, 18: 305-320.PubMed CentralPubMedView ArticleGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bembem LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer MLI, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.PubMed CentralPubMedGoogle Scholar
- Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92: 255-264.PubMedView ArticleGoogle Scholar
- Rudd S: Expressed sequence tags: alternative or complement to whole genome sequences?. Trends Plant Sci. 2003, 8: 321-329.PubMedView ArticleGoogle Scholar
- Zhou Y, Gao F, Liu R, Feng J, Li H: De novo sequencing and analysis of root transcriptome using 454 pyrosequencing to discover putative genes associated with drought tolerance in Ammopiptanthus mongolicus. BMC Genomics. 2012, 13: 266-PubMed CentralPubMedView ArticleGoogle Scholar
- Sun C, Li Y, Wu Q, Luo H, Sun Y, Song J, Lui EM, Chen S: De novo sequencing and analysis of the American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis. BMC Genomics. 2010, 11: 262-PubMed CentralPubMedView ArticleGoogle Scholar
- Weber AP, Weber KL, Carr K, Wilkerson C, Ohlrogge JB: Sampling the arabidopsis transcriptome with massive parallel pyrosequencing. Plant Physiol. 2007, 144: 32-42.PubMed CentralPubMedView ArticleGoogle Scholar
- Vega-Arreguin JC, Ibarra-Laclette E, Jimenez-Moraila B, Martinez O, Vielle-Calzada JP, Herrera-Estrella L, Herrera-Estrella A: Deep sampling of the palomero maize transcriptome by a high throughput strategy of pyrosequencing. BMC Genomics. 2009, 10: 299-PubMed CentralPubMedView ArticleGoogle Scholar
- Alagna F, D’Agostino N, Torchia L, Servili M, Rao R, Pietrella M, Giuliano G, Chiusano ML, Baldoni L, Perrotta G: Comparative 454 pyrosequencing of transcripts from two olive genotypes during fruit development. BMC Genomics. 2009, 10: 399-PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards CE, Parchman TL, Weekley CW: Assembly, gene annotation and marker development using 454 floral transcriptome sequences in Ziziphus celata (rhamnaceae), a highly endangered, Florida endemic plant. DNA Res. 2012, 19: 1-9.PubMed CentralPubMedView ArticleGoogle Scholar
- Hendre PS, Aggarwal RK: DNA markers: development and application for genetic improvement of coffee. Genomics Assisted Crop Improvement. Genomics Applications in Crops. Edited by: Varshney RK, Tuberosa R. 2007, Netherlands: Springer, 399-434. Volume 2View ArticleGoogle Scholar
- Ashrafi H, Hill T, Stoffel K, Kozik A, Yao J, Reyes S, Wo C, Van Deynze A: De novo assembly of the pepper transcriptome (Capsicum annuum): a benchmark for in discovery of SNPs, SSRs and candidate silico genes. BMC Genomics. 2012, 13: 571-PubMed CentralPubMedView ArticleGoogle Scholar
- Wang W, Huang S, Liu Y, Fang Z, Yang L, Hua W, Yuan S, Liu S, Sun J, Zhuang M, Zhang Y, Zeng A: Construction and analysis of a high-density genetic linkage map in cabbage (Brassica oleracea L. var. capitata). BMC Genomics. 2012, 13: 523-PubMed CentralPubMedView ArticleGoogle Scholar
- Smith LB, King GJ: The distribution of BoCAL-a alleles in Brassica oleracea is consistent with a genetic model for curd development and domestication of the cauliflower. Mol Breed. 2000, 6: 603-613.View ArticleGoogle Scholar
- Lowe AJ, Moule C, Trick M, Edwards KJ: Efficient large-scale development of microsatellites for marker and mapping applications in Brassica crop species. Theor Appl Genet. 2004, 108: 1103-1112.PubMedView ArticleGoogle Scholar
- Piquemal J, Cinquin E, Couton F, Rondeau C, Seignoret E, Doucet I, Perret D, Villeger M-J, Vincourt P, Blanchard P: Construction of an oilseed rape (Brassica napus L.) genetic map with SSR markers. Theor Appl Genet. 2005, 111: 1514-1523.PubMedView ArticleGoogle Scholar
- Long Y, Shi J, Qiu D, Li R, Zhang C, Wang J, Hou J, Zhao J, Shi L, Park BS, Choi SR, Lim YP, Meng J: Flowering time quantitative trait loci analysis of oilseed Brassica in multiple environments and genome wide alignment with arabidopsis. Genetics. 2007, 177 (4): 2433-2444.PubMed CentralPubMedGoogle Scholar
- Radoev M, Becker HC, Ecke W: Genetic analysis of heterosis for yield and yield components in rapeseed (Brassica napus L.) by quantitative trait locus mapping. Genetics. 2008, 179: 1547-1558.PubMed CentralPubMedView ArticleGoogle Scholar
- Cheng X, Xu J, Xia S, Gu J, Yang Y, Fu J, Qian X, Zhang S, Wu J, Liu K: Development and genetic mapping of microsatellite markers from genome survey sequences in Brassica napus. Theor Appl Genet. 2009, 118 (6): 1121-1131.PubMedView ArticleGoogle Scholar
- Kim H, Choi SR, Bae J, Hong CP, Lee SY, Hossain MJ, Nguyen DV, Jin M, Park BS, Bang JW, Bancroft I, Lim YP: Sequenced BAC anchored reference genetic map that reconciles the ten individual chromosomes of Brassica rapa. BMC Genomics. 2009, 10: 432-PubMed CentralPubMedView ArticleGoogle Scholar
- Mun J-H, Kwon S-J, Yang T-J, Seol Y-J, Jin M, Kim J-A, Lim M-H, Kim JS, Baek S, Choi B-S, Yu H-J, Kim D-S, Kim N, Lim K-B, Lee S-I, Hahn J-H, Lim YP, Bancroft I, Park B-S: Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication. Genome Biol. 2009, 10: R111-PubMed CentralPubMedView ArticleGoogle Scholar
- Uzunova MI, Ecke W: Abundance, polymorphism and genetic mapping of microsatellites in oilseed rape (Brassica napus L). Plant Breed. 1999, 118: 323-326.View ArticleGoogle Scholar
- Burgess B, Mountford H, Hopkins CJ, Love C, Ling AE, Spangenberg GC, Edwards D, Batley J: Identification and characterization of simple sequence repeat (SSR) markers derived in silico from Brassica oleracea genome shotgun sequences. Mol Ecol Notes. 2006, 6: 1191-1194.View ArticleGoogle Scholar
- Suwabe K, Iketani H, Nunome T, Kage T, Hirai M: Isolation and characterization of microsatellites in Brassica rapa L. Theor Appl Genet. 2002, 104: 1092-1098.PubMedView ArticleGoogle Scholar
- Louarn S, Torp AM, Holme IB, Andersen SB, Jensen BD: Database derived microsatellite markers (SSRs) for cultivar differentiation in Brassica oleracea. Genet Resour Crop Evol. 2007, 54: 1717-1725.View ArticleGoogle Scholar
- Sandler L, Hiraizum Y, Sandler I: Meiotic drive in natural populations of Drosophila melanogaster. I. The cytogenetic basis of segregation distortion. Genetics. 1959, 44: 233-250.PubMed CentralPubMedGoogle Scholar
- Uzunova M, Ecke W, Weissleder K, Uzunova M, Ecke W, Weissleder K, Röbbelen G: Mapping the genome of rapeseed (Brassica napus L.). I. Construction of an RFLP linkage map and localization of QTLs for seed glucosinolate content. Theor Appl Genet. 1995, 90: 194-204.PubMedView ArticleGoogle Scholar
- Varghese JP, Rudolph B, Uzunova MI, Ecke W: Use of 5’-anchored primers for the enhanced recovery of specific microsatellite markers in Brassica napus L. Theor Appl Genet. 2000, 101: 115-119.View ArticleGoogle Scholar
- Choi SR, Teakle GR, Plaha P, Kim JH, Allender CJ, Beynon E, Piao ZY, Soengas P, Han TH, King GJ, Barker GC, Hand P, Lydiate DJ, Batley J, Edwards D, Koo DH, Bang JW, Park BS, Lim YP: The reference genetic linkage map for the multinational Brassica rapa genome sequencing project. Theor Appl Genet. 2007, 115: 777-792.PubMedView ArticleGoogle Scholar
- Batley J, Hopkins CJ, Cogan NOI, Hand M, Jewell E, Kaur J, Kaur S, Li X, Ling AE, Love C, Mountford H, Todorovic M, Vardy M, Walkiewicz M, Spangenberg GC, Edwards D: Identification and characterization of simple sequence repeat markers from Brassica napus expressed sequences. Mol Ecol Notes. 2007, 7: 886-889.View ArticleGoogle Scholar
- Iniguez-Luy FL, Voort AV, Osborn TC: Development of a set of public SSR markers derived from genomic sequence of a rapid cycling Brassica oleracea L. Genotype. Theor Appl Genet. 2008, 117 (6): 977-985.PubMedView ArticleGoogle Scholar
- Kresovich S, Szewc-McFadden AK, Bliek SM, McFerson JR: Abundance and characterization of simple-sequence repeats (SSRs) isolated from a size-fractionated genomic library of Brassica napus L. (rapeseed). Theor Appl Genet. 1995, 91: 206-211.PubMedView ArticleGoogle Scholar
- Szewc-McFadden AK, Kresovich S, Bliek SM, Mitchell SE, McFerson JR: Identification of polymorphic, conserved simple sequence repeats (SSRs) in cultivated Brassica species. Theor Appl Genet. 1996, 93: 534-538.PubMedView ArticleGoogle Scholar
- Bell CJ, Ecker JR: Assignment of 30 microsatellite loci to the linkage map of arabidopsis. Genomics. 1994, 19: 137-144.PubMedView ArticleGoogle Scholar
- Lagercrantz U, Ellegren H, Andersson L: The abundance of various polymorphic microsatellite motifs differs between plants and vertebrates. Nucleic Acids Res. 1993, 21 (5): 1111-1115.PubMed CentralPubMedView ArticleGoogle Scholar
- Sebastian RL, Howell EC, King GJ, Marshall DF, Kearsey MJ: An integrated AFLP and RFLP Brassica oleracea linkage map from two morphologically distinct doubled-haploid mapping populations. Theor Appl Genet. 2000, 100: 75-81.View ArticleGoogle Scholar
- Saito M, Kubo N, Matsumoto S, Suwabe K, Tsukada M, Hirai M: Fine mapping of the clubroot resistance gene, Crr3, in Brassica rapa. Theor Appl Genet. 2006, 114: 81-91.PubMedView ArticleGoogle Scholar
- Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a non-model organism using 454 pyrosequencing. Mol Ecol. 2008, 17: 1636-1647.PubMedView ArticleGoogle Scholar
- Kumar S, Blaxter ML: Comparing de novo assemblers for 454 transcriptome data. BMC Genomics. 2010, 11: 571-PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Z, Fang B, Chen J, Zhang X, Luo Z, Huang F, Chen X, Li Y: De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweetpotato (Ipomoea batatas). BMC Genomics. 2010, 11: 726-PubMed CentralPubMedView ArticleGoogle Scholar
- Fu N, Wang Q, Shen H-L: De novo assembly, gene annotation and marker development using illumina paired-end transcriptome sequences in celery (Apium graveolens L). PLoS ONE. 2013, 8 (2): e57686-PubMed CentralPubMedView ArticleGoogle Scholar
- Wong CE, Bhalla PL, Ottenhof H, Singh MB: Transcriptional profiling of the pea shoot apical meristem reveals processes underlying its function and maintenance. BMC Plant Biol. 2008, 8: 73-PubMed CentralPubMedView ArticleGoogle Scholar
- Triwitayakorn K, Chatkulkawin P, Kanjanawattanawong S, Sraphet S, Yoocha T, Sangsrakru D, Chanpraset J, Ngamphiw C, Jomchai N, Therawattanasuk K, Tangphatsornruang S: Transcriptome sequencing of Hevea brasiliensis for development of microsatellite markers and construction of a genetic linkage map. DNA Res. 2011, 18: 471-482.PubMed CentralPubMedView ArticleGoogle Scholar
- Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005, 23: 48-55.PubMedView ArticleGoogle Scholar
- Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200.PubMedView ArticleGoogle Scholar
- Gao L, Tang J, Li H, Jia J: Analysis of microsatellites in major crops assessed by computational and experimental approaches. Mol Breed. 2003, 12: 245-261.View ArticleGoogle Scholar
- Kumpatla SP, Mukhopadhyay S: Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species. Genome. 2005, 48: 985-998.PubMedView ArticleGoogle Scholar
- Cardle L, Ramsay L, Milbourne D, Macaulay M, Marshall D, Waugh R: Computational and experimental characterization of physically clustered simple sequence repeats in plants. Genetics. 2000, 156: 847-854.PubMed CentralPubMedGoogle Scholar
- Cloutier S, Niu Z, Datla R, Duguid S: Development and analysis of EST-SSRs for flax (Linum usitatissimum L). Theor Appl Genet. 2009, 119: 53-63.PubMedView ArticleGoogle Scholar
- Wang Z, Li J, Luo Z, Huang L, Chen X, Fang B, Li Y, Chen J, Zhang X: Characterization and development of EST-derived SSR markers in cultivated sweetpotato (Ipomoea batatas). BMC Plant Biol. 2011, 11: 139-PubMed CentralPubMedView ArticleGoogle Scholar
- Metzgar D, Bytof J, Wills C: Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000, 10: 72-80.PubMed CentralPubMedGoogle Scholar
- Yi G, Lee JM, Lee S, Choi D, Kim B-D: Exploitation of pepper EST–SSRs and an SSR-based linkage map. Theor Appl Genet. 2006, 114: 113-130.PubMedView ArticleGoogle Scholar
- Shirasawa K, Oyama M, Hirakawa H, Sato S, Tabata S, Fujioka T, Kimizuka-Takagi C, Sasamoto S, Watanabe A, Kato M, Kishida Y, Kohara M, Takahashi C, Tsuruoka H, Wada T, Sakai T, Isobe S: An EST-SSR linkage Map of Raphanus sativus and comparative genomics of the brassicaceae. DNA Res. 2011, 18: 221-232.PubMed CentralPubMedView ArticleGoogle Scholar
- He C, Poysa V, Yu K: Development and characterization of simple sequence repeat (SSR) markers and their use in determining relationships among Lycopersicon esculentum cultivars. Theor Appl Genet. 2003, 106: 363-373.PubMedGoogle Scholar
- Eujayl I, Sorrells ME, Baum M, Wolters P, Powell W: Isolation of EST-derived microsatellite markers for genotyping the A and B genomes of wheat. Theor Appl Genet. 2002, 104: 399-407.PubMedView ArticleGoogle Scholar
- Gupta PK, Rustgi S: Molecular markers from the transcribed/expressed region of the genome in higher plants. Funct Integr Genom. 2004, 4: 139-162.View ArticleGoogle Scholar
- Cho YG, Ishii T, Temnykh S, Chen X, Lipovich L, McCouch SR, Park WD, Ayres N, Cartinhour S: Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L). Theor Appl Genet. 2000, 100: 713-722.View ArticleGoogle Scholar
- Neff M, Neff J, Chory J, Pepper A: dCAPS, a simple technique for the genetic analysis of single nucleotide polymorphisms: experimental applications in Arabidopsis thaliana genetics. Plant J. 1998, 14 (3): 387-392.PubMedView ArticleGoogle Scholar
- Lukens LN, Quijada PA, Udall J, Pires JC, Schranz ME, Osborn TC: Genome redundancy and plasticity within ancient and recent brassica crop species. Biol J Linn Soc. 2004, 82: 665-674.View ArticleGoogle Scholar
- Yang T-J, Kim JS, Kwon SJ, Lim KB, Kim JA, Choi BS, Jin M, Park JY, Lim MH, Kim HI, Lim YP, Kang JJ, Hong J-H, Kim C-B, Bhak J, Bancroft I, Park B-S: Sequence-level analysis of the diploidization process in the triplicated FLOWERING LOCUS C region of Brassica rapa. Plant Cell. 2006, 18: 1339-1347.PubMed CentralPubMedView ArticleGoogle Scholar
- Song QJ, Marek LF, Shoemaker RC, Lark KG, Concibido VC, Delannay X, Specht JE, Cregan PB: A new integrated genetic linkage map of the soybean. Theor Appl Genet. 2004, 109: 122-128.PubMedView ArticleGoogle Scholar
- Cregan PB, Jarvik T, Bush AL, Shoemaker RC, Lark KG, Kahler AL, Kaya N, VanToai TT, Lohnes DG, Chung J, Specht JE: An integrated genetic linkage map of the soybean. Crop Sci. 1999, 39: 1464-1490.View ArticleGoogle Scholar
- Lyttle TW: Segregation distorters. Annu Rev Genet. 1991, 25: 511-557.PubMedView ArticleGoogle Scholar
- Harushima Y, Nakagahra M, Yano M, Sasaki T, Kurata N: Diverse variation of reproductive barriers in three intraspecific rice crosses. Genetics. 2002, 160: 313-322.PubMed CentralPubMedGoogle Scholar
- Mangelsdorf PC, Jones DF: The expression of mendelian factors in the gametophyte of maize. Genetics. 1926, 11: 423-455.PubMed CentralPubMedGoogle Scholar
- Xu Y, Zhu L, Xiao J, Huang N, McCouch SR: Chromosomal regions associated with segregation distortion of molecular markers in F2, backcross, doubled haploid, and recombinant inbred populations in rice (Oryza sativa L). Mol Gen Genet. 1997, 253 (5): 535-545.PubMedView ArticleGoogle Scholar
- Faris JD, Laddomada B, Gill BS: Molecular mapping of segregation distortion loci in Aegilops tauschii. Genetics. 1998, 149: 319-327.PubMed CentralPubMedGoogle Scholar
- Torjek O, Witucka-Wall H, Meyer RC, von Korff M, Kusterer B, Rautengarten C, Altmann T: Segregation distortion in Arabidopsis C24/Co1-0 and Co1-0/C24 recombinant inbred line populations is due to reduced fertility caused by epistatic interaction of two loci. Theor Appl Genet. 2006, 113: 1551-1561.PubMedView ArticleGoogle Scholar
- Alheit KV, Reif JC, Maurer HP, Hahn V, Weissmann EA, Miedaner T, Würschum T: Detection of segregation distortion loci in triticale (x Triticosecale Wittmack) based on a high-density DArT marker consensus genetic linkage map. BMC Genomics. 2011, 12: 380-PubMed CentralPubMedView ArticleGoogle Scholar
- Liu X, Guo L, You J, Liu X, He Y, Yuan J, Liu G, Feng Z: Progress of segregation distortion in genetic mapping of plants. Res J Agro. 2010, 4 (4): 78-83.View ArticleGoogle Scholar
- Zhang LY, Wang SQ, Li HH, Deng QM, Zheng AP, Li SC, Li P, Li ZL, Wang JK: Effects of missing marker and segregation distortion on QTL mapping in F2 populations. Theor Appl Genet. 2010, 121: 1071-1082.PubMedView ArticleGoogle Scholar
- McDaniel SF, Willis JH, Shaw AJ: A linkage map reveals a complex basis for segregation distortion in an interpopulation cross in the moss Ceratodon purpureus. Genetics. 2007, 176: 2489-2500.PubMed CentralPubMedView ArticleGoogle Scholar
- Plomion C, O’Malley DM, Durel CE: Genomic analysis in maritime pine (Pinus pinaster): comparison of two RAPD maps using selfed and open-pollinated seeds of the same individual. Theor Appl Genet. 1995, 90: 1028-1034.PubMedView ArticleGoogle Scholar
- Perfectti F, Pascual L: Segregation distortion of isozyme loci in cherimoya (Annona cherimola mill). Theor Appl Genet. 1996, 93 (3): 440-446.PubMedView ArticleGoogle Scholar
- Izzah NK, Lee J, Perumal S, Park JY, Ahn K, Fu D, Kim G-B, Nam Y-W, Yang T-J: Microsatellite-based analysis of genetic diversity in 91 commercial Brassica oleracea L. cultivars belonging to six varietal groups. Genet Resour Crop Evol. 2013, 60: 1967-1986.View ArticleGoogle Scholar
- Allen GC, Flores-Vergara MA, Krasynanski S, Kumar S, Thompson WF: A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat Prot. 2006, 1: 2320-2325.View ArticleGoogle Scholar
- Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36: 3420-3435.PubMed CentralPubMedView ArticleGoogle Scholar
- Neff M, Turk E, Kalishman M: Web-based primer design for single nucleotide polymorphism analysis. Trends Genet. 2002, 18 (12): 613-615.PubMedView ArticleGoogle Scholar
- Van Ooijen J, Voorrips R: JoinMap 3.0, Software for the Calculation of Genetic Linkage Maps, Plant Research International. 2001, The Netherlands: WageningenGoogle Scholar
- Kosambi D: The estimation of map distance from recombination values. Ann Hum Genet. 1943, 12: 172-175.Google Scholar
- Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 77-78.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.