Fine mapping of the clubroot resistance gene CRb and development of a useful selectable marker in Brassica rapa

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.

Comparative transcriptome analysis of canola carrying a single vs stacked resistance genes against clubroot

Pyramiding resistance genes may expand the efficacy and scope of a canola variety against clubroot (Plasmodiophora brassicae), a serious threat to canola production in western Canada. However, the mechanism(s) of multigenic resistance, especially the potential interaction among clubroot resistance (CR) genes, are not well understood. In this study, transcriptome was compared over three canola (Brassica napus L.) inbred/hybrid lines carrying a single CR gene in chromosome A03 (CRaM, Line 16) or A08 (Crr1rutb, Line 20), and both genes (CRaM+Crr1rutb, Line 15) inoculated with a field population (L-G2) of P. brassicae pathotype X, a new variant found in western Canada recently. The line16 was susceptible, while lines 15 and 20 were partially resistant. Functional annotation identified differential expression of genes (DEGs) involved in biosynthetic processes responsive to stress and regulation of cellular process; The Venn diagram showed that the partially resistant lines 15 and 20 shared 1,896 differentially expressed genes relative to the susceptible line 16, and many of these DEGs are involved in defense responses, activation of innate immunity, hormone biosynthesis and programmed cell death. The transcription of genes involved in Pathogen-Associated Molecular Pattern (PAMP)-Triggered and Effector-Triggered Immunity (PTI and ETI) was particularly up-regulated, and the transcription level was higher in line 15 (CRaM + Crr1rutb) than in line 20 (Crr1rutb only) for most of the DEGs. These results indicated that the partial resistance to the pathotype X was likely conferred by the CR gene Crr1rutb for both lines 15 and 20 that functioned via the activation of both PTI and ETI signaling pathways. Additionally, these two CR genes might have synergistic effects against the pathotype X, based on the higher transcription levels of defense-related DEGs expressed by inoculated line 15, highlighting the benefit of gene stacking for improved canola resistance as opposed to a single CR gene alone.

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.

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

Canola with Stacked Genes Shows Moderate Resistance and Resilience against a Field Population of Plasmodiophora brassicae (Clubroot) Pathotype X

Genetic resistance is a cornerstone for managing clubroot (Plasmodiophora brassicae). However, when used repeatedly, a clubroot resistance (CR) gene can be broken rapidly. In this study, canola inbred/hybrid lines carrying one or two CR genes (Rcr1/CRa^M and Crr1^rutb) were assessed against P. brassicae pathotype X by repeated exposure to the same inoculum source under a controlled environment. Lines carrying two CR genes, either Rcr1 + Crr1^rutb or CRa^M + Crr1^rutb, showed partial resistance. Selected lines were inoculated with a field pathotype X population (L-G3) at 5 × 10⁶ resting spores/g soil, and all clubs were returned to the soil they came from six weeks after inoculation. The planting was repeated for five cycles, with diseased roots being returned to the soil after each cycle. The soil inoculum was quantified using qPCR before each planting cycle. All lines with a single CR gene were consistently susceptible, maintaining high soil inoculum levels over time. The lines carrying two CR genes showed much lower clubroot severity, resulting in a 10-fold decline in soil inoculum. These results showed that the CR-gene stacking provided moderate resistance against P. brassicae pathotype X, which may also help reduce the pathogen inoculum buildup in soil.

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.

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.

Improvement of Resistance to Clubroot Disease in the Ogura CMS Restorer Line R2163 of Brassica napus

Oilseed rape (Brassica napus) has significant heterosis and Ogura CMS is a major way to use it. Ogura CMS has the advantages of complete and stable male sterility and easy-to-breed maintainers. Therefore, to breed better restorers has become an important goal for this system. Incidentally, clubroot is a soil-borne disease that is difficult to control by fungicidal chemicals, and it has been the main disease of oilseed rape in recent years in China, severely restricting the development of the oilseed rape industry. At present, the most effective method for controlling clubroot disease is to cultivate resistant varieties. One Ogura CMS restorer line (R2163) has shown much better combining ability, but lacks the clubroot disease resistance. This study was carried out to improve R2163 through marker-assisted backcross breeding (MABB). The resistant locus PbBa8.1 was introduced into the restorer R2163, and we then selected R2163R with clubroot disease resistance. Using the new restorer R2163R as the male parent and the sterile lines 116A and Z11A as the female parent, the improved, new resistant hybrids Kenyouza 741R and Huayouza 706R performed well, providing strong resistance and good agronomic traits. This work advances the utilization of heterosis and breeding for clubroot disease resistance in B. napus.

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.

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.

Development of molecular markers based on CRa gene sequencing of different clubroot disease-resistant cultivars of Chinese cabbage

BACKGROUND

CRa is a key gene in Chinese cabbage (Brassica rapa ssp. pekinensis) that confers resistance to Plasmodiophora brassicae. In order to efficiently screen the clubroot resistance (CR) gene CRa in breeding, two functional codominant markers of the CRa gene were developed.

METHODS AND RESULTS

