Genetic control of oil content in oilseed rape (Brassica napus L.)

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

Transcriptomic Analysis of the Reduction in Seed Oil Content through Increased Nitrogen Application Rate in Rapeseed (Brassica napus L.)

Nitrogen is essential for improving the seed oil yield of rapeseed (Brassica napus L.). However, the molecular mechanism by which increased nitrogen rates impact seed oil content is largely unknown. Therefore, a field experiment was conducted to determine how three nitrogen application rates (120, 240, and 360 kg ha-1) regulated seed oil content via transcriptomic analysis. The results showed that the seed yield and the protein and total N contents increased from N1 to N3, with average increases of 57.2%, 16.9%, and 79.5%, respectively. However, the seed oil content significantly decreased from N1 to N3, with an average decrease of 8.6%. These results were repeated over a number of years. The quantity of oil protein bodies observed under a transmission electron microscope was in accordance with the ultimate seed oil and protein contents. As the nitrogen application rate increased, a substantial number of genes involved in the photosynthesis, glycolysis, and phenylpropanoid biosynthesis pathways were up-regulated, as were TF families, such as AP2/ERF, MYB, and NAC. The newly identified genes were mainly involved in carbohydrate, lipid, and amino acid metabolism. Metabolic flux analysis showed that most of the genes involved in glycolysis and fatty acid biosynthesis had higher transcript levels in the early development stages. Our results provide new insights into the molecular regulation of rapeseed seed oil content through increased nitrogen application rates.

A review of sources of resistance to turnip yellows virus (TuYV) in Brassica species

Turnip yellows virus (TuYV; previously known as beet western yellows virus) causes major diseases of Brassica species worldwide resulting in severe yield-losses in arable and vegetable crops. It has also been shown to reduce the quality of vegetables, particularly cabbage where it causes tip burn. Incidences of 100% have been recorded in commercial crops of winter oilseed rape (Brassica napus) and vegetable crops (particularly Brassica oleracea) in Europe. This review summarises the known sources of resistance to TuYV in B. napus (AACC genome), Brassica rapa (AA genome) and B. oleracea (CC genome). It also proposes names for the quantitative trait loci (QTLs) responsible for the resistances, Turnip Yellows virus Resistance (TuYR), that have been mapped to at least the chromosome level in the different Brassica species. There is currently only one known source of resistance deployed commercially (TuYR1). This resistance is said to have originated in B. rapa and was introgressed into the A genome of oilseed rape via hybridisation with B. oleracea to produce allotetraploid (AACC) plants that were then backcrossed into oilseed rape. It has been utilised in the majority of known TuYV-resistant oilseed rape varieties. This has placed significant selection pressure for resistance-breaking mutations arising in TuYV. Further QTLs for resistance to TuYV (TuYR2-TuYR9) have been mapped in the genomes of B. napus, B. rapa and B. oleracea and are described here. QTLs from the latter two species have been introgressed into allotetraploid plants, providing for the first time, combined resistance from both the A and the C genomes for deployment in oilseed rape. Introgression of these new resistances into commercial oilseed rape and vegetable brassicas can be accelerated using the molecular markers that have been developed. The deployment of these resistances should lessen selection pressure for resistance-breaking isolates of TuYV and thereby prolong the effectiveness of each other and extant resistance.

Characterization of novel loci controlling seed oil content in Brassica napus by marker metabolite-based multi-omics analysis

BACKGROUND

Seed oil content is an important agronomic trait of Brassica napus (B. napus), and metabolites are considered as the bridge between genotype and phenotype for physical traits.

RESULTS

Using a widely targeted metabolomics analysis in a natural population of 388 B. napus inbred lines, we quantify 2172 metabolites in mature seeds by liquid chromatography mass spectrometry, in which 131 marker metabolites are identified to be correlated with seed oil content. These metabolites are then selected for further metabolite genome-wide association study and metabolite transcriptome-wide association study. Combined with weighted correlation network analysis, we construct a triple relationship network, which includes 21,000 edges and 4384 nodes among metabolites, metabolite quantitative trait loci, genes, and co-expression modules. We validate the function of BnaA03.TT4, BnaC02.TT4, and BnaC05.UK, three candidate genes predicted by multi-omics analysis, which show significant impacts on seed oil content through regulating flavonoid metabolism in B. napus.

CONCLUSIONS

This study demonstrates the advantage of utilizing marker metabolites integrated with multi-omics analysis to dissect the genetic basis of agronomic traits in crops.

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

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.

Identification of candidate genes regulating seed oil content by QTL mapping and transcriptome sequencing in Brassica napus

Increasing oil production is a major goal in rapeseed (Brassica napus) molecular breeding programs. Identifying seed oil content (SOC)-related candidate genes is an important step towards achieving this goal. We performed quantitative trait locus (QTL) mapping of SOC in B. napus using a high-density SNP genetic map constructed from recombinant inbred lines and the Illumina Infinium^TM 60K SNP array. A total of 26 QTLs were detected in three years on A01, A03, A05, A06, A09, C01, C03 and C05, which accounted for 3.69%~18.47% of the phenotypic variation in SOC. Of these, 13 QTLs are reported here for the first time. 1713 candidate genes in the 26 QTLs confidence interval were obtained. We then identified differentially expressed genes (DEGs) between the high- and low-SOC accessions, to narrow down our focus to 21 candidate genes (Y1-Y21) related to SOC, and we will focus on 11 (Y1-Y11) candidate genes that contribute to the formation of high-SOC. In addition to providing insight into the genetic basis of SOC in B. napus, the loci identified and candidate genes in this study can be used in molecular breeding strategies to increase SOC in this important seed crop.

Transcriptome and Regional Association Analyses Reveal the Effects of Oleosin Genes on the Accumulation of Oil Content in Brassica napus

Rapeseed stores lipids in the form of oil bodies. Oil bodies in the seeds of higher plants are surrounded by oleosins. Adjusting oleosin protein levels can prevent the fusion of oil bodies and maintain oil body size during seed development. However, oil contents are affected by many factors, and studies on the complex molecular regulatory mechanisms underlying the variations in seed oil contents of B. napus are limited. In this study, a total of 53 BnOLEO (B. napus oleosin) genes were identified in the genome of B. napus through a genome-wide analysis. The promoter sequences of oleosin genes consisted of various light-, hormone-, and stress-related cis-acting elements, along with transcription factor (TF) binding sites, for 25 TF families in 53 BnOLEO genes. The differentially expressed oleosin genes between two high- and two low-oil-content accessions were explored. BnOLEO3-C09, BnOLEO4-A02, BnOLEO4-A09, BnOLEO2-C04, BnOLEO1-C01, and BnOLEO7-A03 showed higher expressions in the high-oil-content accessions than in low-oil-content accessions, at 25, 35, and 45 days after pollination (DAP) in two different environments. A regional association analysis of 50 re-sequenced rapeseed accessions was used to further analyze these six BnOLEO genes, and it revealed that the nucleotide variations in the BnOLEO1-C01 and BnOLEO7-A03 gene regions were related to the phenotypic variations in seed oil content. Moreover, a co-expression network analysis revealed that the BnOLEO genes were directly linked to lipid/fatty acid metabolism, TF, lipid transport, and carbohydrate genes, thus forming a molecular network involved in seed oil accumulation. These favorable haplotypes can be utilized in molecular marker-assisted selection in order to further improve seed oil contents in rapeseed.

Comprehensive evaluation of high-oleic rapeseed (Brassica napus) based on quality, resistance, and yield traits: A new method for rapid identification of high-oleic acid rapeseed germplasm

