Fine mapping of the clubroot resistance gene, Crr3, in Brassica rapa

Resilience of Canola to Plasmodiophora brassicae (Clubroot) Pathotype 3H under Different Resistance Genes and Initial Inoculum Levels

In this study, we explored the resilience of a clubroot resistance (CR) stacking model against a field population of Plasmodiophora brassicae pathotype 3H. This contrasts with our earlier work, where stacking CRaM and Crr1rutb proved only moderately resistant to pathotype X. Canola varieties carrying Rcr1/Crr1rutb and Rcr1 + Crr1rutb were repeatedly exposed to 3H at low (1 × 104/g soil) and high (1 × 107/g soil) initial resting spore concentrations over five planting cycles under controlled environments to mimic intensive canola production. Initially, all resistant varieties showed strong resistance. However, there was a gradual decline in resistance over time for varieties carrying only a single CR gene, particularly with Crr1rutb alone and at the high inoculum level, where the disease severity index (DSI) increased from 9% to 39% over five planting cycles. This suggests the presence of virulent pathotypes at initially low levels in the 3H inoculum. In contrast, the variety with stacked CR genes remained resilient, with DSI staying below 3% throughout, even at the high inoculum level. Furthermore, the use of resistant varieties, carrying either a single or stacked CR genes, reduced the total resting spore numbers in soil over time, while the inoculum level either increased or remained high in soils where susceptible Westar was continuously grown. Our study demonstrates greater resistance resilience for stacking Rcr1 and Crr1rutb against the field population of 3H. Additionally, the results suggest that resistance may persist even longer in fields with lower levels of inoculum, highlighting the value of extended crop rotation (reducing inoculum) alongside strategic CR-gene deployment to maximize resistance resilience.

Identification of Clubroot (Plasmodiophora brassicae) Resistance Loci in Chinese Cabbage (Brassica rapa ssp. pekinensis) with Recessive Character

The soil-borne pathogen Plasmodiophora brassicae is the causal agent of clubroot, a major disease in Chinese cabbage (Brassica rapa ssp. pekinensis). The host's resistance genes often confer immunity to only specific pathotypes and may be rapidly overcome. Identification of novel clubroot resistance (CR) from germplasm sources is necessary. In this study, Bap246 was tested by being crossed with different highly susceptible B. rapa materials and showed recessive resistance to clubroot. An F2 population derived from Bap246 × Bac1344 was used to locate the resistance Quantitative Trait Loci (QTL) by Bulk Segregant Analysis Sequencing (BSA-Seq) and QTL mapping methods. Two QTL on chromosomes A01 (4.67-6.06 Mb) and A08 (10.42-11.43 Mb) were found and named Cr4Ba1.1 and Cr4Ba8.1, respectively. Fifteen and eleven SNP/InDel markers were used to narrow the target regions in the larger F2 population to 4.67-5.17 Mb (A01) and 10.70-10.84 Mb (A08), with 85 and 19 candidate genes, respectively. The phenotypic variation explained (PVE) of the two QTL were 30.97% and 8.65%, respectively. Combined with gene annotation, mutation site analysis, and real-time quantitative polymerase chain reaction (qRT-PCR) analysis, one candidate gene in A08 was identified, namely Bra020861. And an insertion and deletion (InDel) marker (co-segregated) named Crr1-196 was developed based on the gene sequence. Bra013275, Bra013299, Bra013336, Bra013339, Bra013341, and Bra013357 in A01 were the candidate genes that may confer clubroot resistance in Chinese cabbage. The resistance resource and the developed marker will be helpful in Brassica breeding programs.

Comparative pangenome analyses provide insights into the evolution of Brassica rapa resistance gene analogues (RGAs)

Brassica rapa is grown worldwide as economically important vegetable and oilseed crop. However, its production is challenged by yield-limiting pathogens. The sustainable control of these pathogens mainly relies on the deployment of genetic resistance primarily driven by resistance gene analogues (RGAs). While several studies have identified RGAs in B. rapa, these were mainly based on a single genome reference and do not represent the full range of RGA diversity in B. rapa. In this study, we utilized the B. rapa pangenome, constructed from 71 lines encompassing 12 morphotypes, to describe a comprehensive repertoire of RGAs in B. rapa. We show that 309 RGAs were affected by presence-absence variation (PAV) and 223 RGAs were missing from the reference genome. The transmembrane leucine-rich repeat (TM-LRR) RGA class had more core gene types than variable genes, while the opposite was observed for nucleotide-binding site leucine-rich repeats (NLRs). Comparative analysis with the B. napus pangenome revealed significant RGA conservation (93%) between the two species. We identified 138 candidate RGAs located within known B. rapa disease resistance QTL, of which the majority were under negative selection. Using blackleg gene homologues, we demonstrated how these genes in B. napus were derived from B. rapa. This further clarifies the genetic relationship of these loci, which may be useful in narrowing-down candidate blackleg resistance genes. This study provides a novel genomic resource towards the identification of candidate genes for breeding disease resistance in B. rapa and its relatives.

Analysis of the role of BrRPP1 gene in Chinese cabbage infected by Plasmodiophora brassicae

INTRODUCTION

The clubroot disease caused by Plasmodiophora brassicae (P. brassicae) poses a serious threat to the economic value of cruciferous crops, which is a serious problem to be solved worldwide. Some resistance genes to clubroot disease in Brassica rapa L. ssp pekinensis cause by P. brassicae have been located on different chromosomes. Among them, Rcr1 and Rcr2 were mapped to the common candidate gene Bra019410, but its resistance mechanism is not clear yet.

METHODS

In this experiment, the differences of BrRPP1 between the resistant and susceptible material of Chinese cabbage were analyzed by gene cloning and qRT-PCR. The gene function was verified by Arabidopsis homologous mutants. The expression site of BrRPP1 gene in cells was analyzed by subcellular localization. Finally, the candidate interaction protein of BrRPP1 was screened by yeast two-hybrid library.

RESULTS

The results showed that the cDNA sequence, upstream promoter sequence and expression level of BrRPP1 were quite different between the resistant and susceptible material. The resistance investigation found that the Arabidopsis mutant rpp1 was more susceptible to clubroot disease than the wild type, which suggested that the deletion of rpp1 reduces resistance of plant to clubroot disease. Subcellular location analysis confirmed that BrRPP1 was located in the nucleus. The interaction proteins of BrRPP1 screened from cDNA Yeast Library by yeast two-hybrid are mainly related to photosynthesis, cell wall modification, jasmonic acid signal transduction and programmed cell death.

DISCUSSION

BrRPP1 gene contains TIR-NBS-LRR domain and belongs to R gene. The cDNA and promoter sequence of BrRPP1 in resistant varieties was different from that in susceptible varieties led to the significant difference of the gene expression of BrRPP1 between the resistant varieties and the susceptible varieties. The high expression of BrRPP1 gene in resistant varieties enhanced the resistance of Chinese cabbage to P. brassicae, and the interaction proteins of BrRPP1 are mainly related to photosynthesis, cell wall modification, jasmonic acid signal transduction and programmed cell death. These results provide important clues for understanding the mechanism of BrRPP1 in the resistance of B. rapa to P. brassicae.

Identification of YABBY Transcription Factors and Their Function in ABA and Salinity Response in Nelumbo nucifera

The plant-specific transcription factor family YABBY plays important roles in plant responses to biotic and abiotic stresses. Although the function of YABBY has been identified in many species, systematic analysis in lotus (Nelumbo nucifera) is still relatively lacking. The present study aimed to characterize all of the YABBY genes in lotus and obtain better insights into NnYABBYs in response to salt stress by depending on ABA signaling. Here, we identified nine YABBY genes by searching the whole lotus genome based on the conserved YABBY domain. Further analysis showed that these members were distributed on six different chromosomes and named from YABBY1 to YABBY9, which were divided into five subgroups, including YAB1, YAB2, YAB5, INO, and CRC. The analysis of cis-elements in promotors revealed that NnYABBYs could be involved in plant hormone signaling and plant responses to abiotic stresses. Quantitative real-time PCR (qRT-PCR) showed that NnYABBYs could be up-regulated or down-regulated by ABA, fluridone, and salt treatment. Subcellular localization indicated that NnYABBY4, NnYABBY5, and NnYABBY6 were mainly localized in the cell membrane and cytoplasm. In addition, the intrinsic trans-activity of NnYABBY was tested by a Y2H assay, which revealed that NnYABBY4, NnYABBY5, and NnYABBY6 are deprived of such a property. This study provided a theoretical basis and reference for the functional research of YABBY for the molecular breeding of lotus.

