Comparative genomic analyses reveal the genetic basis of the yellow-seed trait in Brassica napus.

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47aZY821, but not the yellow-seed allele BnA09MYB47aGH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus.

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

Metabolite Characteristics Analysis of Siliques and Effects of Lights on the Accumulation of Glucosinolates in Siliques of Rapeseed.

Glucosinolates (GSLs) are naturally occurring secondary metabolites found in the Brassicaceae family, which mainly synthesize in the siliques with a wide range of functions. In this study, we investigated the effects of lights on metabolites in siliques of rapeseed through ultra high-performance liquid chromatography (UPLC)-heated electrospray ionization (HESI)-tandem mass spectrometry (MS/MS). A total of 249 metabolites, including 29 phenolic acids, 38 flavonoids, 22 GSLs, 93 uncalculated and 67 unknown compounds, were identified in siliques of rapeseed. Meanwhile, 62 metabolites showed significant differences after shading treatment, which were mainly GSLs and unknown compounds. Interestingly, the amounts of 10 GSLs had high accumulation levels in siliques, while the expression levels of their corresponding biosynthetic genes (AOP, GSL-OH, IGMT, and ST5a) were obviously reduced after shading treatment. Further evidence showed that the amounts of GSLs were significantly reduced in seeds, in accordance with the expression profiles of transporter genes (BnaGTRs). Our findings indicated that lights could affect the accumulation and transportation of GSLs from siliques to seeds in rapeseed. Therefore, this study facilitates a better understanding of metabolic characteristics of siliques and provides insight into the importance of light for GSLs accumulation and transportation in siliques and seeds of rapeseed.

Genome-Wide Analysis of the YABBY Transcription Factor Family in Rapeseed (Brassica napus L.).

The YABBY family of plant-specific transcription factors play important regulatory roles during the development of leaves and floral organs, but their functions in Brassica species are incompletely understood. Here, we identified 79 YABBY genes from Arabidopsis thaliana and five Brassica species (B. rapa, B. nigra, B. oleracea, B. juncea, and B. napus). A phylogenetic analysis of YABBY proteins separated them into five clusters (YAB1-YAB5) with representatives from all five Brassica species, suggesting a high degree of conservation and similar functions within each subfamily. We determined the gene structure, chromosomal location, and expression patterns of the 21 BnaYAB genes identified, revealing extensive duplication events and gene loss following polyploidization. Changes in exon-intron structure during evolution may have driven differentiation in expression patterns and functions, combined with purifying selection, as evidenced by Ka/Ks values below 1. Based on transcriptome sequencing data, we selected nine genes with high expression at the flowering stage. qRT-PCR analysis further indicated that most BnaYAB family members are tissue-specific and exhibit different expression patterns in various tissues and organs of B. napus. This preliminary study of the characteristics of the YABBY gene family in the Brassica napus genome provides theoretical support and reference for the later functional identification of the family genes.

Genome-wide identification AINTEGUMENTA-like (AIL) genes in Brassica species and expression patterns during reproductive development in Brassica napus L.

The AINTEGUMENTA-like (AIL) proteins, which belong to the AP2 family, play important roles in regulating the growth and development of plant organs. The AIL family has not yet been comprehensively studied in rapeseed (Brassica napus), an allotetraploid and model organism for the study of polyploid evolution. In the present study, 99 AIL family genes were identified and characterized from B. rapa, B. oleracea, B. napus, B. juncea, and B. nigra using a comprehensive genome-wide study, including analyses of phylogeny, gene structure, chromosomal localization, and expression pattern. Using a phylogenetic analysis, the AIL genes were divided into eight groups, which were closely related to the eight AtAIL genes, and which shared highly conserved structural features within the same subfamily. The non-synonymous/synonymous substitution ratios of the paralogs and orthologs were less than 1, suggesting that the AIL genes mainly experienced purifying selection during evolution. In addition, the RNA sequencing data and qRT-PCR analysis revealed that the B. napus AIL genes exhibited organ- and developmental stage-specific expression patterns. Certain genes were highly expressed in the developing seeds (BnaAIL1, BnaAIL2, BnaAIL5, and BnaAIL6), the roots (BnaANT, BnaAIL5, and BnaAIL6), and the stem (BnaAIL7B). Our results provide valuable information for further functional analysis of the AIL family in B. napus and related Brassica species.