To scientifically evaluate and utilize high-oleic acid rape germplasm resources and cultivate new varieties suitable for planting in the Hunan Province, 30 local high-oleic acid rape germplasms from Hunan were used as materials. The 12 personality indices of quality, yield, and resistance were comprehensively evaluated by variability, correlation, principal component, and cluster analyses. The results of variability showed that except for oleic acid, the lowest coefficient of variation was oil content, which was 0.06. Correlation analysis showed that oil content was positively correlated with main traits such as yield per plant and oleic acid, which could be used in the early screening of high-oleic rape germplasm. The results of principal component analysis showed that the 12 personality indicators were integrated into four principal components, and the cumulative contribution rate was 62.487%. The value of comprehensive coefficient ‘F’ was positively correlated with the first, second, and fourth principal components and negatively correlated with the third principal component. Cluster analysis showed that 30 high-oleic rape germplasms could be divided into four categories consisting of 9 (30%), 6 (20%), 7 (23%), and 8 (27%) high-oleic rape germplasms, each with the characteristics of "high disease resistance", "high yield", "high protein", and "more stability". This study not only provides a reference basis for high-oleic rape breeding but also provides a theoretical basis for their early screening.

Genome-Wide Association Analysis Combined With Quantitative Trait Loci Mapping and Dynamic Transcriptome Unveil the Genetic Control of Seed Oil Content in Brassica napus L

Rapeseed, an allotetraploid oil crop, provides vegetable oil for human consumption. The growing demand for oilseeds has necessitated the development of rapeseed varieties with improved quality. Therefore, a clear understanding of the genetic basis underlying the seed oil content (SOC) is required. In this study, a natural population comprising 204 diverse accessions and recombinant inbred lines (RILs) derived from Brassica napus and Sinapis alba via distant hybridization were collected for genome-wide association analysis (GWAS) and quantitative trait loci (QTL) mapping of the SOC trait, respectively. The variable coefficient of the RIL and natural populations ranged from 7.43 to 10.43% and 8.40 to 10.91%. Then, a high-density linkage map was constructed based on whole genome re-sequencing (WGS); the map harbored 2,799 bin markers and covered a total distance of 1,835.21 cM, with an average marker interval of 0.66 cM. The QTLs for SOC on chromosome A07 were stably detected in both single and multiple environments. Finally, a novel locus qA07.SOC was identified as the major QTL for SOC based on the GWAS and RIL populations. In addition, the RNA-seq results showed that photosynthesis, lipid biosynthesis proteins, fatty acid metabolism, and unsaturated fatty acid biosynthesis were significantly different between the developed seeds of the two parents of the RIL population. By comparing the variation information and expression levels of the syntenic genes within qA07.SOC and its syntenic genomic regions, as well as through haplotype analysis via GWAS, BnaA07.STR18, BnaA07.NRT1, and BnaA07g12880D were predicted as candidate genes in the qA07.SOC interval. These stable QTLs containing candidate genes and haplotypes can potentially provide a reliable basis for marker-assisted selection in B. napus breeding for SOC.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

Chemical Compositions of Walnut (Juglans Spp.) Oil: Combined Effects of Genetic and Climatic Factors

Walnut oil is a high-value oil product. Investigation of the variation and the main climatic factors affecting the oil’s chemical composition is vital for breeding and oil quality improvement. Therefore, the fatty acid, micronutrients, and secondary metabolites compositions and contents in walnut oil were determined in three species: Juglans regia L. (common walnut), J. sigillata Dode (iron walnut), and their hybrids (Juglans sigillata Dode × J. regia L.), which were cultivated aat different sites. The major fatty acids were linoleic (51.39–63.12%), oleic (18.40–33.56%), and linolenic acid (6.52–11.69%). High variation in the contents of fatty acids, micronutrients, and secondary metabolites was found between both species and sites. Interestingly, myristic, margaric, and margaroleic acid were only detected in the hybrids’ walnut oil, yet α-tocopherol was only detected in common and iron walnut oil. Climatic factors significantly affected the composition and content of fatty acid, whereas δ-tocopherol was mostly dependent on the genetic factors. The average relative humidity explained the most variation in the fatty acids, micronutrients, and secondary metabolites, which showed a significant positive and negative effect on the monounsaturated fatty acids and polyunsaturated fatty acids, respectively. These findings contribute to the provision of better guidance in matching sites with walnut trees, and improvement of the nutritional value of walnut oil.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

Development of B. carinata with super-high erucic acid content through interspecific hybridization

KEY MESSAGE

Disomic alien chromosome addition Brassica carinata lines with super-high erucic acid content were developed through interspecific hybridization with B. juncea and characterized using molecular, cytological and biochemical techniques. Brassica carinata [A.] Braun (BBCC, 2n = 34) is a climate-resilient oilseed. Its seed oil is high in erucic acid (> 40%), rendering it well suited for the production of biofuel and other bio-based applications. To enhance the competitiveness of B. carinata with high erucic B. napus (HEAR), lines with super-high erucic acid content were developed through interspecific hybridization. To this end, a fad2B null allele from Brassica juncea (AABB, 2n = 36) was introgressed into B. carinata, resulting in a B. carinata fad2B mutant with erucic acid levels of over 50%. Subsequently, the FAE allele from B. rapa spp. yellow sarson (AA, 2n = 20) was transferred to the fad2B B. carinata line, yielding lines with erucic acid contents of up to 57.9%. Molecular analysis using the Brassica 90 K Illumina Infinium™ SNP genotyping array identified these lines as disomic alien chromosome addition lines, with two extra A08 chromosomes containing the BrFAE gene. The alien chromosomes from B. rapa were clearly distinguished by molecular cytogenetics in one of the addition lines. Analysis of microspore-derived offspring and hybrids from crosses with a CMS B. carinata line showed that the transfer rate of the A08 chromosome into male gametes was over 98%, resulting in almost completely stable transmission of an A08 chromosome copy into the progeny. The increase in erucic acid levels was accompanied by changes in the proportions of other fatty acids depending on the genetic changes that were introduced in the interspecific hybrids, providing valuable insights into erucic acid metabolism in Brassica.

Mapping of seed quality traits in the C genome of Brassica napus by using a population carrying genome content of B. oleracea and their effect on other traits

Increasing seed oil and protein contents and reducing the content of seed glucosinolates (GSLs) in Brassica oilseed crops are important objectives in breeding. By using an oilseed rape (B. napus L.) doubled-haploid (DH) population carrying genome content introgressed from Chinese kale (B. oleracea L.), we mapped quantitative trait loci (QTL) for these seed quality traits and investigated their effect on other traits including seed yield. A stable QTL for seed oil content was identified on chromosome C5 at 40-42 Mb position and a QTL for seed GSL content was identified on C9 at 7-8 Mb position. The C5 and C9 QTL alleles for high oil and GSL contents were derived from Chinese kale, demonstrating that high-oil QTL allele can be found in the parental species of oilseed rape. The low-GSL QTL allele of C9 exerted a significant positive effect on seed protein content, demonstrating that selection for this QTL allele contributed to higher protein content in canola seed. These two QTL were not affected by field environmental conditions and did not exert a significant effect on days to flowering and seed yield. Thus, the genomic regions and the molecular markers identified in this study should be useful in molecular breeding of the seed quality traits in oilseed rape.

Regional association analysis coupled with transcriptome analyses reveal candidate genes affecting seed oil accumulation in Brassica napus

