Mapping of a novel clubroot resistance QTL using ddRAD-seq in Chinese cabbage (Brassica rapa L.)

Multiomics analysis of a resistant European turnip ECD04 during clubroot infection reveals key hub genes underlying resistance mechanism

The clubroot disease has become a worldwide threat for crucifer crop production, due to its soil-borne nature and difficulty to eradicate completely from contaminated field. In this study we used an elite resistant European fodder turnip ECD04 and investigated its resistance mechanism using transcriptome, sRNA-seq, degradome and gene editing. A total of 1751 DEGs were identified from three time points after infection, among which 7 hub genes including XTH23 for cell wall assembly and two CPK28 genes in PTI pathways. On microRNA, we identified 17 DEMs and predicted 15 miRNA-target pairs (DEM-DEG). We validated two pairs (miR395-APS4 and miR160-ARF) by degradome sequencing. We investigated the miR395-APS4 pair by CRISPR-Cas9 mediated gene editing, the result showed that knocking-out APS4 could lead to elevated clubroot resistance in B. napus. In summary, the data acquired on transcriptional response and microRNA as well as target genes provide future direction especially gene candidates for genetic improvement of clubroot resistance on Brassica species.

Fine mapping and candidate gene analysis of CRA8.1.6, which confers clubroot resistance in turnip (Brassica rapa ssp. rapa)

Clubroot disease poses a significant threat to Brassica crops, necessitating ongoing updates on resistance gene sources. In F2 segregants of the clubroot-resistant inbred line BrT18-6-4-3 and susceptible DH line Y510, the genetic analysis identified a single dominant gene responsible for clubroot resistance. Through bulk segregant sequencing analysis and kompetitive allele-specific polymerase chain reaction assays, CRA8.1.6 was mapped within 110 kb (12,255-12,365 Mb) between markers L-CR11 and L-CR12 on chromosome A08. We identified B raA08g015220.3.5C as the candidate gene of CRA8.1.6. Upon comparison with the sequence of disease-resistant material BrT18-6-4-3, we found 249 single-nucleotide polymorphisms, seven insertions, six deletions, and a long terminal repeat (LTR) retrotransposon (5,310 bp) at 909 bp of the first intron. However, the LTR retrotransposon was absent in the coding sequence of the susceptible DH line Y510. Given the presence of a non-functional LTR insertion in other materials, it showed that the LTR insertion might not be associated with susceptibility. Sequence alignment analysis revealed that the fourth exon of the susceptible line harbored two deletions and an insertion, resulting in a frameshift mutation at 8,551 bp, leading to translation termination at the leucine-rich repeat domain's C-terminal in susceptible material. Sequence alignment of the CDS revealed a 99.4% similarity to Crr1a, which indicate that CRA8.1.6 is likely an allele of the Crr1a gene. Two functional markers, CRA08-InDel and CRA08-KASP1, have been developed for marker-assisted selection in CR turnip cultivars. Our findings could facilitate the development of clubroot-resistance turnip cultivars through marker-assisted selection.

Proteome- and metabolome-level changes during early stages of clubroot infection in Brassica napus canola

Clubroot is a destructive root disease of canola (Brassica napus L.) caused by Plasmodiophora brassicae Woronin. Despite extensive research into the molecular responses of B. napus to P. brassicae, there is limited information on proteome- and metabolome-level changes in response to the pathogen, especially during the initial stages of infection. In this study, we have investigated the proteome- and metabolome- level changes in the roots of clubroot-resistant (CR) and -susceptible (CS) doubled-haploid (DH) B. napus lines, in response to P. brassicae pathotype 3H at 1-, 4-, and 7-days post-inoculation (DPI). Root proteomes were analyzed using nanoflow liquid chromatography coupled with tandem mass spectrometry (nano LC-MS/MS). Comparisons of pathogen-inoculated and uninoculated root proteomes revealed 2515 and 1556 differentially abundant proteins at one or more time points (1-, 4-, and 7-DPI) in the CR and CS genotypes, respectively. Several proteins related to primary metabolites (e.g., amino acids, fatty acids, and lipids), secondary metabolites (e.g., glucosinolates), and cell wall reinforcement-related proteins [e.g., laccase, peroxidases, and plant invertase/pectin methylesterase inhibitors (PInv/PMEI)] were identified. Eleven nucleotides and nucleoside-related metabolites, and eight fatty acids and sphingolipid-related metabolites were identified in the metabolomics study. To our knowledge, this is the first report of root proteome-level changes and associated alterations in metabolites during the early stages of P. brassicae infection in B. napus.

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.

A CRISPR/Cas9-based vector system enables the fast breeding of selection-marker-free canola with Rcr1-rendered clubroot resistance

