In order to make the identification and utilization of nutritional quality components in Chinese cabbage more predictive, to obtain ideal raw materials, and to help screen functional Chinese cabbage varieties that have high nutritional value, we conducted a comprehensive evaluation of the nutritional quality of different Chinese cabbage varieties. In this study, 17 nutritional quality indexes of 35 Chinese cabbage varieties, including crude fiber (CF), crude protein (CP), vitamin C (VC), glucose (Glc), fructose (Fru), malic acid (MA), citric acid (CA), oxalic acid (OA), total amino acid (TAA) and CA, K, Mg, P, Cu, Fe, Mn and Zn, were analyzed using diversity analysis, correlation analysis, principal component analysis, membership function analysis and cluster analysis. The results showed that there were different degrees of variation in the 17 nutritional quality indexes, and the coefficients of variation ranged from 11.45% to 91.47%. The correlation analysis found that there were significant or extremely significant correlations between different nutrient elements of Chinese cabbage, which indicated that principal component analysis could be carried out, and the comprehensive score (D value) of different materials could be obtained using principal component analysis and the membership function method. The nutritional quality of Chinese cabbage was classified into five grades by cluster analysis. Finally, a mathematical model for evaluating the nutritional quality of Chinese cabbage was established using the D value and multiple stepwise regression methods, and 10 key indexes were selected from the 17 indexes, which could be used for the rapid identification of the nutritional quality of Chinese cabbage. This study provided a theoretical basis for the nutritional quality evaluation and variety breeding of Chinese cabbage.

The swollen root is an important agronomic trait and is a determinant of yield for turnips, which are cultivated as both vegetables and fodder. However, the genetic mechanism of swollen root formation is poorly understood. In this study, we analyzed the F₂ and BC₁P₂ populations derived from a cross between “10601” (European turnip with swollen root, Brassica rapa ssp. rapifera, AA, 2n = 2× = 20) and “10603” (Chinese cabbage with normal root, Brassica rapa ssp. pekinensis, AA, 2n = 2× = 20), and suggested that the swollen root is a quantitative trait. Two major quantitative trait loci (QTLs), FR1.1 (Fleshy root 1.1) and FR7.1 (Fleshy root 7.1), were identified by QTL-seq analysis and further confirmed by QTL mapping in F₂ and BC₁P₂ populations. The QTL FR1.1 with a likelihood of odd (LOD) of 7.01 explained 17.2% of the total phenotypic variations for root diameter and the QTL FR7.1 explained 23.0% (LOD = 9.38) and 31.0% (LOD = 13.27) of the total phenotypic variations in root diameter and root weight, respectively. After a recombinant screening, the major QTL FR7.1 was further narrowed down to a 220 kb region containing 47 putative genes. A candidate gene, Bra003652, which is a homolog of AT1G78240 that plays an essential role in cell adhesion and disorganized tumor-like formation in Arabidopsis thaliana, was identified in this region. In addition, expression and parental allele analysis supported that Bra003652 was a possible candidate gene of QTL FR7.1 for swollen root formation in turnip. Our research may provide new insight into the molecular mechanism of swollen root formation in root crops.

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.

Brassica rapa is an important Chinese vegetable crop that is beneficial to human health. The primary factor affecting B. rapa yield is low temperature, which promotes bolting and flowering, thereby lowering its commercial value. However, quickened bolting and flowering can be used for rapid breeding. Therefore, studying the underlying molecular mechanism of vernalization in B. rapa is crucial for solving production-related problems. Here, the transcriptome of two B. rapa accessions were comprehensively analyzed during different vernalization periods. During vernalization, a total of 974,584,022 clean reads and 291.28 Gb of clean data were obtained. Compared to the reference genome of B. rapa, 44,799 known genes and 2280 new genes were identified. A self-organizing feature map analysis of 21,035 differentially expressed genes was screened in two B. rapa accessions, ‘Jin Wawa’ and ‘Xiao Baojian’. The analysis indicated that transcripts related to the plant hormone signal transduction, starch and sucrose metabolism, photoperiod and circadian clock, and vernalization pathways changed notably at different vernalization periods. Moreover, different expression patterns of TPS, UGP, CDF, VIN1, and seven hormone pathway genes were observed during vernalization between the two accessions. The transcriptome results of this study provide a new perspective on the changes that occur during B. rapa vernalization, as well as serve as an excellent reference for B. rapa breeding.

Allopolyploidy is an evolutionary and mechanistically intriguing process involving the reconciliation of two or more sets of diverged genomes and regulatory interactions, resulting in new phenotypes. In this study, we explored the gene expression patterns of eight F2 synthetic Brassica napus using RNA sequencing. We found that B. napus allopolyploid formation was accompanied by extensive changes in gene expression. A comparison between F2 and the parent shows a certain proportion of differentially expressed genes (DEG) and activation\silent gene, and the two genomes (female parent (AA)\male parent (CC) genomes) showed significant differences in response to whole-genome duplication (WGD); non-additively expressed genes represented a small portion, while Gene Ontology (GO) enrichment analysis showed that it played an important role in responding to WGD. Besides, genome-wide expression level dominance (ELD) was biased toward the AA genome, and the parental expression pattern of most genes showed a high degree of conservation. Moreover, gene expression showed differences among eight individuals and was consistent with the results of a cluster analysis of traits. Furthermore, the differential expression of waxy synthetic pathways and flowering pathway genes could explain the performance of traits. Collectively, gene expression of the newly formed allopolyploid changed dramatically, and this was different among the selfing offspring, which could be a prominent cause of the trait separation. Our data provide novel insights into the relationship between the expression of differentially expressed genes and trait segregation and provide clues into the evolution of allopolyploids.

Turnip mosaic virus (TuMV), which is distributed almost all over the world and has a wide range of hosts, mainly brassica crops, was first described in Brassica rapa in the USA. Plant volatile compounds play an important role in the host searching behavior of natural enemies of herbivorous insects. In this study, TuMV-inoculated resistant and susceptible B. rapa lines were tested using volatile metabolome and transcriptome analyses. In volatile metabolome analysis, the volatile organic compounds (VOCs) were different after inoculation with TuMV in resistant B80124 and susceptible B80461, and the degree of downregulation of differentially expressed metabolites was more obvious than the degree of upregulation. Through transcriptome analysis, 70% of differentially expressed genes were in biological process, especially focusing on defense response, flavonoid biosynthetic process, and toxin metabolic process, which indicates that TuMV stress maybe accelerate the increase of VOCs. Integrating the metabolome and transcriptome analyses, after inoculating with TuMV, auxin regulation was upregulated, and ARF, IAA and GH3 were also upregulated, which accelerated cell enlargement and plant growth in tryptophan metabolism. The different genes in zeatin biosynthesis pathways were downregulated, which reduced cell division and shoot initiation. However, the metabolite pathways showed upregulation in brassinosteroid biosynthesis and α-linolenic acid metabolism, which could cause cell enlargement and a stress response. This study determined the difference in volatiles between normal plants and infected plants and may lay a foundation for anti-TuMV research in B. rapa.