Regional association analysis of 50 re-sequenced Chinese semi-winter rapeseed accessions in combination with co-expression analysis reveal candidate genes affecting oil accumulation in Brassica napus. One of the breeding goals in rapeseed production is to enhance the seed oil content to cater to the increased demand for vegetable oils due to a growing global population. To investigate the genetic basis of variation in seed oil content, we used 60 K Brassica Infinium SNP array along with phenotype data of 203 Chinese semi-winter rapeseed accessions to perform a genome-wide analysis of haplotype blocks associated with the oil content. Nine haplotype regions harbouring lipid synthesis/transport-, carbohydrate metabolism- and photosynthesis-related genes were identified as significantly associated with the oil content and were mapped to chromosomes A02, A04, A05, A07, C03, C04, C05, C08 and C09, respectively. Regional association analysis of 50 re-sequenced Chinese semi-winter rapeseed accessions combined with transcriptome datasets from 13 accessions was further performed on these nine haplotype regions. This revealed natural variation in the BnTGD3-A02 and BnSSE1-A05 gene regions correlated with the phenotypic variation of the oil content within the A02 and A04 chromosome haplotype regions, respectively. Moreover, co-expression network analysis revealed that BnTGD3-A02 and BnSSE1-A05 were directly linked with fatty acid beta-oxidation-related gene BnKAT2-C04, thus forming a molecular network involved in the potential regulation of seed oil accumulation. The results of this study could be used to combine favourable haplotype alleles for further improvement of the seed oil content in rapeseed.

Genome- and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus

Seed oil content (SOC) is a highly important and complex trait in oil crops. Here, we decipher the genetic basis of natural variation in SOC of Brassica napus by genome- and transcriptome-wide association studies using 505 inbred lines. We mapped reliable quantitative trait loci (QTLs) that control SOC in eight environments, evaluated the effect of each QTL on SOC, and analyzed selection in QTL regions during breeding. Six-hundred and ninety-two genes and four gene modules significantly associated with SOC were identified by analyzing population transcriptomes from seeds. A gene prioritization framework, POCKET (prioritizing the candidate genes by incorporating information on knowledge-based gene sets, effects of variants, genome-wide association studies, and transcriptome-wide association studies), was implemented to determine the causal genes in the QTL regions based on multi-omic datasets. A pair of homologous genes, BnPMT6s, in two QTLs were identified and experimentally demonstrated to negatively regulate SOC. This study provides rich genetic resources for improving SOC and valuable insights toward understanding the complex machinery that directs oil accumulation in the seeds of B. napus and other oil crops.

Genomic Prediction and Genetic Correlation of Agronomic, Blackleg Disease, and Seed Quality Traits in Canola (Brassica napus L.)

Genomic selection accelerates genetic progress in crop breeding through the prediction of future phenotypes of selection candidates based on only their genomic information. Here we report genetic correlations and genomic prediction accuracies in 22 agronomic, disease, and seed quality traits measured across multiple years (2015-2017) in replicated trials under rain-fed and irrigated conditions in Victoria, Australia. Two hundred and two spring canola lines were genotyped for 62,082 Single Nucleotide Polymorphisms (SNPs) using transcriptomic genotype-by-sequencing (GBSt). Traits were evaluated in single trait and bivariate genomic best linear unbiased prediction (GBLUP) models and cross-validation. GBLUP were also expanded to include genotype-by-environment G × E interactions. Genomic heritability varied from 0.31to 0.66. Genetic correlations were highly positive within traits across locations and years. Oil content was positively correlated with most agronomic traits. Strong, not previously documented, negative correlations were observed between average internal infection (a measure of blackleg disease) and arachidic and stearic acids. The genetic correlations between fatty acid traits followed the expected patterns based on oil biosynthesis pathways. Genomic prediction accuracy ranged from 0.29 for emergence count to 0.69 for seed yield. The incorporation of G × E translates into improved prediction accuracy by up to 6%. The genomic prediction accuracies achieved indicate that genomic selection is ready for application in canola breeding.

Marker assisted selection of new high oleic and low linolenic winter oilseed rape (Brassica napus L.) inbred lines revealing good agricultural value

Development of oilseed rape (Brassica napus L.) breeding lines producing oil characterized by high oleic and low linolenic acid content is an important goal of rapeseed breeding programs worldwide. Such kind of oil is ideal for deep frying and can also be used as a raw material for biodiesel production. By performing chemical mutagenesis using ethyl methanesulfonate, we obtained mutant winter rapeseed breeding lines that can produce oil with a high content of oleic acid (C18:1, more than 75%) and a low content of linolenic acid (C18:3, less than 3%). However, the mutant lines revealed low agricultural value as they were characterized by low seed yield, low wintering, and high content of glucosinolates in seed meal. The aim of this work was to improve the mutant lines and develop high-oleic and low-linolenic recombinants exhibiting both good oil quality and high agronomic value. The plant materials used in this study included high-oleic and low-linolenic mutant breeding lines and high-yielding domestic canola-type breeding lines of good agricultural value with high oleic acid content and extremely low glucosinolates content. Field trials were conducted in four environments, in a randomized complete block design. Phenotyping was performed for wintering, yield of seed and oil, and seed quality traits. Genotype × environment interaction was investigated with respect to the content of C18:1 and C18:3 acids in seed oil. Genotyping was done for the selection of homozygous high oleic and low linolenic lines using allele-specific CAPS markers and SNaPshot assay, respectively. Finally, new high oleic and low linolenic winter rapeseed recombinant lines were obtained for use as a starting material for the development of new varieties that may be of high value on the oil crop market.

Identification and QTL mapping of resistance to Turnip yellows virus (TuYV) in oilseed rape, Brassica napus

Partially dominant resistance to Turnip yellows virus associated with one major QTL was identified in the natural allotetraploid oilseed rape cultivar Yudal. Turnip yellows virus (TuYV) is transmitted by the peach-potato aphid (Myzus persicae) and causes severe yield losses in commercial oilseed rape crops (Brassica napus). There is currently only one genetic resource for resistance to TuYV available in brassica, which was identified in the re-synthesised B. napus line 'R54'. In our study, 27 mostly homozygous B. napus accessions, either doubled-haploid (DH) or inbred lines, representing a diverse subset of the B. napus genepool, were screened for TuYV resistance/susceptibility. Partial resistance to TuYV was identified in the Korean spring oilseed rape, B. napus variety Yudal, whilst the dwarf French winter oilseed rape line Darmor-bzh was susceptible. QTL mapping using the established Darmor-bzh × Yudal DH mapping population (DYDH) revealed one major QTL explaining 36% and 18% of the phenotypic variation in two independent experiments. A DYDH line was crossed to Yudal, and reciprocal backcross (BC₁) populations from the F₁ with either the susceptible or resistant parent revealed the dominant inheritance of the TuYV resistance. The QTL on ChrA04 was verified in the segregating BC₁ population. A second minor QTL on ChrC05 was identified in one of the two DYDH experiments, and it was not observed in the BC₁ population. The TuYV resistance QTL in 'R54' is within the QTL interval on Chr A04 of Yudal; however, the markers co-segregating with the 'R54' resistance are not conserved in Yudal, suggesting an independent origin of the TuYV resistances. This is the first report of the QTL mapping of TuYV resistance in natural B. napus.

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.

Resolution of quantitative resistance to clubroot into QTL-specific metabolic modules