Breeding for disease resistance in major crops is of crucial importance for global food security and sustainability. However, common biotechnologies such as traditional transgenesis or genome editing do not provide an ideal solution, whereas transgenic crops free of selection markers such as cisgenic/intragenic crops might be suitable. In this study, after cloning and functional verification of the Rcr1 gene for resistance to clubroot (Plasmodiophora brassicae), we confirmed that the genes Rcr1, Rcr2, Rcr4, and CRa from Brassica rapa crops and the resistance gene from B. napus oilseed rape cv. 'Mendel' on chromosome A03 were identical in their coding regions. We also determined that Rcr1 has a wide distribution in Brassica breeding materials and renders potent resistance against multiple representative clubroot strains in Canada. We then modified a CRISPR/Cas9-based cisgenic vector system and found that it enabled the fast breeding of selection-marker-free transgenic crops with add-on traits, with selection-marker-free canola (B. napus) germplasms with Rcr1-rendered stable resistance to clubroot disease being successfully developed within 2 years. In the B. napus background, the intragenic vector system was able to remove unwanted residue sequences from the final product with high editing efficiency, and off-target mutations were not detected. Our study demonstrates the potential of applying this breeding strategy to other crops that can be transformed by Agrobacterium. Following the streamlined working procedure, intragenic germplasms can be developed within two generations, which could significantly reduce the breeding time and labor compared to traditional introgression whilst still achieving comparable or even better breeding results.

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

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

Development of SNP markers linked to purple blotch resistance for marker-assisted selection in onion (Allium cepa L.) breeding

Purple blotch (PB), caused by Alternaria porri (Ellis) Cifferi, is one of the most destructive diseases of onion worldwide. Rapid development and deployment of resistant onion varieties is the most effective approach to control this disease. A single dominant gene, ApR1 was previously linked to PB resistance in onion cultivar 'Arka Kalyan'. In this study, an advanced RIL population derived from a cross between the resistant (Arka Kalyan) and susceptible (Agrifound Rose) cultivar of onion was used to fine map the resistant locus with SNP markers. Twenty plants from the RIL population, ten each with disease resistance and susceptibility trait, were subjected to restriction site-associated DNA sequencing (RAD-Seq) and generated 7388 single nucleotide polymorphisms (SNPs). Correlation analysis between marker genotypes and PB disease phenotype on the 20 plants identified 27 SNPs as candidate markers linked to ApR1 gene for PB resistance. Six candidate SNPs were converted to Kompetitive Allele-Specific PCR (KASP) markers designated as ApRsnip5, ApRsnip8, ApRsnip14, ApRsnip21, ApRsnip23 and ApRsnip25. Marker-trait association based on disease phenotyping and KASP genotyping data on 153 RILs confirmed that all six KASP markers were tightly associated with ApR1 gene within the genetic distance of 1.3 CentiMorgan (cM). ApRsnip14 co-segregated with the ApR1 locus. Further, the six KASP markers were tested on 27 onion lines with different genetic backgrounds. ApRsnip14, ApRsnip21, ApRsnip5 and ApRsnip23 not only showed the correct resistance allele in 3 resistance genotypes, but also clustered together in the remaining 24 susceptible lines. Alternatively, ApRsnip8 and ApRsnip25 exhibited false positives in two onion lines which do not have the R-gene. Overall, our results suggest that ApRsnip14 and ApRsnip23 with their close linkage to ApR1 locus and greater applicability on breeding germplasm are recommended in marker-assisted selection for PB resistance in onion breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03562-7.

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.

R gene triplication confers European fodder turnip with improved clubroot resistance

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

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.

Identification of Two Major QTLs in Brassica napus Lines With Introgressed Clubroot Resistance From Turnip Cultivar ECD01