Advances in Biological Control and Resistance Genes of Brassicaceae Clubroot Disease-The Study Case of China

Clubroot disease is a soil-borne disease caused by Plasmodiophora brassicae. It occurs in cruciferous crops exclusively, and causes serious damage to the economic value of cruciferous crops worldwide. Although different measures have been taken to prevent the spread of clubroot disease, the most fundamental and effective way is to explore and use disease-resistance genes to breed resistant varieties. However, the resistance level of plant hosts is influenced both by environment and pathogen race. In this work, we described clubroot disease in terms of discovery and current distribution, life cycle, and race identification systems; in particular, we summarized recent progress on clubroot control methods and breeding practices for resistant cultivars. With the knowledge of these identified resistance loci and R genes, we discussed feasible strategies for disease-resistance breeding in the future.

Early-stage responses to Plasmodiophora brassicae at the transcriptome and metabolome levels in clubroot resistant and susceptible oilseed Brassica napus

Clubroot, a devastating soil-borne root disease, in Brassicaceae is caused by Plasmodiophora brassicae Woronin (P. brassicae W.), an obligate biotrophic protist. Plant growth and development, as well as seed yield of Brassica crops, are severely affected due to this disease. Several reports described the molecular responses of B. napus to P. brassicae; however, information on the early stages of pathogenesis is limited. In this study, we have used transcriptomics and metabolomics to characterize P. brassicae pathogenesis at 1-, 4-, and 7-days post-inoculation (DPI) in clubroot resistant (CR) and susceptible (CS) doubled-haploid (DH) canola lines. When we compared between inoculated and uninoculated groups, a total of 214 and 324 putative genes exhibited differential expression (q-value < 0.05) at one or more time-points in the CR and CS genotypes, respectively. When the inoculated CR and inoculated CS genotypes were compared, 4765 DEGs were differentially expressed (q-value < 0.05) at one or more time-points. Several metabolites related to organic acids (e.g., citrate, pyruvate), amino acids (e.g., proline, aspartate), sugars, and mannitol, were differentially accumulated in roots in response to pathogen infection when the CR and CS genotypes were compared. Several DEGs also corresponded to differentially accumulated metabolites, including pyrroline-5-carboxylate reductase (BnaC04g11450D), citrate synthase (BnaC02g39080D), and pyruvate kinase (BnaC04g23180D) as detected by transcriptome analysis. Our results suggest important roles for these genes in mediating resistance to clubroot disease. To our knowledge, this is the first report of an integrated transcriptome and metabolome analysis aimed at characterizing the molecular basis of resistance to clubroot in canola.

Mapping of a novel clubroot disease resistance locus in Brassica napus and related functional identification

Clubroot disease, caused by Plasmodiophora brassicae, is a devastating disease that results in substantial yield loss in Brassicaceae crops worldwide. In this study, we identified a clubroot disease resistance (CR) Brassica napus, "Kc84R," which was obtained by mutation breeding. Genetic analysis revealed that the CR trait of "Kc84R" was controlled by a single dominant locus. We used the bulked segregant analysis sequencing (BSA-seq) approach, combined with genetic mapping based on single nucleotide polymorphism (SNP) markers to identify CR loci from the F₂ population derived from crossing CR "Kc84R" and clubroot susceptible "855S." The CR locus was mapped to a region between markers BnSNP14198336 and BnSNP14462201 on the A03 chromosome, and this fragment of 267 kb contained 68 annotated candidate genes. Furthermore, we performed the CR relation screening of candidate genes with the model species Arabidopsis. An ERF family transcriptional activator, BnERF034, was identified to be associated with the CR, and the corresponding Arabidopsis homozygous knockout mutants exhibited more pronounced resistance compared with the wild-type Col-0 and the transgenic lines of BnERF034 in response to P. brassicae infection. Additionally, the expression analysis between resistant and susceptible materials indicated that BnERF034 was identified to be the most likely CR candidate for the resistance in Kc84R. To conclude, this study reveals a novel gene responsible for CR. Further analysis of BnERF034 may reveal the molecular mechanisms underlying the CR of plants and provide a theoretical basis for Brassicaceae resistance breeding.

Germplasm Enhancement and Identification of Loci Conferring Resistance against Plasmodiophora brassicae in Broccoli

In order to breed broccoli and other Brassica materials to be highly resistant to clubroot disease, 41 Brassicaceae varieties were developed and identified between 2020 and 2021. Seven known clubroot genes were used for screening these materials. In addition, the resistant and susceptible broccoli cultivars were designed for observing their differences in the infection process with Plasmodiophora brassicae. The results showed that 90% of total materials had carried more than two clubroot resistance genes: one material carried two disease resistance genes, four materials carried seven genes for clubroot resistance, two materials carried six genes for clubroot resistance, and in total 32% of these materials carried five genes for clubroot resistance. As a result, several new genotypes of Brassicaceae germplasm were firstly created and obtained based on distant hybridization and identification of loci conferring resistance against Plasmodiophora brassicae in this study. We found and revealed that similar infection models of Plasmodiophora brassicae occurred in susceptible and resistant cultivars of broccoli, but differences in infection efficiency of Plasmodiophora brassicae also existed in both materials. For resistant broccoli plants, a small number of conidia formed in the root hair, and only a few spores could enter the cortex without forming sporangia while sporangia could form in susceptible plants. Our study could provide critical Brassica materials for breeding resistant varieties and new insight into understanding the mechanism of plant resistance.

Introgression of clubroot resistant gene into L. from based on homoeologous exchange

Clubroot is a soil-borne disease in cabbage (Brassica oleracea L. var. capitata L.) caused by Plasmodiophora brassicae, which poses a great threat to cabbage production. However, clubroot resistance (CR) genes in Brassica rapa could be introduced into the cabbage via breeding to make it clubroot resistant. In this study, CR genes from B. rapa were introduced into the cabbage genome and the mechanism of gene introgression was explored. Two methods were used to create CR materials: (i) The fertility of CR Ogura CMS cabbage germplasms containing CRa was restored by using an Ogura CMS restorer. After cytoplasmic replacement and microspore culture, CRa-positive microspore individuals were obtained. (ii) Distant hybridization was performed between cabbage and B. rapa, which contained three CR genes (CRa, CRb, and Pb8.1). Finally, BC2 individuals containing all three CR genes were obtained. Inoculation results showed that both CRa-positive microspore individuals and BC2 individuals containing three CR genes were resistant to race 4 of P. brassicae. Sequencing results from CRa-positive microspore individuals with specific molecular markers and genome-wide association study (GWAS) showed penetration at the homologous position of the cabbage genome by a 3.42 Mb CRa containing a fragment from B. rapa; indicating homoeologous exchange (HE) as the theoretical basis for the introgression of CR resistance. The successful introduction of CR into the cabbage genome in the present study can provide useful clues for creating introgression lines within other species of interest.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

Identification of Novel Locus RsCr6 Related to Clubroot Resistance in Radish (Raphanus sativus L.)

Clubroot is a devastating disease that causes substantial yield loss worldwide. However, the inheritance and molecular mechanisms of clubroot resistance during pathogen infection in radish remain largely unclear. In this study, we investigated the inheritance of clubroot resistance in the F₂ population derived from crossing clubroot-resistant (CR) and clubroot-susceptible inbred lines "GLX" and "XNQ," respectively. Genetic analysis revealed that a single dominant gene controlled the clubroot resistance of "GLX" with a Mendelian ratio of resistance and susceptibility of nearly 3:1. Bulked segregant analysis combined with whole-genome resequencing (BSA-seq) was performed to detect the target region of RsCr6 on chromosome Rs8. Linkage analysis revealed that the RsCr6 locus was located between two markers, HB321 and HB331, with an interval of approximately 92 kb. Based on the outcomes of transcriptome analysis, in the RsCr6 locus, the R120263140 and R120263070 genes with a possible relation to clubroot resistance were considered candidate genes. In addition, three core breeding materials containing the two reported quantitative trait loci (QTLs) and our novel locus RsCr6 targeting clubroot resistance were obtained using marker-assisted selection (MAS) technology. This study reveals a novel locus responsible for clubroot resistance in radishes. Further analysis of new genes may reveal the molecular mechanisms underlying the clubroot resistance of plants and provide a theoretical basis for radish resistance breeding.

Allelic variation of a clubroot resistance gene (Crr1a) in Japanese cultivars of Chinese cabbage (Brassica rapa L.)

