Design of new genome- and gene-sourced primers and identification of QTL for seed oil content in a specially high-oil Brassica napus cultivar

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.

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.

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.

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.

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.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil 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.

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.

Genetic Properties of a Nested Association Mapping Population Constructed With Semi-Winter and Spring Oilseed Rapes

Nested association mapping (NAM) populations have been widely applied to dissect the genetic basis of complex quantitative traits in a variety of crops. In this study, we developed a Brassica napus NAM (BN-NAM) population consisting of 15 recombination inbred line (RIL) families with 2,425 immortal genotypes. Fifteen high-density genetic linkage maps were constructed by genotyping by sequencing (GBS) based on all RIL families, with further integration into a joint linkage map (JLM) having 30,209 unique markers in common with multiple linkage maps. Furthermore, an ultra-density whole-genome variation map was constructed by projecting 4,444,309 high-quality variants onto the JLM. The NAM population captured a total of 88,542 recombination events (REs). The uneven distribution of recombination rate along chromosomes is positively correlated with the densities of genes and markers, but negatively correlated with the density of transposable elements and linkage disequilibrium (LD). Analyses of population structure and principal components revealed that the BN-NAM population could be divided into three groups with weak stratification. The LD decay distance across genome varied between 170 and 2,400 Kb, with LD decay more rapid in the A than in the C sub-genome. The pericentromeric regions contained large LD blocks, especially in the C sub-genome. This NAM population provides a valuable resource for dissecting the genetic basis of important traits in rapeseed, especially in semi-winter oilseed rape.

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.

High-density ddRAD linkage and yield-related QTL mapping delimits a chromosomal region responsible for oil content in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is one of the most important oil crops almost all over the world. Seed-related traits, including oil content (OC), silique length (SL), seeds per silique (SS), and seed weight (SW), are primary targets for oil yield improvement. To dissect the genetic basis of these traits, 192 recombinant inbred lines (RILs) were derived from two parents with distinct oil content and silique length. High-density linkage map with a total length of 1610.4 cM were constructed using 1,329 double-digestion restriction site associated DNA (ddRAD) markers, 107 insertion/deletions (INDELs), and 90 well-distributed simple sequence repeats (SSRs) markers. A total of 37 consensus quantitative trait loci (QTLs) were detected for the four traits, with individual QTL explained 3.1-12.8% of the phenotypic variations. Interestingly, one OC consensus QTL (cqOCA10b) on chromosome A10 was consistently detected in all three environments, and explained 9.8% to 12.8% of the OC variation. The locus was further delimited into an approximately 614 kb genomic region, in which the flanking markers could be further evaluated for marker-assisted selection in rapeseed OC improvement and the candidate genes targeted for map-based cloning and genetic manipulation.

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.

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.

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.

Arabidopsis Seed Content QTL Mapping Using High-Throughput Phenotyping: The Assets of Near Infrared Spectroscopy

Seed storage compounds are of crucial importance for human diet, feed and industrial uses. In oleo-proteaginous species like rapeseed, seed oil and protein are the qualitative determinants that conferred economic value to the harvested seed. To date, although the biosynthesis pathways of oil and storage protein are rather well-known, the factors that determine how these types of reserves are partitioned in seeds have to be identified. With the aim of implementing a quantitative genetics approach, requiring phenotyping of 100s of plants, our first objective was to establish near-infrared reflectance spectroscopic (NIRS) predictive equations in order to estimate oil, protein, carbon, and nitrogen content in Arabidopsis seed with high-throughput level. Our results demonstrated that NIRS is a powerful non-destructive, high-throughput method to assess the content of these four major components studied in Arabidopsis seed. With this tool in hand, we analyzed Arabidopsis natural variation for these four components and illustrated that they all displayed a wide range of variation. Finally, NIRS was used in order to map QTL for these four traits using seeds from the Arabidopsis thaliana Ct-1 × Col-0 recombinant inbred line population. Some QTL co-localized with QTL previously identified, but others mapped to chromosomal regions never identified so far for such traits. This paper illustrates the usefulness of NIRS predictive equations to perform accurate high-throughput phenotyping of Arabidopsis seed content, opening new perspectives in gene identification following QTL mapping and genome wide association studies.

