Data in support of genetic architecture of glucosinolate variations in Brassica napus

BnIR: A multi-omics database with various tools for Brassica napus research and breeding

In the post-genome-wide association study era, multi-omics techniques have shown great power and potential for candidate gene mining and functional genomics research. However, due to the lack of effective data integration and multi-omics analysis platforms, such techniques have not still been applied widely in rapeseed, an important oil crop worldwide. Here, we report a rapeseed multi-omics database (BnIR; http://yanglab.hzau.edu.cn/BnIR), which provides datasets of six omics including genomics, transcriptomics, variomics, epigenetics, phenomics, and metabolomics, as well as numerous "variation-gene expression-phenotype" associations by using multiple statistical methods. In addition, a series of multi-omics search and analysis tools are integrated to facilitate the browsing and application of these datasets. BnIR is the most comprehensive multi-omics database for rapeseed so far, and two case studies demonstrated its power to mine candidate genes associated with specific traits and analyze their potential regulatory mechanisms.

Fine mapping and candidate gene analysis of a seed glucosinolate content QTL, qGSL-C2, in rapeseed (Brassica napus L.)

QTL mapping and candidate gene analysis indicate that allelic variations in BnaC2.MYB28 resulted from homeologous exchange and determine difference in seed glucosinolate content. A low seed glucosinolate content has long been an important breeding objective in rapeseed improvement. However, the molecular mechanisms underlying seed GSL content variations remain to be elucidated in allotetraploid Brassica napus. Here, we developed a double haploid population from a cross between two B. napus accessions that possess relatively low, but significantly different seed GSL contents and identified a major QTL, qGSL-C2, on chromosome C02 that explains 30.88-72.87% of the phenotypic variation observed in five environments. Using near-isogenic lines, we further delimited qGSL-C2 to a physical region of 49 kb on the B. rapa chromosome A02 which is highly homologous to the target C02 interval. Among five candidate genes, BnaC2.MYB28, a homologue of the Arabidopsis MYB28 encoding a putative R2R3-MYB-type transcription factor functioning in aliphatic methionine-derived GSL synthesis, was most likely to be the target gene underlying the QTL. Sequence analysis revealed multiple insertion/deletion and SNP variations in the genomic region between the alleles of the NILs. Furthermore, the allelic variations in BnaC2.MYB28 in the natural B. napus population were significantly associated with seed GSL content. Remarkably, the phylogenetic analysis and sequence comparison suggested that while the BnaC2.MYB28 allele from the parental line G120 was inherited from B. oleracea BolC2.MYB28, its counterpart from the other parent, 9172, most likely evolved from B. rapa BraA2.MYB28 via possible homeologous exchange. Our study promotes greater understanding of the molecular regulation of seed GSL content and provides useful molecular markers for seed GSL improvement in B. napus.

Genetic architecture of glucosinolate variation in Brassica napus

The diverse biological activities of glucosinolate (GSL) hydrolysis products play significant biological and economical roles in the defense system and nutritional qualities of Brassica napus (oilseed rape). Yet, genomic-based study of the B. napus GSL regulatory mechanisms are scarce due to the complexity of working with polyploid species. To address these challenges, we used transcriptome-based GWAS approach, Associative Transcriptomics (AT), across a diversity panel of 288 B. napus genotypes to uncover the underlying genetic basis controlling quantitative variation of GSLs in B. napus vegetative tissues. Single nucleotide polymorphism (SNP) markers and gene expression markers (GEMs) associations identify orthologues of MYB28/HAG1 (AT5G61420), specifically the copies on chromosome A9 and C2, to be the key regulators of aliphatic GSL variation in leaves. We show that the positive correlation observed between aliphatic GSLs in seed and leaf is due to the amount synthesized, as controlled by Bna.HAG1.A9 and Bna.HAG1.C2, rather than by variation in the transport processes. In addition, AT and differential expression analysis in root tissues implicate an orthologue of MYB29/HAG3 (AT5G07690), Bna.HAG3.A3, as controlling root aromatic GSL variation. Based on the root expression data we also propose Bna.MAM3.A3 to have a role in controlling phenylalanine chain elongation for aromatic GSL biosynthesis. This work uncovers a regulator of homophenylalanine-derived aromatic GSLs and implicates the shared biosynthetic pathways between aliphatic and aromatic GSLs.