In Chinese cabbage development the interplay between shoot apex activity and vernalization is pivotal for flowering timing. The intricate relationship between various cell types in the shoot apex meristem and their roles in regulating flowering gene expression in Chinese cabbage is not yet fully understood. A thorough analysis of single-cell types in the Chinese cabbage shoot apex and their influence on flowering genes and vernalization is essential for deeper insight. Our study first established a single-cell transcriptomic atlas of Chinese cabbage after 25 days of non-vernalization. Analyzing 19 602 single cells, we differentiated them into 15 distinct cell clusters using established marker genes. We found that key genes in shoot apex development and flowering were primarily present in shoot meristematic cells (SMCs), companion cells (CCs), and mesophyll cells (MCs). MADS-box protein FLOWERING LOCUS C 2 (BrFLC2), a gene suppressing flowering, was observed in CCs, mirroring patterns found in Arabidopsis. By mapping developmental trajectories of SMCs, CCs, and MCs, we elucidated the evolutionary pathways of crucial genes in shoot apex development and flowering. The creation of a single-cell transcriptional atlas of the Chinese cabbage shoot apex under vernalization revealed distinct alterations in the expression of known flowering genes, such as VERNALIZATION INSENSITIVE 3 (VIN3), VERNALIZATION 1 (VRN1), VERNALIZATION 2 (VRN2), BrFLC, and FLOWERING LOCUS T (FT), which varied by cell type. Our study underscores the transformative impact of single-cell RNA sequencing (scRNA-seq) for unraveling the complex differentiation and vernalization processes in the Chinese cabbage shoot apex. These insights are pivotal for enhancing breeding strategies and cultivation management of this vital vegetable.

In plants, the accumulation of carotenoids can maintain the balance of the photosystem and improve crop nutritional quality. Therefore, the molecular mechanisms underlying carotenoid synthesis and accumulation should be further explored. In this study, carotenoid accumulation differed significantly among parental Brassica rapa. Genetic analysis was carried out using the golden inner leaf ‘1900264′ line and the light−yellow inner leaf ‘1900262′ line, showing that the golden inner leaf phenotype was controlled by a single dominant gene. Using bulked−segregant analysis sequencing, BraA09g007080.3C encoding the ORANGE protein was selected as a candidate gene. Sequence alignment revealed that a 4.67 kb long terminal repeat insertion in the third exon of the BrGOLDEN resulted in three alternatively spliced transcripts. The spatiotemporal expression results indicated that BrGOLDEN might regulate the expression levels of carotenoid−synthesis−related genes. After transforming BrGOLDEN into Arabidopsis thaliana, the seed−derived callus showed that BrGOLDEN_Ins and BrGOLDEN_Del lines presented a yellow color and the BrGOLDEN_Ldel line presented a transparent phenotype. In addition, using the yeast two−hybrid assay, BrGOLDEN_Ins, BrGOLDEN_Ldel, and Brgolden_wt exhibited strong interactions with BrPSY1, but BrGOLDEN_Del did not interact with BrPSY1 in the split−ubiquitin membrane system. In the secondary and 3D structure analysis, BrGOLDEN_Del was shown to have lost the PNFPSFIPFLPPL sequences at the 125 amino acid position, which resulted in the α−helices of BrGOLDEN_Del being disrupted, restricting the formation of the 3D structure and affecting the functions of the protein. These findings may provide new insights into the regulation of carotenoid synthesis in B. rapa.

Chinese cabbage, which is a cold season crop, can still be damaged at an overly low temperature. It is crucial to study the mechanism of the resistance to low temperature of Chinese cabbage. In this study, the Chinese cabbage ‘XBJ’ was used as the material, and nine different low temperatures and control samples were treated. Using RNA-seq and lignin content determination, we analyzed 27 samples, and the stained sections of them were observed. A total of 8845 genes were screened for the WGCNA analysis, yielding 17 modules. The GO and KEGG analyses of the modules was highly associated with a low-temperature treatment. The pathways such as ‘starch and sucrose metabolism’ and ‘plant hormone signal transduction’ were enriched in modules related to low temperature. Interestingly, L-15DAT-associated MEcoral2 was found to have 14 genes related to the ‘lignin biosynthetic process’ in the GO annotation. The combination of the determination of the lignin content and the treatment of the stained sections showed that the lignin content of the low-temperatures samples were indeed higher than that of the control. We further explored the expression changes of the lignin synthesis pathway and various genes and found that low temperature affects the expression changes of most genes in the lignin synthesis pathway, leading to the speculation that the lignin changes at low temperature are a defense mechanism against low temperatures. The 29 BrCOMT gene sequence derived from the RNA-seq was non-conserved, and eight BrCOMT genes were differentially expressed. This study provides a new insight into how lignin is affected by low temperature.

