Pyramiding of triple Clubroot resistance loci conferred superior resistance without negative effects on agronomic traits in Brassica napus.

Clubroot disease caused by Plasmodiophora brassicae is becoming a serious threat to rapeseed (Brassica napus) production worldwide. Breeding resistant varieties using CR (clubroot resistance) loci is the most promising solution. Using marker-assisted selection and speed-breeding technologies, we generated Brassica napus materials in homozygous or heterozygous states using CRA3.7, CRA08.1, and CRA3.2 loci in the elite parental line of the Zhongshuang11 background. We developed three elite lines with two CR loci in different combinations and one line with three CR loci at the homozygous state. In our study, we used six different clubroot strains (Xinmin, Lincang, Yuxi, Chengdu, Chongqing, and Jixi) which are categorized into three groups based on our screening results. The newly pyramided lines with two or more CR loci displayed better disease resistance than the parental lines carrying single CR loci. There is an obvious gene dosage effect between CR loci and disease resistance levels. For example, pyramided lines with triple CR loci in the homozygous state showed superior resistance for all pathogens tested. Moreover, CR loci in the homozygous state are better on disease resistance than the heterozygous state. More importantly, no negative effect was observed on agronomic traits for the presence of multiple CR loci in the same background. Overall, these data suggest that the pyramiding of triple clubroot resistance loci conferred superior resistance with no negative effects on agronomic traits in Brassica napus.

Pathotype Characterization of Plasmodiophora brassicae by European Clubroot Differential and Williams Sets in China.

Clubroot disease caused by the soil-borne Plasmodiophora brassicae is devastating to Brassicaceae crops and spreading rapidly in China in recent years, resulting in great yield losses annually. Virulence of P. brassicae populations specializes and is in dynamic change in the fields. Information on the pathotypes and their distributions is crucial to control the clubroot disease. Presently, the pathotypes of P. brassicae prevalent in China, however, are not well determined. In this study, we used 16 Brassica hosts, including the European Clubroot Differential (ECD) and Williams sets, to designate the pathotypes of 33 P. brassicae populations from 13 provinces. The 33 P. brassicae populations could be divided into 26 pathotypes by the ECD set or seven pathotypes by the Williams set, revealing ECD16/15/31 and ECD16/31/31 or P4 and P2 as the predominant pathotypes. We found that the Brassica rapa differentials ECD01 to ECD04 showed stable and high levels of resistance to most pathotypes of P. brassicae in China, thereby providing valuable resources for clubroot-resistance breeding of Brassicaceae crops. The ECD set exhibited much higher discernibility and further divided the isolates that belonged to the P4 pathotype into 10 ECD pathotypes. Isolates of ECD16/23/31 and ECD16/15/31 were strongly virulent on Huashuang 5R, the first and widely used clubroot-resistant cultivar of oilseed rape in China. As we learn, 26 pathotypes are the most diverse populations of P. brassicae characterized until now in China. Our study provides new insights into virulence specialization of P. brassicae and their geographical distributions, contributing to exploitation of clubroot-resistant resources and the field layout of the present resistant Brassica crops in China.

RNA-Seq Bulked Segregant Analysis of an Exotic B. napus ssp. napobrassica (Rutabaga) F2 Population Reveals Novel QTLs for Breeding Clubroot-Resistant Canola.

In this study, a rutabaga (Brassica napus ssp. napobrassica) donor parent FGRA106, which exhibited broad-spectrum resistance to 17 isolates representing 16 pathotypes of Plasmodiophora brassicae, was used in genetic crosses with the susceptible spring-type canola (B. napus ssp. napus) accession FG769. The F2 plants derived from a clubroot-resistant F1 plant were screened against three P. brassicae isolates representing pathotypes 3A, 3D, and 3H. Chi-square (χ2) goodness-of-fit tests indicated that the F2 plants inherited two major clubroot resistance genes from the CR donor FGRA106. The total RNA from plants resistant (R) and susceptible (S) to each pathotype were pooled and subjected to bulked segregant RNA-sequencing (BSR-Seq). The analysis of gene expression profiles identified 431, 67, and 98 differentially expressed genes (DEGs) between the R and S bulks. The variant calling method indicated a total of 12 (7 major + 5 minor) QTLs across seven chromosomes. The seven major QTLs included: BnaA5P3A.CRX1.1, BnaC1P3H.CRX1.2, and BnaC7P3A.CRX1.1 on chromosomes A05, C01, and C07, respectively; and BnaA8P3D.CRX1.1, BnaA8P3D.RCr91.2/BnaA8P3H.RCr91.2, BnaA8P3H.Crr11.3/BnaA8P3D.Crr11.3, and BnaA8P3D.qBrCR381.4 on chromosome A08. A total of 16 of the DEGs were located in the major QTL regions, 13 of which were on chromosome C07. The molecular data suggested that clubroot resistance in FGRA106 may be controlled by major and minor genes on both the A and C genomes, which are deployed in different combinations to confer resistance to the different isolates. This study provides valuable germplasm for the breeding of clubroot-resistant B. napus cultivars in Western Canada.

