Impact Factor 5.753 | CiteScore 8.2
More on impact ›


Front. Plant Sci., 18 June 2018 |

Identification of Quantitative Trait Loci for Clubroot Resistance in Brassica oleracea With the Use of Brassica SNP Microarray

Lisha Peng1,2,3,4, Lili Zhou1,2,3, Qinfei Li1,2,3, Dayong Wei1,2,3, Xuesong Ren1,2,3, Hongyuan Song1,2,3, Jiaqin Mei5,6*, Jun Si1,2,3* and Wei Qian5,6
  • 1College of Horticulture and Landscape Architecture, Southwest University, Chongqing, China
  • 2Key Laboratory of Horticulture Science for Southern Mountains Regions, Ministry of Education, Chongqing, China
  • 3Chongqing Key Laboratory of Olericulture, Chongqing, China
  • 4Chongqing Yudongnan Academy of Agricultural Sciences, Chongqing, China
  • 5College of Agronomy and Biotechnology, Southwest University, Chongqing, China
  • 6Academy of Agricultural Sciences, Southwest University, Chongqing, China

Increasing clubroot resistance (CR) of Brassica oleracea by ascertaining the molecular mechanisms has been the key focus in modern B. oleracea breeding. In order to identify the quantitative trait loci (QTLs) associated with CR in B. oleracea, 94 F2 vegetative lines which were developed by tissue culture of selfed seeds from the F1 generation between a clubroot-resistant B. oleracea inbred line and a susceptible line, were identified for disease incidence and six CR-associated traits under a lab inoculation by Plasmodiophora brassicae and were genotyped with the 60K Brassica SNP array. Significant correlations were detected for numbers of fibrous roots and P. brassicae content in roots with disease incidence. Nine linkage groups were constructed from 565 bins which covered around 3,000 SNPs, spanning 1,028 cM of the B. oleracea genome with an average distance of 1.82 cM between adjacent bins. A total of 23 QTLs were identified for disease incidence and the other two correlated traits, individually explaining 6.1–17.8% of the phenotypic variation. Several overlaps were detected among traits, including one three-traits-overlapped locus on linkage group C08 and two important overlapped regions between the two CR-associated traits on C06. The QTLs were compared with known CR loci/genes and the novelty of our QTLs was discussed.


Clubroot, caused by the soil-borne obligate Plasmodiophora brassicae, is a devastating disease in Brassica crops including cabbage (Brassica oleracea L. var. capitata) which is one of the most important vegetable crops in the world (Hirai, 2006; Dixon, 2009). It causes serious yield loss in cabbage since the pathogen always induces galls on the plant root, thus hinders the uptake of water and nutrients and finally leads to abnormal growth (Dixon and Robinson, 2010). It is hard to control clubroot by cultural managements or chemical fungicides due to the long period of the pathogen surviving in soil. Therefore, developing resistance cultivars is the most effective way to control this disease.

Identification of clubroot resistance (CR) quantitative trait loci (QTLs) or genes is of great importance in the resistance breeding. Complete resistant accessions against specific pathogen isolates were found in European fodder turnips (B. rapa) and at least eight resistant loci were identified in B. rapa including Crr1, Crr2, Crr3, Crr4, CRa, CRb, CRc, and CRk (Matsumoto et al., 1998; Piao et al., 2002, 2004; Hirai et al., 2004; Saito et al., 2006; Suwabe et al., 2006; Sakamoto et al., 2008; Hatakeyama et al., 2013; Kato et al., 2013; Huang et al., 2017; Yu et al., 2017). These loci were reported to control CR in a qualitative plus quantitative manner (Piao et al., 2002; Suwabe et al., 2006; Sakamoto et al., 2008). Likewise, a few CR QTLs, such as CR2a, CR2b, pb-3, pb-4, pb-Bol and PbBo(Anju)14, were identified from B. oleracea (Landry et al., 1992; Voorrips et al., 1997; Rocherieux et al., 2004; Moriguchi et al., 2010; Nagaoka et al., 2010; Lee et al., 2016) which might be partially resistant to clubroot (Crisp et al., 1989; Piao et al., 2009). Being differently, the CR trait is possibly controlled by many quantitative loci in B. oleracea (Piao et al., 2009; Nagaoka et al., 2010; Tomita et al., 2013). Therefore, identify more CR loci from different B. oleracea resources with wide genetic basis will be benefit to the breeding of clubroot-resistant cabbages.