Plant disease resistance is often under quantitative genetic control. Thus, in a given interaction, plant cellular responses to infection are influenced by resistance or susceptibility alleles at different loci. In this study, a genetic linkage analysis was used to address the complexity of the metabolic responses of Brassica napus roots to infection by Plasmodiophora brassicae. Metabolome profiling and pathogen quantification in a segregating progeny allowed a comparative mapping of quantitative trait loci (QTLs) involved in resistance and in metabolic adjustments. Distinct metabolic modules were associated with each resistance QTL, suggesting the involvement of different underlying cellular mechanisms. This approach highlighted the possible role of gluconasturtiin and two unknown metabolites in the resistance conferred by two QTLs on chromosomes C03 and C09, respectively. Only two susceptibility biomarkers (glycine and glutathione) were simultaneously linked to the three main resistance QTLs, suggesting the central role of these compounds in the interaction. By contrast, several genotype-specific metabolic responses to infection were genetically unconnected to resistance or susceptibility. Likewise, variations of root sugar profiles, which might have influenced pathogen nutrition, were not found to be related to resistance QTLs. This work illustrates how genetic metabolomics can help to understand plant stress responses and their possible links with disease.

Identification of candidate genes controlling oil content by combination of genome-wide association and transcriptome analysis in the oilseed crop Brassica napus

BACKGROUND

Increasing seed oil content is one of the most important targets for rapeseed (Brassica napus) breeding. However, genetic mechanisms of mature seed oil content in Brassica napus (B. napus) remain little known. To identify oil content-related genes, a genome-wide association study (GWAS) was performed using 588 accessions.

RESULTS

High-throughput genome resequencing resulted in 385,692 high-quality single nucleotide polymorphism (SNPs) with a minor allele frequency (MAF) > 0.05. We identified 17 loci that were significantly associated with seed oil content, among which 12 SNPs were distributed on the A3 (11 loci) and A1 (one loci) chromosomes, and five novel significant SNPs on the C5 (one loci) and C7 (four loci) chromosomes, respectively. Subsequently, we characterized differentially expressed genes (DEGs) between the seeds and silique pericarps on main florescences and primary branches of extremely high- and low-oil content accessions (HO and LO). A total of 64 lipid metabolism-related DEGs were identified, 14 of which are involved in triacylglycerols (TAGs) biosynthesis and assembly. Additionally, we analyzed differences in transcription levels of key genes involved in de novo fatty acid biosynthesis in the plastid, TAGs assembly and lipid droplet packaging in the endoplasmic reticulum (ER) between high- and low-oil content B. napus accessions.

CONCLUSIONS

The combination of GWAS and transcriptome analyses revealed seven candidate genes located within the confidence intervals of significant SNPs. Current findings provide valuable information for facilitating marker-based breeding for higher seed oil content in B. napus.

Stable and novel QTL identification and new insights into the genetic networks affecting seed fiber traits in Brassica napus

QTL mapping for fiber-related traits and elucidation of a stable and novel QTL affecting seed lignin content, cellulose content and seed oil content. Dissection of the genetic networks for fiber biosynthesis is important for improving the seed oil content and meal value of Brassica napus. In this study, the genetic basis of seed fiber biosynthesis in B. napus was investigated via quantitative trait locus (QTL) analysis using a doubled haploid population derived from 'KenC-8' crossed with 'N53-2.' Seed lignin content (LC), cellulose content (CC) and hemicellulose content (HC) were significantly negatively correlated with seed oil content (OC). Co-localization QTLs among LC, CC, HC and OC on A09 were found with contributions ranging from 9.87 to 48.50%. Seven co-localization QTLs involved in the fiber component and OC were further verified by bulked segregant analysis (BSA). The unique QTL uqA9-12 might be a real and new QTL that was commonly identified by QTL mapping and BSA and simultaneously affected LC, CC and OC with opposite additive effects. A potential regulatory network controlling seed fiber biosynthesis was constructed to dissect the complex mechanism of seed fiber and oil accumulation, and numerous candidate genes were identified in the fiber-related QTL regions. These results provided an enrichment of QTLs and potential candidates for fiber biosynthesis, as well as useful new information for understanding the complex genetic mechanism underlying rapeseed seed fiber accumulation.

Stable Quantitative Resistance Loci to Blackleg Disease in Canola (Brassica napus L.) Over Continents

The hemibiotrophic fungus, Leptosphaeria maculans is the most devastating pathogen, causing blackleg disease in canola (Brassica napus L). To study the genomic regions involved in quantitative resistance (QR), 259-276 DH lines from Darmor-bzh/Yudal (DYDH) population were assessed for resistance to blackleg under shade house and field conditions across 3 years. In different experiments, the broad sense heritability varied from 43 to 95%. A total of 27 significant quantitative trait loci (QTL) for QR were detected on 12 chromosomes and explained between 2.14 and 10.13% of the genotypic variance. Of the significant QTL, at least seven were repeatedly detected across different experiments on chromosomes A02, A07, A09, A10, C01, and C09. Resistance alleles were mainly contributed by 'Darmor-bzh' but 'Yudal' also contributed few of them. Our results suggest that plant maturity and plant height may have a pleiotropic effect on QR in our conditions. We confirmed that Rlm9 which is present in 'Darmor-bzh' is not effective to confer resistance in our Australian field conditions. Comparative mapping showed that several R genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors map in close proximity (within 200 Kb) of the significant trait-marker associations on the reference 'Darmor-bzh' genome assembly. More importantly, eight significant QTL regions were detected across diverse growing environments: Australia, France, and United Kingdom. These stable QTL identified herein can be utilized for enhancing QR in elite canola germplasm via marker- assisted or genomic selection strategies.

Drought stress has transgenerational effects on seeds and seedlings in winter oilseed rape (Brassica napus L.)

BACKGROUND

Drought stress has a negative effect on both seed yield and seed quality in Brassica napus (oilseed rape, canola). Here we show that while drought impairs the maternal plant performance, it also increases the vigour of progeny of stressed maternal plants. We investigated the transgenerational influence of abiotic stress by detailed analysis of yield, seed quality, and seedling performance on a growth-related and metabolic level. Seeds of eight diverse winter oilseed rape genotypes were generated under well-watered and drought stress conditions under controlled-environment conditions in large plant containers.

RESULTS

We found a decrease in seed quality in seeds derived from mother plants that were exposed to drought stress. At the same time, the seeds that developed under stress conditions showed higher seedling vigour compared to non-stressed controls.This effect on seed quality and seedling vigour was found to be independent of maternal plant yield performance.

CONCLUSIONS

Drought stress has a positive transgenerational effect on seedling vigour. Three potential causes for stress-induced improvement of seedling vigour are discussed: (1) Heterotic effects caused by a tendency towards a higher outcrossing rate in response to stress; (2) an altered reservoir of seed storage metabolites to which the seedling resorts during early growth, and (3) inter-generational stress memory, formed by stress-induced changes in the epigenome of the seedling.

A major QTL on chromosome C05 significantly reduces acid detergent lignin (ADL) content and increases seed oil and protein content in oilseed rape (Brassica napus L.)

