Comparative transcriptomic and metabolomic analyses of carotenoid biosynthesis reveal the basis of white petal color in Brassica napus.

MAIN CONCLUSION: The molecular mechanism underlying white petal color in Brassica napus was revealed by transcriptomic and metabolomic analyses. Rapeseed (Brassica napus L.) is one of the most important oilseed crops worldwide, but the mechanisms underlying flower color in this crop are known less. Here, we performed metabolomic and transcriptomic analyses of the yellow-flowered rapeseed cultivar 'Zhongshuang 11' (ZS11) and the white-flowered inbred line 'White Petal' (WP). The total carotenoid contents were 1.778-fold and 1.969-fold higher in ZS11 vs. WP petals at stages S2 and S4, respectively. Our findings suggest that white petal color in WP flowers is primarily due to decreased lutein and zeaxanthin contents. Transcriptome analysis revealed 10,116 differentially expressed genes with a fourfold or greater change in expression (P-value less than 0.001) in WP vs. ZS11 petals, including 1,209 genes that were differentially expressed at four different stages and 20 genes in the carotenoid metabolism pathway. BnNCED4b, encoding a protein involved in carotenoid degradation, was expressed at abnormally high levels in WP petals, suggesting it might play a key role in white petal formation. The results of qRT-PCR were consistent with the transcriptome data. The results of this study provide important insights into the molecular mechanisms of the carotenoid metabolic pathway in rapeseed petals, and the candidate genes identified in this study provide a resource for the creation of new B. napus germplasms with different petal colors.

Genome-wide identification and expression profiling of the carotenoid cleavage dioxygenase (CCD) gene family in Brassica napus L.

Carotenoid cleavage dioxygenase (CCD), a key enzyme in carotenoid metabolism, cleaves carotenoids to form apo-carotenoids, which play a major role in plant growth and stress responses. CCD genes had not previously been systematically characterized in Brassica napus (rapeseed), an important oil crop worldwide. In this study, we identified 30 BnCCD genes and classified them into nine subgroups based on a phylogenetic analysis. We identified the chromosomal locations, gene structures, and cis-promoter elements of each of these genes and performed a selection pressure analysis to identify residues under selection. Furthermore, we determined the subcellular localization, physicochemical properties, and conserved protein motifs of the encoded proteins. All the CCD proteins contained a retinal pigment epithelial membrane protein (RPE65) domain. qRT-PCR analysis of expression of 20 representative BnCCD genes in 16 tissues of the B. napus cultivar Zhong Shuang 11 ('ZS11') revealed that members of the BnCCD gene family possess a broad range of expression patterns. This work lays the foundation for functional studies of the BnCCD gene family.

Identification and Characterization of Major Constituents in Different-Colored Rapeseed Petals by UPLC-HESI-MS/MS.

Oilseed rape (Brassica napus L.) is the second highest yielding oil crop worldwide. In addition to being used as an edible oil and a feed for livestock, rapeseed has high ornamental value. In this study, we identified and characterized the main floral major constituents, including phenolic acids and flavonoids components, in rapeseed accessions with different-colored petals. A total of 144 constituents were identified using ultrahigh-performance liquid chromatography-HESI-mass spectrometry (UPLC-HESI-MS/MS), 57 of which were confirmed and quantified using known standards and mainly contained phenolic acids, flavonoids, and glucosinolates compounds. Most of the epicatechin, quercetin, and isorhamnetin derivates were found in red and pink petals of B. napus, while kaempferol derivates were in yellow and pale white petals. Moreover, petal-specific compounds, including a putative hydroxycinnamic acid derivative, sinapoyl malate, 1-O-sinapoyl-ß-d-glucose, feruloyl glucose, naringenin-7-O-glucoside, cyanidin-3-glucoside, cyanidin-3,5-di-O-glucoside, petunidin-3-O-ß-glucopyranoside, isorhamnetin-3-O-glucoside, kaempferol-3-O-glucoside-7-O-glucoside, quercetin-3,4'-O-di-ß-glucopyranoside, quercetin-3-O-glucoside, and delphinidin-3-O-glucoside, might contribute to a variety of petal colors in B. napus. In addition, bound phenolics were tentatively identified and contained three abundant compounds (p-coumaric acid, ferulic acid, and 8-O-4'-diferulic acid). These results provide insight into the molecular mechanisms underlying petal color and suggest strategies for breeding rapeseed with a specific petal color in the future.

Association Mapping Analysis of Fatty Acid Content in Different Ecotypic Rapeseed Using mrMLM.