In our previous identification, a cabbage inbred line was found to be highly resistant to clubroot disease, thus an F2 segregating population was developed from the cross between this line and a susceptible B. oleracea. On the other hand, a 60K Brassica SNP microarray was successfully released (Clarke et al., 2016) and used for the construction of genetic linkage map with high density of markers in B. oleracea (Mei et al., 2017). In the present study, disease incidence and several CR related traits were investigated in this F2 population in lab tests, and QTLs for these traits were identified by using the linkage group built by data from the SNP microarray. Our study will be helpful for clubroot resistance breeding in B. oleracea.

Materials and Methods

Plant Materials and Resistance Evaluation

An F2 segregating population, comprising of 94 cloned-lines developed by tissue culture, was developed from hybridization between two inbred lines of B. oleracea, ‘263’ and ‘GZ87’ with diverse resistance levels against P. brassicae. The development of clones for each F2 genotype was conducted according to Luo et al. (2000). The vegetative plants were transplanted into 90 mm × 80 mm pots after rooting and kept in a climate chamber for 1 week (16/8 h light/dark cycle under 26/20°C). Then each plant was inoculated by watering 2 mL P. brassicae (the 4th race) resting spores suspension (2 × 108 spores/mL) at the stem base using a pipette method (Carlsson et al., 2004; Nagaoka et al., 2010) and kept in the climate chamber with the same condition as before. Phenotypic data were collected at 6 weeks post-inoculation, including disease incidence (DIC, percentage of diseased plants in total plants), fresh weight per plant (FW), fibrous root weight/root weight (FR/R), length of root (LR), number of fibrous roots (NFRs), ratio of root surface covered with fibrous roots (RFR), and P. brassicae content in roots (PCR). Four biological replicates of resistance test were conducted in two rounds. Ten to twenty two plants per line were tested in each replicate for the calculation of DIC and five plants for other six indexes. Pearson’s simple correlations were calculated between traits of interest via SAS software (SAS Institute, 1999).

Genotyping, Map Construction, and QTL Analysis

Genomic DNA was extracted from young leaf of plant using the CTAB method (Allen et al., 2006). DNA samples of two parental lines and the 94 F2 lines were genotyped by the Brassica 60K Bead Chip Array (Infinium®, Illumina, Inc., San Diego, CA, United States) (Clarke et al., 2016). SNPs were aligned to the reference genomes of B. oleracea (version 1.1)1 using a local BLAST search. Bins were developed for each chromosome using Perl language by combining SNPs with identical genotypes across the F2 population (Zhang et al., 2014) and then used for genetic map construction in IciMapping version 4.1 with default parameters (Wang et al., 2016). QTLs were analyzed using the inclusive composite interval mapping (ICIM) model in IciMapping with optimized parameters (Step = 1 cM; PIN = 0.005). A permutation test with 1,000 permutations was performed for each trait to calculate the threshold of LOD score at the significance level of P = 0.05.


Phenotypic Performance of Parents and F2 Lines

According to the field identification for three successive years at Fuling, Chongqing, China (E107.6459, N29.5709) where happens serious clubroot disease every year caused mainly by the 4th race of P. brassicae (identified according to Williams host system, data not shown), ‘GZ87’ exhibited complete resistance with a disease index of 0, while ‘263’ which is a founder parent in our cabbage breeding program showed stable susceptibility to P. brassicae (with average DIC of 72.6% and disease index of 33.7). Under the lab identification, similarly, ‘GZ87’ exhibited obviously higher resistance (DIC = 0) than ‘263’ (DIC = 80.8) (Figure 1). Wide variations were detected among the F2 lines for all the seven traits, of which FW, LR, and RFR showed normal distributions, while DIC, FR/R, NFR, and PCR exhibited skew distributions (Table 1). Three traits, particularly PCR and NFR, were found to be significantly correlated with DIC (r= -0.881 and -0.640, respectively) (Table 2). Therefore, QTLs were subsequently screened for DIC, PCR, and NFR.


FIGURE 1. The root morphology of ‘GZ87’ and ‘263’ after infection by Plasmodiophora brassicae under field (A) and lab (B) assays.


TABLE 1. The phenotypic performance of clubroot resistance associated traits in the parental lines and the F2 population of Brassica oleracea.


TABLE 2. Correlation among the clubroot resistance associated traits in Brassica oleracea.

