Deciphering the transcriptional regulatory networks that control size, color, and oil content in Brassica rapa seeds

Abstract Background Brassica rapa is an important oilseed and vegetable crop species and is the A subgenome donor of two important oilseed Brassica crops, Brassica napus and Brassica juncea. Although seed size (SZ), seed color (SC), and oil content (OC) substantially affect seed yield and quality, the mechanisms regulating these traits in Brassica crops remain unclear. Results We collected seeds from a pair of B. rapa accessions with significantly different SZ, SC, and OC at seven seed developmental stages (every 7 days from 7 to 49 days after pollination), and identified 28,954 differentially expressed genes (DEGs) from seven pairwise comparisons between accessions at each developmental stage. K-means clustering identified a group of cell cycle-related genes closely connected to variation in SZ of B. rapa. A weighted correlation analysis using the WGCNA package in R revealed two important co-expression modules comprising genes whose expression was positively correlated with SZ increase and negatively correlated with seed yellowness, respectively. Upregulated expression of cell cycle-related genes in one module was important for the G2/M cell cycle transition, and the transcription factor Bra.A05TSO1 seemed to positively stimulate the expression of two CYCB1;2 genes to promote seed development. In the second module, a conserved complex regulated by the transcription factor TT8 appear to determine SC through downregulation of TT8 and its target genes TT3, TT18, and ANR. In the third module, WRI1 and FUS3 were conserved to increase the seed OC, and Bra.A03GRF5 was revealed as a key transcription factor on lipid biosynthesis. Further, upregulation of genes involved in triacylglycerol biosynthesis and storage in the seed oil body may increase OC. We further validated the accuracy of the transcriptome data by quantitative real-time PCR of 15 DEGs. Finally, we used our results to construct detailed models to clarify the regulatory mechanisms underlying variations in SZ, SC, and OC in B. rapa. Conclusions This study provides insight into the regulatory mechanisms underlying the variations of SZ, SC, and OC in plants based on transcriptome comparison. The findings hold great promise for improving seed yield, quality and OC through genetic engineering of critical genes in future molecular breeding.