Plasmodiophora brassicae, an obligate biotrophic pathogen-causing clubroot disease, can seriously affect Brassica crops worldwide, especially Chinese cabbage. Understanding the transcriptome and metabolome profiling changes during the infection of P. brassicae will provide key insights in understanding the defense mechanism in Brassica crops. In this study, we estimated the phytohormones using targeted metabolome assays and transcriptomic changes using RNA sequencing (RNA-seq) in the roots of resistant (BrT24) and susceptible (Y510-9) plants at 0, 3, 9, and 20 days after inoculation (DAI) with P. brassicae. Differentially expressed genes (DEGs) in resistant vs. susceptible lines across different time points were identified. The weighted gene co-expression network analysis of the DEGs revealed six pathways including “Plant–pathogen interaction” and “Plant hormone signal transduction” and 15 hub genes including pathogenic type III effector avirulence factor gene (RIN4) and auxin-responsive protein (IAA16) to be involved in plants immune response. Inhibition of Indoleacetic acid, cytokinin, jasmonate acid, and salicylic acid contents and changes in related gene expression in R-line may play important roles in regulation of clubroot resistance (CR). Based on the combined metabolome profiling and hormone-related transcriptomic responses, we propose a general model of hormone-mediated defense mechanism. This study definitely enhances our current understanding and paves the way for improving CR in Brassica rapa.

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.

Chinese cabbage (Brassica rapa ssp. pekinensis) is one of the most important vegetables in Asia and is cultivated across the world. Ogura-type cytoplasmic male sterility (Ogura-CMS) has been widely used in the hybrid breeding industry for Chinese cabbage and many other cruciferous vegetables. Although, the cause of Ogura-CMS has been localized to the orf138 locus in the mitochondrial genome, however, the mechanism by which nuclear genes respond to the mutation of the mitochondrial orf138 locus is unclear. In this study, a series of whole genome small RNA, degradome and transcriptome analyses were performed on both Ogura-CMS and its maintainer Chinese cabbage buds using deep sequencing technology. A total of 289 known miRNAs derived from 69 families (including 23 new families first reported in B. rapa) and 426 novel miRNAs were identified. Among these novel miRNAs, both 3-p and 5-p miRNAs were detected on the hairpin arms of 138 precursors. Ten known and 49 novel miRNAs were down-regulated, while one known and 27 novel miRNAs were up-regulated in Ogura-CMS buds compared to the fertile plants. Using degradome analysis, a total of 376 mRNAs were identified as targets of 30 known miRNA families and 100 novel miRNAs. A large fraction of the targets were annotated as reproductive development related. Our transcriptome profiling revealed that the expression of the targets was finely tuned by the miRNAs. Two novel miRNAs were identified that were specifically highly expressed in Ogura-CMS buds and sufficiently suppressed two pollen development essential genes: sucrose transporter SUC1 and H+-ATPase 6. These findings provide clues for the contribution of a potential miRNA regulatory network to bud development and pollen engenderation. This study contributes new insights to the communication between the mitochondria and chromosome and takes one step toward filling the gap in the regulatory network from the orf138 locus to pollen abortion in Ogura-CMS plants from a miRNA perspective.

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34Brassica rapaaccessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies inBrassica rapa.

Flower color is an important trait in plants. However, genes responsible for the white flower trait in Chinese cabbage are rarely reported. In this study, we constructed an F₂ population derived from the Y640-288 (white flower) and Y641-87 (yellow flower) lines for the fine mapping of the white flower gene BrWF3 in Chinese cabbage. Genetic analysis indicated that BrWF3 was controlled by a single recessive gene. Using BSA-seq and KASP assays, BrWF3 was fine-mapped to an interval of 105.6 kb. Functional annotation, expression profiling, and sequence variation analyses confirmed that the AtPES2 homolog, Bra032957, was the most likely candidate gene for BrWF3. Carotenoid profiles and transmission electron microscopy analysis suggested that BrWF3 might participate in the production of xanthophyll esters (particularly violaxanthin esters), which in turn disrupt chromoplast development and the formation of plastoglobules (PGs). A SNP deletion in the third exon of BrWF3 caused the loss of protein function, and interfered with the normal assembly of PGs, which was associated with reduced expression levels of genes involved in carotenoid metabolism. Furthermore, we developed and validated the functional marker TXBH83 for BrWF3. Our results provide insight into the molecular mechanism underlying flower color pigmentation and reveal valuable information for marker-assisted selection (MAS) breeding in Chinese cabbage.