The Peanut (Arachis hypogaea L.) Gene AhLPAT2 Increases the Lipid Content of Transgenic Arabidopsis Seeds

Lysophosphatidic acid acyltransferase (LPAT), which converts lysophosphatidic acid (LPA) to phosphatidic acid (PA), catalyzes the addition of fatty acyl moieties to the sn-2 position of the LPA glycerol backbone in triacylglycerol (TAG) biosynthesis. We recently reported the cloning and temporal-spatial expression of a peanut (Arachis hypogaea) AhLPAT2gene, showing that an increase in AhLPAT2 transcript levels was closely correlated with an increase in seed oil levels. However, the function of the enzyme encoded by the AhLPAT2 gene remains unclear. Here, we report that AhLPAT2 transcript levels were consistently higher in the seeds of a high-oil cultivar than in those of a low-oil cultivar across different seed developmental stages. Seed-specific overexpression of AhLPAT2 in Arabidopsis results in a higher percentage of oil in the seeds and greater-than-average seed weight in the transgenic plants compared with the wild-type plants, leading to a significant increase in total oil yield per plant. The total fatty acid (FA) content and the proportion of unsaturated FAs also increased. In the developing siliques of AhLPAT2-overexpressing plants, the expression levels of genes encoding crucial enzymes involved in de novo FA synthesis, acetyl-CoA subunit (AtBCCP2) and acyl carrier protein 1 (AtACP1) were elevated. AhLPAT2 overexpression also promoted the expression of several key genes related to TAG assembly, sucrose metabolism, and glycolysis. These results demonstrate that the expression of AhLPAT2 plays an important role in glycerolipid production in peanuts.

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.

Microarray expression analysis of the main inflorescence in Brassica napus

The effect of the number of pods on the main inflorescence (NPMI) on seed yield in Brassica napus plants grown at high density is a topic of great economic and scientific interest. Here, we sought to identify patterns of gene expression that determine the NPMI during inflorescence differentiation. We monitored gene expression profiles in the main inflorescence of two B. napus F₆ RIL pools, each composed of nine lines with a low or high NPMI, and their parental lines, Zhongshuang 11 (ZS11) and 73290, using a Brassica 90K elements oligonucleotide array. We identified 4,805 genes that were differentially expressed (≥1.5 fold-change) between the low- and high-NPMI samples. Of these, 82.8% had been annotated and 17.2% shared no significant homology with any known genes. About 31 enriched GO clusters were identified amongst the differentially expressed genes (DEGs), including those involved in hormone responses, development regulation, carbohydrate metabolism, signal transduction, and transcription regulation. Furthermore, 92.8% of the DEGs mapped to chromosomes that originated from B. rapa and B. oleracea, and 1.6% of the DEGs co-localized with two QTL intervals (PMI10 and PMI11) known to be associated with the NPMI. Overexpression of BnTPI, which co-localized with PMI10, in Arabidopsis suggested that this gene increases the NPMI. This study provides insight into the molecular factors underlying inflorescence architecture, NPMI determination and, consequently, seed yield in B. napus.

A combined linkage and regional association mapping validation and fine mapping of two major pleiotropic QTLs for seed weight and silique length in rapeseed (Brassica napus L.)

BACKGROUND

Seed weight (SW) and silique length (SL) are important determinants of the yield potential in rapeseed (Brassica napus L.). However, the genetic basis of both traits is poorly understood. The main objectives of this study were to dissect the genetic basis of SW and SL in rapeseed through the preliminary mapping of quantitative trait locus (QTL) by linkage analysis and fine mapping of the target major QTL by regional association analysis.

RESULTS