SNP Detection and Genetic Map Construction

Among 52,157 SNPs in the Brassica 60K SNP array, 21,646 (41.5%) could be aligned to the B. oleracea genome. These SNPs distributed unevenly on the nine chromosomes of B. oleracea, with the highest number on chromosome C03 (3,558) and the fewest on chromosome C05 (1,264), being in accordant with the chromosome length of B. oleracea. Among the C-genome SNPs, 3,218 (14.9%) were polymorphic between the two parents, forming 615 bins. After removing the serious missing bins, the skewed bins (P≤ 0.01) and the redundant bins, 565 bins were allocated into nine linkage groups (LGs), with the fewest on C01 (16 bins) and the most number on C03 (130 bins) (Table 3). The lengths of LGs ranged from 17.2 to 187.7 cM, summed in a total length of 1,027.98 cM, with an average distance of 1.82 cM between neighboring bins (Figure 2 and Table 3). The bins showed general coincidence between their genetic and physical positions for all LGs except C01 which has a limited number of bins.


TABLE 3. Overview of bins on the nine linkage groups of Brassica oleracea.


FIGURE 2. The Brassica oleracea genetic linkage map and quantitative trait loci (QTLs) for disease incidence (DIC), the number of fibrous roots (NFRs), and the P. brassicae content in roots (PCR).

QTLs for CR-Associated Traits

With significant LOD thresholds of 4.4, 4.5, and 4.4 which were set the by 1,000-permutation test in IciMapping, 2, 13, and 8 QTLs for DIC, NFR, and PCR, were detected respectively (Figures 2, 3 and Table 4), individually explaining 6.1–17.7% of the phenotypic variation. Overlaps were detected between two rounds, such as DIC.I-1 and DIC.II-1, and PCR.I-2 and PCR.II-5. Overlaps were further detected among traits, including one three-traits-overlapped QTL on LG C08 and four overlapped QTLs between NFR and PCR on LG C03, C04, and C06. The QTLs on LG C08 (DIC.I-1, DIC.II-1, NFR.II-8, and PCR.II-6) exhibited 10.2–17.8% of the total phenotypic variation, with an overlapped genetic range of 64.5–67.5 cM which corresponding to a 0.6 Mb interval (32.7–33.3 Mb) on chromosome C08. It is worth mentioning that two overlapped QTL regions on LG C06 between NFR and PCR also exhibited high R2 values. In details, NFR.I-5, PCR.I-2, and PCR.II-5 overlapped on 86.5–92.5 cM which corresponding to 1.98–3.54 Mb on chromosome C06, and NFR.II-5 and PCR.II-4 overlapped on 77.5–83.5 cM, being corresponded to 3.26–3.94 Mb of C06.


FIGURE 3. Quantitative trait loci scanning for clubroot resistance associated traits in B. oleracea in the first (red) and the second replication (green).


TABLE 4. Quantitative trait loci (QTLs) for the number of disease incidence (DI), fibrous roots (NFR), the root surface covered with fibrous roots (RFR), and the Plasmodiophora brassicae content in roots (PCR) in Brassica oleracea.

Comparative Analysis Among CR Loci

According to Bolbase2 which provides information on syntenic regions between B. oleracea and B. rapa (Yu et al., 2013), a QTL interval on C03 (6.28–6.32 Mb, overlapped by NFR.II-3 and PCR.II-2) in the present study was partially syntenic to 0.59–6.22 Mb on A02 of B. rapa where located the CRc- and Pb(Anju)2- linked marker ‘m6R’ (from 2,112,653 to 2,113,153 bp) (Sakamoto et al., 2008; Nagaoka et al., 2010; Lee et al., 2016). A QTL on C07 (NFR.II-7) was found to be syntenic to 11,077,112–11,414,470 bp on A08 of B. rapa, being partially overlapped with a reported CR loci Crr1 which was fine mapped to 10,692,602–11,617,700 bp on A08 in B. rapa (Hasan and Rahman, 2016). The three important regions on C08 and C06 identified no accordance with reported CR loci.