Plasmodiophora brassicae causes clubroot disease in brassica crops worldwide. Brassica rapa, a progenitor of Brassica napus (canola), possesses important sources for resistance to clubroot. A doubled haploid (DH) population consisting of 84 DH lines were developed from a Backcross2 (BC₂) plant through an interspecific cross of B. rapa turnip cv. ECD01 (resistant, R) with canola line DH16516 (susceptible, S) and then backcrossed with DH16516 as the recurrent parent. The DH lines and their parental lines were tested for resistance to four major pathotypes (3A, 3D, 3H, and 5X) of P. brassicae identified from canola. The R:S segregation ratio for pathotype 3A was 1:3, and 3:1 for pathotypes 3D, 3H, and 5X. From genotyping by sequencing (GBS), a total of 355.3 M short reads were obtained from the 84 DH lines, ranging from 0.81 to 11.67 M sequences per line. The short reads were aligned into the A-genome of B. napus "Darmor-bzh" version 4.1 with a total of 260 non-redundant single-nucleotide polymorphism (SNP) sites. Two quantitative trait loci (QTLs), Rcr10 ^ECD01 and Rcr9 ^ECD01 , were detected for the pathotypes in chromosomes A03 and A08, respectively. Rcr10 ^ECD01 and Rcr9 ^ECD01 were responsible for resistance to 3A, 3D, and 3H, while only one QTL, Rcr9 ^ECD01 , was responsible for resistance to pathotype 5X. The logarithm of the odds (LOD) values, phenotypic variation explained (PVE), additive (Add) values, and confidence interval (CI) from the estimated QTL position varied with QTL, with a range of 5.2-12.2 for LOD, 16.2-43.3% for PVE, 14.3-25.4 for Add, and 1.5-12.0 cM for CI. The presence of the QTLs on the chromosomes was confirmed through the identification of the percentage of polymorphic variants using bulked-segregant analysis. There was one gene encoding a disease resistance protein and 24 genes encoding proteins with function related to plant defense response in the Rcr10 ^ECD01 target region. In the Rcr9 ^ECD01 region, two genes encoded disease resistance proteins and 10 genes encoded with defense-related function. The target regions for Rcr10 ^ECD01 and Rcr9 ^ECD01 in B. napus were homologous to the 11.0-16.0 Mb interval of chromosome A03 and the 12.0-14.5 Mb interval of A08 in B. rapa "Chiifu" reference genome, respectively.

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

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

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

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

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.

Transcriptome and Coexpression Network Analyses Reveal Hub Genes in Chinese Cabbage (Brassica rapa L. ssp. pekinensis) During Different Stages of Plasmodiophora brassicae Infection

Clubroot, caused by the soil-borne protist Plasmodiophora brassicae, is one of the most destructive diseases of Chinese cabbage worldwide. However, the clubroot resistance mechanisms remain unclear. In this study, in both clubroot-resistant (DH40R) and clubroot-susceptible (DH199S) Chinese cabbage lines, the primary (root hair infection) and secondary (cortical infection) infection stages started 2 and 5 days after inoculation (dai), respectively. With the extension of the infection time, cortical infection was blocked and complete P. brassica resistance was observed in DH40R, while disease scales of 1, 2, and 3 were observed at 8, 13, and 22 dai in DH199S. Transcriptome analysis at 0, 2, 5, 8, 13, and 22 dai identified 5,750 relative DEGs (rDEGs) between DH40R and DH199S. The results indicated that genes associated with auxin, PR, disease resistance proteins, oxidative stress, and WRKY and MYB transcription factors were involved in clubroot resistance regulation. In addition, weighted gene coexpression network analysis (WGCNA) identified three of the modules whose functions were highly associated with clubroot-resistant, including ten hub genes related to clubroot resistance (ARF2, EDR1, LOX4, NHL3, NHL13, NAC29, two AOP1, EARLI 1, and POD56). These results provide valuable information for better understanding the molecular regulatory mechanism of Chinese cabbage clubroot resistance.

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

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

Genetic Mapping by Sequencing More Precisely Detects Loci Responsible for Anaerobic Germination Tolerance in Rice

Direct seeded rice (DSR) is a mainstay for planting rice in the Americas, and it is rapidly becoming more popular in Asia. It is essential to develop rice varieties that are suitable for this type of production system. ASD1, a landrace from India, possesses several traits desirable for direct-seeded fields, including tolerance to anaerobic germination (AG). To map the genetic basis of its tolerance, we examined a population of 200 F_2:3 families derived from a cross between IR64 and ASD1 using the restriction site-associated DNA sequencing (RAD-seq) technology. This genotyping platform enabled the identification of 1921 single nucleotide polymorphism (SNP) markers to construct a high-resolution genetic linkage map with an average interval of 0.9 cM. Two significant quantitative trait loci (QTLs) were detected on chromosomes 7 and 9, qAG7 and qAG9, with LOD scores of 7.1 and 15.0 and R² values of 15.1 and 29.4, respectively. Here, we obtained more precise locations of the QTLs than traditional simple sequence repeat and low-density SNP genotyping methods and may help further dissect the genetic factors of these QTLs.

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

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

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

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

Genome-wide SNP discovery from Darjeeling tea cultivars - their functional impacts and application toward population structure and trait associations

Genotyping by sequencing and identification of functionally relevant nucleotide variations in crop accessions are the key steps to unravel genetic control of desirable traits. Elite cultivars of Darjeeling tea were undergone SNP genotyping by double-digest restriction-site associated DNA sequencing method. This study reports a set of 54,206 high-quality SNP markers discovered from ~10.4 GB sequence data, encompassing 15 chromosomes of the reference tea genome. Genetic relatedness among the accessions conforms to the analyses of Bayesian clustering, UPGMA, and PCoA methods. Genomic positions of the discovered SNPs and their putative effect on annotated genes designated a thoughtful understanding of their functional aspects in tea system biology. A group of 95 genes was identified to be affected by high impact variants. Genome-wide association analyses of 21 agronomic and biochemical phenotypes resulted in trait-linked polymorphic loci with strong confidence (p < 0.05 and 0.001).

