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 F1 plants were resistant to clubroot. There was a 1:1 ratio of R:S in the BC1 and 3R:1S in the F2, 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

Studies on the temporal, structural, and interacting features of the clubroot resistance gene Rcr1 using CRISPR/Cas9-based systems

Clubroot disease is a severe threat to Brassica crops globally, particularly in western Canada. Genetic resistance, achieved through pyramiding clubroot resistance (CR) genes with different modes of action, is the most important strategy for managing the disease. However, studies on the CR gene functions are quite limited. In this study, we have conducted investigations into the temporal, structural, and interacting features of a newly cloned CR gene, Rcr1, using CRISPR/Cas9 technology. For temporal functionality, we developed a novel CRISPR/Cas9-based binary vector, pHHIGR-Hsp18.2, to deliver Rcr1 into a susceptible canola line (DH12075) and observed that early expression of Rcr1 is critical for conferring resistance. For structural functionality, several independent mutations in specific domains of Rcr1 resulted in loss-of-function, highlighting their importance for CR phenotype. In the study of the interacting features of Rcr1, a cysteine protease gene and its homologous allele in canola were successfully disrupted via CRISPR/Cas9 as an interacting component with Rcr1 protein, resulting in the conversion from clubroot resistant to susceptible in plants carrying intact Rcr1. These results indicated an indispensable role of these two cysteine proteases in Rcr1-mediated resistance response. This study, the first of its kind, provides valuable insights into the functionality of Rcr1. Further, the new vector pHHIGR-Hsp18.2 demonstrated an inducible feature on the removal of add-on traits, which should be useful for functional genomics and other similar research in brassica crops

A CRISPR/Cas9-Based System with Controllable Auto-Excision Feature Serving Cisgenic Plant Breeding and Beyond

Transgenic or genetically modified crops have great potential in modern agriculture but still suffer from heavy regulations worldwide due to biosafety concerns. As a promising alternative route, cisgenic crops have received higher public acceptance and better reviews by governing authorities. To serve the purpose of cisgenic plant breeding, we have developed a CRISPR/Cas9-based vector system, which is capable of delivering target gene-of-interest (GOI) into recipient plants while removing undesired genetic traces in the plants. The new system features a controllable auto-excision feature, which is realized by a core design of embedded multi-clonal sequence and the use of inducible promoters controlling the expression of Cas9 nuclease. In the current proof-of-concept study in Arabidopsis thaliana (L.) Heynh., we have successfully incorporated a GOI into the plant and removed the selection marker and CRISPR/Cas9 components from the final product. Following the designed workflow, we have demonstrated that novel cisgenic plant germplasms with desired traits could be developed within one to two generations. Further characterizations of the vector system have shown that heat treatment at 37 °C could significantly improve the editing efficiency (up to 100%), and no off-target mutations were identified in the Arabidopsis background. This novel vector system is the first CRISPR/Cas9-based genome editing tool for cisgenic plant breeding and should prove powerful for other similar applications in the bright future of precision molecular breeding

Genome-wide transcriptome reveals mechanisms underlying Rlm1-mediated blackleg resistance on canola