Brassica rapa is one of the most important leafy vegetables worldwide, and has a long history of cultivation. However, it has not been possible to completely control the damage of turnip mosaic virus (TuMV), a serious virus in B. rapa, to production. In this study, the genome-wide identification and expression detection of eIF family genes from B. rapa in response to TuMV resistance were analyzed, including the identification of eIF family genes, chromosomal distribution, three-dimensional (3D) structure and sequence logo analyses, and the expression characterization as well as differential metabolite analysis of eIF family genes in resistant/susceptible lines, which may further prove the whole-genome tripling (WGT) event in B. rapa evolution and provide evidence for the functional redundancy and functional loss of multicopy eIF genes in evolution. A qRT-PCR analysis revealed that the relative expressions of eIF genes in a susceptible line (80461) were higher than those in a resistant line (80124), which may prove that, when TuMV infects host plants, the eIF genes can combine with the virus mRNA 5′ end cap structure and promote the initiation of virus mRNA translation in the susceptible B. rapa line. In addition, the metabolite substances were detected, the differences in metabolites between disease-resistant and disease-susceptible plants were mainly manifested by altered compounds such as flavonoids, jasmonic acid, salicylic acid, ketones, esters, etc., which inferred that the different metabolite regulations of eIF family genes and reveal the resistance mechanisms of eIF genes against TuMV in brassica crops. This study may lay a new theoretical foundation for revealing eIF family gene resistance to TuMV in B. rapa, as well as advancing our understanding of virus–host interactions.

The trichomes of mustard leaves have significance due to their ability to combat unfavorable external conditions and enhance disease resistance. It was demonstrated that the MYB-bHLH-WD40 (MBW) ternary complex consists of MYB, basic Helix-Loop-Helix (bHLH), and WD40-repeat (WD40) family proteins and plays a key role in regulating trichome formation and density. The bHLH gene family, particularly the Myelocytomatosis (MYC) proteins that possess the structural bHLH domain (termed bHLH-MYC), are crucial to the formation and development of leaf trichomes in plants. bHLH constitutes one of the largest families of transcription factors in eukaryotes, of which MYC is a subfamily member. However, studies on bHLH-MYC transcription factors in mustard have yet to be reported. In this study, a total of 45 bHLH-MYC transcription factors were identified within the Brassica juncea genome, and a comprehensive series of bioinformatic analyses were conducted on their structures and properties: an examination of protein physicochemical properties, an exploration of conserved structural domains, an assessment of chromosomal positional distributions, an analysis of the conserved motifs, an evaluation of the gene structures, microsynteny analyses, three-dimensional structure prediction, and an analysis of sequence signatures. Finally, transcriptome analyses and a subcellular localization examination were performed. The results revealed that these transcription factors were unevenly distributed across 18 chromosomes, showing relatively consistent conserved motifs and gene structures and high homology. The final results of the transcriptome analysis and gene annotation showed a high degree of variability in the expression of bHLH-MYC transcription factors. Five genes that may be associated with trichome development (BjuVA09G22490, BjuVA09G13750, BjuVB04G14560, BjuVA05G24810, and BjuVA06G44820) were identified. The subcellular localization results indicated that the transcription and translation products of these five genes were expressed in the same organelle: the nucleus. This finding provides a basis for elucidating the roles of bHLH-MYC family members in plant growth and development, and the molecular mechanisms underlying trichome development in mustard leaves.

Bolting and flowering time are critical agronomic traits affecting the commercial value and breeding efficiency of Chinese cabbage (Brassica rapa L. ssp. pekinensis). Although vernalization is a key environmental signal promoting flowering, its regulatory mechanisms remain poorly understood in this crop. Here, we identify the flowering repressor gene BrMAF5 and its antisense long non-coding RNA BrMAF5L as negative regulators of vernalization-induced flowering. During vernalization, transcript levels of both genes showed a decreasing trend as the vernalization period extended. Functional assays in Arabidopsis thaliana demonstrated that ectopic expression of BrMAF5 or BrMAF5L significantly delayed flowering, accompanied by increased expression of floral repressors (AtFLC, AtTEM1) and reduced expression of floral activators (AtFT, AtSOC1). Moreover, protein interaction analysis revealed that BrMAF5 associates with BrACP4 and BrRAB1A, linking it to fatty acid metabolism and membrane trafficking pathways. Collectively, our results reveal a novel regulatory module in vernalization-mediated flowering. These findings pave the way for developing bolting-resistant Brassicaceae crops by identifying promising molecular targets.