Characterization of cold stress responses in different rapeseed ecotypes based on metabolomics and transcriptomics analyses.

The winter oilseed ecotype is more tolerant to low temperature than the spring ecotype. Transcriptome and metabolome analyses of leaf samples of five spring Brassica napus L. (B. napus) ecotype lines and five winter B. napus ecotype lines treated at 4 °C and 28 °C were performed. A total of 25,460 differentially expressed genes (DEGs) of the spring oilseed ecotype and 28,512 DEGs of the winter oilseed ecotype were identified after cold stress; there were 41 differentially expressed metabolites (DEMs) in the spring and 47 in the winter oilseed ecotypes. Moreover, more than 46.2% DEGs were commonly detected in both ecotypes, and the extent of the changes were much more pronounced in the winter than spring ecotype. By contrast, only six DEMs were detected in both the spring and winter oilseed ecotypes. Eighty-one DEMs mainly belonged to primary metabolites, including amino acids, organic acids and sugars. The large number of specific genes and metabolites emphasizes the complex regulatory mechanisms involved in the cold stress response in oilseed rape. Furthermore, these data suggest that lipid, ABA, secondary metabolism, signal transduction and transcription factors may play distinct roles in the spring and winter ecotypes in response to cold stress. Differences in gene expression and metabolite levels after cold stress treatment may have contributed to the cold tolerance of the different oilseed ecotypes.

Comparative Analysis of the Metabolic Profiles of Yellow- versus Black-Seeded Rapeseed Using UPLC-HESI-MS/MS and Transcriptome Analysis.

The high levels of secondary metabolites in rapeseed play important roles in determining the oil quality and feeding value. Here, we characterized the metabolic profiles in seeds of various yellow- and black-seeded rapeseed accessions. Two hundred and forty-eight features were characterized, including 31 phenolic acids, 54 flavonoids, 24 glucosinolates, 65 lipid compounds, and 74 other polar compounds. The most abundant phenolic acids and various flavonoids (epicatechin, isorhamnetin, kaempferol, quercetin, and their derivatives) were widely detected and showed significant differences in distribution between the yellow- and black-seeded rapeseed. Furthermore, the related genes (e.g., BnTT3, BnTT18, BnTT10, BnTT12, and BnBAN) involved in the proanthocyanidin pathway had lower expression levels in yellow-seeded rapeseed, strongly suggesting that the seed coat color could be mainly determined by the levels of epicatechin and their derivatives. These results improve our understanding of the primary constituents of rapeseed and lay the foundation for breeding novel varieties with a high nutritional value.

Genome-wide exploration and characterization of miR172/euAP2 genes in Brassica napus L. for likely role in flower organ development.

BACKGROUND: APETALA2-like genes encode plant-specific transcription factors, some of which possess one microRNA172 (miR172) binding site. The miR172 and its target euAP2 genes are involved in the process of phase transformation and flower organ development in many plants. However, the roles of miR172 and its target AP2 genes remain largely unknown in Brassica napus (B. napus). RESULTS: In this study, 19 euAP2 and four miR172 genes were identified in the B. napus genome. A sequence analysis suggested that 17 euAP2 genes were targeted by Bna-miR172 in the 3' coding region. EuAP2s were classified into five major groups in B.napus. This classification was consistent with the exon-intron structure and motif organization. An analysis of the nonsynonymous and synonymous substitution rates revealed that the euAP2 genes had gone through purifying selection. Whole genome duplication (WGD) or segmental duplication events played a major role in the expansion of the euAP2 gene family. A cis-regulatory element (CRE) analysis suggested that the euAP2s were involved in the response to light, hormones, stress, and developmental processes including circadian control, endosperm and meristem expression. Expression analysis of the miR172-targeted euAP2s in nine different tissues showed diverse spatiotemporal expression patterns. Most euAP2 genes were highly expressed in the floral organs, suggesting their specific functions in flower development. BnaAP2-1, BnaAP2-5 and BnaTOE1-2 had higher expression levels in late-flowering material than early-flowering material based on RNA-seq and qRT-PCR, indicating that they may act as floral suppressors. CONCLUSIONS: Overall, analyses of the evolution, structure, tissue specificity and expression of the euAP2 genes were peformed in B.napus. Based on the RNA-seq and experimental data, euAP2 may be involved in flower development. Three euAP2 genes (BnaAP2-1, BnaAP2-5 and BnaTOE1-2) might be regarded as floral suppressors. The results of this study provide insights for further functional characterization of the miR172 /euAP2 module in B.napus.