In the past, QTLs were mainly identified based on genetic linkage map constructed by DNA molecular markers such as RAPD, RFLP, AFLP, and SSR markers (Landry et al., 1992; Voorrips et al., 1997; Piao et al., 2002, 2004; Suwabe et al., 2006), however, this approach is high-input but low-output in identifying candidate genes. Comparatively, the Brassica SNP microarray which has been widely applied in B. napus could provide high-density genetic maps and thus greatly narrowed the QTL regions for interested traits (Liu et al., 2013; Wei et al., 2016). Since no B. oleracea-specific SNP microarray was available before this study, the 60K Brassica SNP microarray which consisted of probes from both A and C genomes (Clarke et al., 2016) was used in the present study. Although nearly 60% (30,511) SNPs were filtered out (could not be aligned to C genome), the 3,218 polymorphic sites from the remaining 40% (21,646) SNPs still enabled us to construct a high-density genetic map of B. oleracea which spanned 1028 cM of the B. oleracea genome with an average of 1.82 cM (0.67 Mb) between neighboring bins. Given the high-throughput and time-saving nature of SNP array, this approach will dramatically accelerate the process of QTL fine-mapping in B. oleracea.

In the studies associated with CR in Brassica, disease index was the most widely used indicator for CR level (Crisp et al., 1989; Landry et al., 1992; Matsumoto et al., 1998; Carlsson et al., 2004; Hirai et al., 2004; Hirai, 2006; Dixon and Robinson, 2010; Hatakeyama et al., 2013; Kato et al., 2013; Liu et al., 2013; Hasan and Rahman, 2016; Lee et al., 2016; Huang et al., 2017; Yu et al., 2017). However, it was hard to judge the disease index of vegetative propagates in the present study since the disease grade of each copy was difficult to class due to the lack of main roots in vegetative plants and the difficulties on distinguishing clubroot tubers from exogenous-hormone-induced calluses in tissue culture. Nevertheless, it was relatively easier to discriminate diseased plants from healthy ones, thus disease incidence was recorded, though its accuracy to indicate the resistance level of host might be not as high as that by using disease index. As supplements of disease incidence, we measured several traits which might be potentially correlated with the resistance of vegetative plants, and we found the P. brassicae content in roots was highly correlated with DIC, being consistent with that in B. rapa (r = 0.95) (Zhu et al., 2015). In addition, the NFR was found to be moderately correlated with DIC and high correlated with PCR. Therefore, the two CR-associated traits particular PCR could be used to indicate the resistance level of plants when needed, for example, in the plants without main roots.

Clubroot-resistant genes or loci were widely identified in B. rapa (Hirai et al., 2004; Piao et al., 2004; Suwabe et al., 2006; Sakamoto et al., 2008; Huang et al., 2017; Yu et al., 2017) and B. oleracea (Landry et al., 1992; Voorrips et al., 1997; Rocherieux et al., 2004; Moriguchi et al., 2010; Nagaoka et al., 2010; Lee et al., 2016). It was hypothesized that the CR of B. rapa was possibly controlled by qualitative plus quantitative loci (Piao et al., 2002; Suwabe et al., 2006; Sakamoto et al., 2008), while CR trait in B. oleracea might be controlled by QTLs in a quantitative manner (Crisp et al., 1989; Voorrips et al., 1997; Piao et al., 2009; Tomita et al., 2013). In the present study, although only one genetic region on C08 was identified for DIC, many QTLs were identified for the other two CR-associated traits. These QTLs including the two QTLs for DIC explained limited phenotypic variations (6.1–17.8%), suggesting the quantitative nature of these traits and indicated that the pyramiding of these loci may confer durable CR in B. oleracea.

Two possible overlaps of our QTLs were detected with previous reported loci after deep comparison of our regions with other studies which provided direct or indirect (primer sequences of linked markers were then blasted to reference genome to get the physical positions) genomic information of CR loci (Matsumoto et al., 1998; Saito et al., 2006; Suwabe et al., 2006; Sakamoto et al., 2008; Nagaoka et al., 2010; Hatakeyama et al., 2013; Kato et al., 2013; Hasan and Rahman, 2016; Lee et al., 2016), however, the QTLs in these two regions only explained 6.7–10.3% of the total phenotypic variation. It is interesting that no overlap was detected for reported CR loci with the three important regions found from C08 and C06 in our study. Although A08 of B. rapa was reported to carry important CR QTLs such as Crr1 and Rcr9 (Hasan and Rahman, 2016; Yu et al., 2017), they are not syntenic to our region on C08. These data suggest that the CR in our resistant B. oleracea GZ87 is possibly controlled by some novel loci. In each of the three genetic regions, we found several pathogen-responsive genes, including genes encoding receptor like proteins (R proteins), receptor binding proteins, auxin efflux transporters, auxin-responsive proteins, RHD3s (root hair defective 3, required for cell expansion and normal root hair development), pathogen-responsive factors (alpha-dioxygenase, calcium-binding, thaumatins, et al.) and transcriptional factors (MYB, ERF, C2H2, HD-ZIP, et al.). Further study will be carried out to identify the candidate genes and develop linked markers for a pyramiding purpose, since accumulation of QTLs could confer broad-spectrum CR in B. oleracea (Tomita et al., 2013).

