Skip to main content

New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus

This article has been updated

Abstract

Background

Increasing seed oil content is the most important breeding goal in Brassica napus, and phenotyping is crucial to dissect its genetic basis in crops. To date, QTL mapping for oil content has been based on whole seeds, and the lipid distribution is far from uniform in different tissues of seeds in B. napus. In this case, the phenotype based on whole seeds was unable to sufficiently reveal the complex genetic characteristics of seed oil content.

Results

Here, the three-dimensional (3D) distribution of lipid was determined for B. napus seeds by magnetic resonance imaging (MRI) and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues, the outer cotyledon (OC), inner cotyledon (IC), radicle (R) and seed coat (SC), which explained up to 13.76% of the phenotypic variation. Notably, 14 tissue-specific QTLs were reported for the first time, 7 of which were novel. Moreover, haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Furthermore, tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in the IC, OC and R than in the SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in the QTLs for OC and IC.

Conclusions

The present study provides further insight into the genetic basis of seed oil content at the tissue-specific level.

Highlight

Three-dimensional phenotyping, QTL identification with a high-density genetic linkage map and transcriptomic analysis were combined to dissect the genetic architecture of oil content in different seed tissues of Brassica napus.

Background

Rapeseed (Brassica napus) is one of the most important oilseed crops in the world. Increasing the oil content is the most important goal for rapeseed breeding at present [1, 2]. Comprehensively understanding the oil content and lipid distribution in seeds of B. napus is important for breeding. Although some genes that encode lipid biosynthesis in plants have been identified [3,4,5,6,7,8], the genetic basis and molecular regulation of lipid variation in different seed tissues in B. napus are poorly understood.

The analytical determination of oil content in rapeseed has been traditionally accomplished by the Soxhlet-based extraction method and residual method with organic solvents as the extraction medium [9], and supercritical fluid extraction was performed with carbon dioxide as the extraction solvent [10]. The seed lipid could also be analyzed using gas chromatography [11]. The abovementioned oil content measurement methods are all carried out after oil dissection or extraction [9,10,11]. However, the destruction of seeds in these methods is inappropriate if these materials need to be cultivated for other studies after oil content measurement [12]. Thus, some rapid, simple, and nondestructive methods, such as pulsed nuclear magnetic resonance [13] and near-infrared spectroscopy [14] have been developed for the determination of oil content and have been widely used for the evaluation of the oil content in seeds [12, 15,16,17,18]. However, the measurement and visualization of oil for specific seed tissues cannot be accomplished based on the two above methods. The conventional procedures used for lipid visualization, such as Nile red and Sudan black B [19], are based on dyes in two dimensions and are destructive [19]. Mass spectrometry imaging (MSI) has been used to visualize the lipid composition in both plant and animal material [20,21,22]. However, quantification in MSI is still a major hurdle, and the methods are destructive [23]. Recently, an approach that allows the rapid, noninvasive and quantitative visualization of lipids in living seeds using frequency-selected MRI has been reported [24, 25], and this method could provide quantitative lipid maps with a resolution close to the cellular level [26]. Based on MRI, each radiofrequency (RF) pulse excites a narrow slice in the two-dimensional (2D) sequence, whereas each RF pulse excites the entire imaging volume, and encoding is used to discriminate spatially in a three-dimensional (3D) sequence, which is helpful for fully revealing the lipid distribution in seeds and providing more detailed 3D lipid quantitative information.

Oil content is a typical quantitative trait that is controlled by a large number of genes with mainly additive and few epistatic gene actions [27,28,29,30]. QTL mapping is an effective approach to dissect the genetic mechanisms of such complex quantitative traits [31], and numerous QTLs for oil content have been identified in B. napus [27, 28, 32,33,34,35,36,37,38,39,40,41,42]. Lipid biosynthesis is complex and regulated by various biosynthetic enzymes and transcription factors [43], which vary dramatically from tissue to tissue within seeds [44]. Therefore, dissecting the genetic basis for oil content differences in different tissues of seeds is of great significance for further increasing the oil content in rapeseed. Recently, Lu et al. [45] revealed the complex spatial distribution patterns of lipids and transcripts by integrating the results from matrix-assisted laser desorption/ionization-mass spectrometry imaging (MALDI-MSI) of lipids in situ, lipidome profiling of extracts from seed tissues, and tissue-specific transcriptome analysis [45]. However, little is known about the spatial distribution patterns and genetic mechanisms of lipids in different seed tissues in B. napus at the population level.

Phenotyping is crucial to reveal genetic variation accompanied by genomics in crops. With the development of genome sequencing technology, precise assessment of plant phenotypes is falling behind the ability to characterize genotypes [46,47,48]. To date, the phenotypic data of oil content for QTL mapping in rapeseed have all been based on whole seeds [28, 39,40,41,42, 49]. As a result, few new loci have been identified and utilized. If the oil content of different tissues is separated as an independent phenotype for QTL mapping, it would reveal the genetic mechanism of lipid accumulation differences in different tissues of seeds and provide new ideas for breeding to improve oil content. Recently, linkage mapping combined with high-throughput sequencing has been applied to dissect genetic mechanisms and is considered to be an effective way to identify candidate genes. For example, 24 genes were identified by combining SNP-trait association analysis and transcriptome sequencing for resistance to Sclerotinia stem rot in B. napus [50]. DEGs were integrated with association mapping and linkage analysis to confirm their roles in the growth process, and 12 candidate genes associated with growth period traits were found [51, 52]. BnaC07. CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of seed coat content by modulating lignin biosynthesis by combining transcriptome-wide association studies and correlation networks of seed coat content-related gene modules.

The present study visualized the 3D distribution of lipids and obtained the actual lipid content in different tissues of seeds in a double haploid mapping population of B. napus using the noninvasive MRI technique in a 30 μm × 30 μm × 30 μm voxel with a 3D-FISP sequence. A total of 10 traits related to oil content in four subdivided tissues (seed coat (SC), outer cotyledon (OC), inner cotyledon (IC), and radicle (R)) of seeds were obtained. Then, QTL mapping for oil content-related traits was implemented for the four subdivided tissues of seeds. Furthermore, tissue-specific transcriptome sequencing from two developmental stages (24 and 33 days after flowering (DAF)) of two parents with high oil content differences was employed to detect the gene expression characteristics of IC, OC, R and SC and assist in candidate gene identification. Finally, by combining QTL mapping, bulked segregant analysis (BSA) and tissue-specific transcriptomic analysis, common and specific candidates involved in lipid biosynthesis in different seed tissues were identified. This study provides new insights for further improving the oil content of seeds in B. napus.

Results

Using MRI and 3D restructuring to dissect the tissue-specific distribution of lipid in B. napus seeds

Mature B. napus seeds with high (51.47%), middle (46.90%) and low (40.10%) oil content were selected to analyze the lipid distribution using MRI. Images of lipids with a voxel resolution of 30 µm × 30 µm × 30 µm were obtained (Fig. 1a–c and d–f and n–p), and the SC, OC, IC and R tissues could be distinguished clearly (Fig. 1h, i). The results showed that the lipid signal intensity in the four subdivided tissues was significantly different, and a clear reduction in the direction from the OC to IC, R and SC was observed (Fig. 1a–c and d–f).

Fig. 1
figure 1

Quantitative imaging of lipids in rapeseed seeds with different oil content. a, b and c represent the final lipid distribution shown as a gray value (signal intensity), and d, e and f represent the final lipid distribution of seeds shown as a color coded with high oil content seeds, middle oil content seeds and low oil content seeds, respectively. g, h and i represent the original rapeseed seed and the reconstruction of the rapeseed seed. j, k, l and m represent the reconstruction of the seed coat, radicle, inner cotyledon and outer cotyledon. n, o and p represent the distribution of lipids in different tissues of seeds with high, middle and low oil content; 1, 2, 3, 4, 5 and 6 represent seed coat, outer cotyledon, inner cotyledon, radicle, whole seed and seed section, respectively. sc seed coat, oc outer cotyledon, ic inner cotyledon, r radicle