Abstract Genetic resistance to blackleg (Leptosphaeria maculans, Lm) of canola (Brassica napus, Bn) has been extensively studied, but the mechanisms underlying the host–pathogen interaction are still not well understood. Here, a comparative transcriptome analysis was performed on a resistant doubled haploid Bn line carrying the resistance gene Rlm1 following inoculation with a virulent (avrLm1) or avirulent (AvrLm1) Lm isolate on cotyledons. A total of 6999 and 3015 differentially expressed genes (DEGs) were identified, respectively, in inoculated local tissues with compatible (susceptible) and incompatible (resistant) interactions. Functional enrichment analysis found several biological processes, including protein targeting to membrane, ribosome and negative regulation of programmed cell death, were over-represented exclusively among up-regulated DEGs in the resistant reaction, whereas significant enrichment of salicylic acid (SA) and jasmonic acid (JA) pathways observed for down-regulated DEGs occurred only in the susceptible reaction. A heat-map analysis showed that both biosynthesis and signaling of SA and JA were induced more significantly in the resistant reaction, implying that a threshold level of SA and JA signaling is required for the activation of Rlm1-mediated resistance. Co-expression network analysis revealed close correlation of a gene module with the resistance, involving DEGs regulating pathogen-associated molecular pattern recognition, JA signaling and transcriptional reprogramming. Substantially fewer DEGs were identified in mock-inoculated (control) cotyledons, relative to those in inoculated local tissues, including those involved in SA pathways potentially contributing to systemic acquired resistance (SAR). Pre-inoculation of cotyledon with either an avirulent or virulent Lm isolate, however, failed to induce SAR on remote tissues of same plant despite elevated SA and PR1 protein. This study provides insights into the molecular mechanism of Rlm1-mediated resistance to blackleg

Single R Gene Introgression Lines for Accurate Dissection of the Brassica - Leptosphaeria Pathosystem

Seven blackleg resistance (R) genes (Rlm1, Rlm2, Rlm3, Rlm4, LepR1, LepR2 & LepR3) were each introgressed into a common susceptible B. napus doubled-haploid (DH) line through reciprocal back-crossing, producing single-R gene introgression lines (ILs) for use in the pathological and molecular study of Brassica – Leptosphaeria interactions. The genomic positions of the R genes were defined through molecular mapping and analysis with transgenic L. maculans isolates was used to confirm the identity of the introgressed genes where possible. Using L. maculans isolates of contrasting avirulence gene (Avr) profiles, we preformed extensive differential pathology for phenotypic comparison of the ILs to other B. napus varieties, demonstrating the ILs can provide for the accurate assessment of Avr - R gene interactions by avoiding non-Avr dependant alterations to resistance responses which can occur in some commonly used B. napus varieties. Whole-genome SNP-based assessment allowed us to define the donor parent introgressions in each IL and provide a strong basis for comparative molecular dissection of the pathosystem

Evaluating Changes in Cell-Wall Components Associated with Clubroot Resistance Using Fourier Transform Infrared Spectroscopy and RT-PCR

Clubroot disease is a serious threat to canola production in western Canada and many parts of the world. Rcr1 is a clubroot resistance (CR) gene identified recently and its molecular mechanisms in mediating CR have been studied using several omics approaches. The current study aimed to characterize the biochemical changes in the cell wall of canola roots connecting to key molecular mechanisms of this CR gene identified in prior studies using Fourier transform infrared (FTIR) spectroscopy. The expression of nine genes involved in phenylpropanoid metabolism was also studied using qPCR. Between susceptible (S) and resistance (R) samples, the most notable biochemical changes were related to an increased biosynthesis of lignin and phenolics. These results were supported by the transcription data on higher expression of BrPAL1. The up-regulation of PAL is indicative of an inducible defence response conferred by Rcr1; the activation of this basal defence gene via the phenylpropanoid pathway may contribute to clubroot resistance conferred by Rcr1. The data indicate that several cell-wall components, including lignin and pectin, may play a role in defence responses against clubroot. Principal components analysis of FTIR data separated non-inoculated samples from inoculated samples, but not so much between inoculated S and inoculated R samples. It is also shown that FTIR spectroscopy can be a useful tool in studying plant-pathogen interaction at cellular levels

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/CRaM and Crr1rutb) 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 + Crr1rutb or CRaM + Crr1rutb, showed partial resistance. Selected lines were inoculated with a field pathotype X population (L-G3) at 5 × 106 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

Canola with Stacked Genes Shows Moderate Resistance and Resilience against a Field Population of <i>Plasmodiophora brassicae</i> (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/CRaM and Crr1rutb) 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 + Crr1rutb or CRaM + Crr1rutb, showed partial resistance. Selected lines were inoculated with a field pathotype X population (L-G3) at 5 × 106 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

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 &gt;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