Clubroot resistance (CR) is an important trait in Chinese cabbage breeding worldwide. Although Crr1a, the gene responsible for clubroot-resistance, has been cloned and shown to encode the NLR protein, its allelic variation and molecular function remain unknown. Here, we investigated the sequence variation and function of three Crr1a alleles cloned from six CR F₁ cultivars of Chinese cabbage. Gain-of-function analysis revealed that Crr1a^Kinami90_a isolated from the cv. 'Kinami 90' conferred clubroot resistance as observed for Crr1a^G004 . Because two susceptible alleles commonly lacked 172 amino acids in the C-terminal region, we investigated clubroot resistance in transgenic Arabidopsis harboring the chimeric Crr1a, in which 172 amino acids of the functional alleles were fused to the susceptible alleles. The fusion of the C-terminal region to the susceptible alleles restored resistance, indicating that their susceptibility was caused by the lack of the C-terminus. We developed DNA markers to detect the two functional Crr1a alleles, and demonstrated that the functional Crr1a alleles were frequently found in European fodder turnips, whereas they were rarely introduced into Japanese CR cultivars of Chinese cabbage. These results would contribute to CR breeding via marker-assisted selection and help our understanding of the molecular mechanisms underlying clubroot resistance.

Construction of an Intragenic SSR-Based Linkage Map and QTL Mapping for Agronomic Traits in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Chinese cabbage (Brassica rapa L. ssp. pekinensis) is one of the most widely cultivated and economically important vegetables in China. Constructing an effective genetic linkage map and mapping quantitative trait loci (QTLs) related to yield and leafy head morphology is of great importance for molecular breeding of Chinese cabbage. Using two diverse Chinese cabbage inbred lines, ZHB and G291, as parents, an F2 segregating population consisting of 240 individuals was prepared for genetic map construction and phenotype investigation in this study. The two parents are significantly different in both shape and size. Sixteen important agronomic traits of F2 individuals were investigated. A genetic map of 105 intragenic simple sequence repeat (SSR) markers distributed across 10 linkage groups (LGs) was constructed, which was 2034.1 cM in length and had an average inter-locus distance of 21.75 cM. We identified 48 QTLs for the tested important agronomic traits on the studied LGs, with LOD scores of 2.51–12.49, which explained the phenotypic variance of 3.41–26.66%. The QTLs identified in this study will facilitate further genetic analysis and marker-assisted genetic improvement of Chinese cabbage.

Marker-Assisted Pyramiding of Genes for Multilocular Ovaries, Self-Compatibility, and Clubroot Resistance in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Molecular marker-assisted gene pyramiding combined with backcrossing has been widely applied for crop variety improvement. Molecular marker identification could be used in the early stage of breeding to achieve the rapid and effective pyramiding of multiple genes. To create high-quality germplasm for Chinese cabbage breeding, multi-gene pyramiding for self-compatibility, multilocular, and clubroot resistance was performed through molecular marker-assisted selection. The results showed that self-compatibility and multilocular traits were controlled by a pair of recessive genes. Two flanking markers, sau_um190 and cun_246a, and marker Teo-1, based on the gene sequence related to multilocular ovaries, were used for multilocular ovary trait selection. Two flanking markers, SCF-6 and SC-12, and marker Sal-SLGI /PK1+PK4, based on the gene sequence, were used for self-compatibility selection. Two flanking markers, TCR74 and TCR79, closely linked to clubroot resistance gene CRb, were used as foreground selection markers. Based on Chinese cabbage genomic information, 111 SSR markers covering 10 chromosomes were applied for background selection. After multiple generations of selection, a multi-gene pyramided line from a BC4F2 population with self-compatibility, multilocular ovaries, and clubroot resistance was obtained with a high genomic background recovery rate. The improved pyramided line is expected to be utilized as a potential material in further breeding programs.

Histopathology of the Plasmodiophora brassicae-Chinese Cabbage Interaction in Hosts Carrying Different Sources of Resistance

Clubroot is a serious soil-borne disease of crucifers caused by the obligate parasite Plasmodiophora brassicae. The genetic basis and histopathology of clubroot resistance in two Chinese cabbage (Brassica rapa ssp. pekinensis) inbred lines Bap055 and Bap246, challenged with pathotype 4 of P. brassicae, was evaluated. The Chinese cabbage cultivar "Juxin" served as a susceptible check. The resistance in Bap055 was found to be controlled by the CRa gene, while resistance in Bap246 fit a model of control by unknown recessive gene. Infection of the roots by P. brassicae was examined by inverted microscopy. Despite their resistance, primary and secondary infection were observed to occur in Bap055 and Bap246. Primary infection was detected at 2 days post-inoculation (DPI) in "Juxin," at 4 DPI in Bap055, and at 6 DPI in Bap246. Infection occurred most quickly on "Juxin," with 60% of the root hairs infected at 10 DPI, followed by Bap055 (31% of the root hairs infected at 12 DPI) and Bap246 (20% of the root hairs infected at 14 DPI). Secondary infection of "Juxin" was first observed at 8 DPI, while in Bap055 and Bap246, secondary infection was first observed at 10 DPI. At 14 DPI, the percentage of cortical infection in "Juxin," Bap055 and Bap246 was 93.3, 20.0, and 11.1%, respectively. Although cortical infection was more widespread in Bap055 than in Bap246, secondary infection in both of these hosts was restricted relative to the susceptible check, and the vascular system remained intact. A large number of binucleate secondary plasmodia were observed in "Juxin" and the vascular system was disrupted at 16 DPI; in Bap055 and Bap246, only a few secondary plasmodia were visible, with no binucleate secondary plasmodia. The defense mechanisms and expression of resistance appears to differ between Chinese cabbage cultivars carrying different sources of resistance.

Identification and Fine-Mapping of Clubroot (Plasmodiophora brassicae) Resistant QTL in Brassica rapa

European fodder turnips (Brassica rapa ssp. rapifera) were identified as sources of clubroot resistance (CR) and have been widely used in Brassica resistance breeding. An F2 population derived from a cross between a resistant turnip and a susceptible Chinese cabbage was used to determine the inheritance and locating the resistance Quantitative Trait Loci (QTLs). The parents showed to be very resistant/susceptible to the field isolates (pathotype 4) of clubroot from Henan in China. After inoculation, 27 very resistant or susceptible individuals were selected to construct bulks, respectively. Next-generation-sequencing-based Bulk Segregant Analysis Sequencing (BSA-Seq) was used and located resistance QTL on chromosome A03 (3.3–7.5 Mb) and A08 (0.01–6.5 Mb), named Bcr1 and Bcr2, respectively. Furthermore, an F3 population including 180 families derived from F2 individuals was phenotyped and used to verify and narrow candidate regions. Ten and seven Kompetitive Allele-Specific PCR (KASP) markers narrowed the target regions to 4.3–4.78 Mb (A03) and 0.02–0.79 Mb (A08), respectively. The phenotypic variation explained (PVE) of the two QTLs were 33.3% and 13.3% respectively. The two candidate regions contained 99 and 109 genes. In the A03 candidate region, there were three candidate R genes, namely Bra006630, Bra006631 and Bra006632. In the A08 candidate region, there were two candidate R genes, namely Bra030815 and Bra030846.

A Proteome-Level Investigation Into Plasmodiophora brassicae Resistance in Brassica napus Canola

Clubroot of Brassicaceae, an economically important soil borne disease, is caused by Plasmodiophora brassicae Woronin, an obligate, biotrophic protist. This disease poses a serious threat to canola and related crops in Canada and around the globe causing significant losses. The pathogen is continuously evolving and new pathotypes are emerging, which necessitates the development of novel resistant canola cultivars to manage the disease. Proteins play a crucial role in many biological functions and the identification of differentially abundant proteins (DAP) using proteomics is a suitable approach to understand plant-pathogen interactions to assist in the development of gene specific markers for developing clubroot resistant (CR) cultivars. In this study, P. brassicae pathotype 3 (P3H) was used to challenge CR and clubroot susceptible (CS) canola lines. Root samples were collected at three distinct stages of pathogenesis, 7-, 14-, and 21-days post inoculation (DPI), protein samples were isolated, digested with trypsin and subjected to liquid chromatography with tandem mass spectrometry (LC-MS/MS) analysis. A total of 937 proteins demonstrated a significant (q-value < 0.05) change in abundance in at least in one of the time points when compared between control and inoculated CR-parent, CR-progeny, CS-parent, CS-progeny and 784 proteins were significantly (q < 0.05) changed in abundance in at least in one of the time points when compared between the inoculated- CR and CS root proteomes of parent and progeny across the three time points tested. Functional annotation of differentially abundant proteins (DAPs) revealed several proteins related to calcium dependent signaling pathways. In addition, proteins related to reactive oxygen species (ROS) biochemistry, dehydrins, lignin, thaumatin, and phytohormones were identified. Among the DAPs, 73 putative proteins orthologous to CR proteins and quantitative trait loci (QTL) associated with eight CR loci in different chromosomes including chromosomes A3 and A8 were identified. Proteins including BnaA02T0335400WE, BnaA03T0374600WE, BnaA03T0262200WE, and BnaA03T0464700WE are orthologous to identified CR loci with possible roles in mediating clubroot responses. In conclusion, these results have contributed to an improved understanding of the mechanisms involved in mediating response to P. brassicae in canola at the protein level.

