RenSeq and whole genome sequencing uncover allelic diversity of clubroot resistance genes in commercial breeding canola lines

Abstract

Clubroot disease, caused by the obligate biotrophic pathogen Plasmodiophora brassicae, is a major threat to canola (Brassica napus) production worldwide. Clubroot-resistant (CR) cultivars remain the most effective disease-management strategy, but the genetic basis of resistance in commercial canola remains poorly understood because many resistance sources are proprietary and associated genotypic information is rarely accessible. Although nucleotide-binding leucine-rich repeat (NLR) immune receptors account for most cloned CR genes, no pan-NLRome has incorporated CR lines used in commercial canola breeding. Here, we combined whole-genome sequencing and resistance gene enrichment sequencing (RenSeq) to assemble and annotate the NLR repertoires of five homozygous CR inbred lines (IH1–IH5) used for commercial breeding and displaying contrasting resistance profiles against predominant Canadian P. brassicae pathotypes. We integrated these NLRomes with the susceptible cultivar Westar to construct a comparative pan-NLRome for canola. Across the five CR lines, total NLR content was highly conserved, ranging from 504 to 517 genes, with TIR-NLRs representing the predominant class. C-JID-containing TIR-NLRs accounted for more than 30% of each NLR repertoire, and integrated-domain analysis identified conserved and genotype-specific NLR-IDs, including previously unreported domains in IH4. Pan-NLRome analysis resolved 366 NLR orthogroups (OGs), 60.7% of which were core, and identified resistant-line-enriched OGs absent from Westar as candidate CR-associated loci. Unexpectedly, a homolog of the functionally characterized CR gene, CRa, was detected in five CR lines. Moreover, a homolog of another CR gene, Crr1a, was detected in both resistant and susceptible lines, indicating that the presence/absence of a gene alone does not predict resistance. Instead, structural variation affecting LRR and C-JID regions suggests that allele-level diversity within conserved NLR loci contributes to CR-associated variation, with implications for allele-specific marker development and durable CR deployment.