P.brassicae could infect most cruciferous crops and Chinese cabbage is reported highly susceptible to clubroot disease (Dixon, 2009). Mapping of clubroot resistant genes and breeding of clubroot-resistant (CR) cultivars is the most effective and environment-friendly strategy for clubroot management (Pang et al., 2020). Total 26 clubroot resistant genes/loci have been reported from European turnips such as Siloga, ECD1, ECD2, ECD3, ECD4 etc. in Brassica rapa (Suwabe et al., 2003; Chen et al., 2013; Hirani et al., 2018; Wang et al., 2022). However, only CRa and Crr1a have been cloned and these two genes all encode Toll-interleukin-1 receptor-like domain-nucleotide binding site-leucine-rich repeat (TIR-NBS-LRR) proteins (Ueno et al., 2012; Hatakeyama et al., 2013). CR genes appeared to confer isolate-specific resistance to P. brassicae (Chen et al., 2013). In addition, resistant genes encoding TIR-NBS-LRR proteins are known as race-specific resistance in higher plants (Larkan et al. 2020). The pathogenic variability of P. brassicae leading to a complex pathotypes classification (Pang et al., 2020) and causing the breakdown of the resistance of CR cultivars in Chinese cabbage and canola (Kuginuki et al., 1999; Strelkov et al., 2016). Therefore, continues exploring isolate specific clubroot resistance gene are important to the breeding of CR cultivars in Brassica.
In the previously study, we identified a clubroot resistant gene QS_B3.1 which was detected from a European turnip ‘Siloga’ (Pang et al., 2014). QS_B3.1 is a major clubroot resistant locus which explained 70.55% of the phenotypic variation. Then, near-isogenic lines (NILs) of Chinese cabbage containing homozygous QS_B3.1 were selected from BC4F3 population and used for the validation of QS_B3.1 (Ning et al., 2015). In this study, we fine mapped the CR gene QS_B3.1 to an interval of approximately 386 kb between marker syau-InDel3024 and syau-InDel3008 in chromosome A03 using a large F2 population. QS_B3.1 was mapped in an R gene hot zone of A03 where five major CR genes have been reported including CRa, CRb, Rcr1, Rcr2 and Rcr4. CRa has been successfully cloned, encoding TIR-NBS-LRR protein and CRb was identified that is the same gene as CRa (Hatakeyama et al., 2017). While, the other CR genes Rcr1, Rcr2 and Rcr4 require confirmation. Four resistance genes including Bra019409, Bra019410, Bra019412 and Bra019413 which encode the TIR-NBS-LRR proteins were identified from 54 genes in this fine mapped region. Moreover, five unknown genes were detected in the candidate genes. Gene expression data indicates Bra019376, Bra019401, Bra019403 and Bra019410 are contribute to the clubroot resistance (Fig. 2). Sequence variation between the CRa and candidate gene Bra019410 revealed QS_B3.1 is distinct from CRa. Therefore, Bra019376, Bra019401, Bra019403 and Bra019410 are candidates for further gene function study.
The genetic variation within the pathogen and Brassica species is the main challenge in clubroot resistance breeding. Breeding for broaden clubroot resistance cultivars are by far the best effective and economic strategy for dealing with this devastating disease in Brassica crops (Pang et al., 2018). Many CR cultivars have been released in Japan, China, Korea and Canada including Chinese cabbage and canola (Diederichsen et al., 2009; Drury et al., 2021). CRa is a single dominant resistance gene against several pathotypes of P. brassicae (Ueno et al., 2012; Pang et al., 2020). Therefore, CRa was the most widely distributed in CR cultivars of Chinese cabbage in Japan (Aruga et al., 2013). During previous work, CR genes and their closely linked markers exploration have been greatly improved the pyramiding multiple CR genes through marker-assisted selection in CR breading. Different CR genes rotation or pyramiding could confer great resistance to multiple pathotypes of P. brassicae in field (Dolatabadian et al., 2021). In the present work, molecular cnu_ssr316, syau_InDel3008 and sau_um026 showed excellent clubroot resistant individuals selection ability in F2 population which indicates that cnu_ssr316, syau_InDel3008 and sau_um026 are closely linked to clubroot resistance gene. The genetic map shows syau_InDel3008 is the closest molecular market to QS_B3.1. Our research provides a valuable germplasm resource against P. brassicae and a co-segregate molecular marker in future breeding programs for Chinese cabbage and oil seed rape.
In conclusion, the major QTL locus QS_B3.1 conferring to clubroot resistance was fine mapped to an interval of 386 kb on chromosome A3 in B. rapa, and 9 putative candidate genes expression pattern were analyzed. Gene expression results indicated that Bra019376, Bra019401, Bra019403 and Bra019410 are highly possible to be the candidate genes for QS_B3.1. In addition, the candidate gene Bra019410 from ‘CR510’ showed sequence variation from CRa that indicated QS_B3.1 is possible to be a novel resistance gene to clubroot. Further candidate gene’s function analysis will be necessary for QS_B3.1 in brasscia rapa. Our research provides a valuable germplasm resource against P. brassicae and QS_B3.1 closely linked marker syau-InDel3008 for marker-assisted selection in the CR cultivars breeding.