QTL mapping for Fusarium wilt resistance based on the whole-genome resequencing and their association with functional genes in Raphanus sativus

Two major QTL associated with resistance to Fusarium wilt (FW) were identified using whole-genome resequencing. Sequence variations and gene expression level differences suggest that TIR-NBS and LRR-RLK are candidate genes associated with FW-resistance. Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. raphani is an important disease in radish, leading to severe decrease in yield and quality. YR4 as a novel genetic source to resistant to FW was confirmed through screening with five pathogen isolates. We have generated F₂ and F_2:3 populations segregated with FW resistance using YR4 and YR18 inbred lines. The disease symptom was evaluated in F_2:3 population (n = 180) in three independent studies over two years. We identified 4 QTL including the two major QTL (FoRsR7.1^59A and FoRsR9.3^59A). FoRsR7.1^59A and FoRsR9.3^59A were detected in three replicated experiments. FoRsR7.1^59A was delimited to the 2.18-Mb physical interval on chromosome R07, with a high LOD value (5.17-12.84) and explained phenotypic variation (9.34%-27.97%). The FoRsR9.3^59A represented relatively low LOD value (3.38-4.52) and explained phenotypic variation (6.24%-8.82%). On the basis of the re-sequencing data for the parental lines, we identified five putative resistance-related genes and 13 unknown genes with sequence variations at the gene and protein levels. A semi-quantitative RT-PCR analysis revealed that Rs382940 (TIR-NBS) and Rs382200 (RLK) were expressed only in 'YR4' from 0 to 6 days after the inoculation. Moreover, Rs382950 (TIR-NBS-LRR) was more highly expressed in 'YR4' from 3 to 6 days after the inoculation. These three genes might be important for FW-resistance in radish. We identified several markers based on these potential candidate genes. The marker set should be useful for breeding system to introduce the FW resistance loci from 'YR4' to improve tolerance to FW.

Clubroot resistance derived from the European Brassica napus cv. 'Tosca' is not effective against virulent Plasmodiophora brassicae isolates from Alberta, Canada

In this study, clubroot resistance in the resynthesized European winter Brassica napus cv. 'Tosca' was introgressed into a Canadian spring canola line '11SR0099', which was then crossed with the clubroot susceptible spring line '12DH0001' to produce F₁ seeds. The F₁ plants were used to develop a doubled haploid (DH) mapping population. The parents and the DH lines were screened against 'old' pathotypes 2F, 3H, 5I, 6M and 8N of the clubroot pathogen, Plasmodiophora brassicae, as well as against the 'new' pathotypes 5X, 5L, 2B, 3A, 3D, 5G, 8E, 5C, 8J, 5K, 3O and 8P. Genotyping was conducted using a Brassica 15K SNP array. The clubroot screening showed that 'Tosca, '11SR0099' and the resistant DH lines were resistant to three (2F, 3H and 5I) of the five 'old' pathotypes and four (2B, 3O, 8E and 8P) of the 12 'new' pathotypes, while being moderately resistant to the 'old' pathotype 8N and the 'new' pathotypes 3D and 5G. 'Tosca' was susceptible to isolates representing pathotype 3A (the most common among the 'new' pathotypes) as well as pathotypes 6M, 5X, 5L, 5K and 8J. Linkage analysis and QTL mapping identified a ca. 0.88-0.95 Mb genomic region on the A03 chromosome of 'Tosca' as conferring resistance to pathotypes 2F, 3H, 5I, 2B, 3D, 5G, 8E, 3O and 8P. The identified QTL genomic region housed the CRk, Crr3 and CRd gene(s). However, the susceptibility of 'Tosca' to most of the common virulent pathotypes makes it unattractive as a sole CR donor in the breeding of commercial canola varieties in western Canada.

Local Duplication of TIR-NBS-LRR Gene Marks Clubroot Resistance in Brassica napus cv. Tosca

Clubroot, caused by Plasmodiophora brassicae infection, is a disease of growing importance in cruciferous crops, including oilseed rape (Brassica napus). The affected plants exhibit prominent galling of the roots that impairs their capacity for water and nutrient uptake, which leads to growth retardation, wilting, premature ripening, or death. Due to the scarcity of effective means of protection against the pathogen, breeding of resistant varieties remains a crucial component of disease control measures. The key aspect of the breeding process is the identification of genetic factors associated with variable response to the pathogen exposure. Although numerous clubroot resistance loci have been described in Brassica crops, continuous updates on the sources of resistance are necessary. Many of the resistance genes are pathotype-specific, moreover, resistance breakdowns have been reported. In this study, we characterize the clubroot resistance locus in the winter oilseed rape cultivar "Tosca." In a series of greenhouse experiments, we evaluate the disease severity of P. brassicae-challenged "Tosca"-derived population of doubled haploids, which we genotype with Brassica 60 K array and a selection of SSR/SCAR markers. We then construct a genetic map and narrow down the resistance locus to the 0.4 cM fragment on the A03 chromosome, corresponding to the region previously described as Crr3. Using Oxford Nanopore long-read genome resequencing and RNA-seq we review the composition of the locus and describe a duplication of TIR-NBS-LRR gene. Further, we explore the transcriptomic differences of the local genes between the clubroot resistant and susceptible, inoculated and control DH lines. We conclude that the duplicated TNL gene is a promising candidate for the resistance factor. This study provides valuable resources for clubroot resistance breeding programs and lays a foundation for further functional studies on clubroot resistance.

Identification of resistance loci against new pathotypes of Plasmodiophora brassicae in Brassica napus based on genome-wide association mapping

Genetic resistance is a successful strategy for management of clubroot (Plasmodiophora brassicae) of brassica crops, but resistance can break down quickly. Identification of novel sources of resistance is especially important when new pathotypes arise. In the current study, the reaction of 177 accessions of Brassica napus to four new, virulent pathotypes of P. brassicae was assessed. Each accession was genotyped using genotyping by sequencing to identify and map novel sources of clubroot resistance using mixed linear model (MLM) analysis. The majority of accessions were highly susceptible (70–100 DSI), but a few accessions exhibited strong resistance (0–20 DSI) to pathotypes 5X (21 accessions), 3A (8), 2B (7), and 3D (15), based on the Canadian Clubroot Differential system. In total, 301,753 SNPs were mapped to 19 chromosomes. Population structure analysis indicated that the 177 accessions belong to seven major populations. SNPs were associated with resistance to each pathotype using MLM. In total, 13 important SNP loci were identified, with 9 SNPs mapped to the A-genome and 4 to the C-genome. The SNPs were associated with resistance to pathotypes 5X (2 SNPs), 3A (4), 2B (5) and 3D (6). A Blast search of 1.6 Mb upstream and downstream from each SNP identified 13 disease-resistance genes or domains. The distance between a SNP locus and the nearest resistance gene ranged from 0.04 to 0.74 Mb. The resistant lines and SNP markers identified in this study can be used to breed for resistance to the most prevalent new pathotypes of P. brassicae in Canada.

Clubroot in Brassica: recent advances in genomics, breeding, and disease management

Clubroot disease, caused by Plasmodiophora brassicae, affects Brassica oilseed and vegetable production worldwide. This review is focused on various aspects of clubroot disease and its management, including understanding the pathogen and resistance in the host plants. Advances in genetics, molecular biology techniques, and omics research have helped to identify several major loci, QTL, and genes from the Brassica genomes involved in the control of clubroot resistance. Transcriptomic studies have helped to extend our understanding of the mechanism of infection by the pathogen and the molecular basis of resistance/susceptibility in the host plants. A comprehensive understanding of the clubroot disease and host resistance would allow developing a better strategy by integrating the genetic resistance with cultural practices to manage this disease from a long-term perspective.

Clubroot Symptoms and Resting Spore Production in a Doubled Haploid Population of Oilseed Rape (Brassica napus) Are Controlled by Four Main QTLs