Preliminary linkage mapping identified thirteen and nine consensus QTLs for SW and SL, respectively. These QTLs explained 0.7-67.1% and 2.1-54.4% of the phenotypic variance for SW and SL, respectively. Of these QTLs, three pairs of SW and SL QTLs were co-localized and integrated into three unique QTLs. In addition, the significance level and genetic effect of the three co-localized QTLs for both SW and SL showed great variation before and after the conditional analysis. Moreover, the allelic effects of the three QTLs for SW were highly consistent with those for SL. Two of the three co-localized QTLs, uq.A09-1 (mean R(2) = 20.1% and 19.0% for SW and SL, respectively) and uq.A09-3 (mean R(2) = 13.5% and 13.2% for SW and SL, respectively), were detected in all four environments and showed the opposite additive-effect direction. These QTLs were validated and fine mapped (their confidence intervals were narrowed down from 5.3 cM to 1 cM for uq.A09-1 and 13.2 cM to 2.5 cM for uq.A09-3) by regional association analysis with a panel of 576 inbred lines, which has a relatively rapid linkage disequilibrium decay (0.3 Mb) in the target QTL region.

CONCLUSIONS

A few QTLs with major effects and several QTLs with moderate effects might contribute to the natural variation of SW and SL in rapeseed. The meta-, conditional and allelic effect analyses suggested that pleiotropy, rather than tight linkage, was the genetic basis of the three pairs of co-localized of SW and SL QTLs. Regional association analysis was an effective and highly efficient strategy for the direct fine mapping of target major QTL identified by preliminary linkage mapping.

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

Increasing seed oil content is one of the most important goals in breeding of rapeseed (B. napus L.). To dissect the genetic basis of oil content in B. napus, a large and new double haploid (DH) population containing 348 lines was obtained from a cross between 'KenC-8' and 'N53-2', two varieties with >10% difference in seed oil content, and this population was named the KN DH population. A genetic linkage map consisting of 403 markers was constructed, which covered a total length of 1783.9 cM with an average marker interval of 4.4 cM. The KN DH population was phenotyped in eight natural environments and subjected to quantitative trait loci (QTL) analysis for oil content. A total of 63 identified QTLs explaining 2.64-17.88% of the phenotypic variation were identified, and these QTLs were further integrated into 24 consensus QTLs located on 11 chromosomes using meta-analysis. A high-density consensus map with 1335 marker loci was constructed by combining the KN DH map with seven other published maps based on the common markers. Of the 24 consensus QTLs in the KN DH population, 14 were new QTLs including five new QTLs in A genome and nine in C genome. The analysis revealed that a larger population with significant differences in oil content gave a higher power detecting new QTLs for oil content, and the construction of the consensus map provided a new clue for comparing the QTLs detected in different populations. These findings enriched our knowledge of QTLs for oil content and should be a potential in marker-assisted breeding of B. napus.

Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.

Characterization of the quantitative trait locus OilA1 for oil content in Brassica napus

Increasing seed oil content has become one of the most important breeding criteria in rapeseed (Brassica napus). However, oil content is a complex quantitative trait. QTL mapping in a double haploid population (SG population) emerging from a cross between a German (Sollux) and Chinese (Gaoyou) cultivars revealed one QTL for oil content on linkage group A1 (OilA1), which was mapped to a 17 cM genetic interval. To further validate and characterize the OilA1, we constructed a high-resolution map using B. rapa sequence resources and developed a set of near-isogenic lines (NILs) by employing a DH line SG-DH267 as donor and Chinese parent Gaoyou as recurrent background. The results showed highly conserved synteny order between B. rapa and B. napus within the linkage group A1 and revealed a possible centromere region between two markers ZAASA1-38 and NTP3 (2.5 cM). OilA1 was firstly validated by 250 BC5F2 plants and was confirmed in a 10.6 cM interval between the markers ZAASA1-47 and ZAASA1-77. Further substitution mapping was conducted by using two generations of QTL-NILs, 283 lines from eight BC5F3:4 families and 428 plants from six BC5F4 sub-NILs and thus narrowed the OilA1 interval to 6.9 cM and 4.3 cM (1.4 Mb), respectively. Field investigations with two replications using homozygous BC5F3:4 sister sub-NILs indicated that NILs, which carry a Sollux chromosome segment across the target region showed significant higher oil content (1.26 %, p < 0.001) than their sister NILs containing Gaoyou chromosome. The OilA1 locus is of particular interest for breeding purpose in China because 80 % of Chinese cultivars do not carry this desirable allele.