A reduction in acid detergent lignin content in oilseed rape resulted in an increase in seed oil and protein content. Worldwide increasing demand for vegetable oil and protein requires continuous breeding efforts to enhance the yield of oil and protein crop species. The oil-extracted meal of oilseed rape is currently mainly used for feeding livestock, but efforts are undertaken to use the oilseed rape protein in food production. One limiting factor is the high lignin content of black-seeded oilseed rape that negatively affects digestibility and sensory quality of food products compared to soybean. Breeding attempts to develop yellow seeded oilseed rape with reduced lignin content have not yet resulted in competitive cultivars. The objective of this work was to investigate the inheritance of seed quality in a DH population derived from the cross of the high oil lines SGDH14 and cv. Express. The DH population of 139 lines was tested in field experiments in 14 environments in north-west Europe. Seeds harvested from open pollinated plants were used for extensive seed quality analysis. A molecular marker map based on the Illumina Infinium 60 K Brassica SNP chip was used to map QTL. Amongst others, one major QTL for acid detergent lignin content, explaining 81% of the phenotypic variance, was identified on chromosome C05. Lines with reduced lignin content nevertheless did not show a yellowish appearance, but showed a reduced seed hull content. The position of the QTL co-located with QTL for oil and protein content of the defatted meal with opposite additive effects, suggesting that the reduction in lignin content resulted in an increase in oil and protein content.

Dissection of the genetic architecture of three seed-quality traits and consequences for breeding in Brassica napus

Genome-wide association studies (GWASs) combining high-throughput genome resequencing and phenotyping can accelerate the dissection of genetic architecture and identification of genes for plant complex traits. In this study, we developed a rapeseed genomic variation map consisting of 4 542 011 SNPs and 628 666 INDELs. GWAS was performed for three seed-quality traits, including erucic acid content (EAC), glucosinolate content (GSC) and seed oil content (SOC) using 3.82 million polymorphisms in an association panel. Six, 49 and 17 loci were detected to be associated with EAC, GSC and SOC in multiple environments, respectively. The mean total contribution of these loci in each environment was 94.1% for EAC and 87.9% for GSC, notably higher than that for SOC (40.1%). A high correlation was observed between phenotypic variance and number of favourable alleles for associated loci, which will contribute to breeding improvement by pyramiding these loci. Furthermore, candidate genes were detected underlying associated loci, based on functional polymorphisms in gene regions where sequence variation was found to correlate with phenotypic variation. Our approach was validated by detection of well-characterized FAE1 genes at each of two major loci for EAC on chromosomes A8 and C3, along with MYB28 genes at each of three major loci for GSC on chromosomes A9, C2 and C9. Four novel candidate genes were detected by correlation between GSC and SOC and observed sequence variation, respectively. This study provides insights into the genetic architecture of three seed-quality traits, which would be useful for genetic improvement of B. napus.

Effect of high night temperature on storage lipids and transcriptome changes in developing seeds of oilseed rape

Global warming causes a faster increase of night temperature than of day temperature in tropical and subtropical zones. Little is known about the effect of high night temperature on storage lipids and transcriptome changes in oilseed rape. This study compared the total fatty acids and fatty acid compositions in seeds of two oilseed rape cultivars between high and low night temperatures. Their transcriptome profiles were also analyzed. High night temperature significantly affected the total fatty acids and fatty acid compositions in seeds of both low and high oil content cultivars, namely Jiuer-13 and Zheyou-50, thereby resulting in 18.9% and 13.7% total fatty acid reductions, respectively. In particular, high night temperature decreased the relative proportions of C18:0 and C18:1 but increased the proportions of C18:2 and C18:3 in both cultivars. In-depth analysis of transcriptome profiles revealed that high night temperature up-regulated gibberellin signaling during the night-time. This up-regulation was associated with the active expression of genes involved in fatty acid catabolism, such as those in β-oxidation and glyoxylate metabolism pathways. Although the effect of temperature on plant lipids has been previously examined, the present study is the first to focus on night temperature and its effect on the fatty acid composition in seeds.

QTL Landscape for Oil Content in Brassica juncea: Analysis in Multiple Bi-Parental Populations in High and "0" Erucic Background

Increasing oil content in oilseed mustard (Brassica juncea) is a major breeding objective-more so, in the lines that have "0" erucic acid content (< 2% of the seed oil) as earlier studies have shown negative pleiotropic effect of erucic acid loci on the oil content, both in oilseed mustard and rapeseed. We report here QTL analysis of oil content in eight different mapping populations involving seven different parents-including a high oil content line J8 (~49%). The parental lines of the mapping populations contained wide variation in oil content and erucic acid content. The eight mapping populations were categorized into two sets-five populations with individuals segregating for erucic acid (SE populations) and the remaining three with zero erucic acid segregants (ZE populations). Meta-analysis of QTL mapped in individual SE populations identified nine significant C-QTL, with two of these merging most of the major oil QTL that colocalized with the erucic acid loci on the linkage groups A08 and B07. QTL analysis of oil content in ZE populations revealed a change in the landscape of the oil QTL compared to the SE populations, in terms of altered allelic effects and phenotypic variance explained by ZE QTL at the "common" QTL and observation of "novel" QTL in the ZE background. The important loci contributing to oil content variation, identified in the present study could be used in the breeding programmes for increasing the oil content in high erucic and "0" erucic backgrounds.

Genomic and Transcriptomic Analysis Identified Gene Clusters and Candidate Genes for Oil Content in Peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea), a major source of vegetable oil in many Asian countries, has become an integral part of human diet globally due to its high nutritional properties and option to consume in different forms. In order to meet the demand of vegetable oil, many peanut breeding programs of China have intensified their efforts in increasing oil content in newly bred varieties for reducing the import of edible oils in China. In this context, transcriptome sequencing data generated on 49 peanut cultivars were analyzed to identify candidate genes and develop molecular markers for seed oil content across multiple environments. Transcriptome analysis identified 5458 differentially expressed genes (DEGs) including 2243 positive DEGs and 3215 negative DEGs involved in oil synthesis process. Genome-wide association study identified 48 significant insertion/deletion (InDel) markers associated with seed oil content across five environments. A comparative genomics and transcriptomics analysis detected a total of 147 common gene clusters located in 17 chromosomes. Interestingly, an InDel cluster associated with seed oil content on A03 chromosome was detected in three different environments. Candidate genes identified on A03 form a haplotype, in which variable alleles were found to be different in oil content in an independent population. This locus is important for understanding the genetic control of peanut oil content and may be useful for marker-assisted selection in peanut breeding programs.

Synteny analysis of genes and distribution of loci controlling oil content and fatty acid profile based on QTL alignment map in Brassica napus

Background

Deciphering the genetic architecture of a species is a good way to understand its evolutionary history, but also to tailor its profile for breeding elite cultivars with desirable traits. Aligning QTLs from diverse population in one map and utilizing it for comparison, but also as a basis for multiple analyses assure a stronger evidence to understand the genetic system related to a given phenotype.

Results

In this study, 439 genes involved in fatty acid (FA) and triacylglycerol (TAG) biosyntheses were identified in Brassica napus. B. napus genome showed mixed gene loss and insertion compared to B. rapa and B. oleracea, and C genome had more inserted genes. Identified QTLs for oil (OC-QTLs) and fatty acids (FA-QTLs) from nine reported populations were projected on the physical map of the reference genome “Darmor-bzh” to generate a map. Thus, 335 FA-QTLs and OC-QTLs could be highlighted and 82 QTLs were overlapping. Chromosome C3 contained 22 overlapping QTLs with all trait studied except for C18:3. In total, 218 candidate genes which were potentially involved in FA and TAG were identified in 162 QTLs confidence intervals and some of them might affect many traits. Also, 76 among these candidate genes were found inside 57 overlapping QTLs, and candidate genes for oil content were in majority (61/76 genes). Then, sixteen genes were found in overlapping QTLs involving three populations, and the remaining 60 genes were found in overlapping QTLs of two populations. Interaction network and pathway analysis of these candidate genes indicated ten genes that might have strong influence over the other genes that control fatty acids and oil formation.

Conclusion