Clubroot, caused by Plasmodiophora brassicae Woronin, is one of the most important diseases of oilseed rape (Brassica napus L.). The rapid erosion of monogenic resistance in clubroot-resistant (CR) varieties underscores the need to diversify resistance sources controlling disease severity and traits related to pathogen fitness, such as resting spore production. The genetic control of disease index (DI) and resting spores per plant (RSP) was evaluated in a doubled haploid (DH) population consisting of 114 winter oilseed rape lines, obtained from the cross 'Aviso' × 'Montego,' inoculated with P. brassicae isolate "eH." Linkage analysis allowed the identification of three quantitative trait loci (QTLs) controlling DI (PbBn_di_A02, PbBn_di_A04, and PbBn_di_C03). A significant decrease in DI was observed when combining effects of the three resistance alleles at these QTLs. Only one QTL, PbBn_rsp_C03, was found to control RSP, reducing resting spore production by 40%. PbBn_rsp_C03 partially overlapped with PbBn_di_C03 in a nucleotide-binding leucine-rich repeat (NLR) gene-containing region. Consideration of both DI and RSP in breeding for clubroot resistance is recommended for the long-term management of this disease.

Two Clubroot-Resistance Genes, Rcr3 and Rcr9wa, Mapped in Brassica rapa Using Bulk Segregant RNA Sequencing

Genetic resistance is widely used to manage clubroot (Plasmodiophora brassicae) in brassica crops, but new pathotypes have recently been identified on canola (Brassica napus) on the Canadian prairies. Resistance effective against both the most prevalent pathotype (3H, based on the Canadian Clubroot Differential system) and the new pathotypes is needed. BC₁ plants of Brassica rapa from a cross of line 96-6990-2 (clubroot resistance originating from turnip cultivar ‘Waaslander’) and a susceptible doubled-haploid line, ACDC, exhibited a 1:1 segregation for resistance against pathotypes 3H and 5X. A resistance gene designated as Rcr3 was mapped initially based on the percentage of polymorphic variants using bulked segregant RNA sequencing (BSR-Seq) and further mapped using Kompetitive Allele Specific PCR. DNA variants were identified by assembling short reads against a reference genome of B. rapa. Rcr3 was mapped into chromosome A08. It was flanked by single nucleotide polymorphisms (SNP) markers (A90_A08_SNP_M12 and M16) between 10.00 and 10.23 Mb, in an interval of 231.6 Kb. There were 32 genes in the Rcr3 interval. Three genes (Bra020951, Bra020974, and Bra020979) were annotated with disease resistance mechanisms, which are potential candidates for Rcr3. Another resistance gene, designated as Rcr9^wa, for resistance to pathotype 5X was mapped, with the flanking markers (A90_A08_SNP_M28 and M79) between 10.85 and 11.17 Mb using the SNP sites identified through BSR-Seq for Rcr3. There were 44 genes in the Rcr9^wa interval, three of which (Bra020827, Bra020828, Bra020814) were annotated as immune-system-process related genes, which are potential candidates for Rcr9^wa.

QTL Mapping and Inheritance of Clubroot Resistance Genes Derived From Brassica rapa subsp. rapifera (ECD 02) Reveals Resistance Loci and Distorted Segregation Ratios in Two F 2 Populations of Different Crosses

In this study, Brassica rapa subsp. rapifera (ECD 02) which exhibits broad-spectrum resistance to many Canadian Plasmodiophora brassicae isolates was crossed with two clubroot-susceptible B. rapa accessions to produce two F ₂ populations. The F ₂ plants were screened against P. brassicae pathotypes 3H, 5X, and 5G. The Chi-square goodness of fit test showed that the vast majority (≈75%) of the crosses that produced the F ₂ populations showed segregation ratios of 9R:7S, 7R:9S, 13R:3S, 3R:13S, 5R:11S, 11R:5S, and 1R:15S. These were modifications of the 15R:1S ratio expected for the inheritance of two dominant major clubroot resistance (CR) genes from ECD 02. The distorted segregation ratios suggest that the two resistance genes are on different chromosomes and that two genes interact in an epistatic manner to confer resistance. Genotyping was conducted with 144 PCR-based markers in the two F ₂ populations. Linkage and QTL analysis with the polymorphic markers identified two QTLs on chromosome A03 to be associated with resistance to P. brassicae pathotypes 5X and 5G in Popl#1 while only the second QTL on chromosome A03 was associated with resistance to pathotypes 5X and 5G in Popl#2. The QTLs clustered in genomic regions on the A03 chromosome of B. rapa where the CRa/CRb ^Kato gene(s) are mapped. In addition, the Crr1 gene on the A08 chromosome of B. rapa was detected in the two F ₂ populations. Therefore, the phenotypic and molecular data confirm the existence of two CR genes in ECD 02. This is the first study that shows that major dominant genes in Brassica interact in a non-additive manner to confer resistance to different P. brassicae pathotypes. Key Message: This study provides knowledge on the inheritance and type of gene action for clubroot resistance derived from Brassica rapa subsp. rapifera (ECD 02). The results indicated that duplicate recessive and recessive suppression epistatic interactions, digenic additivity and complementary gene action between the CRa/CRb ^Kato gene(s) on the A03 and the Crr1 gene on the A08 chromosome of B. rapa controlled clubroot resistance to P. brassicae pathotypes 3H, 5X and 5G.

Quantitative Trait Locus Mapping of Clubroot Resistance and Plasmodiophora brassicae Pathotype Banglim-Specific Marker Development in Brassica rapa

Clubroot resistance is an economically important trait in Brassicaceae crops. Although many quantitative trait loci (QTLs) for clubroot resistance have been identified in Brassica, disease-related damage continues to occur owing to differences in host variety and constant pathogen variation. Here, we investigated the inheritance of clubroot resistance in a double haploid population developed by crossing clubroot resistant and susceptible lines "09CR500" and "09CR501", respectively. The resistance of "09CR500" to Plasmodiophora brassicae pathotype "Banglim" was controlled as a single dominant gene, with the segregation of resistance and susceptibility being nearly 1:1. PbBrA08^Banglim was identified as having a logarithm of odds value of 7.9-74.8, and a phenotypic variance of 26.0-97.1% with flanking marker "09CR.11390652" in A08. After aligning QTL regions to the B. rapa reference genome, 11 genes were selected as candidates. PbBrA08^Banglim was located near Crr1, CRs, and Rcr9 loci, but differences were validated by marker analysis, gene structural variations, and gene expression levels, as well as phenotypic responses to the pathotype. Genotyping using the "09CR.11390652" marker accurately distinguished the Banglim-resistance phenotypes in the double haploid population. Thus, the developed marker will be useful in Brassica breeding programs, marker-assisted selection, and gene pyramiding to identify and develop resistant cultivars.

Genetics of Clubroot and Fusarium Wilt Disease Resistance in Brassica Vegetables: The Application of Marker Assisted Breeding for Disease Resistance

The genus Brassica contains important vegetable crops, which serve as a source of oil seed, condiments, and forages. However, their production is hampered by various diseases such as clubroot and Fusarium wilt, especially in Brassica vegetables. Soil-borne diseases are difficult to manage by traditional methods. Host resistance is an important tool for minimizing disease and many types of resistance (R) genes have been identified. More than 20 major clubroot (CR) disease-related loci have been identified in Brassica vegetables and several CR-resistant genes have been isolated by map-based cloning. Fusarium wilt resistant genes in Brassica vegetables have also been isolated. These isolated R genes encode the toll-interleukin-1 receptor/nucleotide-binding site/leucine-rice-repeat (TIR-NBS-LRR) protein. DNA markers that are linked with disease resistance allele have been successfully applied to improve disease resistance through marker-assisted selection (MAS). In this review, we focused on the recent status of identifying clubroot and Fusarium wilt R genes and the feasibility of using MAS for developing disease resistance cultivars in Brassica vegetables.

Utilization of Ogura CMS germplasm with the clubroot resistance gene by fertility restoration and cytoplasm replacement in Brassica oleracea L