Author Contributions

LP, LZ, QL, and DW conducted the lab experiments and data analysis. XR and HS conducted the field experiments. LP, JM, JS, and WQ wrote the manuscript. JM, JS, and WQ directed the project.


This study was financially supported by the National Key Research and Development Program of China (2017YFD0101804), the Science and Technology Innovation Program for the Social Undertakings and the People’s Livelihood in Chongqing (cstc2015shms-ztzx80005; cstc2015shms-ztzx80007; cstc2015shms-ztzx80009; and cstc2016shmszx0674), and the Fundamental Research Funds for the Central Universities (XDJK2018B022).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


  1. ^
  2. ^


Allen, G. C., FloresVergara, M. A., Krasynanski, S., Kumar, S., and Thompson, W. F. (2006). A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat. Protoc. 1, 2320–2325. doi: 10.1038/nprot.2006.384

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlsson, M., Bothmer, R. V., and Merker, A. (2004). Screening and evaluation of resistance to downy mildew (Peronospora parasitica) and clubroot (Plasmodiophora brassicae) in genetic resources of Brassica oleracea. Hereditas 141, 293–300. doi: 10.1111/j.1601-5223.2004.01818.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, W. E., Higgins, E. E., Plieske, J., Wieseke, R., Sidebottom, C., Khedikar, Y., et al. (2016). A high-density SNP genotyping array for Brassica napus and its ancestral diploid species based on optimised selection of single-locus markers in the allotetraploid genome. Theor. Appl. Genet. 129, 1–13. doi: 10.1007/s00122-016-2746-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Crisp, P., Crute, I. R., Sutherland, R. A., Angell, S. M., Bloor, K., Burgess, H., et al. (1989). The exploitation of genetic resources of Brassica oleracea in breeding for resistance to clubroot ( Plasmodiophora brassicae ). Euphytica 42, 215–226.

Google Scholar

Dixon, G. R. (2009). The occurrence and economic impact of Plasmodiophora brassicae and clubroot disease. J. Plant Growth Regul. 28, 194–202. doi: 10.1007/s00344-009-9090-y

CrossRef Full Text | Google Scholar

Dixon, G. R., and Robinson, D. L. (2010). The susceptibility of Brassica oleracea cultivars to Plasmodiophora brassicae (clubroot). Plant Pathol. 35, 101–107. doi: 10.1111/j.1365-3059.1986.tb01987.x

CrossRef Full Text | Google Scholar

Hasan, M. J., and Rahman, H. (2016). Genetics and molecular mapping of resistance to Plasmodiophora brassicae pathotypes 2,3,5,6, and 8 in rutabaga (Brassica napus var. napobrassica). Genome 59, 809–815. doi: 10.1139/gen-2016-0034

PubMed Abstract | CrossRef Full Text | Google Scholar

Hatakeyama, K., Suwabe, K., Tomita, R. N., Kato, T., Nunome, T., Fukuoka, H., et al. (2013). Identification and characterization of Crr1a, a gene for resistance to clubroot disease (Plasmodiophora brassicae Woronin) in Brassica rapa L. PLoS One 8:e54745. doi: 10.1371/journal.pone.0054745

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirai, M. (2006). Genetic analysis of clubroot resistance in Brassica crops. Breed. Sci. 56, 223–229. doi: 10.1270/jsbbs.56.223

CrossRef Full Text | Google Scholar