To better understand the distribution of lipid and calculate oil content in different tissues in seeds, 3D reconstruction was implemented combining MRI images and Amira 5.4.0 software. The whole seed (Fig. 1g, h) could be clearly separated into SC (Fig. 1j), OC (Fig. 1m), IC (Fig. 1l) and R (Fig. 1k) at the 3D level, and the volume of each tissue could be calculated. The total oil content (TO) and mean oil content (MO) of different tissues could be calculated based on the calibration curve (r2 = 0.99473) established based on the lipid signal intensity and lipid content measured by gas chromatography. The lipid signal intensity in different tissues was different (Fig. 1a–f and Fig. 1n–p), which indicated that there was different oil content in different types of seed tissues. These results indicated that through MRI and 3D reconstruction, lipid density and content in four fine tissues of B. napus seed could be quickly and accurately quantified without destroying the seed structure.

Noninvasive 3D lipid phenotyping for four subdivided tissues in the KN segregating population

The abovementioned method was used for noninvasive phenotyping analysis for 200 DH lines in the KN population [49]. Finally, 10 novel traits related to oil content for different tissues of B. napus seeds, including OCTO, ICTO, RTO, SCTO, OCMO, ICMO, RMO, SCMO, WSTO and WSMO, were obtained. The detailed operating procedures of the MRI and Amira 5.4.0 are provided in the Methods section and supplementary video (Additional file 3: Video S1; Fig. 2). Based on the novel phenotypic data, a wide range of variation and transgressive segregation for 10 traits were observed in the KN population (Additional file 1: Fig. S2; Additional file 2: Table S1; Fig. 3). WSTO ranged from 0.79 mg to 3.66 mg with an average of 2.07 mg. For the four subdivided tissues, OC had the highest oil content (average of 1.28 mg, ranging from 0.39 to 2.27 mg), followed by IC (average of 0.60 mg, ranging from 0.18 to 1.12 mg), R (average of 0.14 mg, ranging from 0.04 to 0.27 mg) and SC (average of 0.058 mg, ranging from 0 to 0.19 mg). For the MO in the KN population, OC had the highest mean oil content, with an average of 0.63 mg/mm3 (ranging from 0.21 mg/mm3 to 1.11 mg/mm3), followed by IC (an average of 0.52 mg/mm3, ranging from 0.19 to 0.92 mg/mm3), R (an average of 0.254 mg/mm3, ranging from 0.06 to 0.71 mg/mm3) and SC (an average of 0.113 mg/mm3, ranging from 0 to 0.39 mg/mm3). The frequency distribution of 10 novel traits showed continuous variation in the KN population (Fig. 3), suggesting a polygenic effect of oil content in different seed tissues. The correlations among these 10 traits in the KN population were also calculated (Additional file 2: Table S2; Additional file 1: Fig. S3). Regardless of the TO or MO, the correlation between OC, IC and the whole seed was the strongest. In particular, OCTO and ICTO had high correlations of 0.962 and 0.832 with WSTO and WSMO, respectively.

Fig. 2
figure 2

Quantitative imaging of lipids in rapeseed seeds of the KN DH population. a, b, c, and d represent the imaging process of lipids in seeds, and e and f represent the three-dimensional reconstruction of seeds based on lipid imaging. g Three-dimensional reconstruction of different seed tissues, including the seed coat, outer cotyledon, inner cotyledon, and radicle. h represents the whole seed and segmentation of different issues of seeds based on lipid imaging

Fig. 3
figure 3

Frequency distribution of oil content in different tissues of seeds in the KN population. SCMO seed coat mean oil content, mg/mm3; ICMO inner cotyledon mean oil content, mg/mm3; OCMO outer cotyledon mean oil content, mg/mm3; RMO radicle mean oil content, mg/mm3; SCTO seed coat total oil content, mg; ICTO inner cotyledon total oil content, mg; OCTO outer cotyledon total oil content, mg; RTO radicle total oil content, mg; WSTO whole seed oil content, mg. The distribution density is overlaid

QTL mapping for oil content-related traits of subdivided tissues in B. napus seeds based on 3D phenotyping

QTL mapping for the ten oil content-related traits (OCTO, ICTO, RTO, SCTO, OCMO, ICMO, RMO, SCMO, WSTO and WSMO) was implemented based on the above phenotyping data and by combining a high-density genetic map containing 3207 markers [62]. A total of 35 QTLs, which could explain 5.22–13.76% of the phenotypic variation (PV), were identified and mainly distributed on chromosomes A1, A8, A9, A10, C1, C3 and C9 (Table 1; Fig. 4a; Additional file 1: Fig. S4). Thirteen QTLs were identified for TO, of which 2 QTLs for WSTO (qWSTO-A9-1 and qWSTO-A9-2) on the A9 chromosome had a PV of more than 10% and had a relatively large additive effect (> 0.2) (Table 1; Fig. 4a). In addition, 4, 2, 4 and 1 for OCTO, ICTO, RTO and SCTO, respectively, were located on A9 with relatively lower PV and additive effect (Table 1; Fig. 4b). Five, five, seven, three and two QTLs were identified for MO of the whole seed and subdivided tissues (WSMO, OCMO, ICMO, RMO and SCMO) and were distributed on the chromosomes of A1, A8, A9 and C3. Most MO-QTLs (12 of 22) were distributed on A9 (Fig. 4b).

Table 1 Identified QTLs and unique QTLs identified by high-density linkage map for oil content in different tissues of seeds in the KN DH population
Fig. 4
figure 4

QTLs for oil content in different tissues of seeds in the KN DH population. a Distribution of identified QTLs for oil content in different tissues of seeds in each linkage group. b, c, and d represent the distribution of common and specific QTLs for the oil content of different seed tissues in the A9, A10 and C3 linkage groups, respectively. WOC refers to identified QTLs for relative oil content detected by near-infrared spectroscopy

The above 35 QTLs were integrated into 23 unique QTLs by meta-analysis (Table 1), and 8 unique QTLs had pleiotropic effects. Four (uqA9-5, uqA9-7, uqA9-10 and uqA9-11) were responsible for MO and TO simultaneously and were all located on A9, and four unique QTLs (uqA1-2, uqA8-1, uqA9-2, and uqA9-12) were responsible for MO or TO of the same trait and were located on A1, A8, A9 and A9, respectively (Table 1). uqA1-2, uqA8-1 and uqA9-11 affected the OCMO, ICMO and WSMO simultaneously, which was consistent with the phenotypic performances that OC and IC had the largest effect on WSMO. Further analysis revealed that the pleiotropic QTLs uqA8-1, uqA9-7, uqA9-10 and uqA9-12 that responded to IC, OC and WS were also identified in our previous studies based on the oil content of whole seeds [62] (Table 1). Furthermore, uqA9-11, which was identified to control RTO, OCMO and ICMO simultaneously, was a novel locus with pleiotropic effects on the oil content of different tissues.

Fourteen unique QTLs that were identified to only control MO and/or TO in one tissue were considered to be tissue-specific QTLs and would provide new potential loci or genes for breeding to improve seed oil content (Table 1; Fig. 4a). Furthermore, uqA9-1 and uqA9-9, which were two adjacent QTLs, were identified for MO and TO of the R, respectively, while uqA9-3 and uqA9-4 were identified for TO and MO of the SC, respectively (Fig. 4b); uqA1-3 and uqA9-6 were identified for MO of the IC, while uqA8-2 and uqA9-8 were identified for MO of the OC. In addition, both uqC1-1 and uqC1-2 controlled the TO of the IC, while uqA10-1 and uqA10-2 were identified to be specific for TO of the OC (Fig. 4c); uqC3 and uqC9 were specific for RMO and RTO, respectively (Fig. 4a and d). All these QTLs were thought to have some effect on oil accumulation in corresponding tissues of the seeds, and these tissue-specific QTLs were the first reported in B. napus. Moreover, combined with BSA results for oil content reported in our previous study [62], it was revealed that many common significant variation loci revealed by BSA analysis were also identified in the present study (15 QTLs, most of them for OC and IC) (Table 2). For instance, uqA8-1 controlling the OCMO, ICMO and WSMO was identified in BSA analysis and colocalized with a QTL (cqOC-A8-3) identified for oil content of the whole seed (Fig. 5; Table 2).

Table 2 AGRs overlapping with QTLs for oil content from the QTL mapping method based on a high-density genetic linkage map
Fig. 5
figure 5