Clubroot disease, a major plant root disease caused by Plasmodiophora brassicae, has become one of the most destructive diseases among cultivated cruciferous vegetables. However, clubroot-resistant Brassica oleracea materials are rare. A few clubroot-resistant cabbage varieties are available on the market, but all are Ogura cytoplasmic male sterile (CMS) types. Therefore, in this study, to reutilize the clubroot-resistant Ogura CMS germplasm of cabbage, a new fertility-restored Ogura CMS material, 16Q2-11, was used as a bridge to transfer the clubroot resistance (CR) gene from the Ogura CMS cytoplasm to the normal cytoplasm by a two-step method (a fertility restoration and cytoplasm replacement method). In the first cross for fertility restoration of Ogura CMS clubroot-resistant cabbage (FRCRC), 16Q2-11 was used as a restorer to cross with Ogura CMS materials containing the CR gene CRb2. Eleven Rfo-positive progenies were generated, of which four contained CRb2: F8-514, F8-620, F8-732 and F8-839. After inoculation with race 4 of P. brassicae, these four CRb2-positive individuals showed resistance. Furthermore, F8-514 and F8-839 were then used as male parents in the second cross of FRCRC to cross with cabbage inbred lines, resulting in the successful introgression of the CRb2 gene into the inbred lines. All offspring produced from this step of cross, which had a normal cytoplasm, showed a high resistance to race 4 of P. brassicae and could be utilized for the breeding of clubroot-resistant cabbage varieties in the future. This is the first time that the Ogura CMS restorer has been used to restore the fertility of Ogura CMS clubroot-resistant cabbages, which could improve germplasm diversity in cabbage and provide a reference method for using CMS germplasm in Brassica crops.

An update on the arsenal: mining resistance genes for disease management of Brassica crops in the genomic era

Brassica species include many economically important crops that provide nutrition and health-promoting substances to humans worldwide. However, as with all crops, their production is constantly threatened by emerging viral, bacterial, and fungal diseases, whose incidence has increased in recent years. Traditional methods of control are often costly, present limited effectiveness, and cause environmental damage; instead, the ideal approach is to mine and utilize the resistance genes of the Brassica crop hosts themselves. Fortunately, the development of genomics, molecular genetics, and biological techniques enables us to rapidly discover and apply resistance (R) genes. Herein, the R genes identified in Brassica crops are summarized, including their mapping and cloning, possible molecular mechanisms, and application in resistance breeding. Future perspectives concerning how to accurately discover additional R gene resources and efficiently utilize these genes in the genomic era are also discussed.

Development of a Sinitic Clubroot Differential Set for the Pathotype Classification of Plasmodiophora brassicae

Plasmodiophora brassicae, which is known for its broad genetic diversity for virulence, is the causal agent of clubroot disease of Brassica crops worldwide. Studies on pathotype characterization with four differential hosts according to Williams' classification system showed the predominance of pathotype 4 in China. However, the genetic variability within pathotype 4 complicates the breeding of durable clubroot-resistant (CR) cultivars. Herein, a Sinitic clubroot differential (SCD) set was developed using a set of eight differential inbred lines of Chinese cabbage with known or novel CR genes. The presence of immense diversity within pathotype 4 of Williams' system was verified, and 11 pathotypes were characterized using the developed SCD system. The scalability and practicability of the system was further confirmed with a subset of 95 field isolates from different Brassica crops and different regions of China and Korea. Sixteen pathotypes were detected from 132 field isolates, named Pb1 to Pb16, respectively. Among them, Pb1 and Pb4 were prevalent in diverse Brassica crops in the southern and northern regions of China. Pb12, Pb13, Pb14, and Pb16 showed area-specific distribution. The SCD set developed herein will provide important genetic resources for pathogenicity studies of P. brassicae and for CR breeding in Chinese cabbage and other Brassica crops.

Genome-Wide Mapping of Loci Associated With Resistance to Clubroot in Brassica napus ssp. napobrassica (Rutabaga) Accessions From Nordic Countries

Rutabaga [Brassica napus ssp. napobrassica (L.) Hanelt] is reported to be an excellent source of clubroot (Plasmodiophora brassicae) resistance genes. In this study, 124 rutabaga accessions from the Nordic countries (Norway, Sweden, Finland, Denmark, and Iceland) were evaluated for their reaction to five single-spore isolates representing P. brassicae pathotypes 2F, 3H, 5I, 6M, and 8N and 12 field isolates representing pathotypes 2B, 3A, 3O, 5C, 5G, 5K, 5L, 5X (two isolates, L-G2 and L-G3), 8E, 8J, and 8P. The accessions were also genotyped using a 15K Brassica SNP array and 60 PCR-based primers linked to previously identified clubroot resistance genes. Six thousand eight hundred sixty-one SNP markers were retained after filtering with TASSEL 5.0, and used to evaluate four general linear models (GLM) and four mixed linear models (MLM). The PCA + K and Q + K MLM models gave the minimal deviance of the observed from the expected distribution in quantile-quantile plots, and hence were used for SNP-clubroot association analyses. In addition, 108 alleles derived from the PCR-based markers and the phenotypic data were analyzed with the PCA + K model. Forty-five SNPs and four PCR-based markers were identified to be associated strongly with resistance to isolates representing 13 pathotypes (2F, 3H, 5I, 6M, 8N, 2B, 3A, 3O, 5C, 5G, 5K, 5L, and 8P). These markers revealed the top and bottom segments of rutabaga chromosome A03 and the middle segment of chromosome A08 as genomic hotspots associated with resistance to the different P. brassicae pathotypes.

Integrating long noncoding RNAs and mRNAs expression profiles of response to Plasmodiophora brassicae infection in Pakchoi (Brassica campestris ssp. chinensis Makino)

The biotrophic protist Plasmodiophora brassicae causes serious damage to Brassicaceae crops grown worldwide. However, the molecular mechanism of the Brassica rapa response remains has not been determined. Long noncoding RNA and mRNA expression profiles in response to Plasmodiophora brassicae infection were investigated using RNA-seq on the Chinese cabbage inbred line C22 infected with P. brassicae. Approximately 5,193 mRNAs were significantly differentially expressed, among which 1,345 were upregulated and 3,848 were downregulated. The GO enrichment analysis shows that most of these mRNAs are related to the defense response. Meanwhile, 114 significantly differentially expressed lncRNAs were identified, including 31 upregulated and 83 downregulated. Furthermore, a total of 2,344 interaction relationships were detected between 1,725 mRNAs and 103 lncRNAs with a correlation coefficient greater than 0.8. We also found 15 P. brassicaerelated mRNAs and 16 lncRNA interactions within the correlation network. The functional annotation showed that 15 mRNAs belong to defense response proteins (66.67%), protein phosphorylation (13.33%), root hair cell differentiation (13.33%) and regulation of salicylic acid biosynthetic process (6.67%). KEGG annotation showed that the vast majority of these genes are involved in the biosynthesis of secondary metabolism pathways and plant-pathogen interactions. These results provide a new perspective on lncRNA-mRNA network function and help to elucidate the molecular mechanism of P. brassicae infection.

iTRAQ-based quantitative analysis reveals proteomic changes in Chinese cabbage (Brassica rapa L.) in response to Plasmodiophora brassicae infection

Clubroot disease is one of the major diseases affecting Brassica crops, especially Chinese cabbage (Brassica rapa L. ssp. pekinensis), which is known to be highly susceptible to the disease. In this study, the obligate biotrophic protist Plasmodiophora brassicae Woronin was used to infect the roots of Chinese cabbage seedlings. The disease symptoms were noticeable at 28 and 35 days after inoculation (DAI) in the susceptible (CM) line. Using isobaric tags for relative and absolute quantitation (iTRAQ) analysis, a total of 5,003 proteins of differential abundance were identified in the resistant/susceptible lines, which could be quantitated by dipeptide or polypeptide segments. Gene ontology (GO) analysis indicated that the differentially expressed proteins (DEPs) between the susceptible (CM) and resistant (CCR) lines were associated with the glutathione transferase activity pathway, which could catalyze the combination of glutathione and other electrophilic compounds to protect plants from disease. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the DEPs may be significantly enriched cytokinin signaling or arginine biosynthesis pathways, both of which are responses to stimuli and are plant defense reactions. The cytokinins may facilitate cell division in the shoot, resulting in the hypertrophy and formation of galls and the presentation of typical clubroot symptoms. In this study, the proteomic results provide a new perspective for creating germplasm resistance to P. brassicae, as well as a genetic basis for breeding to improve Chinese cabbage.

Analysing the genetic architecture of clubroot resistance variation in Brassica napus by associative transcriptomics

Clubroot is a destructive soil-borne pathogen of Brassicaceae that causes significant recurrent reductions in yield of cruciferous crops. Although there is some resistance in oilseed rape (a crop type of the species Brassica napus), the genetic basis of that resistance is poorly understood. In this study, we used an associative transcriptomics approach to elucidate the genetic basis of resistance to clubroot pathotype ECD 17/31/31 across a genetic diversity panel of 245 accessions of B. napus. A single nucleotide polymorphism (SNP) association analysis was performed with 256,397 SNPs distributed across the genome of B. napus and combined with transcript abundance data of 53,889 coding DNA sequence (CDS) gene models. The SNP association analysis identified two major loci (on chromosomes A2 and A3) controlling resistance and seven minor loci. Within these were a total of 86 SNP markers. Altogether, 392 genes were found in these regions. Another 21 genes were implicated as potentially involved in resistance using gene expression marker (GEM) analysis. After GO enrichment analysis and InterPro functional analysis of the identified genes, 82 candidate genes were identified as having roles in clubroot resistance. These results provide useful information for marker-assisted breeding which could lead to acceleration of pyramiding of multiple clubroot resistance genes in new varieties.