Genome-wide identification and characterization of Gretchen Hagen3 (GH3) family genes in Brassica napus.

The hormone auxin is involved in many biological processes throughout a plant's lifecycle. However, genes in the GH3 (Gretchen Hagen3) family, one of the three major auxin-responsive gene families, have not yet been identified in oilseed rape (Brassica napus). In this study, we identified 63 BnaGH3 genes in oilseed rape using homology searches. We analyzed the chromosome locations, gene structures, and phylogenetic relationships of the BnaGH3 genes, as well as the cis-elements in their promoters. Most BnaGH3 genes are located on chromosomes A03, A09, C02, C03, and C09, each with 4-7 members. In addition, we analyzed the expression patterns of BnaGH3 genes in seven tissues by transcriptome sequencing and quantitative RT-PCR analysis of plants under exogenous IAA treatment. The BnaGH3 genes showed different expression patterns in various tissues. BnaA.GH3.2-1 and BnaC.GH3.2-1 were expressed in the seed and seed coat during development and in response to IAA treatment. These results shed light on the possible roles of the GH3 gene family in oilseed rape.

Genome-Wide Identification of the LAC Gene Family and Its Expression Analysis Under Stress in Brassica napus.

Lignin is an important biological polymer in plants that is necessary for plant secondary cell wall ontogenesis. The laccase (LAC) gene family catalyzes lignification and has been suggested to play a vital role in the plant kingdom. In this study, we identified 45 LAC genes from the Brassica napus genome (BnLACs), 25 LAC genes from the Brassica rapa genome (BrLACs) and 8 LAC genes from the Brassica oleracea genome (BoLACs). These LAC genes could be divided into five groups in a cladogram and members in same group had similar structures and conserved motifs. All BnLACs contained hormone- and stress- related elements determined by cis-element analysis. The expression of BnLACs was relatively higher in the root, seed coat and stem than in other tissues. Furthermore, BnLAC4 and its predicted downstream genes showed earlier expression in the silique pericarps of short silique lines than long silique lines. Three miRNAs (miR397a, miR397b and miR6034) target 11 BnLACs were also predicted. The expression changes of BnLACs under series of stresses were further investigated by RNA sequencing (RNA-seq) and quantitative real-time polymerase chain reaction (qRT-PCR). The study will give a deeper understanding of the LAC gene family evolution and functions in B. napus.

Joint QTL mapping and transcriptome sequencing analysis reveal candidate flowering time genes in Brassica napus L.

BACKGROUND: Optimum flowering time is a key agronomic trait in Brassica napus. To investigate the genetic architecture and genetic regulation of flowering time in this important crop, we conducted quantitative trait loci (QTL) analysis of flowering time in a recombinant inbred line (RIL) population, including lines with extreme differences in flowering time, in six environments, along with RNA-Seq analysis. RESULTS: We detected 27 QTLs distributed on eight chromosomes among six environments, including one major QTL on chromosome C02 that explained 11-25% of the phenotypic variation and was stably detected in all six environments. RNA-Seq analysis revealed 105 flowering time-related differentially expressed genes (DEGs) that play roles in the circadian clock/photoperiod, autonomous pathway, and hormone and vernalization pathways. We focused on DEGs related to the regulation of flowering time, especially DEGs in QTL regions. CONCLUSIONS: We identified 45 flowering time-related genes in these QTL regions, eight of which are DEGs, including key flowering time genes PSEUDO RESPONSE REGULATOR 7 (PRR7) and FY (located in a major QTL region on C02). These findings provide insights into the genetic architecture of flowering time in B. napus.