The present results provided new information for genetic basis of FA and TAG formation in B. napus. A map including QTLs from numerous populations was built, which could serve as reference to study the genome profile of B. napus, and new potential genes emerged which might affect seed oil. New useful tracks were showed for the selection of population or/and selection of interesting genes for breeding improvement purpose.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4176-6) contains supplementary material, which is available to authorized users.

Gene Introgression in Weeds Depends on Initial Gene Location in the Crop: Brassica napus-Raphanus raphanistrum Model

The effect of gene location within a crop genome on its transfer to a weed genome remains an open question for gene flow assessment. To elucidate this question, we analyzed advanced generations of intergeneric hybrids, derived from an initial pollination of known oilseed rape varieties (Brassica napus, AACC, 2n = 38) by a local population of wild radish (Raphanus raphanistrum, RrRr, 2n = 18). After five generations of recurrent pollination, 307 G5 plants with a chromosome number similar to wild radish were genotyped using 105 B. napus specific markers well distributed along the chromosomes. They revealed that 49.8% of G5 plants carried at least one B. napus genomic region. According to the frequency of B. napus markers (0-28%), four classes were defined: Class 1 (near zero frequency), with 75 markers covering ∼70% of oilseed rape genome; Class 2 (low frequency), with 20 markers located on 11 genomic regions; Class 3 (high frequency), with eight markers on three genomic regions; and Class 4 (higher frequency), with two adjacent markers detected on A10. Therefore, some regions of the oilseed rape genome are more prone than others to be introgressed into wild radish. Inheritance and growth of plant progeny revealed that genomic regions of oilseed rape could be stably introduced into wild radish and variably impact the plant fitness (plant height and seed number). Our results pinpoint that novel technologies enabling the targeted insertion of transgenes should select genomic regions that are less likely to be introgressed into the weed genome, thereby reducing gene flow.

Genetic dissection of seed oil and protein content and identification of networks associated with oil content in Brassica napus

High-density linkage maps can improve the precision of QTL localization. A high-density SNP-based linkage map containing 3207 markers covering 3072.7 cM of the Brassica napus genome was constructed in the KenC-8 × N53-2 (KNDH) population. A total of 67 and 38 QTLs for seed oil and protein content were identified with an average confidence interval of 5.26 and 4.38 cM, which could explain up to 22.24% and 27.48% of the phenotypic variation, respectively. Thirty-eight associated genomic regions from BSA overlapped with and/or narrowed the SOC-QTLs, further confirming the QTL mapping results based on the high-density linkage map. Potential candidates related to acyl-lipid and seed storage underlying SOC and SPC, respectively, were identified and analyzed, among which six were checked and showed expression differences between the two parents during different embryonic developmental periods. A large primary carbohydrate pathway based on potential candidates underlying SOC- and SPC-QTLs, and interaction networks based on potential candidates underlying SOC-QTLs, was constructed to dissect the complex mechanism based on metabolic and gene regulatory features, respectively. Accurate QTL mapping and potential candidates identified based on high-density linkage map and BSA analyses provide new insights into the complex genetic mechanism of oil and protein accumulation in the seeds of rapeseed.

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

BACKGROUND

B. napus (oilseed) is an important source of edible vegetable oil, and its nutritional and economic value is determined by its fatty acid composition and content.

RESULTS

Using the Brassica 60 K SNP array, we performed a genome-wide association study of fatty acid composition in a population of 520 genetically diverse oilseed accessions. Using the PCA + K model in TASSEL 5.2.1, we identified 62 genomic regions that were significantly associated with the composition of seven fatty acids, and five consensus regions that mapped to the A2, A8, A9, C1, and C3 chromosomes, respectively, of the Brassica napus Darmor-bzh genome. We then identified 24 orthologs of the functional candidate genes involved in fatty acid biosynthesis, excluding BnaA.FAE1 and BnaC.FAE1 on the A8 and C3 homologous genome blocks, which are known to have critical roles in the fatty acid biosynthesis pathway, and potential orthologs of these genes (e.g., LACS9, KCR1, FAB1, LPAT4, KCS17, CER4, TT16, and ACBP5).

CONCLUSIONS

Our results demonstrate the power of association mapping in identifying genes of interest in B. napus and provide insight into the genetic basis of fatty acid biosynthesis in B. napus. Furthermore, our findings may facilitate marker-based breeding efforts aimed at improving fatty acid composition and quality in B. napus.

A Major Locus for Manganese Tolerance Maps on Chromosome A09 in a Doubled Haploid Population of Brassica napus L

Soil acidity poses a major threat to productivity of several crops; mainly due to the prevalence of toxic levels of Al³⁺ and Mn²⁺. Crop productivity could be harnessed on acid soils via the development of plant varieties tolerant to phytotoxic levels of these cations. In this study, we investigated the extent of natural variation for Mn²⁺ tolerance among ten parental lines of the Australian and International canola mapping populations. Response to Mn²⁺ toxicity was measured on the bases of cotyledon chlorosis, shoot biomass, and leaf area in nutrient solution under control (9 μM of MnCl₂⋅4H₂O) and Mn treatment (125 μM of MnCl₂⋅4H₂O). Among parental lines, we selected Darmor-bzh and Yudal that showed significant and contrasting variation in Mn²⁺ tolerance to understand genetic control and identify the quantitative trait loci (QTL) underlying Mn²⁺ tolerance. We evaluated parental lines and their doubled haploid (DH) progenies (196 lines) derived from an F₁ cross, Darmor-bzh/Yudal for Mn²⁺ tolerance. Mn²⁺-tolerant genotypes had significantly higher shoot biomass and leaf area compared to Mn²⁺-sensitive genotypes. A genetic linkage map based on 7,805 DArTseq markers corresponding to 2,094 unique loci was constructed and further utilized for QTL identification. A major locus, BnMn²⁺.A09 was further mapped with a SNP marker, Bn-A09-p29012402 (LOD score of 34.6) accounting for most of the variation in Mn²⁺ tolerance on chromosome A09. This is the first report on the genomic localization of a Mn²⁺ tolerance locus in B. napus. Additionally, an ortholog of A. thaliana encoding for cation efflux facilitator transporter was located within 3,991 bp from significant SNP marker associated with BnMn²⁺.A09. A suite of genome sequence based markers (DArTseq and Illumina Infinium SNPs) flanking the BnMn²⁺.A09 locus would provide an invaluable tool for various molecular breeding applications to improve canola production and profitability on Mn²⁺ toxic soils.

Seed Quality Traits Can Be Predicted with High Accuracy in Brassica napus Using Genomic Data

Improving seed oil yield and quality are central targets in rapeseed (Brassica napus) breeding. The primary goal of our study was to examine and compare the potential and the limits of marker-assisted selection and genome-wide prediction of six important seed quality traits of B. napus. Our study is based on a bi-parental population comprising 202 doubled haploid lines and a diverse validation set including 117 B. napus inbred lines derived from interspecific crosses between B. rapa and B. carinata. We used phenotypic data for seed oil, protein, erucic acid, linolenic acid, stearic acid, and glucosinolate content. All lines were genotyped with a 60k SNP array. We performed five-fold cross-validations in combination with linkage mapping and four genome-wide prediction approaches in the bi-parental population. Quantitative trait loci (QTL) with large effects were detected for erucic acid, stearic acid, and glucosinolate content, blazing the trail for marker-assisted selection. Despite substantial differences in the complexity of the genetic architecture of the six traits, genome-wide prediction models had only minor impacts on the prediction accuracies. We evaluated the effects of training population size, marker density and phenotyping intensity on the prediction accuracy. The prediction accuracy in the independent and genetically very distinct validation set still amounted to 0.14 for protein content and 0.17 for oil content reflecting the utility of the developed calibration models even in very diverse backgrounds.

Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus

Seed oil content is an important agricultural trait in rapeseed breeding. Although numerous quantitative trait locus (QTL) have been identified, most of them cannot be applied in practical breeding mainly due to environmental instability or large confidence intervals. The purpose of this study was to identify and validate high quality and more stable QTLs by combining linkage mapping and genome-wide association study (GWAS). For linkage mapping, we constructed two F₂ populations from crosses of high-oil content (∼50%) lines 6F313 and 61616 with a low-oil content (∼40%) line 51070. Two high density linkage maps spanned 1987cM (1659 bins) and 1856cM (1746 bins), respectively. For GWAS, we developed more than 34,000 high-quality SNP markers based on 227 accessions. Finally, 40 QTLs and 29 associations were established by linkage and association mapping in different environments. After merging the results, 32 consensus QTLs were obtained and 7 of them were identified by both mapping methods. Seven overlapping QTLs covered an average confidence interval of 183kb and explained the phenotypic variation of 10.23 to 24.45%. We further developed allele-specific PCR primers to identify each of the seven QTLs. These stable QTLs should be useful in gene cloning and practical breeding application.

Genetic basis of nitrogen use efficiency and yield stability across environments in winter rapeseed

Background

Nitrogen use efficiency is an important breeding trait that can be modified to improve the sustainability of many crop species used in agriculture. Rapeseed is a major oil crop with low nitrogen use efficiency, making its production highly dependent on nitrogen input. This complex trait is suspected to be sensitive to genotype × environment interactions, especially genotype × nitrogen interactions. Therefore, phenotyping diverse rapeseed populations under a dense network of trials is a powerful approach to study nitrogen use efficiency in this crop. The present study aimed to determine the quantitative trait loci (QTL) associated with yield in winter oilseed rape and to assess the stability of these regions under contrasting nitrogen conditions for the purpose of increasing nitrogen use efficiency.

Results

Genome-wide association studies and linkage analyses were performed on two diversity sets and two doubled-haploid populations. These populations were densely genotyped, and yield-related traits were scored in a multi-environment design including seven French locations, six growing seasons (2009 to 2014) and two nitrogen nutrition levels (optimal versus limited). Very few genotype × nitrogen interactions were detected, and a large proportion of the QTL were stable across nitrogen nutrition conditions. In contrast, strong genotype × trial interactions in which most of the QTL were specific to a single trial were found. To obtain further insight into the QTL × environment interactions, genetic analyses of ecovalence were performed to identify the genomic regions contributing to the genotype × nitrogen and genotype × trial interactions. Fifty-one critical genomic regions contributing to the additive genetic control of yield-associated traits were identified, and the structural organization of these regions in the genome was investigated.

Conclusions

Our results demonstrated that the effect of the trial was greater than the effect of nitrogen nutrition levels on seed yield-related traits under our experimental conditions. Nevertheless, critical genomic regions associated with yield that were stable across environments were identified in rapeseed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0432-z) contains supplementary material, which is available to authorized users.

A genome-wide association study reveals novel elite allelic variations in seed oil content of Brassica napus

A set of additive loci for seed oil content were identified using association mapping and one of the novel loci on the chromosome A5 was validated by linkage mapping. Increasing seed oil content is one of the most important goals in the breeding of oilseed crops including Brassica napus, yet the genetic basis for variations in this important trait remains unclear. By genome-wide association study of seed oil content using 521 B. napus accessions genotyped with the Brassica 60K SNP array, we identified 50 loci significantly associated with seed oil content using three statistical models, the general linear model, the mixed linear model and the Anderson-Darling test. Together, the identified loci could explain approximately 80 % of the total phenotypic variance, and 29 of these loci have not been reported previously. Furthermore, a novel locus on the chromosome A5 that could increase 1.5-1.7 % of seed oil content was validated in an independent bi-parental linkage population. Haplotype analysis showed that the favorable alleles for seed oil content exhibit cumulative effects. Our results thus provide valuable information for understanding the genetic control of seed oil content in B. napus and may facilitate marker-based breeding for a higher seed oil content in this important oil crop.

Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus

BACKGROUND

Yield of rapeseed is determined by three components: silique number, seed number per silique and thousand seed weight. Seed number per silique and thousand seed weight are influenced by silique length, seed density, silique breadth, silique thickness and silique volume. Some QTLs for silique traits have been reported in B. napus, however, no studies have focused on the six agronomic traits (seed number per silique, silique length, silique breadth, silique thickness, seed density and silique volume) simultaneously, and the genetic determinism of such complex traits have not been fully elucidated.

RESULTS

In this study, the six silique traits were evaluated using 348 lines of a doubled haploid population, the KN population. The results showed that 2, 4, 1, 1 and 2 QTLs explaining > 10 % of phenotypic variation were obtained for silique length, silique breadth, silique thickness, seed number per silique and silique volume, respectively. Notably, three major effect QTLs (cqSB-C6-1, cqSB-C6-2 and cqSV-C6-3) were identified in at least three environments, and 17 unique QTLs controlling at least two traits were obtained. A high-density consensus map containing 1225 markers was constructed for QTL comparison by combining the KN map with other five published maps. The comparative results revealed that 14, 13 and 11 QTLs for silique breadth, silique thickness and silique volume might be the potential new QTLs because few QTLs for these traits were reported in B. napus. In addition, potential new QTLs for silique length (11), seed number per silique (6) and seed density (5) were also identified. Twenty-five candidate genes underlying 27 QTLs for silique related traits were obtained.

CONCLUSIONS

This study constructed QTL analysis in B. napus, and obtained 60 consensus QTLs for six silique related traits. The potential new QTLs will enhance our understanding of the genetic control of silique traits, and the stable QTLs provided the targets for improving seed yield in future. These findings provided comprehensive insights into the genetic network affecting silique traits at QTL level in B. napus.

Genomic Prediction of Testcross Performance in Canola (Brassica napus)

Genomic selection (GS) is a modern breeding approach where genome-wide single-nucleotide polymorphism (SNP) marker profiles are simultaneously used to estimate performance of untested genotypes. In this study, the potential of genomic selection methods to predict testcross performance for hybrid canola breeding was applied for various agronomic traits based on genome-wide marker profiles. A total of 475 genetically diverse spring-type canola pollinator lines were genotyped at 24,403 single-copy, genome-wide SNP loci. In parallel, the 950 F1 testcross combinations between the pollinators and two representative testers were evaluated for a number of important agronomic traits including seedling emergence, days to flowering, lodging, oil yield and seed yield along with essential seed quality characters including seed oil content and seed glucosinolate content. A ridge-regression best linear unbiased prediction (RR-BLUP) model was applied in combination with 500 cross-validations for each trait to predict testcross performance, both across the whole population as well as within individual subpopulations or clusters, based solely on SNP profiles. Subpopulations were determined using multidimensional scaling and K-means clustering. Genomic prediction accuracy across the whole population was highest for seed oil content (0.81) followed by oil yield (0.75) and lowest for seedling emergence (0.29). For seed yieId, seed glucosinolate, lodging resistance and days to onset of flowering (DTF), prediction accuracies were 0.45, 0.61, 0.39 and 0.56, respectively. Prediction accuracies could be increased for some traits by treating subpopulations separately; a strategy which only led to moderate improvements for some traits with low heritability, like seedling emergence. No useful or consistent increase in accuracy was obtained by inclusion of a population substructure covariate in the model. Testcross performance prediction using genome-wide SNP markers shows considerable potential for pre-selection of promising hybrid combinations prior to resource-intensive field testing over multiple locations and years.