Two QTLs controlling Clubroot resistance identified from Bulked Segregant Sequencing in Pakchoi (Brassica campestris ssp. chinensis Makino)

Clubroot, caused by Plasmodiophora Brassicae, is a serious soil-borne disease in worldwide. In recent years, progression of clubroot is rapid and serious in Shanghai, China. In this study, The inheritance of clubroot resistance (CR) were determined in pakchoi using F₂ segregation population that were developed by crossing highly resistant line ‘CR38’ and susceptible line ‘CS22’. Two novel QTLs, qBrCR38-1 and qBrCR38-2, was identified by BSA-seq (Bulked Segregant Sequencing) resistant to P. brassicae physiological race 7. Two significant peak qBrCR38-1 and qBrCR38-2 were observed by three statistical methods between interval of 19.7–20.6 Mb in chromosome A07 and 20.0–20.6 Mb in chromosome A08, respectively. In addition, Polymorphic SNPs identified within target regions were converted to kompetitive allele-specific PCR (KASP) assays. In target regions of qBrCR38-1 and qBrCR38-2, there were twenty SNP sites identified, eleven KASP markers of which are significantly associated to CR (P < 0.05). Seven candidate genes were identified and found to be involved in disease resistance (TIR-NBS-LRR proteins), defense responses of bacterium and fungi and biotic/abiotic stress response in the target regions harboring the two QTLs. Two novel QTLs and candidate genes identified from the present study provide insights into the genetic mechanism of CR in B.rapa, and the associated SNPs can be effectively used for marker-assisted breeding.

Clubroot resistance gene Rcr6 in Brassica nigra resides in a genomic region homologous to chromosome A08 in B. rapa

BACKGROUND

Clubroot, caused by Plasmodiophora brassicae Woronin, is a very important disease of Brassica species. Management of clubroot relies heavily on genetic resistance. In a cross of Brassica nigra lines PI 219576 (highly resistant, R) × CR2748 (highly susceptible, S) to clubroot, all F₁ plants were resistant to clubroot. There was a 1:1 ratio of R:S in the BC₁ and 3R:1S in the F₂, which indicated that a single dominant gene controlled clubroot resistance in PI 219576. This gene was designated Rcr6. Mapping of Rcr6 was performed using genome sequencing information from A-genome of B. rapa and B-genome of B. nigra though bulked segregant RNA sequencing (BSR-Seq) and further mapping with Kompetitive Allele Specific PCR (KASP) analysis.

RESULTS

Reads of R and S bulks from BSR-Seq were initially aligned onto B. rapa (A-genome; B. nigra has the B-genome) where Rcr6 was associated with chromosome A08. KASP analysis showed that Rcr6 was flanked by SNP markers homologous to the region of 14.8-15.4 Mb of chromosome A08. There were 190 genes annotated in this region, with five genes (Bra010552, Bra010588, Bra010589, Bra010590 and Bra010663) identified as encoding the toll-interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat (TIR-NBS-LRR; TNL) class of proteins. The reads from BSR-Seq were then aligned into a draft B-genome of B. nigra, where Rcr6 was mapped on chromosome B3. KASP analysis indicated that Rcr6 was located on chromosome B3 in a 0.5 Mb region from 6.1-6.6 Mb. Only one TNL gene homologous to the B. rapa gene Bra010663 was identified in the target region. This gene is a likely candidate for Rcr6. Subsequent analysis of the Rcr6 equivalent region based on a published B. nigra genome was performed. This gene is located into chromosome B7 of the published B-genome, homologous to BniB015819.

CONCLUSION

Rcr6 was the first gene identified and mapped in the B-genome of Brassica species. It resides in a genomic region homologous to chromosome A08 of A-genome. Based on this finding, it could possibly integrate into A08 of B. napus using marker assisted selection with SNP markers tightly linked to Rcr6 developed in this study.

Analysis of genome-wide variants through bulked segregant RNA sequencing reveals a major gene for resistance to Plasmodiophora brassicae in Brassica oleracea

Two cabbage (Brassica oleracea) cultivars 'Tekila' and 'Kilaherb' were identified as resistant to several pathotypes of Plasmodiophora brassicae. In this study, we identified a clubroot resistance gene (Rcr7) in 'Tekila' for resistance to pathotype 3 of P. brassicae from a segregating population derived from 'Tekila' crossed with the susceptible line T010000DH3. Genetic mapping was performed by identifying the percentage of polymorphic variants (PPV), a new method proposed in this study, through bulked segregant RNA sequencing. Chromosome C7 carried the highest PPV (42%) compared to the 30-34% in the remaining chromosomes. A peak with PPV (56-73%) was found within the physical interval 41-44 Mb, which indicated that Rcr7 might be located in this region. Kompetitive Allele-Specific PCR was used to confirm the association of Rcr7 with SNPs in the region. Rcr7 was flanked by two SNP markers and co-segregated with three SNP markers in the segregating population of 465 plants. Seven genes encoding TIR-NBS-LRR disease resistance proteins were identified in the target region, but only two genes, Bo7g108760 and Bo7g109000, were expressed. Resistance to pathotype 5X was also mapped to the same region as Rcr7. B. oleracea lines including 'Kilaherb' were tested with five SNP markers for Rcr7 and for resistance to pathotype 3; 11 of 25 lines were resistant, but 'Kilaherb' was the only line that carried the SNP alleles associated with Rcr7. The presence of Rcr7 in 'Kilaherb' for resistance to both pathotypes 3 and 5X was confirmed through linkage analysis.

Combinations of Independent Dominant Loci Conferring Clubroot Resistance in All Four Turnip Accessions (Brassica rapa) From the European Clubroot Differential Set

Clubroot disease is devastating to Brassica crop production when susceptible cultivars are planted in infected fields. European turnips are the most resistant sources and their resistance genes have been introduced into other crops such oilseed rape (Brassica napus L.), Chinese cabbage and other Brassica vegetables. The European clubroot differential (ECD) set contains four turnip accessions (ECD1-4). These ECD turnips exhibited high levels of resistance to clubroot when they were tested under controlled environmental conditions with Canadian field isolates. Gene mapping of the clubroot resistance genes in ECD1-4 were performed and three independent dominant resistance loci were identified. Two resistance loci were mapped on chromosome A03 and the third on chromosome A08. Each ECD turnip accession contained two of these three resistance loci. Some resistance loci were homozygous in ECD accessions while others showed heterozygosity based on the segregation of clubroot resistance in 20 BC₁ families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

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

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.

Identification and Mapping of the Clubroot Resistance Gene CRd in Chinese Cabbage (Brassica rapa ssp. pekinensis)

The rapid spread of clubroot disease, which is caused by Plasmodiophora brassicae, threatens Brassicaceae crop production worldwide. Breeding plants that have broad-spectrum disease resistance is one of the best ways to prevent clubroot. In the present study, eight Chinese cabbage germplasms were screened using published clubroot-resistant (CR) loci-/gene-linked markers. A CR gene Crr3 potential carrier "85-74" was detected which linked to marker BRSTS61; however, "85-74" shows different responses to local pathogens "LAB-19," "LNND-2," and "LAB-10" from "CR-73" which harbors Crr3. We used a next-generation sequencing-based bulked segregant analysis approach combined with genetic mapping to detect CR genes in an F₂ segregant population generated from a cross between the Chinese cabbage inbred lines "85-74" (CR) and "BJN3-1" (clubroot susceptible). The "85-74" line showed resistance to a local pathogen "LAB-19" which was identified as race 4; a genetic analysis revealed that the resistance was conferred by a single dominant gene. The CR gene which we named CRd was mapped to a 60 kb (1 cM) region between markers yau389 and yau376 on chromosome A03. CRd is located upstream of Crr3 which was confirmed based on the physical positions of Crr3 linked markers. The identification of CRd linked markers can be applied to marker-assisted selection in the breeding of new CR cultivars of Chinese cabbage and other Brassica crops.

Identification and Mapping of the Clubroot Resistance Gene CRd in Chinese Cabbage (Brassica rapa ssp. pekinensis)

