Abstract
Clubroot, a soil-borne disease caused by Plasmodiophora brassicae, is one of the most destructive diseases of Brassica oleracea all over the world. However, the mechanism of clubroot resistance remains unclear. In this research, transcriptome sequencing was conducted on root samples from both resistant (R) and susceptible (S) B. oleracea plants infected by P. brassicae. Then the comparative analysis was carried out between the R and S samples at different time points during the infection stages to reveal clubroot resistance related pathways and candidate genes. Compared with 0 days after inoculation, a total of 4991 differential expressed genes were detected from the S pool, while only 2133 were found from the R pool. Gene function enrichment analysis found that the effector-triggered immunity played a major role in the R pool, while the pathogen-associated molecular pattern triggered immune response was stronger in the S pool. Simultaneously, candidate genes were identified through weighted gene co-expression network analysis, with Bol010786 (CNGC13) and Bol017921 (SD2-5) showing potential for conferring resistance to clubroot. The findings of this research provide valuable insights into the molecular mechanisms underlying clubroot resistance and present new avenues for further research aimed at enhancing the clubroot resistance of B. oleracea through breeding.