Reddish purple Chinese cabbage (RPCC) is a popular variety of Brassica rapa (AA = 20). It is rich in anthocyanins, which have many health benefits. We detected novel anthocyanins including cyanidin 3-(feruloyl) diglucoside-5-(malonoyl) glucoside and pelargonidin 3-(caffeoyl) diglucoside-5-(malonoyl) glucoside in RPCC. Analyses of transcriptome data revealed 32,395 genes including 3345 differentially expressed genes (DEGs) between 3-week-old RPCC and green Chinese cabbage (GCC). The DEGs included 218 transcription factor (TF) genes and some functionally uncharacterized genes. Sixty DEGs identified from the transcriptome data were analyzed in 3-, 6- and 9-week old seedlings by RT-qPCR, and 35 of them had higher transcript levels in RPCC than in GCC. We detected cis-regulatory motifs of MYB, bHLH, WRKY, bZIP and AP2/ERF TFs in anthocyanin biosynthetic gene promoters. A network analysis revealed that MYB75, MYB90, and MYBL2 strongly interact with anthocyanin biosynthetic genes. Our results show that the late biosynthesis genes BrDFR, BrLDOX, BrUF3GT, BrUGT75c1-1, Br5MAT, BrAT-1, BrAT-2, BrTT19-1, and BrTT19-2 and the regulatory MYB genes BrMYB90, BrMYB75, and BrMYBL2-1 are highly expressed in RPCC, indicative of their important roles in anthocyanin biosynthesis, modification, and accumulation. Finally, we propose a model anthocyanin biosynthesis pathway that includes the unique anthocyanin pigments and genes specific to RPCC.

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

The B-box zinc-finger transcription factors are important for plant growth, development, and various physiological processes such as photomorphogenesis, light signaling, and flowering, as well as for several biotic and abiotic stress responses. However, there is relatively little information available regarding Brassica B-box genes and their expression. In this study, we identified 51, 52, and 101 non-redundant genes encoding B-box proteins in Brassica rapa (BrBBX genes), B. oleracea (BoBBX genes), and B. napus (BnBBX genes), respectively. A whole-genome identification, characterization, and evolutionary analysis (synteny and orthology) of the B-box gene families in the diploid species B. rapa (A genome) and B. oleracea (C genome) and in the allotetraploid species B. napus (AC genome) revealed segmental duplications were the major contributors to the expansion of the BrassicaBBX gene families. The BrassicaBBX genes were classified into five subgroups according to phylogenetic relationships, gene structures, and conserved domains. Light-responsive cis-regulatory elements were detected in many of the BBX gene promoters. Additionally, BrBBX expression profiles in different tissues and in response to various abiotic stresses (heat, cold, salt, and drought) or hormones (abscisic acid, methyl jasmonate, and gibberellic acid) were analyzed by qRT-PCR. The data indicated that many B-box genes (e.g., BrBBX13, BrBBX15, and BrBBX17) may contribute to plant development and growth as well as abiotic stress tolerance. Overall, the identified BBX genes may be useful as functional genetic markers for multiple stress responses and plant developmental processes.

Radish flower color, bolting time, and flowering time are important traits for attracting certain pollinators and affect fleshy root quality. In this study, an analysis of the anthocyanidins in radish flowers by high-performance liquid chromatography revealed that differences in the cyanidin content are likely to be associated with the variability in radish flower colors (i.e., purple and white petals). A quantitative trait loci (QTL) analysis identified nine QTLs on three Raphanus sativus linkage groups. Three QTLs—qRFC1, qRBT1, and qRFT1—which were consistently detected and explained a high proportion of the observed variation (10.30% to 34.57%), were considered as the major QTLs responsible for flower color, bolting time, and flowering time, respectively. A total of 16 and 11 candidate genes within the major QTL regions for flower color and bolting/flowering times, respectively, were preliminarily annotated. Six genes (Rs018140, Rs018950, Rs019220, Rs020080, Rs020590, and Rs021450) related to flower color were differentially expressed in the parental lines. On the basis of nucleotide and amino acid sequence diversity between the parental lines, Rs314940, Rs315000, Rs315310, and Rs315960 were identified as candidate genes mediating the radish bolting and flowering times. This study revealed the genetic complexity of the radish flower color, bolting time, and flowering time traits. The identified candidate genes in the QTL regions may be useful for radish breeding programs and also for functional characterization in radish.

Clubroot resistance is an economically important trait in Brassicaceae crops. Although many quantitative trait loci (QTLs) for clubroot resistance have been identified in Brassica, disease-related damage continues to occur owing to differences in host variety and constant pathogen variation. Here, we investigated the inheritance of clubroot resistance in a double haploid population developed by crossing clubroot resistant and susceptible lines “09CR500” and “09CR501”, respectively. The resistance of “09CR500” to Plasmodiophora brassicae pathotype “Banglim” was controlled as a single dominant gene, with the segregation of resistance and susceptibility being nearly 1:1. PbBrA08^Banglim was identified as having a logarithm of odds value of 7.9–74.8, and a phenotypic variance of 26.0–97.1% with flanking marker “09CR.11390652” in A08. After aligning QTL regions to the B. rapa reference genome, 11 genes were selected as candidates. PbBrA08^Banglim was located near Crr1, CRs, and Rcr9 loci, but differences were validated by marker analysis, gene structural variations, and gene expression levels, as well as phenotypic responses to the pathotype. Genotyping using the “09CR.11390652” marker accurately distinguished the Banglim-resistance phenotypes in the double haploid population. Thus, the developed marker will be useful in Brassica breeding programs, marker-assisted selection, and gene pyramiding to identify and develop resistant cultivars.

Protein post-translational modification by phosphorylation is essential for the activity and stability of proteins in higher plants and underlies their responses to diverse stimuli. There are more than 300 leucine-rich repeat receptor-like kinases (LRR-RLKs), a major group of receptor-like kinases (RLKs) that plays an important role in growth, development, and biotic stress responses in higher plants. To analyze auto- and transphosphorylation patterns and kinase activities in vitro, 43 full-length complementary DNA (cDNA) sequences were cloned from genes encoding LRR-RLKs. Autophosphorylation activity was found in the cytoplasmic domains (CDs) of 18 LRR-RLKs; 13 of these LRR-RLKs with autophosphorylation activity showed transphosphorylation in Escherichia coli. BRI1-Associated Receptor Kinase (BAK1), which is critically involved in the brassinosteroid and plant innate immunity signal transduction pathways, showed strong auto- and transphosphorylation with multi-specific kinase activity within 2 h of induction of Brassica oleraceae BAK1-CD (BoBAK1-CD) in E. coli; moreover, the carboxy-terminus of LRR-RLKs regulated phosphorylation and kinase activity in Arabidopsis thaliana and vegetative crops.