Brassica napus L. is a widely cultivated oil crop and provides important resources of edible vegetable oil, and its quality is determined by fatty acid composition and content. To explain the genetic basis and identify more minor loci for fatty acid content, the multi-locus random-SNP-effect mixed linear model (mrMLM) was used to identify genomic regions associated with fatty acid content in a genetically diverse population of 435 rapeseed accessions, including 77 winter-type, 55 spring-type, and 303 semi-winter-type accessions grown in different environments. A total of 149 quantitative trait nucleotides (QTNs) were found to be associated with fatty acid content and composition, including 34 QTNs that overlapped with the previously reported loci, and 115 novel QTNs. Of these, 35 novel QTNs, located on chromosome A01, A02, A03, A05, A06, A09, A10, and C02, respectively, were repeatedly detected across different environments. Subsequently, we annotated 95 putative candidate genes by BlastP analysis using sequences from Arabidopsis thaliana homologs of the identified regions. The candidate genes included 34 environmentally-insensitive genes (e.g., CER4, DGK2, KCS17, KCS18, MYB4, and TT16) and 61 environment-sensitive genes (e.g., FAB1, FAD6, FAD7, KCR1, KCS9, KCS12, and TT1) as well as genes invloved in the fatty acid biosynthesis. Among these, BnaA08g08280D and BnaC03g60080D differed in genomic sequence between the high- and low-oleic acid lines, and might thus be the novel alleles regulating oleic acid content. Furthermore, RT-qPCR analysis of these genes showed differential expression levels during seed development. Our results highlight the practical and scientific value of mrMLM or QTN detection and the accuracy of linking specific QTNs to fatty acid content, and suggest a useful strategy to improve the fatty acid content of B. napus seeds by molecular marker-assisted breeding.

Genome-Wide Identification and Structural Analysis of bZIP Transcription Factor Genes in Brassica napus.

The basic region/leucine zipper motif (bZIP) transcription factor family is one of the largest families of transcriptional regulators in plants. bZIP genes have been systematically characterized in some plants, but not in rapeseed (Brassica napus). In this study, we identified 247 BnbZIP genes in the rapeseed genome, which we classified into 10 subfamilies based on phylogenetic analysis of their deduced protein sequences. The BnbZIP genes were grouped into functional clades with Arabidopsis genes with similar putative functions, indicating functional conservation. Genome mapping analysis revealed that the BnbZIPs are distributed unevenly across all 19 chromosomes, and that some of these genes arose through whole-genome duplication and dispersed duplication events. All expression profiles of 247 bZIP genes were extracted from RNA-sequencing data obtained from 17 different B. napus ZS11 tissues with 42 various developmental stages. These genes exhibited different expression patterns in various tissues, revealing that these genes are differentially regulated. Our results provide a valuable foundation for functional dissection of the different BnbZIP homologs in B. napus and its parental lines and for molecular breeding studies of bZIP genes in B. napus.

Genome-wide association mapping and Identification of candidate genes for fatty acid composition in Brassica napus L. using SNP markers.

BACKGROUND: B. napus (oilseed) is an important source of edible vegetable oil, and its nutritional and economic value is determined by its fatty acid composition and content. RESULTS: Using the Brassica 60 K SNP array, we performed a genome-wide association study of fatty acid composition in a population of 520 genetically diverse oilseed accessions. Using the PCA + K model in TASSEL 5.2.1, we identified 62 genomic regions that were significantly associated with the composition of seven fatty acids, and five consensus regions that mapped to the A2, A8, A9, C1, and C3 chromosomes, respectively, of the Brassica napus Darmor-bzh genome. We then identified 24 orthologs of the functional candidate genes involved in fatty acid biosynthesis, excluding BnaA.FAE1 and BnaC.FAE1 on the A8 and C3 homologous genome blocks, which are known to have critical roles in the fatty acid biosynthesis pathway, and potential orthologs of these genes (e.g., LACS9, KCR1, FAB1, LPAT4, KCS17, CER4, TT16, and ACBP5). CONCLUSIONS: Our results demonstrate the power of association mapping in identifying genes of interest in B. napus and provide insight into the genetic basis of fatty acid biosynthesis in B. napus. Furthermore, our findings may facilitate marker-based breeding efforts aimed at improving fatty acid composition and quality in B. napus.

Identification of Candidate Genes for Seed Glucosinolate Content Using Association Mapping in Brassica napus L.

Rapeseed contains glucosinolates, a toxic group of sulfur-containing glucosides, which play critical roles in defense against herbivores and microbes. However, the presence of glucosinolates in rapeseed reduces the value of the meal as feed for livestock. We performed association mapping of seed glucosinolate (GS) content using the 60K Brassica Infinium single nucleotide polymorphism (SNP) array in 520 oilseed rape accessions. A total of 11 peak SNPs significantly associated with GS content were detected in growing seasons of 2013 and 2014 and were located on B. napus chromosomes A08, A09, C03, and C09, respectively. Two associated regions of GS content covered by these markers were further verified, and three B. napus homologous genes involved in the biosynthesis and accumulation of GS were identified. These genes were multigene family members and were distributed on different chromosomes. Moreover, two genes (BnGRT2 and BnMYB28) associated with GS content were validated by the qRT-PCR analysis of their expression profiles. The further identification and functionalization of these genes will provide useful insight into the mechanism underlying GS biosynthesis and allocation in B. napus, and the associated SNPs markers could be helpful for molecular maker-assisted breeding for low seed GS in B. napus.