In this study, through comparing the CRa allele sequences in resistant and susceptible cultivars of Chinese cabbage, we found two insertion and deletion of sequence variations in the fourth exon between resistant and susceptible cultivars. Two functional codominant markers for CRa gene were obtained based on the variations, namely, CRaEX04-1 and CRaEX04-3. The lengths of the extended fragment of CRaEX04-1 marker were 321 bp and 186 bp in resistant and susceptible cultivars, respectively. In contrast, those of CRaEX04-3 were 704 bp and 413 bp, respectively. We verified the genetic stability between the developed markers and CRa gene using 57 Chinese cabbage cultivars with known resistance and two genetic populations. The results showed that the marker identification was completely consistent with the known phenotypes in 57 cultivars. The marker identification results followed the 3:1 of Mendel's first law in the F₂ population, and the 1:1 of Mendel's first law in the BC₁.

CONCLUSIONS

CRaEX04-1 and CRaEX04-3 can be used as a practical molecular marker for breeding and germplasm resource creation of clubroot disease-resistant Chinese cabbage.

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.

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.

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.

Development of novel clubroot resistant rapeseed lines (Brassica napus L.) effective against Japanese field isolates by marker assisted selection

Clubroot is an important disease infectible to cruciferous plants and a major threat to rapeseed production in Japan. However, no clubroot resistant rapeseed cultivars have been released. We surveyed pathotype variation of six isolates collected from rapeseed fields and found they were classified as pathotype groups 2 and 4 using Japanese F₁ Chinese cabbage cultivars. We produced the resynthesized clubroot resistant Brassica napus harboring two resistant loci, Crr1 and Crr2, by interspecific crossing and developed resistant rapeseed lines for southern and northern regions by marker-assisted selection and backcrossing. We improved the DNA marker for erucic acid content to remove linkage drag between Crr1 and high erucic acid content and successfully selected lines with clubroot resistance and zero erucic acid for northern regions. A novel line, 'Tohoku No. 106', suitable for southern regions showed stable resistance against all six isolates and high performance in infested fields. We conclude that Crr1 and Crr2 are important genes for CR rapeseed breeding and marker-assisted selection is effective in improving clubroot resistance.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Molecular Breeding Strategy and Challenges Towards Improvement of Downy Mildew Resistance in Cauliflower (Brassica oleracea var. botrytis L.)

Cauliflower (Brassica oleracea var. botrytis L.) is one of the important, nutritious and healthy vegetable crops grown and consumed worldwide. But its production is constrained by several destructive fungal diseases and most importantly, downy mildew leading to severe yield and quality losses. For sustainable cauliflower production, developing resistant varieties/hybrids with durable resistance against broad-spectrum of pathogens is the best strategy for a long term and reliable solution. Identification of novel resistant resources, knowledge of the genetics of resistance, mapping and cloning of resistance QTLs and identification of candidate genes would facilitate molecular breeding for disease resistance in cauliflower. Advent of next-generation sequencing technologies (NGS) and publishing of draft genome sequence of cauliflower has opened the flood gate for new possibilities to develop enormous amount of genomic resources leading to mapping and cloning of resistance QTLs. In cauliflower, several molecular breeding approaches such as QTL mapping, marker-assisted backcrossing, gene pyramiding have been carried out to develop new resistant cultivars. Marker-assisted selection (MAS) would be beneficial in improving the precision in the selection of improved cultivars against multiple pathogens. This comprehensive review emphasizes the fascinating recent advances made in the application of molecular breeding approach for resistance against an important pathogen; Downy Mildew (Hyaloperonospora parasitica) affecting cauliflower and Brassica oleracea crops and highlights the QTLs identified imparting resistance against this pathogen. We have also emphasized the critical research areas as future perspectives to bridge the gap between availability of genomic resources and its utility in identifying resistance genes/QTLs to breed downy mildew resistant cultivars. Additionally, we have also discussed the challenges and the way forward to realize the full potential of molecular breeding for downy mildew resistance by integrating marker technology with conventional breeding in the post-genomics era. All this information will undoubtedly provide new insights to the researchers in formulating future breeding strategies in cauliflower to develop durable resistant cultivars against the major pathogens in general and downy mildew in particular.

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.

Comparison and evolutionary analysis of Brassica nucleotide binding site leucine rich repeat (NLR) genes and importance for disease resistance breeding

The Brassica genus contains many agriculturally significant oilseed and vegetable crops, however the crop yield is threatened by a range of fungal and bacterial pathogens. Nucleotide Binding Site Leucine Rich Repeat (NLR) genes play important roles in plant innate immunity. The evolution of NLR genes is influenced by genomic processes and pathogen selection. At the whole genome level, whole genome duplications (WGDs) generate abundant gene copies, most of which are lost during genome fractionation. At sub-genomic levels, some retained copies undergo duplication forming clusters which facilitate rapid evolution through recombination. The number, distribution and genetic variations of the NLR genes vary among Brassica species and within populations suggesting differential selection pressure exerted by pathogen populations throughout the evolutionary history. A study of the evolution of disease resistance genes in agriculturally important plants such as Brassicas helps gain insights into their function and inform the identification of resistance genes for breeding of resistant lines.

