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.