Frontiers in Dissecting and Managing Brassica Diseases: From Reference-Based RGA Candidate Identification to Building Pan-RGAomes

The Brassica genus contains abundant economically important vegetable and oilseed crops, which are under threat of diseases caused by fungal, bacterial and viral pathogens. Resistance gene analogues (RGAs) are associated with quantitative and qualitative disease resistance and the identification of candidate RGAs associated with disease resistance is crucial for understanding the mechanism and management of diseases through breeding. The availability of Brassica genome assemblies has greatly facilitated reference-based quantitative trait loci (QTL) mapping for disease resistance. In addition, pangenomes, which characterise both core and variable genes, have been constructed for B. rapa, B. oleracea and B. napus. Genome-wide characterisation of RGAs using conserved domains and motifs in reference genomes and pangenomes reveals their clustered arrangements and presence of structural variations. Here, we comprehensively review RGA identification in important Brassica genome and pangenome assemblies. Comparison of the RGAs in QTL between resistant and susceptible individuals allows for efficient identification of candidate disease resistance genes. However, the reference-based QTL mapping and RGA candidate identification approach is restricted by the under-represented RGA diversity characterised in the limited number of Brassica assemblies. The species-wide repertoire of RGAs make up the pan-resistance gene analogue genome (pan-RGAome). Building a pan-RGAome, through either whole genome resequencing or resistance gene enrichment sequencing, would effectively capture RGA diversity, greatly expanding breeding resources that can be utilised for crop improvement.

Revised Taxonomy and Expanded Biodiversity of the Phytomyxea (Rhizaria, Endomyxa)

Phytomyxea (phytomyxids) is a group of obligate biotrophic pathogens belonging to the Rhizaria. Some phytomyxids are well studied and include known plant pathogens such as Plasmodiophora brassicae, the causal agent of clubroot disease. Despite this economic importance, the taxonomy and biodiversity of this group are largely cryptic, with many species described in the premolecular area. Some of these species were key for establishing the morphotaxonomic concepts that define most genera to this day, but systematic efforts to include and integrate those species into molecular studies are still lacking. The aim of this study was to expand our understanding of phytomyxid biodiversity in terrestrial environments. Thirty-eight environmental samples from habitats in which novel and known diversity of Phytomyxea was expected were analysed. We were able to generate 18S rRNA sequences from Ligniera verrucosa, a species which is well defined based on ultrastructure. Phylogenetic analyses of the collected sequences rendered the genera Lignera, Plasmodiophora and Spongospora polyphyletic, and identified two novel and apparently diverse lineages (clade 17, clade 18). Based on these findings and on data from previous studies, we formally establish the new genera Pseudoligniera n. gen. for L. verrucosa,Hillenburgia n. gen. for Spongospora nasturtii and revert Plasmodiophora diplantherae to its original name Ostenfeldiella diplantherae.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

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

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

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.

Expression and Role of Response Regulating, Biosynthetic and Degrading Genes for Cytokinin Signaling during Clubroot Disease Development

The obligate biotroph Plasmodiophora brassicae causes clubroot disease in oilseeds and vegetables of the Brassicaceae family, and cytokinins play a vital role in clubroot formation. In this study, we examined the expression patterns of 17 cytokinin-related genes involved in the biosynthesis, signaling, and degradation in Chinese cabbage inoculated with the Korean pathotype group 4 isolate of P. brassicae, Seosan. This isolate produced the most severe clubroot symptoms in Chinese cabbage cultivar “Bullam-3-ho” compared to three other Korean geographical isolates investigated. BrIPT1, a cytokinin biosynthesis gene, was induced on Day 1 and Day 28 in infected root tissues and the upregulation of this biosynthetic gene coincided with the higher expression of the response regulators BrRR1, on both Days and BrRR6 on Day 1 and 3. BrRR3 and 4 genes were also induced during gall enlargement on Day 35 in leaf tissues. The BrRR4 gene, which positively interact with phytochrome B, was consistently induced in leaf tissues on Day 1, 3, and 14 in the inoculated plants. The cytokinin degrading gene BrCKX3-6 were induced on Day 14, before gall initiation. BrCKX2,3,6 were induced until Day 28 and their expression was downregulated on Day 35. This insight improves our current understanding of the role of cytokinin signaling genes in clubroot disease development.

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.

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

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

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

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