Genome-Wide Characterization and Analysis of Metallothionein Family Genes That Function in Metal Stress Tolerance in Brassica napus L.

Brassica plants exhibit both high biomass productivity and high rates of heavy metal absorption. Metallothionein (MT) proteins are low molecular weight, cysteine-rich, metal-binding proteins that play crucial roles in protecting plants from heavy metal toxicity. However, to date, MT proteins have not been systematically characterized in Brassica. In this study, we identified 60 MTs from Arabidopsis thaliana and five Brassica species. All the MT family genes from Brassica are closely related to Arabidopsis MTs, encoding putative proteins that share similar functions within the same clades. Genome mapping analysis revealed high levels of synteny throughout the genome due to whole genome duplication and segmental duplication events. We analyzed the expression levels of 16 Brassica napus MTs (BnaMTs) by RNA-sequencing and real-time RT-PCR (RT-qPCR) analysis in plants under As3+ stress. These genes exhibited different expression patterns in various tissues. Our results suggest that BnaMT3C plays a key role in the response to As3+ stress in B. napus. This study provides insight into the phylogeny, origin, and evolution of MT family members in Brassica, laying the foundation for further studies of the roles of MT proteins in these important crops.

Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (Lib.) infection.

Sclerotinia stem rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is a devastating disease resulting in yield losses and decreases in seed quality in oilseed rape (Brassica napus) worldwide. However, the molecular mechanisms underlying the response of oilseed rape to S. sclerotiorum infection at the transcriptional and post-transcriptional levels are poorly understood. Here, we used an integrated omics approach (transcriptome, sRNAome, and degradome sequencing) on the Illumina platform to compare the RNA expression and post-transcriptional profiles of oilseed rape plants inoculated or not with S. sclerotiorum. In total, 7,065 differentially expressed genes (DEGs) compared with the mock-inoculated control at 48 hours post inoculation were identified. These DEGs were associated with protein kinases, signal transduction, transcription factors, hormones, pathogenesis-related proteins, secondary metabolism, and transport. In the sRNA-Seq analysis, 77 known and 176 novel miRNAs were identified; however, only 10 known and 41 novel miRNAs were differentially expressed between the samples inoculated or not with S. sclerotiorum. Degradome sequencing predicted 80 cleavage sites with 64 miRNAs. Integrated mRNA, sRNA and degradome sequencing analysis reveal oilseed rape complex responses to S. sclerotiorum infection. This study provides a global view of miRNA and mRNA expression profiles in oilseed rape following S. sclerotiorum infection.

Genome-Wide Identification of MicroRNAs in Response to Cadmium Stress in Oilseed Rape (Brassica napus L.) Using High-Throughput Sequencing.

MicroRNAs (miRNAs) have important roles in regulating stress-response genes in plants. However, identification of miRNAs and the corresponding target genes that are induced in response to cadmium (Cd) stress in Brassica napus remains limited. In the current study, we sequenced three small-RNA libraries from B. napus after 0 days, 1 days, and 3 days of Cd treatment. In total, 44 known miRNAs (belonging to 27 families) and 103 novel miRNAs were identified. A comprehensive analysis of miRNA expression profiles found 39 differentially expressed miRNAs between control and Cd-treated plants; 13 differentially expressed miRNAs were confirmed by qRT-PCR. Characterization of the corresponding target genes indicated functions in processes including transcription factor regulation, biotic stress response, ion homeostasis, and secondary metabolism. Furthermore, we propose a hypothetical model of the Cd-response mechanism in B. napus. Combined with qRT-PCR confirmation, our data suggested that miRNAs were involved in the regulations of TFs, biotic stress defense, ion homeostasis and secondary metabolism synthesis to respond Cd stress in B. napus.

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.

qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms.

Real-time quantitative polymerase chain reaction (qPCR) is one of the most important methods for analyzing the expression patterns of target genes. However, successful qPCR experiments rely heavily on the use of high-quality primers. Various qPCR primer databases have been developed to address this issue, but these databases target only a few important organisms. Here, we developed the qPrimerDB database, founded on an automatic gene-specific qPCR primer design and thermodynamics-based validation workflow. The qPrimerDB database is the most comprehensive qPCR primer database available to date, with a web front-end providing gene-specific and pre-computed primer pairs across 147 important organisms, including human, mouse, zebrafish, yeast, thale cress, rice and maize. In this database, we provide 3331426 of the best primer pairs for each gene, based on primer pair coverage, as well as 47760359 alternative gene-specific primer pairs, which can be conveniently batch downloaded. The specificity and efficiency was validated for qPCR primer pairs for 66 randomly selected genes, in six different organisms, through qPCR assays and gel electrophoresis. The qPrimerDB database represents a valuable, timesaving resource for gene expression analysis. This resource, which will be routinely updated, is publically accessible at http://biodb.swu.edu.cn/qprimerdb.

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.

Genetic and transcriptomic analyses of lignin- and lodging-related traits in Brassica napus.

KEY MESSAGE: Candidate genes associated with lignin and lodging traits were identified by combining phenotypic, genotypic, and gene expression data in B. napus. Brassica napus is one of the world's most important oilseed crops, but its yield can be dramatically reduced by lodging, bending, and falling of its vertical stems. Lignin has been shown to contribute to stem mechanical strength. In this study, we found that the syringyl/guaiacyl (S/G) monolignol ratio exhibits a significant negative correlation with disease and lodging resistance. A total of 92 and 50 SNP and SSR loci, respectively, were found to be significantly associated with five traits, breaking force, breaking strength, lodging coefficient, acid detergent lignin content, and the S/G monolignol ratio using GWAS. To identify novel genes involved in lignin biosynthesis, transcriptome sequencing of high- (H) and low (L)-ADL content accessions was performed. The up-regulated genes were mainly involved in glycoside catabolic processes (especially glucosinolate catabolism) and cell wall biogenesis, while down-regulated genes were involved in glucosinolate biosynthesis, indicating that crosstalk exists between glucosinolate metabolic processes and lignin biosynthesis. Integrating this differential expression with the GWAS analysis, we identified four candidate genes regulating lignin, including glycosyl hydrolase (BnaA01g00480D), CYT1 (BnaA04g22820D), and two encoding transcription factors, SHINE1 (ERF family) and DAR6 (LIM family). This study provides insight into the genetic control of lodging and lignin in B. napus.

Genome-Wide Association Mapping of Seed Coat Color in Brassica napus.

Seed coat color is an extremely important breeding characteristic of Brassica napus. To elucidate the factors affecting the genetic architecture of seed coat color, a genome-wide association study (GWAS) of seed coat color was conducted with a diversity panel comprising 520 B. napus cultivars and inbred lines. In total, 22 single-nucleotide polymorphisms (SNPs) distributed on 7 chromosomes were found to be associated with seed coat color. The most significant SNPs were found in 2014 near Bn-scaff_15763_1-p233999, only 43.42 kb away from BnaC06g17050D, which is orthologous to Arabidopsis thaliana TRANSPARENT TESTA 12 (TT12), an important gene involved in the transportation of proanthocyanidin precursors into the vacuole. Two of eight repeatedly detected SNPs can be identified and digested by restriction enzymes. Candidate gene mining revealed that the relevant regions of significant SNP loci on the A09 and C08 chromosomes are highly homologous. Moreover, a comparison of the GWAS results to those of previous quantitative trait locus (QTL) studies showed that 11 SNPs were located in the confidence intervals of the QTLs identified in previous studies based on linkage analyses or association mapping. Our results provide insights into the genetic basis of seed coat color in B. napus, and the beneficial allele, SNP information, and candidate genes should be useful for selecting yellow seeds in B. napus breeding.