Efficient marker-assisted breeding for clubroot resistance in elite Pol-CMS rapeseed varieties by updating the PbBa8.1 locus.

Clubroot disease poses a severe threat to rapeseed (Brassica napus) production worldwide and has recently been spreading across China at an unprecedented pace. Breeding and cultivation of resistant varieties constitute a promising and environment-friendly approach to mitigating this threat. In this study, the clubroot resistance locus PbBa8.1 was successfully transferred into SC4, a shared paternal line of three elite varieties in five generations by marker-assisted backcross breeding. Kompetitive allele specific PCR (KASP) markers of clubroot resistance gene PbBa8.1 and its linked high erucic acid gene (FAE1) were designed and applied for foreground selection, and 1,000 single-nucleotide polymorphisms (SNPs) were selected and used for the background selection. This breeding strategy produced recombinants with the highest recovery ratio of the recurrent parent genome (> 95%) at BC2F2 while breaking the linkage with FAE1 during the selection. An updated version of the paternal line (SC4R) was generated at BC2F3, showing significantly improved clubroot resistance at the seedling stage via artificial inoculation, and was comparable to that of the donor parent. Field trials of the three elite varieties and their updated versions in five environments indicated similar agronomic appearance and final yield. The introduced breeding strategy precisely pyramids the PbBa8.1 and FAE1 loci with the assistance of technical markers in a shorter period and could be applied to other desirable traits for directional improvement in the future. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01305-9.

Transferring of clubroot-resistant locus CRd from Chinese cabbage (Brassica rapa) to canola (Brassica napus) through interspecific hybridization.

Clubroot, caused by Plasmodiophora brassicae is one of the most severe threats to brassica species in China and worldwide. Breeding for clubroot resistant varieties is one of the best ways to overcome this disease. In this study, we introduced clubroot resistance (CR) gene CRd from Chinese cabbage (85-74) into elite Brassica napus inbred line Zhongshuang 11 through interspecific hybridization and subsequent backcrossing with whole-genome molecular marker-assisted selection (MAS). The resistant test of CRd to P. brassicae isolates was evaluated in the greenhouse as well as in field conditions. Close linkage markers and the whole-chromosome background marker selection approach improved the recovery rate from 78.3% in BC1 to 100% in BC3F1. The improved clubroot-resistant variety, Zhongshuang11R, was successfully selected in the BC3F2 generation. The greenhouse and field resistant tests revealed that Zhongshuang11R was resistant to P. brassicae pathotypes. The agronomic characteristics of Zhongshuang11R were similar to those of its recurrent parental line, including oil content, composition of fatty acid, plant height, primary effective branches, grain yield per plant and thousand-seed weight. In addition, the oil quality could satisfy the quality requirements for commercial rapeseed oil. Our results will enrich the resistant resources of canola and will certainly accelerate clubroot resistance breeding programs in B. napus.

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 F1 cultivars of Chinese cabbage. Gain-of-function analysis revealed that Crr1aKinami90_a isolated from the cv. 'Kinami 90' conferred clubroot resistance as observed for Crr1aG004 . 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.

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. BC1 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 Rcr9wa, 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 Rcr9wa interval, three of which (Bra020827, Bra020828, Bra020814) were annotated as immune-system-process related genes, which are potential candidates for Rcr9wa.

Classification of "nabana" (Brassica rapa) cultivars and landraces based on simple sequence repeat markers.

Brassica rapa or B. napus vegetables for eating as young inflorescences and stalks are called "nabana". Japanese nabana includes "flower-bud type" and "stem-and-leaf type". Chinese and European types are also known (cai-xin, zicaitai, and broccoletto). We classified nabana belonging to B. rapa and other B. rapa vegetables. In a simple sequence repeat-based phylogram, 49 ingroup samples were classified into four groups (I-IV). Flower-bud and stem-and-leaf types were separated into groups I and III, respectively, with a slight overlap in group II. Cai-xin and non-heading Chinese cabbages were included in group IV. Broccoletto was placed in group III, close to turnips. Zicaitai cultivars were included in group II. We tested for clubroot resistance (CR) and its marker genotypes in nabana because of their agronomical importance. Ten cultivars were resistant to group 4 pathogen but not to group 2. Most of the CR cultivars had heterozygous resistance alleles in the CRb and Crr1 loci, consistent with inoculation tests. Our results suggest that Japanese nabana lines and foreign types were differentiated according to their consumption parts and cultivar origins, respectively. This study elucidates the relationships and CR properties of nabana and provides valuable information for the breeding of nabana cultivars.