Allelic Variation of BnaC.TT2.a and Its Association with Seed Coat Color and Fatty Acids in Rapeseed (Brassica napus L.)

Efficient molecular markers for the selection of rapeseed genetic materials with high seed oil content and ideal fatty acid (FA) composition are preferred by rapeseed breeders. Recently, we reported the molecular mechanism of TRANSPARENT TESTA 2 (TT2) in inhibiting seed FA biosynthesis in Arabidopsis. However, evidence showing the association of rapeseed TT2 homologs and seed FA production are still insufficient. In this study, we collected 83 rapeseed (Brassica napus L.) landraces from different geographical backgrounds to conduct association mapping of BnaC.TT2.a in relation to seed coat color and FA biosynthesis. Population background was corrected by 84 pairs of SSR markers that were uniformly distributed among the linkage groups of the Tapidor-Ningyou-7 DH population. A single copy of BnaC.TT2.a for single nucleotide polymorphism (SNP) assay was cloned by a pair of previously reported specific primers. From the analysis of BnaC.TT2.a allelic variations using GLM+Q model, four SNPs on intron 1 of BnaC.TT2.a that were associated with seed FA were discovered. Moreover, an InDel at position 738 on exon 3 of BnaC.TT2.a indicated a change of protein function that was significantly associated with seed coat color, linoleic acid (C18:2), and total FA content. These findings revealed the role of BnaC.TT2.a in regulating the seed color formation and seed FA biosynthesis in rapeseed, thereby suggesting effective molecular markers for rapeseed breeding.

Genetic variation and inheritance of phytosterol and oil content in a doubled haploid population derived from the winter oilseed rape Sansibar × Oase cross

Identification of QTL for phytosterol content, oil content, fatty acids content, protein content of defatted meal, and seed weight by multiple interval mapping in a Brassica napus DH population. Phytosterols are minor seed constituents in oilseed rape which have recently drawn wide-interest from the food and nutrition industry due to their health benefit in lowering LDL cholesterol in humans. To understand the genetic basis of phytosterol content and its relationship with other seed quality traits in oilseed rape, QTL mapping was performed in a segregating DH population derived from the cross of two winter oilseed rape varieties, Sansibar and Oase, termed SODH population. Both parental lines are of canola quality which differ in phytosterol and oil content in seed. A genetic map was constructed for SODH population based on a total of 1638 markers organized in 23 linkage groups and covering a map length of 2350 cM with a mean marker interval of 2.0 cM. The SODH population and the parental lines were cultivated at six environments in Europe and were phenotyped for phytosterol content, oil content, fatty acids content, protein content of the defatted meal, and seed weight. Multiple interval mapping identified between one and six QTL for nine phytosterol traits, between two and six QTL for four fatty acids, five QTL for oil content, four QTL for protein content of defatted meal, and three QTL for seed weight. Colocalizations of QTL for different traits were more frequently observed than individual isolated QTL. Major QTL (R (2) ≥ 25%) were all located in the A genome, and the possible candidate genes were investigated by physical localization of the QTL to the reference genome sequence of Brassica rapa.

Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed (Brassica napus L.)

To facilitate the pseudochromosomes assembly and gene cloning in rapeseed, we developed a reference genetic population/map (named BnaZNF2) from two sequenced cultivars, Zhongshuang11 and No.73290, those exhibit significant differences in many traits, particularly yield components. The BnaZNF2 genetic map exhibited perfect collinearity with the physical map of B. napus, indicating its high quality. Comparative mapping revealed several genomic rearrangements between B. napus and B. rapa or B. oleracea. A total of eight and 16 QTLs were identified for pod number and seed number per pod, respectively, and of which three and five QTLs are identical to previously identified ones, whereas the other five and 11 are novel. Two new major QTL respectively for pod number and seed number per pod, qPN.A06-1 and qSN.A06-1 (R(2 )= 22.8% and 32.1%), were colocalised with opposite effects, and only qPN.A06-1 was confirmed and narrowed by regional association analysis to 180 kb including only 33 annotated genes. Conditional QTL analysis and subsequent NILs test indicated that tight linkage, rather than pleiotropy, was the genetic causation of their colocalisation. Our study demonstrates potential of this reference genetic population/map for precise QTL mapping and as a base for positional gene cloning in rapeseed.

Identification and genetic characterization by high-throughput SNP analysis of intervarietal substitution lines of rapeseed (Brassica napus L.) with enhanced embryogenic potential

Seven intervarietal substitution lines were identified with embryogenic potentials up to 40.4 times that of the recurrent parent, providing an ideal material for further in depth studies of this trait. To identify genomic regions that carry genetic factors controlling embryogenic potential of isolated microspores of rapeseed, marker segregations were analysed in a segregating population of haploid microspore-derived embryos and a BC1 population from a cross between 'Express 617' and 'RS239'. After map construction 15 intervarietal substitution lines from the same cross with 'Express 617' as recurrent parent were selected with donor segments covering five genomic regions that had shown skewed segregations in the population of microspore-derived embryos but not in the BC1 population. By comparing the embryogenic potential of microspores of the 15 substitution lines and 'Express 617', seven lines were identified with significantly enhanced embryogenic potential ranging from 4.1 to 40.4 times that of 'Express 617'. To improve the genetic characterization of the selected lines, they were subjected to a high-throughput SNP analysis using the Illumina Infinium 60K chip for rapeseed. Based on 7,960 mapped SNP markers, one to eight donor segments per line, which cover 0.64-6.79% of the 2,126.1 cM of the SNP map, were found. The SNP analysis also gave evidence that homoeologous exchanges had occurred during the development of the substitution line population, increasing the genetic diversity within this population. By comparing donor segments between lines with significantly enhanced embryogenic potential and non-significant lines, 12 genomic regions were identified that may contain genetic factors controlling embryogenic potential in rapeseed. These regions range in size from 0 (represented by just one marker) to 26.8 cM and cover together just 5.42% of the SNP map.

Interacted QTL mapping in partial NCII design provides evidences for breeding by design

The utilization of heterosis in rice, maize and rapeseed has revolutionized crop production. Although elite hybrid cultivars are mainly derived from the F₁ crosses between two groups of parents, named NCII mating design, little has been known about the methodology of how interacted effects influence quantitative trait performance in the population. To bridge genetic analysis with hybrid breeding, here we integrated an interacted QTL mapping approach with breeding by design in partial NCII mating design. All the potential main and interacted effects were included in one full model. If the number of the effects is huge, bulked segregant analysis were used to test which effects were associated with the trait. All the selected effects were further shrunk by empirical Bayesian, so significant effects could be identified. A series of Monte Carlo simulations was performed to validate the new method. Furthermore, all the significant effects were used to calculate genotypic values of all the missing F₁ hybrids, and all these F₁ phenotypic or genotypic values were used to predict elite parents and parental combinations. Finally, the new method was adopted to dissect the genetic foundation of oil content in 441 rapeseed parents and 284 F₁ hybrids. As a result, 8 main-effect QTL and 37 interacted QTL were found and used to predict 10 elite restorer lines, 10 elite sterile lines and 10 elite parental crosses. Similar results across various methods and in previous studies and a high correlation coefficient (0.76) between the predicted and observed phenotypes validated the proposed method in this study.