uqA8-1 controlling OCMO, ICMO, and WSMO was also identified on BSA. Original QTL identification shown by curves above the line of the linkage group (top). The Δ (SNP index) plot with statistically significantly associated regions (two green significant threshold lines, P < 0.01 and P < 0.05) is drawn at the bottom. In the two coordinate axes, the X-axis represents the position of linkage group A8 and chromosome A08, and the Y-axis represents the LOD value and Δ (SNP index). The additive effect was also shown under linkage Group A8. In the middle, the CI and peak of qOCMO-A8-1, qICMO-A8-1, qWSMO-A8-1, qOCMO-A8-2 and qWOC-A8-2 are shown. WOC represents QTLs identified for relative oil content. The middle connecting line links QTLs and their corresponding AGRs

Favorable alleles for oil content from different tissues exhibit cumulative effects

To understand the effects of QTLs for different tissues on the oil content of the whole seed, allelic genotypes in TO-QTLs and phenotypes for different tissues and the whole seed were assessed based on DH lines from the KN population (Additional file 2: Table S3). The results revealed that the TO of tissues (IC, OC and SC) and WSTO gradually increased with an increasing number of favorable haplotypes (Fig. 6). DH lines with favorable haplotypes in three QTLs, including two ICTO-QTLs (qICTO-C1-1 and qICTO-C1-2) and one OCTO-QTL (qOCTO-A10-2), showed higher OCTO and WSTO than DH lines without favorable haplotypes and DH lines with favorable haplotypes in only two ICTO-QTLs (qICTO-C1-1 and qICTO-C1-2) (Fig. 6b, c). In addition, the SCTO changed insignificantly when QTLs with favorable haplotypes involved in RTO, ICTO and OCTO aggregated, but those DH lines that harbored 11 QTLs with favorable haplotypes involved in SCTO, RTO, ICTO and OCTO showed a higher SCTO (Fig. 6d). Most importantly, more QTL accumulation caused a significant increase in WSTO, especially when OCTO-QTLs were affiliated (Fig. 6e). These results suggested that the genetic control of the seed oil content in B. napus exhibits an additive effect, and accumulation of OCTO played a key role in the improvement of total seed oil content.

Fig. 6
figure 6

Cumulative effects of tissue-specific TO-QTLs. The numbers on the ordinate represent the number of accumulating QTLs with favorable haplotypes. 0 represents the phenotype value of DH lines with undesirable haplotypes at all 12 assessed TO-QTLs. 1 represents DH lines with a favorable haplotype in only one QTL (qRTO-C9-2). 2 represents DH lines with favorable haplotypes in two QTLs (qICTO-C1-1 and qICTO-C1-2). 3 represents DH lines with favorable haplotypes in three QTLs (qICTO-C1-1, qICTO-C1-2 and qOCTO-A10-2). 11 represents those DH lines that harbored 11 QTLs with favorable haplotypes involved in SCTO, RTO, ICTO and OCTO. a, b, c and d show the phenotypic values of RTO, ICTO, OCTO and SCTO variation in DH lines with different types and numbers of favorable alleles, respectively. e shows the WSTO variation along with the accumulation of different types and numbers of favorable alleles

Characteristics of gene expression in different seed tissues

A total of 48 RNA samples from four seed tissues (SC, R, IC and OC) of Ken-C8 (with lower oil content) and N53-2 (with higher oil content) at 24 and 33 DAF (with three biological replicates) were collected for transcriptome sequencing, and 2675 million clean reads, with an average of 55.72 million reads per sample, were generated (Additional file 2: Table S4). The total mapped reads ranged from 35.08 to 66.27 million reads, with 61.00 million reads per sample on average, and the average unique mapping rate was 86.00%. All 48 samples showed high reproducibility (r2 > 0.9) among biological replicates in different tissues of Ken-C8 and N53-2 (Additional file 1: Fig. S5).

By comparing the 4 tissues at the 2 stages between Ken-C8 and N53-2, 15,064, 15,231, 15,286 and 26,429 DEGs were identified in IC, OC, R and SC, respectively. Further analysis revealed that 6637, 6083, 7107 and 10495 DEGs existed at 24 and 33 DAF in IC, OC, R and SC, respectively (Fig. 7a, b). Further analysis revealed that the IC, OC and R were clustered together instead of different lines and developmental stages together for each tissue based on the expression level of DEGs (Additional file 1: Fig. S6), which indicated that the IC, OC and R were highly similar and different from SC from the perspective of gene expression. In addition, the KEGG enrichment analysis for tissue-specific DEGs (Ken-C8 vs. N53-2) at 24 and 33 DAF in the IC, OC, R and SC showed that pyruvate metabolism was significantly enriched in IC, OC and R (Additional file 1: Fig. S7), which indicated a relatively greater tendency toward more pyruvate utilization for the accumulation of major compounds, including oil, protein and other carbohydrates, in these three tissues.

Fig. 7
figure 7

Differentially expressed gene analysis. a Venn diagrams show the common and different DEGs (Ken-C8 vs. N53-2) in the four tissues. b Venn diagrams show the common and different DEGs among different tissues at 24 and/or 33 DAP. c Clustering analysis based on the expression level of the DEGs involved in the TAG biosynthesis pathway, including plastidial fatty acid synthesis, fatty acid elongation, TAG synthesis, TAG degradation and beta-oxidation, among different samples

Furthermore, the DEGs involved in triacylglycerol (TAG) biosynthesis (including plastidial fatty acid (FA) synthesis, fatty acid elongation, TAG synthesis, TAG degradation and beta-oxidation) among different tissues were identified, and cluster analysis showed that the tissue samples could be divided into three distinct clusters: four SC samples from Ken-C8 and N53-2 collected at 24 and 33 DAF were clustered together (I), and samples of IC, OC and R collected at 33 (II) and 24 (III) DAF were clustered together (Fig. 7c). DEGs with relatively higher expression levels involved in plastidial fatty acid synthesis were found in cluster III, which indicated that fatty acid synthesis was more active in IC, OC and R at 24 DAF. For example, BnaC09. BCCP2 (the subunit of heteromeric ACCase) and BnaA01. LPD2 (the E3 component of the pyruvate dehydrogenase complex) showed higher expression levels in IC, OC and R at 24 DAF than at 33 DAF (Additional file 1: Fig. S8). Otherwise, genes involved in TAG synthesis with higher expression levels were also found in cluster II, which indicated that TAG synthesis was more active in IC, OC and R at 33 DAF, such as BnaA09g02110D and BnaC04g32530D (both encode the oil-body oleosin), which were more highly expressed in IC, OC and R at 33 DAF than at 24 DAF (Additional file 1: Fig. S8). In contrast, the β-oxidation of TAG tended to be more active in SC at 24 and 33 DAF (Additional file 1: Fig. S8).

To more intuitively understand the differences in the transcriptional regulation of lipid formation in IC, OC, R and SC between Ken-C8 and N53-2, a TAG biosynthesis- and accumulation-related pathway was reconstructed by integrating differential gene expression in different tissues according to RNA-Seq (Fig. 8). BnaC06. FatA and BnaC06. KASII were differentially expressed in all four tissues, and higher expression in N53-2 with higher oil content was observed. These key genes involved in the FA synthetic pathway showed significantly different expression levels in the four tissues, which suggested that various tissues of seeds in N53-2 have more FA synthesis activity than those of Ken-C8. The genes involved in the TAG synthetic pathway, such as BnaA03. GPDH, BnaA03. LPAAT5 and BnaC06. DGAT3 showed higher expression levels in the four tissues of N53-2, especially in IC and OC, which suggested a difference in TAG synthesis ability in different tissues between the two parents with high and low seed oil content. Interestingly, three copies of two important transcription factors, WRI1 and LEC1, were found to have significantly different expression only in the SC at 33 DAP, which suggested that the SC of N53-2 maintained a higher FA synthesis ability than that of Ken-C8 at the late stage of seed development, which was verified by the higher expression level of BnaC09. CAC2, BnaA02. HAD, BnaC02. HAD and BnaA02. KASII, which is regulated by WRI1 and LEC1 (Fig. 8a). This phenomenon also indicated that the high fatty acid synthesis activity of the tissue SC is still maintained in the late stage of seed development, which may be an important factor for the formation of high oil content.

Fig. 8
figure 8

TAG biosynthesis and accumulation-related pathways regulating the oil content. Gene expression differences between the two parents are shown in the pathway, and homologous copies showing differential expression in any tissue are listed next to the protein. The 24 and 33 DAP were defined by horizontal lines, and the positions of colored squares below the horizontal line are as follows: IC, OC, R and SC. The genes located in the QTL interval are marked with colored dots, and the color annotation is located in the top right corner. The expression difference of DEGs identified in the QTL interval for the four tissues