Genetics and molecular mapping of resistance to Plasmodiophora brassicae pathotypes 2, 3, 5, 6, and 8 in rutabaga (Brassica napus var. napobrassica).

Clubroot disease, caused by Plasmodiophora brassicae, is a threat to the production of Brassica crops including oilseed B. napus. In Canada, several pathotypes of this pathogen, such as pathotypes 2, 3, 5, 6, and 8, were identified, and resistance to these pathotypes was found in a rutabaga (B. napus var. napobrassica) genotype. In this paper, we report the genetic basis and molecular mapping of this resistance by use of F2, backcross (BC1), and doubled haploid (DH) populations generated from crossing of this rutabaga line to a susceptible spring B. napus canola line. The F1, F2, and BC1 populations were evaluated for resistance to pathotype 3, and the DH population was evaluated for resistance to pathotypes 2, 3, 5, 6, and 8. A 3:1 segregation in F2 and a 1:1 segregation in BC1 were found for resistance to pathotype 3, and a 1:1 segregation was found in the DH population for resistance to all pathotypes. Molecular mapping by using the DH population identified a genomic region on chromosome A8 carrying resistance to all five pathotypes. This suggests that a single gene or a cluster of genes, located in this genomic region, is involved in the control of resistance to these pathotypes.

Clubroot-Induced Changes in the Root and Rhizosphere Microbiome of Susceptible and Resistant Canola.

Clubroot is a soilborne disease of canola (Brassica napus) and other crucifers caused by the obligate parasite Plasmodiophora brassicae. In western Canada, clubroot is usually managed by planting-resistant cultivars, but the emergence of resistance-breaking pathotypes of P. brassicae represents a major threat to sustainable canola production. The rhizosphere and root contain beneficial microorganisms that can improve plant health. In this study, we evaluated the effect of two P. brassicae isolates (termed A and B) with different levels of virulence on the root and rhizosphere microbiomes of clubroot-resistant and clubroot-susceptible canola. Additionally, potential biocontrol microorganisms were identified based on taxa antagonistic to clubroot. Although both P. brassicae isolates were classified as pathotype 3A, isolate A caused a higher disease severity index in the resistant canola genotype compared with isolate B. Metabarcoding analysis indicated a shift in the bacterial and fungal communities in response to inoculation with either field isolate. Root endophytic bacterial and fungal communities responded to changes in inoculation, isolate type, sampling time, and canola genotype. In contrast, fungal communities associated with the rhizosphere exhibited significant differences between sampling times, while bacterial communities associated with the rhizosphere exhibited low variability.

Resynthesizing Brassica napus with race specific resistance genes and race non-specific QTLs to multiple races of Plasmodiophora brassicae.

Clubroot disease in canola (Brassica napus) continues to spread across the Canadian prairies. Growing resistant cultivars is considered the most economical means of controlling the disease. However, sources of resistance to clubroot in B. napus are very limited. In this study, we conducted interspecific crosses using a B. rapa line (T19) carrying race-specific resistance genes and two B. oleracea lines, ECD11 and JL04, carrying race non-specific QTLs. Employing embryo rescue and conventional breeding methods, we successfully resynthesized a total of eight B. napus lines, with four derived from T19 × ECD11 and four from T19 × JL04. Additionally, four semi-resynthesized lines were developed through crosses with a canola line (DH16516). Testing for resistance to eight significant races of Plasmodiophora brassicae was conducted on seven resynthesized lines and four semi-resynthesized lines. All lines exhibited high resistance to the strains. Confirmation of the presence of clubroot resistance genes/QTLs was performed in the resynthesized lines using SNP markers linked to race-specific genes in T19 and race non-specific QTLs in ECD11. The developed B. napus germplasms containing clubroot resistance are highly valuable for the development of canola cultivars resistant to clubroot.

Lignin accumulation in cell wall plays a role in clubroot resistance.