Hirai, M., Harada, T., Kubo, N., Tsukada, M., Suwabe, K., and Matsumoto, S. (2004). A novel locus for clubroot resistance in Brassica rapa and its linkage markers. Theor. Appl. Genet. 108, 639–643. doi: 10.1007/s00122-003-1475-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Z., Peng, G., Liu, X., Deora, A., Falk, K. C., Gossen, B. D., et al. (2017). Fine mapping of a clubroot resistance gene in Chinese cabbage using SNP markers identified from bulked segregant RNA sequencing. Front. Plant Sci. 8:1448. doi: 10.3389/fpls.2017.01448

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, T., Hatakeyama, K., Fukino, N., and Matsumoto, S. (2013). Fine mapping of the clubroot resistance gene CRb and development of a useful selectable marker in Brassica rapa. Breed. Sci. 63, 116–124. doi: 10.1270/jsbbs.63.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Landry, B. S., Hubert, N., Crete, R., Chang, M. S., Lincoln, S. E., and Etoh, T. (1992). A genetic map of Brassica oleracea based on RFLP markers detected with expressed DNA sequences and mapping resistance genes to race 2 of Plasmodiophora brassicae (Woronin). Genome 35, 409–420. doi: 10.1139/g92-061

CrossRef Full Text | Google Scholar

Lee, J., Izzah, N. K., Choi, B. S., Joh, H. J., Lee, S. C., Perumal, S., et al. (2016). Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.). DNA Res. 23, 29–41. doi: 10.1093/dnares/dsv034

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Qu, C., Wittkop, B., Yi, B., Xiao, Y., He, Y., et al. (2013). A high-density SNP map for accurate mapping of seed fibre QTL in Brassica napus L. PLoS One 8:e83052. doi: 10.1371/journal.pone.0083052

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, P., Lang, Z., Deng, J., and Wang, Z. (2000). Application of in vitro organ culture in wide-cross breeding of rapeseed. Euphytica 114, 217–221. doi: 10.1023/A:1003911507540

CrossRef Full Text | Google Scholar

Matsumoto, E., Yasui, C., Ohi, M., and Tsukada, M. (1998). Linkage analysis of RFLP markers for clubroot resistance and pigmentation in Chinese cabbage (Brassica rapa ssp. pekinensis). Euphytica 104:79. doi: 10.1023/A:1018370418201

CrossRef Full Text | Google Scholar

Mei, J., Wang, J., Li, Y., Shuai, T., Wei, D., Shao, C., et al. (2017). Mapping of genetic locus for leaf trichome in Brassica oleracea. Theor. Appl. Genet. 130, 1953–1959. doi: 10.1007/s00122-017-2936-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Moriguchi, K., Kimizuka-Takagi, C., Ishii, K., and Nomura, K. (2010). A genetic map aased on RAPD, RFLP, isozyme, morphological markers and QTL analysis for clubroot resistance in Brassica oleracea. Breed. Sci. 49, 257–265. doi: 10.1270/jsbbs.49.257

CrossRef Full Text

Nagaoka, T., Doullah, M. A., Matsumoto, S., Kawasaki, S., Ishikawa, T., Hori, H., et al. (2010). Identification of QTLs that control clubroot resistance in Brassica oleracea and comparative analysis of clubroot resistance genes between B. rapa and B. oleracea. Theor. Appl. Genet. 120, 1335–1346. doi: 10.1007/s00122-010-1259-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Piao, Z. Y., Deng, Y. Q., Choi, S. R., Park, Y. J., and Lim, Y. P. (2004). SCAR and CAPS mapping of CRb, a gene conferring resistance to Plasmodiophora brassicae in Chinese cabbage ( Brassica rapa ssp. pekinensis). Theor. Appl. Genet. 108, 1458–1465. doi: 10.1007/s00122-003-1577-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Piao, Z. Y., Park, Y. J., Ryun Choi, S., Pyo Hong, C., Jee Young, P., Young Sang, C., et al. (2002). Conversion of an AFLP marker linked to clubroot resistance gene in Chinese cabbage into a SCAR marker. Hortic. Environ. Biotechnol. 43, 653–659.

PubMed Abstract | Google Scholar

Piao, Z. Y., Ramchiary, N., Yongpyo, L., and Dixon, G. R. (2009). Genetics of clubroot resistance in Brassica species. J. Plant Growth Regul. 28, 252–264. doi: 10.1007/s00344-009-9093-8

CrossRef Full Text | Google Scholar

