- Research article
- Open Access
Genome-wide SNP identification and QTL mapping for black rot resistance in cabbage
BMC Plant Biology volume 15, Article number: 32 (2015)
Black rot is a destructive bacterial disease causing large yield and quality losses in Brassica oleracea. To detect quantitative trait loci (QTL) for black rot resistance, we performed whole-genome resequencing of two cabbage parental lines and genome-wide SNP identification using the recently published B. oleracea genome sequences as reference.
Approximately 11.5 Gb of sequencing data was produced from each parental line. Reference genome-guided mapping and SNP calling revealed 674,521 SNPs between the two cabbage lines, with an average of one SNP per 662.5 bp. Among 167 dCAPS markers derived from candidate SNPs, 117 (70.1%) were validated as bona fide SNPs showing polymorphism between the parental lines. We then improved the resolution of a previous genetic map by adding 103 markers including 87 SNP-based dCAPS markers. The new map composed of 368 markers and covers 1467.3 cM with an average interval of 3.88 cM between adjacent markers. We evaluated black rot resistance in the mapping population in three independent inoculation tests using F2:3 progenies and identified one major QTL and three minor QTLs.
We report successful utilization of whole-genome resequencing for large-scale SNP identification and development of molecular markers for genetic map construction. In addition, we identified novel QTLs for black rot resistance. The high-density genetic map will promote QTL analysis for other important agricultural traits and marker-assisted breeding of B. oleracea.
Cabbage (Brassica oleracea L.) is one of the most important vegetable crops, and is consumed as a food worldwide due to its healthy compounds for humans. Besides its economic importance, cabbage is considered a valuable plant for the study of genome evolution because it contains a CC genome, which represents one of three basic diploid Brassica species in the U’s triangle . Recently, two draft genome sequences of B. oleracea were reported [2,3], and the availability of this reference genome enhances our understanding of the genome architecture of B. oleracea and the evolution of Brassica species, as well as facilitates identification of genes associated with important traits for breeding.
Black rot is one of the most devastating diseases to crucifers including B. oleracea and is caused by the vascular bacterium Xanthomonas campestris pv. campestris (Pammel) Dowson (Xcc). The disease infects the host plants through hydathodes, wounded tissue, insects and stomata [4,5]. The main disease symptoms are V-shaped chlorotic lesions at the margins of leaves, necrosis and darkening of leaf veins, which lead to serious production losses in vegetable crops . Accordingly, development of cultivars resistant to black rot has been a priority for breeders.
Several methods have been attempted to control black rot disease, including crop diversification and rotation, production of disease-free seed, pre-treatment of seed with bactericide, elimination of potential pathogen sources such as infected crop debris and weeds, and planting of resistant cultivars . Among these, utilization of resistant cultivars is one of the most effective and efficient ways to reduce disease incidence and crop loss. However, the development of commercially acceptable resistant varieties has proven to be extremely difficult due to the lack of studies on genetics and breeding for resistance in cabbage. Two major factors hinder black rot resistance breeding in B. oleracea: multigenic control of resistance and emergence of new races of the pathogen that overcome host resistance . Nine races of Xcc have been identified , among which races 1 and 4 are the major pathogens causing black rot disease in B. oleracea crops . Therefore, obtaining B. oleracea cultivars that have resistance to both races is considered a prerequisite to control black rot disease .
Molecular markers are highly useful for genomic analysis and allow exploration of heritable traits and the corresponding genomic variation . DNA markers are now key components of crop improvement programs, and are applied to identify cultivars, analyze genetic diversity, construct linkage maps and identify quantitative trait loci (QTL) . Advances in molecular markers have facilitated the identification of interesting traits via marker-assisted selection (MAS) in plant improvement. Marker-based approaches represent an effective and rapid strategy for identifying and transferring useful genes in breeding programs . Furthermore, the identification of markers linked to QTL can allow analysis of the consistency of QTL effects across different environments and genetic backgrounds, and increase the frequency of favorable alleles during selection . Several QTLs for black rot resistance in B. oleracea have been reported, including two on linkage groups 1 and 9, and two additional QTLs on linkage group 2 , as well as two other major QTLs on linkage groups 2 and 9, and two minor QTLs on linkage groups 3 and 7 . Moreover, three QTLs analyzed using SNP markers in the F2 mapping population derived from a cross between resistant cabbage and susceptible broccoli were found on linkage groups 2, 4 and 5, and exhibited significant effects in black rot resistance . Recently, three further QTLs for black rot resistance were also detected in linkage groups 5, 8 and 9 . In total, 14 QTLs with major and minor effects have been mapped on eight different B. oleracea chromosomes, suggesting that resistance to black rot disease is complex and quantitatively controlled by multiple genes in B. oleracea.
Successful QTL mapping requires a large number of genetic markers . Markers based on simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) are commonly used due to their advantages over other types of genetic markers. SSR markers are highly reproducible, highly polymorphic, and amenable to automation. However, next-generation sequencing (NGS) technology makes SNP markers preferable to SSR markers . SNPs have proved to be universal as well as the most abundant forms of genetic variation even among individuals of the same species . Therefore, SNP markers exhibit higher polymorphism than SSR markers [20,21].
In this study, we have resequenced two parental cabbage lines up to 20× genome coverage and conducted a genome-wide survey for SNPs. We validated the SNPs and developed derived cleaved amplified polymorphic sequences (dCAPS) markers for resistance against black rot disease. The genome-wide catalog of SNPs, the high-density map derived from a mapping population generated from elite cabbage breeding lines with a narrow genetic background, and the QTLs reported herein all will be valuable for both breeding and genetic research in B. oleracea.
Whole-genome resequencing of two cabbage parental lines and SNP detection
Whole genome sequencing data included about 114 million raw reads for C1184 and 113 million for C1234 (Table 1). The recently assembled B. oleracea genome sequence consists of 488.6 Mb, including 446.9 Mb in 9 pseudo-chromosomes and 41.2 Mb of unanchored scaffolds, and corresponding to almost 75% of the estimated genome size (648 Mb) . Our new sequencing data represented approximately 18-fold genome coverage for both parental lines based on the estimated genome size. We mapped each set of paired reads onto the nine pseudo-chromosomes of reference genome sequence. In total, almost 94 million raw reads (82.1%) and 88 million (77.6%) from C1184 and C1234, respectively, were successfully aligned to the reference genome. The average mapping depth was 21.2- and 20-fold for C1184 and C1234, respectively.
The total number of SNPs relative to the reference sequence and average SNP densities were very similar in both parental lines. Approximately 1.20 and 1.24 million high-quality SNPs are identified in C1184 and C1234, respectively, by comparison to the reference genome. On average, a SNP was detected in each 372.8 bp in C1184, and each 360.0 bp in C1234. Chromosome C03 of both lines had the most SNPs, whereas the fewest SNPs were found on chromosome C06 of C1184 and chromosome C04 of C1234.
These SNPs were merged and used to detect SNPs between the two parental lines (Table 2). As a result, a total of 674,521 SNPs were found throughout nine chromosomes, with an average of 1 SNP per 662.5-bp interval. The highest density of SNPs was found on chromosome C03, with a SNP per 541 bp, while the lowest density was on chromosome C05, with one SNP per 818.9 bp. Analysis of the distribution of SNPs per 100 kb along the nine chromosomes revealed areas of high and low SNP density on each chromosome (Figure 1).
Development of dCAPS markers and construction of genetic map
We used the SNPs between C1184 and C1234 for development of dCAPS markers. Based on the physical positions of all markers used in a previous genetic map for B. oleracea , new dCAPS markers were designed for the regions of low marker density. Among 167 markers amplified, 117 (70.1%) were polymorphic between the two parental lines (Table 2 and Additional file 1: Table S1). Among the 117 polymorphic markers, 26 showed heterozygosity in one of parental lines (Table 2). We used 87 of these polymorphic dCAPS markers for genotyping of each individual in the F2 population (Additional file 1: Table S1). Additionally, 16 other types of polymorphic markers including five EST-based dCAPS markers, five MIP markers, three IBP markers, two genomic SSR markers, and one INDEL marker were also genotyped with the same population. Among 103 newly analyzed markers, 25 markers showed a segregation pattern distorted from the 1:2:1 Mendelian ratio in the F2 population, based on chi-square goodness of fit at the 0.05 probability level (Additional file 2: Table S2). There were six segregation distortion regions (SDRs) in the previous map , and all dCAPS markers designed from the SDRs of C01 and C05 showed the same distortion ratio.
The 103 novel polymorphic marker loci (Additional files 1 and 2: Tables S1 and S2) were added to the previous 265 markers  to develop a higher density genetic map. All 368 markers were placed on the map, and a linkage map was generated with nine linkage groups (LGs) in which each LG had more than 32 markers (Figure 2, Table 3). The improved B. oleracea genetic map spanned 1,467.3 cM, which is 135.4 cM more than the previous map, and the average distance between neighboring loci was reduced to 3.88 from 5.02 cM. Most of the new dCAPS markers were mapped to the originally estimated position of each chromosome sequence. The exceptions included BoRSdcaps1-35, which was designed on chromosome C01 but mapped to chromosome C02, and BoRSdcaps5-18, designed on chromosome C05 but mapped to chromosome C09.
Black rot resistance assays and QTL analysis
We performed three independent inoculation trials over three years. The final disease index for F2 plants was determined by calculating the average value of the black rot disease indices for 10 ~ 15 F2:3 progeny plants for each trial. Although all three inoculation tests were performed under the same conditions, the disease symptoms for each test were not consistent and tended to become more severe in later years (Additional file 3: Figure S1), possibly due to differences in plant growth or storage term for the F3 seeds or to weather differences between years.
QTL analyses were performed for each of three trials. We detected significant QTLs, based on higher LOD scores than the thresholds calculated in the permutation tests; LOD threshold values for the tests in 2012, 2013, and 2014 were 3.063, 2.912, and 2.906, respectively. In the first test performed in 2012, there were three significant QTL regions: BRQTL-C1_1 and BRQTL-C1_2 on chromosome C01, and BRQTL-C3 on chromosome C03 (Figure 2). Among these, BRQTL-C1_2 had the highest LOD score, additive effect, and variance explained (Table 4, Figure 2). The second test identified only a single QTL, which was included within BRQTL-C1_2 detected in the 2012 test, although this QTL had smallest LOD score among all QTLs identified in the three tests. The last test, carried out in 2014, identified BRQTL-C1_1 and 2 as well as a novel QTL in chromosome 6, BRQTL-C6. BRQTL-C1_2 in the 2014 test was identified as a smaller region than in 2012, but had the highest LOD score among all QTLs and accounted for 27.3% of the variation.
NBS-encoding genes in QTL regions
In most plants, disease resistance-related genes (R genes) encode proteins containing nucleotide binding sites (NBS) and a series of leucine-rich repeats (LRRs), termed NBS-LRR proteins. NBS-LRR proteins recognize and correspond to pathogen avirulence factors, and lead to defense responses and hypersensitive reactions . Hence, we compared our genetic map to the pseudo-chromosome sequences  and searched for NBS-LRR genes within the QTL regions (Table 5). BRQTL-C1_1 was found between markers H073E22-3 and BoRSdcaps1-11, and BRQTL-C1_2 was between BoRSdcaps1-13 and BoEdcaps4 (Table 4). We identified eight NBS-LRR-encoding genes between H073E22-3 and BoEdcaps4 showing BRQTL-C1_1 and BRQTL-C1_2 QTLs. Seven NBS-LRR type R genes were detected within 1 Mb of the BoESSR291 marker, which is located near the BRQTL-C3 region. BRQTL-C6 contained five NBS-LRR type R genes.
We compared the sequences of these 21 candidate R genes against the Brassica Database (BRAD; http://brassicadb.org/) . All 21 sequences showed similarity to disease resistance proteins, of which 19 and 11 sequences had syntenic genes in B. rapa and A. thaliana, respectively. According to the gene annotation, two candidate disease resistance genes (Bo1g094680 and Bo1g094710) in BRQTL-C1, and seven genes in BRQTL-C3 were found as gene clusters (Table 5).
Seven of nine NBS-LRR genes in BRQTL-C1 were syntenic with the R genes in the counterpart regions of chromosome A01 in B. rapa (Table 5). Orthologous genes of two NBS-LRR genes, Bo1g094680 and Bo1G094710, located within a 63-Kb portion of BRQTL-C1 appeared as tandem array at the counterpart syntenic region in B. rapa and A. thaliana (Figure 3a). All NBS-LRR genes in BRQTL-C6 also showed highly conserved syntenic relationships with counterpart regions in B. rapa and A. thaliana. However, a 72-Kb region near BRQTL-C3 contained a cluster of seven NBS-LRR genes, whereas the syntenic region in B. rapa contained a cluster of only three such genes, and the corresponding syntenic region in A. thaliana did not have any R genes (Figure 3b).
Frequency and utility of SNPs revealed by whole-genome resequencing
An appropriate reference sequence allows whole-genome sequence data from individuals to be aligned, and thus resequencing can be used to detect genetic variation between different samples as high-confidence sequence differences . Accordingly, we performed whole-genome resequencing for two cabbage lines to detect genome-wide SNPs for marker development to construct a high-density genetic map. There are two available draft genome sequences for B. oleracea; although the total assembled sequence of Liu et al. (539.9 Mb)  is larger than that of Parkin et al. (488.6 Mb) , the size of the nine pseudo-chromosomes of the latter (446.9 Mb) is larger than that of the former (388.8 Mb). Therefore, we chose the genome sequence of Parkin et al.  as a reference for this research because of its advantage for sequence-guided SNP marker development.
Approximately 80% of our newly generated PE sequence reads were successfully aligned to the reference genome. Almost 1.2 million SNPs were found in both lines compared to the reference, and 670,000 SNPs were found between C1184 and C1234. The number and density of SNPs between both parental lines were much lower than those detected by comparison with the reference. This could be related to the fact that the plant material used for the reference genome sequencing was kale-like B. oleracea , whereas the two parental lines used in this work were typical cabbages. In addition, we detected fewer SNPs between these lines than the 1.42 million SNPs (averaging one SNP in every 360 bp) previously reported between two other cabbages . Our plant materials have been used as elite breeding resources by a Korean company, and thus the genetic relationship between these two parental inbred lines is likely much closer than to the reference accession or the relationship between the two cabbages used by Liu et al. ; this close relationship likely underlies the relatively low number of SNPs we identified.
Among 167 dCAPS markers that successfully produced PCR products, 70.1% showed polymorphism. This rate was much higher than that of EST-derived dCAPS markers, in which the polymorphic rate was 58.44% when evaluated with low-sequencing depth 454 RNA-seq reads . Paired-end read data generated by Illumina sequencing could allow more accurate alignment of raw reads compared to single reads from 454 RNA-seq, and thus the SNP calling process would also become more precise. The 29.9% of dCAPS markers showing no polymorphism might reflect false mapping of reads to paralogous regions, as there is high sequence similarity between the triplicated genomes and among recently duplicated chromosome segments [26-28]. This could also be the reason two dCAPS markers, BoRSdcaps1-35 and BoRSdcaps5-18, were mapped to unexpected chromosomes.
Collectively, our results demonstrate that whole-genome resequencing data generated by NGS techniques can be highly useful for large-scale discovery of SNPs and development of SNP-based molecular markers. Further study will enable high-throughput genotyping with SNPs detected here.
Improvement of the genetic map between cabbage breeding lines
By obtaining large numbers of reliable SNPs and utilizing them for development of DNA markers, we were able to improve the genetic map of cabbage. The genetic map now spans a total 1,467.3 cM after our addition of SNP markers developed for the relatively large gaps (greater than 20 cM) in the previous map . Consequently, the 12 gaps in the previous genetic map are now reduced to 6 gaps and the average interval is smaller than before. The 368 markers used for the improved genetic map are promising for general cabbage breeding purposes because the map was built using a mapping population between two elite breeding lines with narrow genetic diversity. By contrast, most of previous genetic map was built using mapping populations derived between lines with wide genetic diversity for academic purposes, for example a cross between double-haploid (DH) lines derived from other subspecies [3,29] or DH lines selected based on simple morphological differences . Therefore, the genetic map in this study will be helpful for molecular breeding associated not only with black rot resistance but also with many other important agricultural traits.
QTL mapping of black rot resistance
We identified four QTL regions that could contribute additively to resistance. The BRQTL-C1_2 QTL region was detected repeatedly in the three independent inoculation tests, had the highest LOD values and also accounted for the highest percentage of the variation in all tests. Accordingly, BRQTL-C1_2 is a strong candidate to be a major QTL for black rot resistance. BRQTL-C1_1, BRQTL-C3, and BRQTL-C6 seem to be minor QTLs, which could be influenced by plant conditions and environmental factors. Although two QTLs identified in chromosome C1 included SDRs, we retained all distorted markers for QTL analysis because distorted markers can also be helpful for QTL mapping when they are addressed properly .
The positions of our black rot resistance QTLs did not coincide with those of the 14 previously reported QTLs [4,5,15,16]. This lack of overlap is probably due to differences in disease resistance sources or inocula used. Some studies did not describe the races used [15,16], while some [4,5] used Xcc race 1. The exact race used in this study has not been classified yet. Integrated and standardized protocols for black rot disease races and testing would facilitate further research. However, even though the same Xcc race was used in our three inoculation tests, disease indices for same F2 lineages were not consistent year to year and thus different QTLs were detected between tests. Resistance to Xcc has been reported to vary depending on accessions of B. oleracea and pathogen races [31,32]. Further, the resistance is likely also affected by complex polygenic control under different environmental conditions. Regardless of the race used (if the same as in previous studies or not), the different QTLs detected here should represent new regions.
Candidate genes for black rot resistance
The genomes of B. oleracea, B. rapa, and A. thaliana share a set of 24 conserved syntenic blocks, A to X, that can be identified among the ancestral karyotype . The complete B. oleracea draft genome also demonstrates generally strong conservation with B. rapa in large segments at the pseudo-molecule level [2,3]. Comparative analysis revealed the presence of conserved R gene orthologs at the syntenic counterparts in B. oleracea, B. rapa and A. thaliana. In particular, the BRQTL-C1 region of C1 in B. oleracea showed large-scale conservation with A01 in B. rapa. Our analysis demonstrated that Bra038144, found in unanchored scaffold000140 of the B. rapa genome, is an ortholog of Bo1g087610 in B. oleracea (Table 5). Based on our finding that AT1G57850, the corresponding orthologous gene in A. thaliana, was also located in a syntenic region, the unanchored B. rapa scaffold000140 is likely derived from chromosome A01.
In plant genomes, hundreds of NBS-LRR genes are distributed as single genes or in tandem arrays as gene clusters, which arise from tandem gene duplications or homologous recombination and homogenization [34,35]. We detected 21 R genes in the four QTL regions, of which 9 were in gene clusters (Table 5). Most of the R genes showed conserved syntenic relationships in Brassica and Arabidopsis (Figure 3a). However, near BRQTL-C3 were NBS-LRR gene clusters that appear to be unique to the Brassica lineage (Figure 3b). This result implied that three-R-gene clusters arose by insertion in the Brassica lineage at BrA03 and subsequently amplified to a seven-R-gene cluster in B. oleracea over the 4.6 million years after divergence of the Brassica species .
Although genomes of Brassica-lineage species underwent whole-genome triplication events, the number of resistance genes was not proportionally increased in the Brassica genome [27,36]. Around 150 ~ 200 R genes were reported in the A. thaliana genome [2,35,37], and 206  ~ 244 (http://brassicadb.org/) and 157 genes  were annotated as R genes in B. rapa and B. oleracea genomes, respectively. The 21 NBS-LRR genes found in the four QTL regions are proportionally higher density compared to other chromosomal regions, supporting the idea that some of these NBS-LRRs could be candidate to control black rot resistance in B. oleracea. Further analysis to reveal the function of these genes will be necessary for identification of the major resistance genes for Xcc.
We performed whole-genome resequencing of two cabbage inbred lines that are parental lines for black rot disease resistance and breeding lines with elite agricultural traits. Based on genome-wide SNP detection and validation with dCAPS markers, we report 670,000 SNPs with 70% accuracy between the parental lines. By combining SNP-based markers into the previous genetic map, we improved the genetic map and identified four QTL regions that contained 21 candidate R genes. We thus demonstrated that whole-genome resequencing can successfully be applied for detection of large-scale SNPs, development of molecular markers, and ultimately construction of a high-density genetic map for QTL analysis and marker-assisted breeding of B. oleracea.
Plant materials and whole-genome resequencing
Two cabbage (Brassica oleracea L. var. capitata) inbred lines, C1184 and C1234, were selected as parents to develop a mapping population. The two lines show different responses to black rot disease; C1184 is susceptible to X. campestris pv. campestris (Xcc), whereas C1234 is resistant. The mapping population consisted of 97 F2 plants generated by crossing between C1184 and C1234, as described previously . Furthermore, the 97 F2 plants were vernalized and self-pollinated to produce seeds of F3 progenies for inoculation tests. All plant materials examined in this study were obtained from Joeun Seeds Co. (Chungcheongbuk-Do, Korea).
Genomic DNAs were extracted from approximately 5 g samples of young leaves from the cabbage parental lines, following the modified cetyltrimethylammonium bromide (CTAB) protocol . The quality and quantity of the DNA were examined using a NanoDrop ND-1000 (NanoDrop Technologies, Inc., USA). More than 5 μg extracted DNA was randomly sheared and quantified using DNA 1000 kit (Agilent Technologies, Inc., USA) according to the manufacturer’s protocol. Sequencing with constructed shotgun libraries of C1184 and C1234 was performed by Illumina Hi-seq 2000. Fragmentation, library construction, and sequencing were carried out by the National Instrumentation Center for Environmental Management (NICEM; Seoul, Korea).
SNP discovery and dCAPS marker design
Overall process of SNP discovery was performed by following the framework described by DePristo et al. . Briefly, Illumina paired reads from the parental lines were aligned to the reference sequence of B. oleracea  using Bowtie2 program . Then, read grouping and removal of PCR duplicates were carried out using Picard (http://picard.sourceforge.net). Misalignments caused by INDELs were corrected by local re-alignment using Genome Analysis Toolkit (GATK) and candidate SNPs were called using Variant Caller, a utility in GATK . To filter variants and avoid false positives, candidate SNPs exhibiting any of the following characteristics were removed: (1) mapping quality score lower than 4; (2) quality less than 30; (3) less than 10× or more than 45× mapping depth.
Initially, SNPs of C1184 and C1234 relative to the reference genome were called separately. All of the identified SNP positions from both parental lines were then merged and compared to each other, and promising SNPs for this research between C1184 and C1234 were identified. The selected SNPs were used to develop dCAPS markers using the dCAPS Finder 2.0 program (http://helix.wustl.edu/dcaps) for design of nearly-matched primers including SNP positions. After designing such mismatched primers for each SNP, the opposite primers were designed using the Primer3 program (http://primer3.wi.mit.edu/). All primers were synthesized by Macrogen (Seoul, Korea).
Molecular marker analysis
The newly developed dCAPS markers were validated by examining polymorphisms between the two parental lines C1184 and C1234. Additional expressed sequence tag (EST)-based dCAPS, intron-based polymorphic (IBP), genomic SSR, and INDEL markers that were not included in the previous genetic map  were also analyzed in this study. Furthermore, five polymorphic markers based on miniature inverted transposable element (MITE) insertion polymorphism (MIP) [42,43] were also used for genotyping the F2 population.
PCR amplifications were performed in a total volume of 25 μL containing 20 ng genomic DNA template, 1 × PCR buffer, 20 pM each primer set, 0.2 mM each dNTP, 1 U Taq DNA polymerase (VIVAGEN, Korea). The amplicons of dCAPS markers were mixed with 3 U appropriate restriction enzymes (New England Biolabs, USA), the corresponding 1 × buffer, and 1 × BSA when necessary, then incubated at 37°C for more than three hours. The digested fragments of dCAPS markers and amplicons of other markers stained by ethidium bromide were visualized on a UV trans-illuminator after electrophoresis using 9% non-denaturing polyacrylamide gels or 1% agarose gels depending on fragment size.
Xanthomonas campestris pv. Campestris KACC 10366, obtained from the Korean Agricultural Culture Collection (KACC; Suwon, Korea) were used for the inoculation tests. Inoculum of the bacterium was scraped and cultured on tryptic soy agar (TSA) plates at 30°C for 48 h. Cultured bacteria were harvested using a spreader and diluted with distilled water to 0.125 OD at 600 nm to prepare bacterial suspension for inoculation.
Inoculation tests, carried out in 2012, 2013, and 2014 under the same conditions at the Korea Research Institute of Chemical Technology (Dae-jeon, Korea), were performed with 10 ~ 15 F3 plants of each individual F2 plants selected for genotyping analysis. The F3 seeds were sown and grown on 5 × 8 plastic pots for 20 d in a greenhouse. Afterwards, 20-d-old plants, usually at a stage with two sufficiently developed true leaves, were inoculated by spraying bacterial suspension until adaxial and abaxial surfaces of leaves were sufficiently wet. Each plastic pot (40 plants) received 80 mL bacterial suspension, and the inoculated plants were moved into a dew chamber with the temperature set at 28°C. After 48 h incubation, all plants were transferred to a room maintained at 25°C and 80% humidity for further 7 d incubation with 12 h light/day, and disease symptoms on two inoculated leaves per each plant were surveyed. The severity of the black rot symptoms were recorded based on infected leaf area, with the following disease indices: (0) less than 15%, (1) 15-30%, (2) 30-55%, (3) 55-75%, (4) more than 75% leaf area showing black rot symptoms (Figure 4).
Map construction and QTL analysis
A total of 103 polymorphic markers were genotyped in the F2 population, and the resulting scores were integrated into genotyping data used for a previous genetic map . Linkage analysis and map construction were performed using JoinMap version 4.1 with the same parameters as in the previous study . The Kosambi mapping function was used to convert recombination frequencies into genetic distances.
A disease index for each F2 individual was calculated as the mean grade of 10 ~ 15 F3 seedlings. QTLs for Xcc resistance were evaluated using composite interval mapping (CIM) analysis with QGene program. CIM was performed with LOD (logarithm of odds) threshold values that were estimated using 1,000 permutation tests at 5% significance with 0.5-cM scan intervals.
Quantitative trait loci
Single nucleotide polymorphisms
Xanthomonas campestris pv. campestris
Simple sequence repeat
Nucleotide binding sites
Genome analysis toolkit
Expressed sequence tag
Derived cleaved amplified polymorphic sequences
Miniature inverted transposable element
MITE insertion polymorphism
Composite interval mapping
UN. Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. J Japanese Bot. 1935;7:389–452.
Liu S, Liu Y, Yang X, Tong C, Edwards D, Parkin IA, et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nature Commun. 2014;5:3930.
Parkin IA, Koh C, Tang H, Robinson SJ, Kagale S, Clarke WE, et al. Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. Genome Biol. 2014;15(6):R77.
Kifuji Y, Hanzawa H, Terasawa Y, Nishio T. QTL analysis of black rot resistance in cabbage using newly developed EST-SNP markers. Euphytica. 2013;190(2):289–95.
Tonu NN, Doullah MA-u, Shimizu M, Karim MM, Kawanabe T, Fujimoto R, et al. Comparison of positions of QTLs conferring resistance to xanthomonas campestris pv. campestris in Brassica oleracea. Am J Plant Sci. 2013;4:11.
Williams P. Black rot: a continuing threat to world crucifers. Plant Dis. 1980;64(8):736–42.
Jensen BD, Massomo SM, Swai IS, Hockenhull J, Andersen SB. Field evaluation for resistance to the black rot pathogen Xanthomonas campestris pv. campestris in cabbage (Brassica oleracea). Eur J Plant Pathol. 2005;113(3):297–308.
Hutton SF, Scott JW, Yang W, Sim S-C, Francis DM, Jones JB. Identification of QTL associated with resistance to bacterial spot race T4 in tomato. Theor Appl Genet. 2010;121(7):1275–87.
Fargier E, Manceau C. Pathogenicity assays restrict the species Xanthomonas campestris into three pathovars and reveal nine races within X. campestris pv. campestris. Plant Pathol. 2007;56(5):805–18.
Vicente J, Conway J, Roberts S, Taylor J. Identification and origin of Xanthomonas campestris pv. campestris races and related pathovars. Phytopathology. 2001;91(5):492–9.
Soengas P, Hand P, Vicente J, Pole J, Pink D. Identification of quantitative trait loci for resistance to Xanthomonas campestris pv. campestris in Brassica rapa. Theor Appl Genet. 2007;114(4):637–45.
Duran C, Appleby N, Edwards D, Batley J. Molecular genetic markers: discovery, applications, data storage and visualisation. Curr Bioinformatics. 2009;4(1):16–27.
Hendre PS, Aggarwal RK. DNA markers: development and application for genetic improvement of coffee. In: Varshney RK, Tuberosa R, editors. Genomics-assisted crop improvement. 2nd ed. Netherlands: Springer; 2007. p. 399–434.
Tanksley SD, McCouch SR. Seed banks and molecular maps: unlocking genetic potential from the wild. Science. 1997;277(5329):1063–6.
Camargo L, Williams P, Osborn T. Mapping of quantitative trait loci controlling resistance of Brassica oleracea to Xanthomonas campestris pv. campestris in the field and greenhouse. Phytopathology. 1995;85(10):1296–300.
Doullah M, Mohsin G, Ishikawa K, Hori H, Okazaki K. Construction of a Linkage Map and QTL analysis for Black Rot Resistance in Brassica oleracea L. Int J Nat Sci. 2011;1(1):1–6.
Yu Z, Guo X. Identification and mapping of disease-resistance QTLs in the eastern oyster, Crassostrea virginica Gmelin. Aquaculture. 2006;254(1):160–70.
Mammadov J, Aggarwal R, Buyyarapu R, Kumpatla S. SNP markers and their impact on plant breeding. Int J Plant Genomics. 2012;2012:11.
Rafalski JA. Novel genetic mapping tools in plants: SNPs and LD-based approaches. Plant Sci. 2002;162(3):329–33.
Wang W, Huang S, Liu Y, Fang Z, Yang L, Hua W, et al. Construction and analysis of a high-density genetic linkage map in cabbage (Brassica oleracea L. var. capitata). BMC Genomics. 2012;13(1):523.
Izzah NK, Lee J, Jayakodi M, Perumal S, Jin M, Park B-S, et al. Transcriptome sequencing of two parental lines of cabbage (Brassica oleracea L. var. capitata L.) and construction of an EST-based genetic map. BMC Genomics. 2014;15(1):149.
Hammond-Kosack KE, Jones JD. Plant disease resistance genes. Annu Rev Plant Biol. 1997;48(1):575–607.
Cheng F, Liu S, Wu J, Fang L, Sun S, Liu B, et al. BRAD, the genetics and genomics database for Brassica plants. BMC Plant Biol. 2011;11(1):136.
Bentley DR. Whole-genome re-sequencing. Curr Opin Genet Dev. 2006;16(6):545–52.
Liu J, Huang S, Sun M, Liu S, Liu Y, Wang W, et al. An improved allele-specific PCR primer design method for SNP marker analysis and its application. Plant Methods. 2012;8(1):34.
Lysak MA, Koch MA, Pecinka A, Schubert I. Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005;15(4):516–25.
Yang T-J, Kim JS, Kwon S-J, Lim K-B, Choi B-S, Kim J-A, et al. Sequence-level analysis of the diploidization process in the triplicated flowering locus C region of Brassica rapa. Plant Cell Online. 2006;18(6):1339–47.
Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, et al. The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011;43(10):1035–9.
Iniguez-Luy FL, Lukens L, Farnham MW, Amasino RM, Osborn TC. Development of public immortal mapping populations, molecular markers and linkage maps for rapid cycling Brassica rapa and B. oleracea. Theor Appl Genet. 2009;120(1):31–43.
Alheit KV, Reif JC, Maurer HP, Hahn V, Weissmann EA, Miedaner T, et al. 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(1):380.
Taylor J, Conway J, Roberts S, Astley D, Vicente J. Sources and origin of resistance to Xanthomonas campestris pv. campestris in Brassica genomes. Phytopathology. 2002;92(1):105–11.
Vicente J, Taylor J, Sharpe A, Parkin I, Lydiate D, King G. Inheritance of race-specific resistance to Xanthomonas campestris pv. campestris in Brassica genomes. Phytopathology. 2002;92(10):1134–41.
Schranz ME, Lysak MA, Mitchell-Olds T. The ABC’s of comparative genomics in the Brassicaceae: building blocks of crucifer genomes. Trends Plant Sci. 2006;11(11):535–42.
Holub EB. The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet. 2001;2(7):516–27.
Leister D. Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance genes. Trends Genet. 2004;20(3):116–22.
Town CD, Cheung F, Maiti R, Crabtree J, Haas BJ, Wortman JR, et al. Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy. Plant Cell. 2006;18(6):1348–59.
Meyers BC, Kozik A, Griego A, Kuang H, Michelmore RW. Genome-wide analysis of NBS-LRR–encoding genes in Arabidopsis. Plant Cell. 2003;15(4):809–34.
Allen G, Flores-Vergara M, Krasynanski S, Kumar S, Thompson W. A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat Protoc. 2006;1(5):2320–5.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Nekrutenko A, Taylor J. Next-generation sequencing data interpretation: enhancing reproducibility and accessibility. Nat Rev Genet. 2012;13(9):667–72.
Sampath P, Lee S-C, Lee J, Izzah NK, Choi B-S, Jin M, et al. Characterization of a new high copy stowaway family MITE, BRAMI-1 in Brassica genome. BMC Plant Biol. 2013;13(1):56.
Sampath P, Murukarthick J, Izzah NK, Lee J, Choi H-I, Shirasawa K, et al. Genome-wide comparative analysis of 20 miniature inverted-repeat transposable element families in Brassica rapa and B. oleracea. PLoS One. 2014;9(4):e94499.
This research was supported by the Golden Seed Project (Center for Horticultural Seed Development, No. 213003-04-1-SB430), Ministry of Agriculture, Food and Rural Affairs (MAFRA), Ministry of Oceans and Fisheries (MOF), Rural Development Administration (RDA) and Korea Forest Service (KFS).
The authors declare that they have no competing interests.
HJL, NKI, HJJ, and JL carried out the molecular experiments. MJ and SP carried out SNP discovery. JS, JY, YS, and KA developed the mapping population and maintained plant materials. JHL and GJC performed inoculation tests and investigated disease symptoms. K-WY, JYP, and S-CL collected plant materials, and provided technical assistance. JL interpreted the results, and wrote the manuscript. YY, HK, I-SN, and T-JY conceived of and managed the research. All authors critically read and approved the final version of the manuscript.
Description of polymorphic markers between C1184 and C1234 used in this study.
Results of the chi-square goodness-of-fit tests for the observed segregation ratios with the genotyped markers among F2 plants.
Disease index distribution of F2 population, evaluated by average scores from inoculated F3 plants.
About this article
Cite this article
Lee, J., Izzah, N.K., Jayakodi, M. et al. Genome-wide SNP identification and QTL mapping for black rot resistance in cabbage. BMC Plant Biol 15, 32 (2015). https://doi.org/10.1186/s12870-015-0424-6
- Whole-genome resequencing
- Genetic linkage map
- Black rot