Clubroot, caused by Plasmodiophora brassicae, is a significant disease affecting brassica crops worldwide and poses a threat to canola (Brassica napus) production in western Canada. Management of this disease heavily relies on the use of resistant cultivars, but resistance erosion is a serious concern due to the highly diverse pathogen populations. Understanding resistance mechanisms may aid in better deployment/rotation of clubroot resistance (CR) genes and improve resistance resilience. In this study, we conducted a comparative analysis using resistant canola varieties carrying either a single (Rcr1) or double CR genes (Rcr1+Crr1rutb ) to decipher the resistance modes associated with these genes. Cell wall (CW) biopolymeric compounds in different root layers were mapped and quantified using Fourier-transform mid-infrared microspectroscopy for changes in CW elements associated with clubroot resistance. Transmission electron and confocal microscopy were used to assess root infection details and relative transcript abundance was analyzed to determine the activation of the lignin-related pathway in relation to resistance. Neither resistant variety affected the primary infection of root hairs/epidermal cells compared to the susceptible "Westar", but both exhibited strong inhibition of cortical infection, effectively 'trapping' the pathogen in the exodermis. The most prominent change observed was increased lignin accumulation associated with resistance. In Westar, the pathogen was able to degrade CW lignin, facilitating access to the root cortex by secondary plasmodia of P. brassicae. In contrast, resistant varieties showed clear lignin accumulation around the penetration site on the exodermis, accompanied by elevated expression of genes involved in the phenylpropanoid pathway. These results suggest that induced lignin accumulation plays a role in clubroot resistance mediated by the CR genes Rcr1 and Crr1rutb in canola, providing cellular and structural evidence that supports the data from earlier transcriptomic studies.

Identification of Micro Ribonucleic Acids and Their Targets in Response to Plasmodiophora brassicae Infection in Brassica napus.

Clubroot disease, which is caused by the soil-borne pathogen Plasmodiophora brassicae War (P. brassicae), is one of the oldest and most destructive diseases of Brassica and cruciferous crops in the world. Plant microRNAs [micro ribonucleic acids (miRNAs)] play important regulatory roles in several developmental processes. Although the role of plant miRNAs in plant-microbe interaction has been extensively studied, there are only few reports on the specific functions of miRNAs in response to P. brassicae. This study investigated the roles of miRNAs and their targets during P. brassicae infection in a pair of Brassica napus near-isogenic lines (NILs), namely clubroot-resistant line 409R and clubroot-susceptible line 409S. Small RNA sequencing (sRNA-seq) and degradome-seq were performed on root samples of 409R and 409S with or without P. brassicae inoculation. sRNA-seq identified a total of 48 conserved and 72 novel miRNAs, among which 18 had a significant differential expression in the root of 409R, while only one miRNA was differentially expressed in the root of 409S after P. brassicae inoculation. The degradome-seq analysis identified 938 miRNA target transcripts, which are transcription factors, enzymes, and proteins involved in multiple biological processes and most significantly enriched in the plant hormone signal transduction pathway. Between 409R and 409S, we found eight different degradation pathways in response to P. brassicae infection, such as those related to fatty acids. By combining published transcriptome data, we identified a total of six antagonistic miRNA-target pairs in 409R that are responsive to P. brassicae infection and involved in pathways associated with root development, hypersensitive cell death, and chloroplast metabolic synthesis. Our results reveal that P. brassicae infection leads to great changes in miRNA pool and target transcripts. More interestingly, these changes are different between 409R and 409S. Clarification of the crosstalk between miRNAs and their targets may shed new light on the possible mechanisms underlying the pathogen resistance against P. brassicae.

Root Transcriptome and Metabolome Profiling Reveal Key Phytohormone-Related Genes and Pathways Involved Clubroot Resistance in Brassica rapa L.

Plasmodiophora brassicae, an obligate biotrophic pathogen-causing clubroot disease, can seriously affect Brassica crops worldwide, especially Chinese cabbage. Understanding the transcriptome and metabolome profiling changes during the infection of P. brassicae will provide key insights in understanding the defense mechanism in Brassica crops. In this study, we estimated the phytohormones using targeted metabolome assays and transcriptomic changes using RNA sequencing (RNA-seq) in the roots of resistant (BrT24) and susceptible (Y510-9) plants at 0, 3, 9, and 20 days after inoculation (DAI) with P. brassicae. Differentially expressed genes (DEGs) in resistant vs. susceptible lines across different time points were identified. The weighted gene co-expression network analysis of the DEGs revealed six pathways including "Plant-pathogen interaction" and "Plant hormone signal transduction" and 15 hub genes including pathogenic type III effector avirulence factor gene (RIN4) and auxin-responsive protein (IAA16) to be involved in plants immune response. Inhibition of Indoleacetic acid, cytokinin, jasmonate acid, and salicylic acid contents and changes in related gene expression in R-line may play important roles in regulation of clubroot resistance (CR). Based on the combined metabolome profiling and hormone-related transcriptomic responses, we propose a general model of hormone-mediated defense mechanism. This study definitely enhances our current understanding and paves the way for improving CR in Brassica rapa.

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 BC1 families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

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.