Rocherieux, J., Glory, P., Giboulot, A., Boury, S., Barbeyron, G., Thomas, G., et al. (2004). Isolate-specific and broad-spectrum QTLs are involved in the control of clubroot in Brassica oleracea. Theor. Appl. Genet. 108, 1555–1563. doi: 10.1007/s00122-003-1580-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Saito, M., Kubo, N., Matsumoto, S., Suwabe, K., Tsukada, M., and Hirai, M. (2006). Fine mapping of the clubroot resistance gene, Crr3, in Brassica rapa. Theor. Appl. Genet. 114, 81–91. doi: 10.1007/s00122-006-0412-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakamoto, K., Saito, A., Hayashida, N., Taguchi, G., and Matsumoto, E. (2008). Mapping of isolate-specific QTLs for clubroot resistance in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor. Appl. Genet. 117, 759–767. doi: 10.1007/s00122-008-0817-0

PubMed Abstract | CrossRef Full Text | Google Scholar

SAS Institute (1999). SAS/STAT User’s Guide : Version 8. Cary, NC: SAS Institute Inc.

Google Scholar

Suwabe, K., Tsukazaki, H., Iketani, H., Hatakeyama, K., Kondo, M., Fujimura, M., et al. (2006). Simple sequence repeat-based comparative genomics between Brassica rapa and Arabidopsis thaliana: the genetic origin of clubroot resistance. Genetics 173, 309–319. doi: 10.1534/genetics.104.038968

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomita, H., Shimizu, M., Asadud, D. M., Fujimoto, R., and Okazaki, K. (2013). Accumulation of quantitative trait loci conferring broad-spectrum clubroot resistance in Brassica oleracea. Mol. Breed. 32, 889–900. doi: 10.1007/s11032-013-9918-9

CrossRef Full Text | Google Scholar

Voorrips, R. E., Jongerius, M. C., and Kanne, H. J. (1997). Mapping of two genes for resistance to clubroot (Plasmodiophora brassicae) in a population of doubled haploid lines of Brassica oleracea by means of RFLP and AFLP markers. Theor. Appl. Genet. 94, 75–82. doi: 10.1007/s001220050384

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Li, H., Zhang, L., and Meng, L. (2016). User’s Manual of QTL IciMapping ver. 4.1. Mexico: Chinese Academy of Agricultural Sciences.

Wei, L., Jian, H., Lu, K., Filardo, F., Yin, N., Liu, L., et al. (2016). Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus. Plant Biotechnol. J. 14, 1368–1380. doi: 10.1111/pbi.12501

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, F., Zhang, X., Peng, G., Falk, K. C., Strelkov, S. E., and Gossen, B. D. (2017). Genotyping-by-sequencing reveals three QTL for clubroot resistance to six pathotypes of Plasmodiophora brassicae in Brassica rapa. Sci. Rep. 7:4516. doi: 10.1038/s41598-017-04903-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J., Zhao, M., Wang, X., Tong, C., Huang, S., Sadia, T., et al. (2013). Bolbase: a comprehensive genomics database for Brassica oleracea. BMC Genomics 14:664. doi: 10.1186/1471-2164-14-664

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Hua, Y., Wang, X., Zhao, H., Shi, L., and Xu, F. (2014). A high-density genetic map identifies a novel major QTL for boron efficiency in oilseed rape (Brassica napus L.). PLoS One 9:e112089. doi: 10.1371/journal.pone.0112089

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H. L., Xiao, C. G., Chen, G. K., Ma, Y. Y., and Qing, L. (2015). Resistance analysis of Chinese cabbages to Plasmodiophora brassicae based on early infection in roots. J. Plant Prot. 42, 510–516.

Google Scholar

Keywords: Brassica oleracea, Plasmodiophora brassicae, clubroot, SNP microarray, quantitative trait loci

Citation: Peng L, Zhou L, Li Q, Wei D, Ren X, Song H, Mei J, Si J and Qian W (2018) Identification of Quantitative Trait Loci for Clubroot Resistance in Brassica oleracea With the Use of Brassica SNP Microarray. Front. Plant Sci. 9:822. doi: 10.3389/fpls.2018.00822

Received: 02 February 2018; Accepted: 28 May 2018;
Published: 18 June 2018.

Edited by:

Jacqueline Batley, The University of Western Australia, Australia

Reviewed by:

Fengqun Yu, Agriculture and Agri-Food Canada (AAFC), Canada
Ryo Fujimoto, Kobe University, Japan

Copyright © 2018 Peng, Zhou, Li, Wei, Ren, Song, Mei, Si and Qian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jiaqin Mei, Jun Si,