The rapid spread of clubroot disease, which is caused by Plasmodiophora brassicae, threatens Brassicaceae crop production worldwide. Breeding plants that have broad-spectrum disease resistance is one of the best ways to prevent clubroot. In the present study, eight Chinese cabbage germplasms were screened using published clubroot-resistant (CR) loci-/gene-linked markers. A CR gene Crr3 potential carrier "85-74" was detected which linked to marker BRSTS61; however, "85-74" shows different responses to local pathogens "LAB-19," "LNND-2," and "LAB-10" from "CR-73" which harbors Crr3. We used a next-generation sequencing-based bulked segregant analysis approach combined with genetic mapping to detect CR genes in an F₂ segregant population generated from a cross between the Chinese cabbage inbred lines "85-74" (CR) and "BJN3-1" (clubroot susceptible). The "85-74" line showed resistance to a local pathogen "LAB-19" which was identified as race 4; a genetic analysis revealed that the resistance was conferred by a single dominant gene. The CR gene which we named CRd was mapped to a 60 kb (1 cM) region between markers yau389 and yau376 on chromosome A03. CRd is located upstream of Crr3 which was confirmed based on the physical positions of Crr3 linked markers. The identification of CRd linked markers can be applied to marker-assisted selection in the breeding of new CR cultivars of Chinese cabbage and other Brassica crops.

Combinations of Independent Dominant Loci Conferring Clubroot Resistance in All Four Turnip Accessions (Brassica rapa) From the European Clubroot Differential Set

Clubroot disease is devastating to Brassica crop production when susceptible cultivars are planted in infected fields. European turnips are the most resistant sources and their resistance genes have been introduced into other crops such oilseed rape (Brassica napus L.), Chinese cabbage and other Brassica vegetables. The European clubroot differential (ECD) set contains four turnip accessions (ECD1-4). These ECD turnips exhibited high levels of resistance to clubroot when they were tested under controlled environmental conditions with Canadian field isolates. Gene mapping of the clubroot resistance genes in ECD1-4 were performed and three independent dominant resistance loci were identified. Two resistance loci were mapped on chromosome A03 and the third on chromosome A08. Each ECD turnip accession contained two of these three resistance loci. Some resistance loci were homozygous in ECD accessions while others showed heterozygosity based on the segregation of clubroot resistance in 20 BC₁ families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Fine Mapping of a Clubroot Resistance Gene in Chinese Cabbage Using SNP Markers Identified from Bulked Segregant RNA Sequencing

Clubroot, caused by Plasmodiophora brassicae, is an important disease of canola (Brassica napus) in western Canada and worldwide. In this study, a clubroot resistance gene (Rcr2) was identified and fine mapped in Chinese cabbage cv. "Jazz" using single-nucleotide polymorphisms (SNP) markers identified from bulked segregant RNA sequencing (BSR-Seq) and molecular markers were developed for use in marker assisted selection. In total, 203.9 million raw reads were generated from one pooled resistant (R) and one pooled susceptible (S) sample, and >173,000 polymorphic SNP sites were identified between the R and S samples. One significant peak was observed between 22 and 26 Mb of chromosome A03, which had been predicted by BSR-Seq to contain the causal gene Rcr2. There were 490 polymorphic SNP sites identified in the region. A segregating population consisting of 675 plants was analyzed with 15 SNP sites in the region using the Kompetitive Allele Specific PCR method, and Rcr2 was fine mapped between two SNP markers, SNP_A03_32 and SNP_A03_67 with 0.1 and 0.3 cM from Rcr2, respectively. Five SNP markers co-segregated with Rcr2 in this region. Variants were identified in 14 of 36 genes annotated in the Rcr2 target region. The numbers of poly variants differed among the genes. Four genes encode TIR-NBS-LRR proteins and two of them Bra019410 and Bra019413, had high numbers of polymorphic variants and so are the most likely candidates of Rcr2.

Genotyping-by-sequencing reveals three QTL for clubroot resistance to six pathotypes of Plasmodiophora brassicae in Brassica rapa

Clubroot, caused by Plasmodiophora brassicae, is an important disease of Brassica crops worldwide. F1 progeny from the Brassica rapa lines T19 (resistant) × ACDC (susceptible) were backcrossed with ACDC, then self-pollinated to produce BC1S1 lines, From genotyping-by-sequencing (GBS) of the parental lines and BC1 plants, about 1.32 M sequences from T19 were aligned into the reference genome of B. rapa with 0.4-fold coverage, and 1.77 M sequences with 0.5-fold coverage in ACDC. The number of aligned short reads per plant in the BC1 ranged from 0.07 to 1.41 M sequences with 0.1-fold coverage. A total of 1584 high quality SNP loci were obtained, distributed on 10 chromosomes. A single co-localized QTL, designated as Rcr4 on chromosome A03, conferred resistance to pathotypes 2, 3, 5, 6 and 8. The peak was at SNP locus A03_23710236, where LOD values were 30.3 to 38.8, with phenotypic variation explained (PVE) of 85-95%. Two QTLs for resistance to a novel P. brassicae pathotype 5x, designated Rcr8 on chromosome A02 and Rcr9 on A08, were detected with 15.0 LOD and 15.8 LOD, and PVE of 36% and 39%, respectively. Bulked segregant analysis was performed to examine TIR-NBS-LRR proteins in the regions harboring the QTL.

Clubroot resistance QTL are modulated by nitrogen input in Brassica napus

Nitrogen levels can modulate the effectiveness of clubroot resistance in an isolate- and host-specific manner. While the same QTL were detected under high and low nitrogen, their effects were altered. Clubroot, caused by Plasmodiophora brassicae, is one of the most damaging diseases of oilseed rape and is known to be affected by nitrogen fertilization. However, the genetic factors involved in clubroot resistance have not been characterized under nitrogen-limiting conditions. This study aimed to assess the variability of clubroot resistance under different nitrogen levels and to characterize the impact of nitrogen supply on genetic resistance factors. Linkage analyses and a genome-wide association study were conducted to detect QTL for clubroot resistance and evaluate their sensitivity to nitrogen. The clubroot response of a set of 92 diverse oilseed rape accessions and 108 lines derived from a cross between 'Darmor-bzh' (resistant) and 'Yudal' (susceptible) was studied in the greenhouse under high- and low-nitrogen conditions, following inoculation with the P. brassicae isolates eH and K92-16. Resistance to each isolate was controlled by a major QTL and a few small-effects QTL. While the same QTL were detected under both high and low nitrogen, their effects were altered. Clubroot resistance to isolate eH, but not K92-16, was greater under a low-N supply versus a high-N supply. New sources of resistance were found among the oilseed rape accessions under both low and high-N conditions. The results are discussed relative to the literature and from a crop improvement perspective.

A Genome-Wide Association Study Reveals New Loci for Resistance to Clubroot Disease in Brassica napus

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. However, the yield and quality of rapeseed were largely decreased by clubroot (Plasmodiophora brassicae Woronin). Therefore, it is of great importance for screening more resistant germplasms or genes and improving the resistance to P. brassicae in rapeseed breeding. In this study, a massive resistant identification for a natural global population was conducted in two environments with race/pathotype 4 of P. brassicae which was the most predominant in China, and a wide range of phenotypic variation was found in the population. In addition, a genome-wide association study of 472 accessions for clubroot resistance (CR) was performed with 60K Brassica Infinium SNP arrays for the first time. In total, nine QTLs were detected, seven of which were novel through integrative analysis. Furthermore, additive effects in genetic control of CR in rapeseed among the above loci were found. By bioinformatic analyses, the candidate genes of these loci were predicted, which indicated that TIR-NBS gene family might play an important role in CR. It is believable that the results presented in our study could provide valuable information for understanding the genetic mechanism and molecular regulation of CR.

A Genome-Wide Association Study Reveals New Loci for Resistance to Clubroot Disease in Brassica napus

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. However, the yield and quality of rapeseed were largely decreased by clubroot (Plasmodiophora brassicae Woronin). Therefore, it is of great importance for screening more resistant germplasms or genes and improving the resistance to P. brassicae in rapeseed breeding. In this study, a massive resistant identification for a natural global population was conducted in two environments with race/pathotype 4 of P. brassicae which was the most predominant in China, and a wide range of phenotypic variation was found in the population. In addition, a genome-wide association study of 472 accessions for clubroot resistance (CR) was performed with 60K Brassica Infinium SNP arrays for the first time. In total, nine QTLs were detected, seven of which were novel through integrative analysis. Furthermore, additive effects in genetic control of CR in rapeseed among the above loci were found. By bioinformatic analyses, the candidate genes of these loci were predicted, which indicated that TIR-NBS gene family might play an important role in CR. It is believable that the results presented in our study could provide valuable information for understanding the genetic mechanism and molecular regulation of CR.