Identifying candidate genes for oil content QTLs in subdivided seed tissues

The candidate genes within QTL intervals were screened according to the relationship between the genetic map and the physical map. Overall, 861 DEGs between Ken-C8 and N53-2 at 24 and 33 DAF were found to be located in the interval of QTLs for different seed tissues. Fourteen of OC, 23 of IC, 19 of R and 8 of SC were located in the genomic region of QTL for IC, OC, R and SC, respectively (Fig. 8b). For example, BnaA09g43600D, which was highly expressed in the IC of Ken-C8 at 24 and 33 DAF, was located in the QTL qICMO-A9-2. BnaA10g09270D, which is located in the OC-specific QTL qOCTO-A10-1, showed higher expression in the OC of N53-2 at 24 and 33 DAF. BnaA09g35850D was differentially expressed in R at 24 and 33 DAF, and it was identified as underlying the R-QTL qRMO-A9-1. BnaA09g39710D (located in qSCTO-A9-2), BnaA09g40190D and BnaA09g40250D (located in qSCTO-A9-3) showed higher expression in the SC of Ken-C8 at 24 and 33 DAF.

Through annotation from the A. thaliana database (TAIR: http://www.arabidopsis.org/), 86 candidate genes within intervals of 19 unique QTLs were identified to be associated with lipid metabolism (Table 3). Three of these related genes (CAC2, KCS17 and FAE1) that are directly involved in FA synthesis and TAG metabolism were identified as pleiotropic genes for OC and IC. BnaA09g48250D, which is within uqA9-10 and affects OCTO and ICMO, was annotated as CAC2 (BnaA09. CAC2), the biotin carboxylase (BC) domain of ACC, which catalyzes the first committed step in fatty acid synthesis. BnaA09. CAC2 was significantly more highly expressed in the OC of N53-2 at 24 and 33 DAF compared to Ken-C8 (Fig. 8). BnaA08g11140D (homologous to AtKCS17) and BnaA08g11130D (homologous to AtFAE1), which mapped within the CI of uqA8-1, simultaneously control ICMO, OCMO and WSMO and encode 3-ketoacyl-CoA synthases involved in the biosynthesis of very long-chain fatty acids. It was exciting that as many as ten candidate genes involved in phospholipase, aliphatic suberin synthesis, lipase, sphingolipid synthesis, fatty acid elongation and cuticular wax synthesis were found to be located in CI of uqA8-1, of which BnaA08g10850D encoding phospholipase showed different expression in IC between N53-2 and Ken-C8. In addition, differential expression of BnaA09g38540D, BnaA09g38630D and BnaA09g39290D (homologous to AtG3, AtCDS5 and AT3G61580, located in the overlapping CI of qRMO-A9-2 and qSCMO-A9-2), which are involved in cuticular wax synthesis, sulfolipid synthesis and sphingolipid synthesis, was observed in all tissues.

Table 3 Identification of candidate genes for oil content QTLs in different seed tissues in the KN DH population

In addition, 54 genes (62.8%) were identified as the underlying genes of unique QTLs having a specific effect on specific seed tissues (Table 3), which revealed the specific function of some genes participating in oil synthesis and accumulation in corresponding tissues. For instance, three oil body-related genes, BnaA01g14480D, BnaA10g09480D and BnaC01g34070D (homologous to AtOLEO1, AtPXG2 and AT3G18570), were mapped to WSMO, OCTO and ICTO, respectively. Homologous genes of AT2G22230, AtACP4 and AtKCS4 involved in fatty acid synthesis and elongation located in QTL CIs were identified as underlying genes related to MO, and BnaC01g33800D and BnaA10g09480D (homologous to AtLPAT5 and ATPXG2 involved in TAG synthesis) were identified to affect ICTO and OCTO, respectively. In addition, 14 candidate genes were mapped to be specifically related to oil accumulation in R seeds, including the TF BnaC03. ABI3 is involved in regulating the transition between embryo maturation and early seedling development. BnaA09g40180D, BnaA09g39950D and BnaA09g40500D, which are annotated as being involved in TAG degradation, lipase and lipid signaling, were identified to be specific for oil content in the SC of seeds.

Discussion

3D phenotyping of oil content in different tissues of seeds provides a novel strategy for genetic dissection for seed oil content

Understanding the genetic basis of oil formation is necessary for increasing the oil content in rapeseed. Research has revealed that the distribution of lipid is far from uniform in different seed tissues [24]. Therefore, dissecting the genetic basis for oil content formation in different seed tissues has great practical significance for rapeseed breeding. Phenotyping is crucial to reveal natural genetic variation accompanied by genomics in crops, so progress in phenotyping technologies is required to accelerate genetic mapping and gene discovery [53]. Most of the current methods for the determination of oil content in rapeseed are either destructive [9,10,11, 45, 54] or difficult to realize quantitative visualization of lipids in seeds [12,13,14, 40,41,42]. MRI has been successfully used to evaluate the relative contribution of different tissues to the overall visualization and quantity of lipid accumulated in various seeds, including B. napus seeds [24, 26, 55], which provides possibilities for QTL mapping based on the oil content of different seed tissues. In this study, by combining lipid MRI and 3D quantitative analysis, we reveal the lipid distribution, which overcomes many drawbacks existing in oil content determination, such as the influence of the gap between the slices on lipid quantitation and the destruction of seeds in the follow-up experiment [11, 24, 26, 45], and precision is compromised when the seed shape is not the same as the model seed for imaging [14, 18]. Here, our study provides a novel method for the acquisition of phenotypic data of oil content for QTL mapping and the dissection of the genetic basis for oil content in different seed tissues.

Conventional breeding for high oil content has always relied on improvements in the capacity of the organs appropriate for lipid storage [56,57,58]. Lipid maps showed that the deposit level of oil was significantly different in different tissues, even within a tissue, which indicated that the breeding ideas can be further expanded. Thus, the currently rendered lipid mapping of the different tissues of seeds allows the identification of more precise ‘topographical targets’ for breeding in the future. Since the cotyledons made the largest contribution to the oil content, the seed oil content could be increased more if the oil content in the dark areas of OC and IC was improved in B. napus (Fig. 1a–c). In addition, the proportion of oil content in SC and R was 18.58% and 10.52% of the volume in the whole seed, respectively; however, the proportion of oil content in SC and R was 2.54% and 5.78% in N53-2 with high seed oil content (51.74%) (Additional file 1: Fig. S1). The results suggest that the oil content could increase by further increasing the oil content of cotyledons and radicles or selecting seeds with thin seed coats and larger cotyledons (especially outer cotyledons) to reach a level equivalent to that achieved in cotyledons, as suggested by Borisjuk et al. [26]. Even though the dark areas of the seeds were exits of the vascular system for directing the flow of water [59], after a preliminary calculation, we deduced that the oil content of rapeseed varieties with 60% might be expected to reach up to 75% when the oil content in the middle dark zone of the cotyledon reached the average level of lipids in the surrounding area of the cotyledon.

New perspective for QTL mapping for oil content in B. napus seeds

In recent years, numerous QTLs for oil content have been identified in B. napus [32,33,34,35,36,37]. These QTLs explained PV ranging from 1.2% to approximately 20% [36, 60, 61], while in previous reports, all QTLs for seed oil content were based on the oil content of the whole seeds. To our knowledge, this genetic dissection of seed oil content based on QTL mapping in the different tissues of seeds is the first reported. More QTLs for oil content were identified when the whole seed was subdivided into four tissues (Table 1). Four QTLs (uqA8-1, uqA9-7, uqA9-10 and uqA9-12) were identified corresponding to QTLs that were identified based on the same population and genetic map [62]. These QTLs were related to MO or TO of WS, OC or IC, in which uqA9-7 and uqA9-12 were the two unique QTLs that control the only two WSTO-QTLs (qWSTO-A9-1 and qWSTO-A9-2), respectively. In addition, the QTL genomic regions of two unique QTLs, uqA9-1 (qRMO-A9-1) and uqA9-2 (qRMO-A9-2 and qSCMO-A9-2), and four unique QTLs, uqA9-6 (qICMO-A9-2), uqA9-7 (qWSMO-A9-2 and qWSTO-A9-1), uqA9-7 (qOCMO-A9-2) and uqA9-7 (qRTO-A9-2), were mapped to be identical to TN-qOC-A9-3 and TN-qOC-A9-4, respectively [28]. Another two QTLs (qOCMO-A8-2 and qOCTO-A10-2, with small effects) were also identified in a previous study [62]. Especially important, 14 tissue-specific QTLs were identified, including some with major effects. In addition, uqA9-6, with up to 12.63% of PV, was a specific QTL for the average oil content of IC and colocalized with TN-qOC-A9-4. Moreover, the BSA results further confirmed the common QTLs identified in this study and Chao et al. [62]. These results indicated that QTL mapping based on MRI was accurate and credible, and it could provide new insights into variation loci of oil content derived from the contribution of which tissues when more detailed information was provided. Here, 11 QTLs were thought to be novel and were not identified in all previous studies, seven of which were tissue-specific QTLs. Among them, uqA9-11 was considered to be a main QTL controlling ICMO, OCMO, RTO and WSMO. In our study, it was revealed that accumulating favorable alleles for TO in different tissues, especially the OC, may be an effective way to increase the WSTO. Adding favorable alleles for RTO and OCTO could significantly increase the WSTO (Fig. 6e). The cumulative effects of QTLs of different tissues thus revealed that seed oil content is partly controlled by tissue-specific quantitative genes with mainly additive effects. Our study indicates that tissue-specific selection for seed oil content should be effective in breeding for oil content enhancement, which further emphasizes the urgent need to develop tissue-specific molecular markers for molecular-assisted selection breeding of high oil content varieties.

Analysis of candidate genes demonstrates the application value of 3D phenotyping strategies

Transcriptome analysis for each seed tissue was performed to offer further support for candidate identification in combination with tissue-specific QTL mapping. In terms of TAG biosynthesis-related genes, IC, OC and R showed the same gene expression characteristics, FA synthesis was more active at 24 DAP, and TAG synthesis was more active at 33 DAF (Fig. 7 and Additional file 1: Fig. S8), which reflected lipid synthesis regularity in the developmental process of seeds. Some important candidate genes related to lipid synthesis and metabolism in the CIs of QTLs were identified in the QTL interval. CAC2, KCS17 and FAE1, which are directly involved in fatty acid (FA) synthesis and elongation, were identified in the CI of pleiotropic QTLs that control OC and IC. The BnaA09. CAC2, which is located in both OC-QTL and IC OC-QTL, showed higher expression levels in OC and IC of N53-2 at 24 and 33 DAF, which will promote the flow of precursors to FA synthesis in OC and IC.

Ten candidate genes involved in phospholipase, aliphatic suberin synthesis, lipase, sphingolipid synthesis, fatty acid elongation and cuticular wax synthesis, such as BnaA08. FAE1, BnaA08. KCS17, BnaA08. CER4 and BnaA08. LPT1, were mapped in CI of uqA8-1, which controls ICMO, OCMO and WSMO. More erucic acid was detected in the seeds of N53-2 than Ken-C8 [63], although BnaA08. FAE1, which has been reported to determine the erucic acid content [64], had not been detected to be differentially expressed in OC and IC between the two parents, and we are more likely to believe that it might be only one causal gene, BnaA08. FAE1 controlled ICMO, OCMO and WSMO, which needs to be further verified by fine mapping. Interestingly, three genes, BnaA09g38540D, BnaA09g38630D and BnaA09g39290D, in the overlapping CI of qRMO-A9-2 and qSCMO-A9-2 were annotated to be involved in cuticular wax synthesis, sulfolipid synthesis and sphingolipid synthesis, respectively, instead of TAG metabolism, which would provide new clues to further explore oil content variation in R and SC. In addition, the genes involved in TAG degradation, lipase and lipid signaling were identified to be specific for oil content in SC, which would provide some evidence for the high accumulation of specific fatty acids in the SC, such as C16:1 [26]. In the annotation of DEGs, SC showed greater differences than IC, OC and SC between the two parents, and important TFs and key structural genes, such as WRI1, LEC1/2, CAC2 and KASII, involved in FA synthesis were differentially expressed, which suggested that FA biosynthesis was significantly different in SC between the high oil content N53-2 and the low oil content Ken-C8. However, the data from 3D phenotyping analysis showed that no oil content difference in SC was found between the two parents (Additional file 2: Table S1). What is the reason for this? We also found that flavonoid and phenylpropanoid biosynthesis pathways were significantly enriched for SC-specific DEGs, which might indicate that the synthesized fatty acids are degraded for the synthesis of flavonoids simultaneously. This discovery provided new insights for candidate identification of oil content QTLs in SC.

In addition, we focused on the candidate genes within CIs of tissue-specific QTLs of IC and OC. BnaA01g14480D, BnaA10g09480D and BnaC01g34070D were mapped in the CIs of qWSMO-A1-1, qOCTO-A10-1 and qICTO-C1-1. They were homologous to three oil-body genes, AtOLEO1, AtPXG2 and AT3G18570, and might be related to the specific growth and building of oil bodies in OC and IC. A large number of candidate genes involved in FA and TAG synthesis were mapped to underlying oil variation of OC or IC, for example, AT2G22230, AtACP4, LPAT5 and ATPXG2, which demonstrates the utility of 3D phenotyping strategies for genetic dissection of the mechanism for oil formation of different tissues of seeds in B. napus.

Conclusions

The phenotype based on whole seeds was unable to sufficiently reflect the complex genetic characteristics of seed oil content. In this study, the 3D distribution of lipids was measured for B. napus seeds by MRI and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds into different tissues. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues of OC, IC, R and SC with up to 13.76% of the phenotypic variation, and 14 tissue-specific QTLs were the first reported, including 7 novel QTLs. Haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in IC, OC and R than in SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in QTLs for OC and IC. The present study provides further insight into the genetic basis for improving seed oil content at the tissue-specific level.

Materials and methods

Plant materials

Three B. napus lines from the KN DH population [49] with different oil content (51.74%, 46.90% and 40.10%) were selected for magnetic resonance imaging (MRI) analysis in different seed tissues. The KN DH population was derived from microspore culture of F1 plants after hybridization between ‘N53-2’ with high oil content (approximately 50%) and male ‘KenC-8’ with low oil content (approximately 40%) [49]. In addition, a random subset of 200 lines from the KN DH population was used for oil content analysis in different seed tissues using MRI. All materials were planted in a field located in Dali of Shaanxi Province (winter rapeseed planting area, E109.93°, N34.80°), China (sown in September 2014 and harvested in May of the following year). The field experiment was performed following standard breeding field protocols as previously described [49], and the seeds of all lines were collected after maturation.

Noninvasive 3D phenotypic analysis and total and mean oil content determination in different seed tissues

MRI experiments were performed as previously described [24] on a vertical 11.75 T wide bore (inner diameter 89 mm) Bruker Avance III spectrometer (Bruker BioSpin GmbH, http://www.bruker-biospin.com/, Rheinstetten, Germany) at 500.13 MHz and equipped with an actively shielded Micro 2.5 gradient system (inner diameter of 40 mm and maximum strength 1 T/m). For MRI analysis, the mature dry seeds were placed into a home-built 5-mm ID solenoid coil, and a global T1 measurement of the lipid signal was conducted to determine the minimum repetition time (TR) necessary to avoid T1 correction during the quantification process. Considering the gap between planes, it could not sufficiently reveal the lipid distribution and lipid quantitative information when using the 2D sequence in MRI [24, 26]. Here, a 3D-FISP sequence with a CHESS suppression module on the water resonance was applied to acquire high-resolution lipid images of the seeds. Eventually, a voxel resolution of up to 30 μm isotropic in the seeds could be achieved (experimental time 10 h 40 min).

To accelerate 3D phenotyping, a special scaffold was produced to hold 12 rapeseed seeds at the same time for imaging. After obtaining data, in-house software was used for data reconstruction. The visualization and manual segmentation for modeling were performed with the 3D visualization software package Amira 5.4.0 (Visage Imaging GmbH, http://www.visageimaging.com/) based on standard manipulation procedures [24]. After segmentation, the average signal and total signal intensity in each tissue was obtained (Additional file 3: Video SI). After imaging, the seeds (approximately 4 mg) were immersed in 1.5 ml of 2.5% H2SO4 methanol solution with 0.01% butylated hydroxytoluene (BHT), 0.4 ml methylbenzene, and 200 μl of 2 mg/ml C17:0 as an internal standard following nitrogen filling. The mixture was incubated at 90 °C for 1 h. Then, 1.8 ml of deionized water and 1 ml of hexane were added to the mixture, and the supernatant was transferred into vials for gas chromatography (GC) analysis. GC was performed as previously described [65], and the actual oil content was calculated by the flame ionization detector response of sample components relative to 17:0 internal standards.

The measured MRI signals of the seeds were plotted vs. the known total lipid content of the seeds through Microsoft Office Excel 2010. The fitted calibration curve was used to quantify the oil content in different tissues of the rapeseed seeds, providing quantitative information on lipid distribution.

Trait evaluation and statistical analysis, QTL mapping and meta-analysis

According to the method established above, the total lipid signal intensity, volume and mean signal intensity for different seed tissues were obtained for each line from a DH population (n = 200) after segmentation was conducted with Amira 5.4.0. After calibration, the total oil content (TO) and the mean oil content (MO) of each seed were determined. In total, 10 oil content-related characteristics in different tissues of seeds in the KN DH population were obtained, which included the whole seed oil content (WSTO), outer cotyledon total oil content (OCTO), inner cotyledon total oil content (ICTO), radicle total oil content (RTO), seed coat total oil content (SCTO), whole seed mean oil content (WSMO), outer cotyledon mean oil content (OCMO), inner cotyledon mean oil content (ICMO), radicle mean oil content (RMO) and seed coat mean oil content (SCMO). The statistical analysis of the phenotypic data of 10 traits was conducted using Microsoft Office Excel 2010 and SPSS 20.0 software.

A high-density genetic linkage map containing 101 non-SNP (SSR, STS) and 3106 SNP-bin markers constructed as previously described [62] was employed for QTL mapping. QTL identification was implemented using the composite interval mapping (CIM) model in Windows QTL Cartographer V2.5 [49]. Parameters were set for the CIM model as described by Chao et al. [62]. The LOD threshold was determined for significant QTLs by a 1,000-permutation test based on a 5% experimental error rate. Significant QTLs identified in each trait were termed ‘identified QTLs’. The method for assigning nomenclature to identified QTLs was described in Wang et al. [49]. The identified QTLs were named by combining the trait and the chromosome number, e.g., qOCTO-A10-1, representing the first QTL identified for OCTO on chromosome A10. QTLs for different traits with overlapping CIs were integrated into unique QTLs using the QTL meta-analysis method in the BioMercator V4.2 program [62].

Identification of genomic regions and candidate genes underlying QTLs for oil content

The collinearity analysis of the linkage map and the B. napus “Darmor-bzh” reference genome [66] was conducted based on SNP probes as described previously [67]. Genomic regions corresponding to the CI of QTLs were defined through both closely linked SNPs within QTL CIs and colinearity between the linkage map and the reference genome. If genes that fell within the QTL genomic region were annotated to be related to lipid metabolism, they were regarded as candidate genes. The genomic regions of QTLs from other populations, DY [36], RNSL [36], and TN [28], were obtained from electronic PCR (e-PCR) [68] for comparison with the QTLs identified in the present study, which was performed with the primer sequences of the molecular markers flanking the QTL confidence intervals using the genomic sequences of ‘Darmor-bzh’ as templates. The orthologs of candidate genes within QTL regions and their annotations were obtained by BLASTn based on the A. thaliana database (TAIR: http://www.arabidopsis.org/). In addition, BSA analysis based on extremely high and low oil content lines in the KN population was from a previous study [62].

RNA extraction, transcriptome sequencing and DEG identification

IC, OC, R and SC tissues of seeds N53-2 and Ken-C8 at 24 and 33 days after flowering (DAF) were separated using a stereomicroscope and then flash-frozen in liquid nitrogen. Three biological replicates were collected, and total RNA was extracted from each seed tissue using an RNAprep pure plant kit (DP432, http://www.tiangen.com/). Extracted RNA samples were sent Novogene Corporation (Beijing, China) for library construction and transcriptome sequencing on an Illumina HiSeq 2500 platform. Approximately, 6 Gb of Illumina cleaned reads were collected for each tissue per sample type. Low-quality reads, connectors, and barcode sequences were eliminated using Trimmomatic-0.39. All downstream analyses were based on clean data with high quality. Then, the clean data were aligned to the reference genome of ‘Darmor-bzh’ [66] using HISAT (v2.2.1, https://www.nature.com/articles/nmeth.3317). Gene expression levels were quantified as count number and FPKM using the programs featurecounts and cuffquant, respectively [69]. Differentially expressed genes (DEGs) were identified using the R package DEseq 2 based on the criteria of false discovery rate (FDR, Benjamini–Hochberg multiple test correction) < 0.01 and absolute fold change > 2 [70].

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Change history

  • 23 June 2023

    The error in the ESM file has been corrected.

References

  1. Wang H. Strategy for rapeseed genetic improvement in China in the coming fifteen years. Chin J Oil Crop Sci. 2004;26(2):98–101.

    Google Scholar 

  2. Hua W, Liu J, Wang H. Molecular regulation and genetic improvement of seed oil content in Brassica napus L. Front Agric Sci Eng. 2016;3(3):186–94.

    Article  Google Scholar 

  3. Ohlrogge J, Browse J. Lipid biosynthesis. Plant Cell. 1995;7(7):957.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Zheng P, Allen WB, Roesler K, Williams ME, Zhang S, Li J, Glassman K, Ranch J, Nubel D, Solawetz W, et al. A phenylalanine in DGAT is a key determinant of oil content and composition in maize. Nat Genet. 2008;40(3):367–72.

    Article  CAS  PubMed  Google Scholar 

  5. Baud S, Lepiniec L. Regulation of de novo fatty acid synthesis in maturing oilseeds of Arabidopsis. Plant Physiol Biochem. 2009;47(6):448–55.

    Article  CAS  PubMed  Google Scholar 

  6. Baud S, Lepiniec L. Physiological and developmental regulation of seed oil production. Prog Lipid Res. 2010;49(3):235–49.

    Article  CAS  PubMed  Google Scholar 

  7. Pouvreau B, Baud S, Vernoud V, Morin V, Py C, Gendrot G, Pichon JP, Rouster J, Paul W, Rogowsky PM. Duplicate maize wrinkled1 transcription factors activate target genes involved in seed oil biosynthesis. Plant Physiol. 2011;156(2):674–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Elhai N, Duncan RW, Stasolla C. Molecular regulation of seed oil accumulation. Adv Nutr Human Metab. 2016;2:e1296.

    Google Scholar 

  9. Hartwig RA, Hurburgh CR. Interlaboratory comparison of soybean protein and oil determinations. J Am Oil Chem Soc. 1991;68(12):949–55.

    Article  CAS  Google Scholar 

  10. Taylor SL, King JW, List GR. Determination of oil content in oilseeds by analytical supercritical fluid extraction. J Am Oil Chem Soc. 1993;70(4):437–9.

    Article  CAS  Google Scholar 

  11. Li Y, Beisson F, Pollard M, Ohlrogge J. Oil content of arabidopsis seeds: the influence of seed anatomy, light and plant-to-plant variation. Phytochemistry. 2006;67(9):904–15.

    Article  CAS  PubMed  Google Scholar 

  12. Rolletschek H, Fuchs J, Friedel S, Börner A, Todt H, Jakob PM, Borisjuk L. A novel noninvasive procedure for high-throughput screening of major seed traits. Plant Biotechnol J. 2015;13(2):188–99.

    Article  CAS  PubMed  Google Scholar 

  13. Tiwari PN, Gambhir PN, Rajan TS. Rapid and nondestructive determination of seed oil by pulsed nuclear magnetic resonance technique. J Am Oil Chem SOC. 1974;51(3):104–9.

    Article  CAS  Google Scholar 

  14. Sato T, Takahata Y, Noda T, Yanagisawa T, Morishita T, Sakai S. Nondestructive determination of fatty acid composition of husked sunflower (Helianthus annua L.) seeds by near-infrared spectroscopy. J Am Oil Chem Soc. 1995;72(10):1177–83.

    Article  CAS  Google Scholar 

  15. Guillén MD, Ruiz A. High resolution 1H nuclear magnetic resonance in the study of edible oils and fats. Trends Food Sci Technol. 2001;12(9):328–38.

    Article  Google Scholar 

  16. Tillman BL, Gorbet DW, Person G. Predicting oleic and linoleic acid content of single peanut seeds using near-infrared reflectance spectroscopy. Crop Sci. 2006;46(5):2121–6.

    Article  CAS  Google Scholar 

  17. Marcone MF, Wang S, Albabish W, Nie S, Somnarain D, Hill A. Diverse food-based applications of nuclear magnetic resonance (NMR) technology. Food Res Int. 2013;51(2):729–47.

    Article  CAS  Google Scholar 

  18. Fassio AS, Restaino EA, Cozzolino D. Determination of oil content in whole corn (Zea mays L.) seeds by means of near infrared reflectance spectroscopy. Comput Electron Agric. 2015;110:171–5.

    Article  Google Scholar 

  19. Titford M. Progress in the development of microscopical techniques for diagnostic pathology. J Histotechnol. 2009;32(1):9–19.

    Article  CAS  Google Scholar 

  20. Horn PJ, Korte AR, Neogi PB, Love E, Chapman KD. Spatial mapping of lipids at cellular resolution in embryos of cotton. Plant Cell. 2012;24(2):622–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Horn PJ, Chapman KD. Lipidomics in situ: insights into plant lipid metabolism from high resolution spatial maps of metabolites. Prog Lipid Res. 2014;54:32–52.

    Article  CAS  PubMed  Google Scholar 

  22. Li B, Dunham SJB, Dong Y, Yoon S, Zeng M, Sweedler JV. Analytical capabilities of mass spectrometry imaging and its potential applications in food science. Trends Food Sci Technol. 2016;47:50–63.

    Article  CAS  Google Scholar 

  23. Sturtevant D, Lee YJ, Chapman KD. Matrix assisted laser desorption/ionization-mass spectrometry imaging (MALDI-MSI) for direct visualization of plant metabolites in situ. Curr Opin Biotechnol. 2016;37:53–60.

    Article  CAS  PubMed  Google Scholar 

  24. Neuberger T, Rolletschek H, Webb A, Borisjuk L. Non-invasive mapping of lipids in plant tissue using magnetic resonance imaging. In: Armstrong D, editor. Lipidomics. Totowa: Springer; 2009. p. 485–96.

    Chapter  Google Scholar 

  25. Borisjuk L, Rolletschek H, Neuberger T. Surveying the plant’ s world by magnetic resonance imaging. Plant J. 2012;70(1):129–46.

    Article  CAS  PubMed  Google Scholar 

  26. Borisjuk L, Rolletschek H, Neuberger T. Nuclear magnetic resonance imaging of lipid in living plants. Prog Lipid Res. 2013;52(4):465–87.

    Article  CAS  PubMed  Google Scholar 

  27. Ecke W, Uzunova M, Weissleder K. Mapping the genome of rapeseed (Brassica napus L.). II. Localization of genes controlling erucic acid synthesis and seed oil content. Theor Appl Genet. 1995;91(6):972–7.

    Article  CAS  PubMed  Google Scholar 

  28. Jiang C, Shi J, Li R, Long Y, Wang H, Li D, Zhao J, Meng J. Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.). Theor Appl Genet. 2014;127(4):957–68.

    Article  CAS  PubMed  Google Scholar 

  29. Si P, Mailer RJ, Galwey N, Turner DW. Influence of genotype and environment on oil and protein concentrations of canola (Brassica napus L.) grown across southern Australia. Aust J Agric Res. 2003;54(4):397–407.

    Article  Google Scholar 

  30. Liu S, Fan C, Li J, Cai G, Yang Q, Wu J, Yi X, Zhang C, Zhou Y. A genome-wide association study reveals novel elite allelic variations in seed oil content of Brassica napus. Theor Appl Genet. 2016;129(6):1203–15.

    Article  CAS  PubMed  Google Scholar 

  31. Mauricio R. Mapping quantitative trait loci in plants: uses and caveats for evolutionary biology. Nat Rev Genet. 2001;2(5):370–81.

    Article  CAS  PubMed  Google Scholar 

  32. Lionneton E, Ravera S, Sanchez L, Aubert G, Delourme R, Ochatt S. Development of an AFLP-based linkage map and localization of QTLs for seed fatty acid content in condiment mustard (Brassica juncea). Genome. 2002;45(6):1203–15.

    Article  CAS  PubMed  Google Scholar 

  33. Delourme R, Falentin C, Huteau V, Clouet V, Horvais R, Gandon B, Specel S, Hanneton L, Dheu JE, Deschamps M, et al. Genetic control of oil content in oilseed rape (Brassica napus L.). Theor Appl Genet. 2006;113(7):1331–45.

    Article  CAS  PubMed  Google Scholar 

  34. Rathke GW, Behrens T, Diepenbrock W. Integrated nitrogen management strategies to improve seed yield, oil content and nitrogen efficiency of winter oilseed rape (Brassica napus L.): a review. Agr Ecosyst Environ. 2006;117(2–3):80–108.

    Article  CAS  Google Scholar 

  35. Zhao J, Becker HC, Zhang D, Zhang Y, Ecke W. Conditional QTL mapping of oil content in rapeseed with respect to protein content and traits related to plant development and grain yield. Theor Appl Genet. 2006;113(1):33–8.

    Article  CAS  PubMed  Google Scholar 

  36. Liu LZ, Li JN. QTL Mapping of oleic acid, linolenic acid and erucic acid content in Brassica napus by using the high density SNP genetic map. Sci Agric Sin. 2014;47:24–32.

    CAS  Google Scholar 

  37. Sun F, Liu J, Hua W, Sun X, Wang X, Wang H. Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus. Plant Sci. 2016;252:388–99.

    Article  CAS  PubMed  Google Scholar 

  38. Teh L, Möllers C. Genetic variation and inheritance of phytosterol and oil content in a doubled haploid population derived from the winter oilseed rape Sansibar×Oase cross. Theor Appl Genet. 2016;129(1):181–99.

    Article  CAS  PubMed  Google Scholar 

  39. Raboanatahiry N, Chao H, Guo L, Gan J, Xiang J, Yan M, Zhang L, Yu L, Li M. Synteny analysis of genes and distribution of loci controlling oil content and fatty acid profile based on QTL alignment map in Brassica napus. BMC Genomics. 2017;18(1):1–15.

    Article  Google Scholar 

  40. Xiao Z, Zhang C, Tang F, Yang B, Zhang L, Liu J, Huo Q, Wang S, Li S, Wei L, et al. Identification of candidate genes controlling oil content by combination of genome-wide association and transcriptome analysis in the oilseed crop Brassica napus. Biotechnol Biofuels. 2019;12(1):1–16.

    Article  Google Scholar 

  41. Tang S, Zhao H, Lu S, Yu L, Zhang G, Zhang Y, Yang Q, Zhou Y, Wang X, Ma W, et al. Genome-and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus. Mol Plant. 2021;14(3):470–87.

    Article  CAS  PubMed  Google Scholar 

  42. Yao M, Guan M, Yang Q, Huang L, Xiong X, Jan HU, Voss-Fels KP, Werner CR, He X, Qian W, et al. Regional association analysis coupled with transcriptome analyses reveal candidate genes affecting seed oil accumulation in Brassica napus. Theor Appl Genet. 2021;134(5):1545–55.

    Article  CAS  PubMed  Google Scholar 

  43. Linder CR. Adaptive evolution of seed oils in plants: accounting for the biogeographic distribution of saturated and unsaturated fatty acids in seed oils. Am Nat. 2000;156(4):442–58.

    Article  PubMed  Google Scholar 

  44. Neuhaus HE, Emes MJ. Nonphotosynthetic metabolism in plastids. Annu Rev Plant Biol. 2000;51:111.

    Article  CAS  Google Scholar 

  45. Lu S, Sturtevant D, Aziz M, Jin C, Li Q, Chapman KD, Guo L. Spatial analysis of lipid metabolites and expressed genes reveals tissue-specific heterogeneity of lipid metabolism in high-and low-oil Brassica napus L. seeds. Plant J. 2018;94(6):915–32.

    Article  CAS  PubMed  Google Scholar 

  46. Houle D, Govindaraju DR, Omholt S. Phenomics: the next challenge. Nat Rev Genet. 2010;11(12):855–66.

    Article  CAS  PubMed  Google Scholar 

  47. Cobb JN, Declerck G, Greenberg A, Clark R, McCouch S. Next-generation phenotyping: requirements and strategies for enhancing our understanding of genotype-phenotype relationships and its relevance to crop improvement. Theor Appl Genet. 2013;126(4):867–87.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Gong H, Xu D, Yuan J, Li X, Guo C, Peng J, Li Y, Schwarz LA, Li A, Hu B, et al. High-throughput dual-colour precision imaging for brain-wide connectome with cytoarchitectonic landmarks at the cellular level. Nat Commun. 2016;7(1):1–12.

    Article  Google Scholar 

  49. Wang X, Wang H, Long Y, Li D, Yin Y, Tian J, Chen L, Liu L, Zhao W, Zhao Y, et al. Identification of QTLs associated with oil content in a high-oil Brassica napus cultivar and construction of a high-density consensus map for QTLs comparison in B. napus. PLoS ONE. 2013;8(12):e80569.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Wei L, Jian H, Lu K, Filardo F, Yin N, Liu L, Qu C, Li W, Du H, Li J. Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus. Plant Biotechnol J. 2016;14(6):1368–80.

    Article  CAS  PubMed  Google Scholar 

  51. Wang T, Wei L, Wang J, Xie L, Li YY, Ran S, Ren L, Lu K, Li J, Timko MP, et al. Integrating GWAS, linkage mapping and gene expression analyses reveals the genetic control of growth period traits in rapeseed (Brassica napus L.). Biotechnol Biofuels. 2020;13(1):1–19.

    Article  Google Scholar 

  52. Zhang Y, Zhang H, Zhao H, Xia Y, Zheng X, Fan R, Tan Z, Duan C, Fu Y, Li L, et al. Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus. Genome Biol. 2022;23(1):1–22.

    Article  Google Scholar 

  53. Yang W, Guo Z, Huang C, Duan L, Chen G, Jiang N, Fang W, Feng H, Xie W, Lian X, et al. Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice. Nat Commun. 2014;5(1):1–9.

    Article  Google Scholar 

  54. Hu Z, Hua W, Zhang L, Deng L, Wang X, Liu G, Hao W, Wang H. Seed structure characteristics to form ultrahigh oil content in rapeseed. PLoS ONE. 2013;8(4):e62099.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Schwender J, Hebbelmann I, Heinzel N, Hildebrandt T, Rogers A, Naik D, Klapperstück M, Braun HP, Schreiber F, Denolf P, et al. Quantitative multilevel analysis of central metabolism in developing oilseeds of oilseed rape during in vitro culture. Plant Physiol. 2015;168(3):828–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Dyer JM, Stymne S, Green AG, Carlsson AS. High-value oils from plants. Plant J. 2008;54(4):640–55.

    Article  CAS  PubMed  Google Scholar 

  57. Godfray HCJ, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, Pretty J, Robinson S, Thomas SM, Toulmin C. Food security: the challenge of feeding 9 billion people. Science. 2010;327(5967):812–8.

    Article  CAS  PubMed  Google Scholar 

  58. Gunstone FD. Non-food use of vegetable oils. Lipid Technol. 2013;25(3):72.

    Article  Google Scholar 

  59. Munz E, Rolletschek H, Oeltze-Jafra S, Fuchs J, Guendel A, Neuberger T, Ortleb S, Jakob PM, Borisjuk L. A functional imaging study of germinating oilseed rape seed. New Phytol. 2017;216(4):1181–90.

    Article  CAS  PubMed  Google Scholar 

  60. Chen G, Geng J, Rahman M, Liu X, Tu J, Fu T, Li G, Mcvetty PBE, Tahir M. Identification of QTL for oil content, seed yield, and flowering time in oilseed rape (Brassica napus). Euphytica. 2010;175(2):161–74.

    Article  CAS  Google Scholar 

  61. Sun M, Hua W, Liu J, Huang S, Wang X, Liu G, Wang H. Design of new genome-and gene-sourced primers and identification of QTL for seed oil content in a specially high-oil Brassica napus cultivar. PLoS ONE. 2012;7: e47037.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Chao H, Wang H, Wang X, Guo L, Gu J, Zhao W, Li B, Chen D, Raboanatahiry N, Li M. Genetic dissection of seed oil and protein content and identification of networks associated with oil content in Brassica napus. Sci Rep. 2017;7(1):1–16.

    Article  Google Scholar 

  63. Bao B, Chao H, Wang H, Zhao W, Zhang L, Raboanatahiry N, Wang X, Wang B, Jia H, Li M. Stable, environmental specific and novel QTL identification as well as genetic dissection of fatty acid metabolism in Brassica napus. Front Plant Sci. 2018;9:1018.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Wang N, Wang Y, Tian F, King GJ, Zhang C, Long Y, Shi L, Meng J. A functional genomics resource for Brassica napus: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING. New Phytol. 2008;180(4):751–65.

    Article  CAS  PubMed  Google Scholar 

  65. Borisjuk L, Nguyen TH, Neuberger T, Rutten T, Tschiersch H, Claus B, Feussner I, Webb AG, Jakob P, Weber H. Gradients of lipid storage, photosynthesis and plastid differentiation in developing soybean seeds. New Phytol. 2005;167(3):761–76.

    Article  CAS  PubMed  Google Scholar 

  66. Chalhoub B, Denoeud F, Liu S, Parkin IAP, Tang H, Wang X, Chiquet J, Belcram H, Tong C, Samans B, et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;345(6199):950–3.

    Article  CAS  PubMed  Google Scholar 

  67. Cai G, Yang Q, Yi B, Fan C, David E, Jacqueline B, Zhou Y, Hector C. A complex recombination pattern in the genome of allotetraploid Brassica napus as revealed by a high-density genetic map. PLoS ONE. 2014;9(10):e109910.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Schuler GD. Sequence mapping by electronic PCR. Genome Res. 1997;7(5):541–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57(1):289–300.

    Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The work was supported by the National Natural Science Foundation of China (32072098 and 31871656) and the National Key Research and Development Program (2022YFD1200402).

Author information

Authors and Affiliations

Authors

Contributions

LG performed research, collection, curation of the date, and date analysis of overall study, and wrote the manuscript. HC participated in QTL analysis, and wrote the manuscript, YY participated in candidate genes analysis. HL and DH helped revised the manuscript. HW and WZ implemented part of the field trials. LZ provided some guides about modified manuscript. ML and CZ designed, led and coordinated the overall study, and provided guidelines for writing the paper. All the authors read and approved the final manuscript.

Corresponding authors

Correspondence to Chunyu Zhang or Maoteng Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Total and mean oil content in different tissues of rapeseed seeds with different oil content.andrepresent the total and mean oil content in different tissues of rapeseed seeds, respectively. Fig. S2. Quantitative imaging of lipids in different tissues of seeds in the KN DH population based on three-dimensional reconstruction. 1, 2, 3, 4, 5 and 6 represent the whole seed, seed coat, inner cotyledon, outer cotyledon, radicle and seed section, respectively. Fig. S3. Pearson correlation coefficients for trait pairs affecting the oil content of rapeseed seeds in the KN DH population. Fig. S4. Distribution of identified QTLs for oil content in different tissues of seeds in the A1, A8, A9, A10, C1, C3 and C9 linkage groups. WOC refers to identified QTLs for relative oil content detected by near-infrared spectroscopy. Fig. S5. The correlations among all 48 samples in different tissues at the two seed sampling stages of Ken-C8 and N53-2. Fig. S6. KEGG enrichment of the tissue-specific DEGsin the four tissues at 24 and 33 DAP. Fig. S7. The expression characteristics of genes involved in fatty acid synthesis, TAG synthesisand β-oxidationin four tissues of Ken-C8 and N53-2 at 24 and 33 DAF. Fig. S8. The expression characteristics of genes involved in fatty acid synthesis, TAG synthesisand β-oxidationin four tissues of Ken-C8 and N53-2 at 24 and 33 DAF.

Additional file 2: Table S1.

Mean values and phenotypic variation for oil content in different tissues of seeds in the KN DH population. Table S2. Pearson correlation coefficients for trait pairs affecting the oil content of seeds in the KN DH population. Table S3. The allele type of each DH line in 12 accessed TO-QTLs. Table S4. Characteristics of sequencing data for RAN-Seq.

Additional file 3: Video S1. 3D lipid phenotyping in different tissues of seeds in the mapping population.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guo, L., Chao, H., Yin, Y. et al. New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus. Biotechnol Biofuels 16, 88 (2023). https://doi.org/10.1186/s13068-023-02324-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-023-02324-0

Keywords