Characterization, identification and expression profiling of genome-wide R-genes in melon and their putative roles in bacterial fruit blotch resistance

Background

Bacterial fruit blotch (BFB), a disease caused by Acidovorax citrulli, results in significant economic losses in melon. The causal QTLs and genes for resistance to this disease have yet to be identified. Resistance (R)-genes play vital roles in resistance to plant diseases. Since the complete genome sequence of melon is available and genome-wide identification of R-genes has been performed for this important crop, comprehensive expression profiling may lead to the identification of putative candidate genes that function in the response to BFB.

Results

We identified melon accessions that are resistant and susceptible to BFB through repeated bioassays and characterized all 70 R-genes in melon, including their gene structures, chromosomal locations, domain organizations, motif distributions, and syntenic relationships. Several disease resistance-related domains were identified, including NBS, TIR, LRR, CC, RLK, and DUF domains, and the genes were categorized based on the domains of their encoded proteins. In addition, we profiled the expression patterns of the genes in melon accessions with contrasting levels of BFB resistance at 12 h, 1 d, 3 d, and 6 d after inoculation with A. citrulli. Six R-genes exhibited consistent expression patterns (MELO3C023441, MELO3C016529, MELO3C022157, MELO3C022146, MELO3C025518, and MELO3C004303), with higher expression levels in the resistant vs. susceptible accession.

Conclusion

We identified six putative candidate R-genes against BFB in melon. Upon functional validation, these genes could be targeted for manipulation via breeding and biotechnological approaches to improve BFB resistance in melon in the future.

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.

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.

Resolution of quantitative resistance to clubroot into QTL-specific metabolic modules

Plant disease resistance is often under quantitative genetic control. Thus, in a given interaction, plant cellular responses to infection are influenced by resistance or susceptibility alleles at different loci. In this study, a genetic linkage analysis was used to address the complexity of the metabolic responses of Brassica napus roots to infection by Plasmodiophora brassicae. Metabolome profiling and pathogen quantification in a segregating progeny allowed a comparative mapping of quantitative trait loci (QTLs) involved in resistance and in metabolic adjustments. Distinct metabolic modules were associated with each resistance QTL, suggesting the involvement of different underlying cellular mechanisms. This approach highlighted the possible role of gluconasturtiin and two unknown metabolites in the resistance conferred by two QTLs on chromosomes C03 and C09, respectively. Only two susceptibility biomarkers (glycine and glutathione) were simultaneously linked to the three main resistance QTLs, suggesting the central role of these compounds in the interaction. By contrast, several genotype-specific metabolic responses to infection were genetically unconnected to resistance or susceptibility. Likewise, variations of root sugar profiles, which might have influenced pathogen nutrition, were not found to be related to resistance QTLs. This work illustrates how genetic metabolomics can help to understand plant stress responses and their possible links with disease.

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.

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.

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.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

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.

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.

Resynthesis of Brassica juncea for resistance to Plasmodiophora brassicae pathotype 3

The oilseed crop Brassica juncea carries many desirable traits; however, resistance to clubroot disease, caused by Plasmodiophora brassicae, is not available in this species. We are the first to report the clubroot resistant resynthesized B. juncea lines, developed through interspecific crosses between a clubroot resistant B. rapa ssp. rapifera and two susceptible B. nigra lines, and the stability of the resistance in self-pollinated generations. The interspecific nature of the resynthesized B. juncea plants was confirmed by using A- and B-genome specific SSR markers, and flow cytometric analysis of nuclear DNA content. Self-pollinated progeny (S₁ and S₂) of the resynthesized B. juncea plants were evaluated for resistance to P. brassicae pathotype 3. The S₁ and S₂ progenies of one of the resynthesized B. juncea lines were resistant to this pathotype. However, resistance was lost in 6 to 13% plants of the S₂ progenies derived from the second resynthesized B. juncea line; this apparently resulted from the loss of the genomic region carrying resistance due to meiotic anomalies.

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.

Genome resequencing and comparative variome analysis in a Brassica rapa and Brassica oleracea collection

The closely related species Brassica rapa and B. oleracea encompass a wide range of vegetable, fodder and oil crops. The release of their reference genomes has facilitated resequencing collections of B. rapa and B. oleracea aiming to build their variome datasets. These data can be used to investigate the evolutionary relationships between and within the different species and the domestication of the crops, hereafter named morphotypes. These data can also be used in genetic studies aiming at the identification of genes that influence agronomic traits. We selected and resequenced 199 B. rapa and 119 B. oleracea accessions representing 12 and nine morphotypes, respectively. Based on these resequencing data, we obtained 2,249,473 and 3,852,169 high quality SNPs (single-nucleotide polymorphisms), as well as 303,617 and 417,004 InDels for the B. rapa and B. oleracea populations, respectively. The variome datasets of B. rapa and B. oleracea represent valuable resources to researchers working on evolution, domestication or breeding of Brassica vegetable crops.

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.