Novel insights into seed fatty acid synthesis and modification pathways from genetic diversity and quantitative trait Loci analysis of the Brassica C genome

Investigation of mutation load and rate in androgenic mutant lines of rapeseed in early generations evaluated by high-density SNP genotyping

Oilseed rape (Brassica napus) is an important oil crop distributed worldwide with a broad adaptation to different climate zones. The cultivation of rapeseed is one of the most commercially viable areas in crop production. Altogether 269,093 ha of rapeseed are cultivated in Kazakhstan. However, all rapeseed cultivars and lines cultivated in Kazakhstan on an industrial scale predominantly belong to the foreign breeding system. Therefore, the formation of a diverse genetic pool for breeding new, highly productive cultivars adopted to the environmental conditions of Kazakhstan is the most important goal in country selection programs. In this work, we have developed ethyl methanesulfonate (EMS) doubled haploid mutant lines from plant material of cultivars 'Galant' and 'Kris' to broad diversity of rapeseed in Kazakhstan. The development of mutant lines was performed via embryo callusogenesis or embryo secondary callusogenesis. Mutants were investigated by Brassica90k SNP array, and we were able to locate 24,657 SNPs from 26,256 SNPs filtered by quality control on the genome assembly (Bra_napus_v2.0). Only 18,831 SNPs were assigned to the available annotated genomic features. The most frequent combination of mutations according to reference controls was adenine with guanine (70%), followed by adenine with cytosine (28.8%), and only minor fractions were cytosine with guanine (0.54%) and adenine with thymine (0.59%). We revealed 5606.27 markers for 'Kris' and 4893.01 markers for 'Galant' by mutation occurrence. Most mutation occurrences were occupied by double mutations where progenitors and offspring were homozygous by different alleles, enabling the selection of appropriate genotypes in a short period of time. Regarding the biological impact of mutations, 861 variants were reported as having a low predicted impact, with 1042 as moderate and 121 as high; all others were reported as belonging to non-coding sequences, intergenic regions, and other features with the effect of modifiers. Protein encoding genes, such as wall-associated receptor kinase-like protein 5, TAO1-like disease resistance protein, receptor-like protein 12, and At5g42460-like F-box protein, contained more than two variable positions, with an impact on their biological activities. Nevertheless, the obtained mutant lines were able to survive and reproduce. Mutant lines, which include moderate and high impact mutations in encoding genes, are a perfect pool not only for MAS but also for the investigation of the fundamental basis of protein functions. For the first time, a collection of mutant lines was developed in our country to improve the selection of local rapeseed cultivars.

Combining quantitative trait locus and co-expression analysis allowed identification of new candidates for oil accumulation in rapeseed

In crops there are quantitative trait loci (QTLs) in which some of the causal quantitative trait genes (QTGs) have not been functionally characterized even in the model plant Arabidopsis. We propose an approach to delineate QTGs in rapeseed by coordinating expression of genes located within QTLs and known orthologs related to traits from Arabidopsis. Using this method in developing siliques 15 d after pollination in 71 lines of rapeseed, we established an acyl-lipid metabolism co-expression network with 21 modules composed of 270 known acyl-lipid genes and 3503 new genes. The core module harbored 76 known genes involved in fatty acid and triacylglycerol biosynthesis and 671 new genes involved in sucrose transport, carbon metabolism, amino acid metabolism, seed storage protein processes, seed maturation, and phytohormone metabolism. Moreover, the core module closely associated with the modules of photosynthesis and carbon metabolism. From the co-expression network, we selected 12 hub genes to identify their putative Arabidopsis orthologs. These putative orthologs were functionally analysed using Arabidopsis knockout and overexpression lines. Four knockout mutants exhibited lower seed oil content, while the seed oil content in 10 overexpression lines was significantly increased. Therefore, combining gene co-expression network analysis and QTL mapping, this study provides new insights into the detection of QTGs and into acyl-lipid metabolism in rapeseed.

Development and Validation of Kompetitive Allele-Specific PCR Assays for Erucic Acid Content in Indian Mustard [Brassica juncea (L.) Czern and Coss.]

Brassica juncea L. is the most widely cultivated oilseed crop in Indian subcontinent. Its seeds contain oil with very high concentration of erucic acid (≈50%). Of late, there is increasing emphasis on the development of low erucic acid varieties because of reported association of the consumption of high erucic acid oil with cardiac lipidosis. Erucic acid is synthesized from oleic acid by an elongation process involving two cycles of four sequential steps. Of which, the first step is catalyzed by β-ketoacyl-CoA synthase (KCS) encoded by the fatty acid elongase 1 (FAE1) gene in Brassica. Mutations in the coding region of the FAE1 lead to the loss of KCS activity and consequently a drastic reduction of erucic acid in the seeds. Molecular markers have been developed on the basis of variation available in the coding or promoter region(s) of the FAE1. However, majority of these markers are not breeder friendly and are rarely used in the breeding programs. Present studies were planned to develop robust kompetitive allele-specific PCR (KASPar) assays with high throughput and economics of scale. We first cloned and sequenced FAE1.1 and FAE1.2 from high and low erucic acid (<2%) genotypes of B. juncea (AABB) and its progenitor species, B. rapa (AA) and B. nigra (BB). Sequence comparisons of FAE1.1 and FAE1.2 genes for low and high erucic acid genotypes revealed single nucleotide polymorphisms (SNPs) at 8 and 3 positions. Of these, three SNPs for FAE1.1 and one SNPs for FAE1.2 produced missense mutations, leading to amino acid modifications and inactivation of KCS enzyme. We used SNPs at positions 735 and 1,476 for genes FAE1.1 and FAE1.2, respectively, to develop KASPar assays. These markers were validated on a collection of diverse genotypes and a segregating backcross progeny. KASPar assays developed in this study will be useful for marker-assisted breeding, as these can track recessive alleles in their heterozygous state with high reproducibility.

Association Mapping of Seed Quality Traits Under Varying Conditions of Nitrogen Application in Brassica juncea L. Czern & Coss

Indian mustard (Brassica juncea) is a major source of vegetable oil in the Indian subcontinent. The seed cake left after the oil extraction is used as livestock feed. We examined the genetic architecture of oil, protein, and glucosinolates by conducting a genome-wide association study (GWAS), using an association panel comprising 92 diverse genotypes. We conducted trait phenotyping over 2 years at two levels of nitrogen (N) application. Genotyping by sequencing was used to identify 66,835 loci, covering 18 chromosomes. Genetic diversity and phenotypic variations were high for the studied traits. Trait performances were stable when averaged over years and N levels. However, individual performances differed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits. Population structure, principal components (PCs) analysis, and discriminant analysis of principal components (DAPCs) were included as covariates to overcome the bias due to the population stratification. We identified 16, 23, and 27 loci associated with oil, protein, and glucosinolates, respectively. We also established LD patterns and haplotype structures for the candidate genes. The average block sizes were larger on A-genome chromosomes as compared to the B- genome chromosomes. Genetic associations differed over N levels. However, meta-analysis of GWAS datasets not only improved the power to recognize associations but also helped to identify common SNPs for oil and protein contents. Annotation of the genomic region around the identified SNPs led to the prediction of 21 orthologs of the functional candidate genes related to the biosynthesis of oil, protein, and glucosinolates. Notable among these are: LACS5 (A09), FAD6 (B05), ASN1 (A06), GTR2 (A06), CYP81G1 (B06), and MYB44 (B06). The identified loci will be very useful for marker-aided breeding for seed quality modifications in B. juncea.

Identification of Climate and Genetic Factors That Control Fat Content and Fatty Acid Composition of Theobroma cacao L. Beans

The main ingredients of chocolate are usually cocoa powder, cocoa butter, and sugar. Both the powder and the butter are extracted from the beans of the cacao tree (Theobroma cacao L.). The cocoa butter represents the fat in the beans and possesses a unique fatty acid profile that results in chocolate's characteristic texture and mouthfeel. Here, we used a linkage mapping population and phenotypic data of 3,292 samples from 420 progeny which led to the identification of 27 quantitative trait loci (QTLs) for fatty acid composition and six QTLs for fat content. Progeny showed extensive variation in fat levels and composition, with the level of palmitic acid negatively correlated to the sum of stearic acid, oleic acid, and linoleic acid. A major QTL explaining 24% of the relative level of palmitic acid was mapped to the distal end of chromosome 4, and those higher levels of palmitic acid were associated with the presence of a haplotype from the "TSH 1188" parent in the progeny. Within this region of chromosome 4 is the Thecc1EG017405 gene, an orthologue and isoform of the stearoyl-acyl carrier protein (ACP) desaturase (SAD) gene in plants, which is involved in fatty acid biosynthesis. Besides allelic differences, we also show that climate factors can change the fatty acid composition in the beans, including a significant positive correlation between higher temperatures and the higher level of palmitic acid. Moreover, we found a significant pollen donor effect from the variety "SIAL 70" which was associated with decreased palmitic acid levels.

Identification of SNP loci and candidate genes related to four important fatty acid composition in Brassica napus using genome wide association study

Rapeseed oil (canola, Brassica napus L.) is an important healthy vegetable oil throughout the world, the nutritional and economical value of which largely depends on its seed fatty acid composition. In this study, based on 201,187 SNP markers developed from the SLAF-seq (specific locus amplified fragment sequencing), a genome wide association study of four important fatty acid content traits (erucic acid, oleic acid, linoleic acid and linolenic acid) in a panel of 300 inbred lines of rapeseed in two environments (JXAU and JXRIS) was carried out. A total of 148 SNP loci significantly associated with these traits were detected by MLM model analysis respectively, and 30 SNP loci on A08 and C03 chromosomes were detected in three traits of erucic acid, oleic acid and linoleic acid contents simultaneously. Furthermore, 108 highly favorable alleles for increasing oleic acid and linoleic acid content, also for decreasing erucic acid content simultaneously were observed. By a basic local alignment search tool (BLAST) search with in a distance of 100 Kb around these significantly SNP-trait associations, we identified 20 orthologs of the functional candidate genes related to fatty acid biosynthesis, including the known vital fatty acid biosynthesis genes of BnaA.FAE1 and BnaC. FAE1 on the A08 and C03 chromosomes, and other potential candidate genes involving in the fatty acid biosynthesis pathway, such as the orthologs genes of FAD2, LACS09, KCS17, CER4, TT16 and ACBP5. This study lays a basis for uncovering the genetic variations and the improvement of fatty acid composition in B. napus.

Increased biomass and lipid production by continuous cultivation of Nannochloropsis salina transformant overexpressing a bHLH transcription factor

Microalgae are promising feedstocks for sustainable and eco-friendly production of biomaterials, which can be improved by genetic engineering. It is also necessary to optimize the processes to produce biomaterials from engineered microalgae. We previously reported that genetic improvements of an industrial microalga Nannochloropsis salina by overexpressing a basic helix-loop-helix transcription factor (NsbHLH2). These transformants showed an improved growth and lipid production particularly during the early phase of culture under batch culture. However, they had faster uptake of nutrients, resulting in earlier starvation and reduced growth during the later stages. We attempted to optimize the growth and lipid production by growing one of the transformants in continuous culture with variable dilution rate and feed nitrogen concentration. Relative to wild-type, NsbHLH2 transformant consumed more nitrate at a high dilution rate (0.5 day ^-1 ), and had greater biomass production. Subsequently, nitrogen limitation at continuous cultivation led to an increased fatty acid methyl ester production by 83.6 mg l ^-1 day ^-1 . To elucidate genetic mechanisms, we identified the genes containing E-boxes, known as binding sites for bHLH transcription factors. Among these, we selected 18 genes involved in the growth and lipid metabolism, and revealed their positive contribution to the phenotypes via quantitative real-time polymerase chain reaction. These results provide proof-of-concept that NsbHLH2 can be used to produce biomass and lipids.

The study of flavonolignan association patterns in fruits of diverging Silybum marianum (L.) Gaertn. chemotypes provides new insights into the silymarin biosynthetic pathway

Silymarin is the phytochemical with medicinal properties extracted from Silybum marianum (L.) Gaertn. fruits. Yet, little information is available about silymarin biosynthesis. Moreover, the generally accepted pathway, formulated thus far, is not in agreement with actual experimental measurements on flavonolignan contents. The present work analyses flavonolignan and taxifolin content in 201 S. marianum samples taking into consideration a wide phenotypic variability. Two stable chemotypes were identified: one characterized by both high silychristin and silybin content (chemotype A) and another by a high silydianin content (chemotype B). Through the correlation analysis of samples divided according to chemotype, it was possible to construct a simplified silymarin biosynthetic pathway that is sufficiently versatile in explaining experimental results responding to the actually unresolved questions about this process. The proposed pathway highlights that three separate and equally sized metabolite pools exist, namely: diastereoisomers A (silybin A plus isosilybin A), diastereoisomers B (silybin B plus isosilybin B) and silychristin. In both A and B diastereoisomers pools, isosilybin A and isosilybin B always represent a given amount of the metabolite flux through the specific metabolite pool suggesting the possible involvement of dirigent protein-like enzymes. We suggest that chemotype B possesses a complete silymarin biosynthetic pathway in which silydianin biosynthesis is enzymatically controlled. On the contrary, chemotype A is probably a natural mutant unable to biosynthesize silydianin. The present simplified pathway for silymarin biosynthesis will constitute an important tool for the further understanding of the reactions that drive flavonolignan biosynthesis in S. marianum.

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.

Genome-Wide Association Study of Genetic Control of Seed Fatty Acid Biosynthesis in Brassica napus

Fatty acids and their composition in seeds determine oil value for nutritional or industrial purposes and also affect seed germination as well as seedling establishment. To better understand the genetic basis of seed fatty acid biosynthesis in oilseed rape (Brassica napus L.) we applied a genome-wide association study, using 91,205 single nucleotide polymorphisms (SNPs) characterized across a mapping population with high-resolution skim genotyping by sequencing (SkimGBS). We identified a cluster of loci on chromosome A05 associated with oleic and linoleic seed fatty acids. The delineated genomic region contained orthologs of the Arabidopsis thaliana genes known to play a role in regulation of seed fatty acid biosynthesis such as Fatty acyl-ACP thioesterase B (FATB) and Fatty Acid Desaturase (FAD5). This approach allowed us to identify potential functional genes regulating fatty acid composition in this important oil producing crop and demonstrates that this approach can be used as a powerful tool for dissecting complex traits for B. napus improvement programs.

Genetic mapping of QTLs controlling fatty acids provided insights into the genetic control of fatty acid synthesis pathway in peanut (Arachis hypogaea L.)

Peanut, a high-oil crop with about 50% oil content, is either crushed for oil or used as edible products. Fatty acid composition determines the oil quality which has high relevance to consumer health, flavor, and shelf life of commercial products. In addition to the major fatty acids, oleic acid (C18:1) and linoleic acid (C18:2) accounting for about 80% of peanut oil, the six other fatty acids namely palmitic acid (C16:0), stearic acid (C18:0), arachidic acid (C20:0), gadoleic acid (C20:1), behenic acid (C22:0), and lignoceric acid (C24:0) are accounted for the rest 20%. To determine the genetic basis and to improve further understanding on effect of FAD2 genes on these fatty acids, two recombinant inbred line (RIL) populations namely S-population (high oleic line 'SunOleic 97R' × low oleic line 'NC94022') and T-population (normal oleic line 'Tifrunner' × low oleic line 'GT-C20') were developed. Genetic maps with 206 and 378 marker loci for the S- and the T-population, respectively were used for quantitative trait locus (QTL) analysis. As a result, a total of 164 main-effect (M-QTLs) and 27 epistatic (E-QTLs) QTLs associated with the minor fatty acids were identified with 0.16% to 40.56% phenotypic variation explained (PVE). Thirty four major QTLs (>10% of PVE) mapped on five linkage groups and 28 clusters containing more than three QTLs were also identified. These results suggest that the major QTLs with large additive effects would play an important role in controlling composition of these minor fatty acids in addition to the oleic and linoleic acids in peanut oil. The interrelationship among these fatty acids should be considered while breeding for improved peanut genotypes with good oil quality and desired fatty acid composition.

New insights into the genetic networks affecting seed fatty acid concentrations in Brassica napus

BACKGROUND

Rapeseed (B. napus, AACC, 2n = 38) is one of the most important oil seed crops in the world, it is also one of the most common oil for production of biodiesel. Its oil is a mixture of various fatty acids and dissection of the genetic network for fatty acids biosynthesis is of great importance for improving seed quality.

RESULTS

The genetic basis of fatty acid biosynthesis in B. napus was investigated via quantitative trail locus (QTL) analysis using a doubled haploid (DH) population with 202 lines. A total of 72 individual QTLs and a large number pairs of epistatic interactions associated with the content of 10 different fatty acids were detected. A total of 234 homologous genes of Arabidopsis thaliana that are involved in fatty acid metabolism were found within the confidence intervals (CIs) of 47 QTLs. Among them, 47 and 15 genes homologous to those of B. rapa and B. oleracea were detected, respectively. After the QTL mapping, the epistatic and the candidate gene interaction analysis, a potential regulatory pathway controlling fatty acid biosynthesis in B. napus was constructed, including 50 enzymes encoded genes and five regulatory factors (LEC1, LEC2, FUS3, WRI1 and ABI3). Subsequently, the interaction between these five regulatory factors and the genes involved in fatty acid metabolism were analyzed.

CONCLUSIONS

In this study, a potential regulatory pathway controlling the fatty acid was constructed by QTL analysis and in silico mapping analysis. These results enriched our knowledge of QTLs for fatty acids metabolism and provided a new clue for genetic engineering fatty acids composition in B. napus.

Development and characterization of low α-linolenic acid Brassica oleracea lines bearing a novel mutation in a 'class a' FATTY ACID DESATURASE 3 gene

Background

Traditional canola (Brassica napus L.; AACC, 2n = 38) cultivars yield seed oil with a relatively high proportion of α-linolenic acid (ALA; C18:3^{cis∆9,12,15}), which is desirable from a health perspective. Unfortunately, due to the instability of this fatty acid, elevated levels also result in oils that exhibit a short shelf life and problems associated with use at high temperatures. As a result, the development of cultivars bearing reduced amounts of ALA in their seeds is becoming a priority. To date, several low ALA B. napus cultivars (~2-3% ALA of total fatty acids) have been developed and molecular analyses have revealed that the low ALA phenotype of lines tested thus far is a result of mutations within two ‘class b’ FATTY ACID DESATURASE 3 (FAD3) genes. Since B. napus possesses six FAD3 genes (two ‘class a’, two ‘class b’ and two ‘class c’) and ALA levels of approximately 2-3% remain in these low ALA lines, it is likely that the mutation of additional FAD3 genes could further decrease the content of this fatty acid.

Results

In this study, we generated low ALA (≤2%) lines of B. oleracea, which is the C genome progenitor species of B. napus, via ethyl methanesulphonate (EMS) mutagenesis. We identified a novel nonsense mutation within the ‘class a’ FAD3 gene (BoFAD3-2) in these lines, which would result in the production of an encoded protein lacking 110 amino acids at its C terminus. When expressed in Saccharomyces cerevisiae, this mutant protein exhibited a drastic decline in its Δ-15 desaturase activity compared to the wild-type (wt) protein. Furthermore, we demonstrated that the expression of the mutant BoFAD3-2 gene was significantly reduced in developing seeds of low ALA lines when compared to expression in wt plants.

Conclusions

Given the additive nature of FAD3 mutations on ALA content and the ease with which B. napus can be re-synthesized from its progenitor species, the mutant isolated here has the potential to be used for the future development of B. napus cultivars exhibiting further reductions in ALA content.

Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.)

This report describes an integrative analysis of seed-oil-content quantitative trait loci (QTL) in Brassica napus , using a high-density genetic map to align QTL among different populations. Rapeseed (Brassica napus) is an important source of edible oil and sustainable energy. Given the challenge involved in using only a few genes to substantially increase the oil content of rapeseed without affecting the fatty acid composition, exploitation of a greater number of genetic loci that regulate the oil content variation among rapeseed germplasm is of fundamental importance. In this study, we investigated variation in the seed-oil content among two related genetic populations of Brassica napus, the TN double-haploid population and its derivative reconstructed-F2 population. Each population was grown in multiple experiments under different environmental conditions. Mapping of quantitative trait loci (QTL) identified 41 QTL in the TN populations. Furthermore, of the 20 pairs of epistatic interaction loci detected, approximately one-third were located within the QTL intervals. The use of common markers on different genetic maps and the TN genetic map as a reference enabled us to project QTL from an additional three genetic populations onto the TN genetic map. In summary, we used the TN genetic map of the B. napus genome to identify 46 distinct QTL regions that control seed-oil content on 16 of the 19 linkage groups of B. napus. Of these, 18 were each detected in multiple populations. The present results are of value for ongoing efforts to breed rapeseed with high oil content, and alignment of the QTL makes an important contribution to the development of an integrative system for genetic studies of rapeseed.

Genetic enhancement of Brassica napus seed quality

The ultimate value of the Brassica napus (canola) seed is derived from the oil fraction, which has long been recognized for its premium dietary attributes, including its low level of saturated fatty acids, high content of monounsaturated fatty acids, and favorable omega-3 fatty acid profile. However, the protein (meal) portion of the seed has also received favorable attention for its essential amino acids, including abundance of sulfur-containing amino acids, such that B. napus protein is being contemplated for large scale use in livestock and fish feed formulations. Efforts to optimize the composition of B. napus oil and protein fractions are well documented; therefore, this article will review research concerned with optimizing secondary metabolites that affect the quality of seed oil and meal, from undesirable anti-nutritional factors to highl value beneficial products. The biological, agronomic, and economic values attributed to secondary metabolites have brought much needed attention to those in Brassica oilseeds and other crops. This review focuses on increasing levels of beneficial endogenous secondary metabolites (such as carotenoids, choline and tochopherols) and decreasing undesirable antinutritional factors (glucosinolates, sinapine and phytate). Molecular genetic approaches are given emphasis relative to classical breeding.

Genotyping by RAD sequencing enables mapping of fatty acid composition traits in perennial ryegrass (Lolium perenne (L.))

Perennial ryegrass (Lolium perenne L.) is the most important forage crop in temperate livestock agriculture. Its nutritional quality has significant impact on the quality of meat and milk for human consumption. Evidence suggests that higher energy content in forage can assist in reducing greenhouse gas emissions from ruminants. Increasing the fatty acid content (especially α-linolenic acid, an omega-3 fatty acid) may thus contribute to better forage, but little is known about the genetic basis of variation for this trait. To this end, quantitative trait loci (QTLs) were identified associated with major fatty acid content in perennial ryegrass using a population derived from a cross between the heterozygous and outbreeding high-sugar grass variety AberMagic and an older variety, Aurora. A genetic map with 434 restriction-associated DNA (RAD) and SSR markers was generated. Significant QTLs for the content of palmitic (C16:0) on linkage groups (LGs) 2 and 7; stearic (C18:0) on LGs 3, 4 and 7; linoleic (C18:2n-6) on LGs 2 and 5; and α-linolenic acids (C18:3n-3) on LG 1 were identified. Two candidate genes (a lipase and a beta-ketoacyl CoA synthase), both associated with C16:0, and separately with C18:2n-6 and C18:0 contents, were identified. The physical positions of these genes in rice and their genetic positions in perennial ryegrass were consistent with established syntenic relationships between these two species. Validation of these associations is required, but the utility of RAD markers for rapid generation of genetic maps and QTL analysis has been demonstrated for fatty acid composition in a global forage crop.

Suppression of the SUGAR-DEPENDENT1 triacylglycerol lipase family during seed development enhances oil yield in oilseed rape (Brassica napus L.)

Increasing the productivity of oilseed crops is an important challenge for plant breeders and biotechnologists. To date, attempts to increase oil production in seeds via metabolic pathway engineering have focused on boosting synthetic capacity. However, in the tissues of many organisms, it is well established that oil levels are determined by both anabolism and catabolism. Indeed, the oil content of rapeseed (Brassica napus L.) has been reported to decline by approximately 10% in the final stage of development, as the seeds desiccate. Here, we show that RNAi suppression of the SUGAR-DEPENDENT1 triacylglycerol lipase gene family during seed development results in up to an 8% gain in oil yield on either a seed, plant or unit area basis in the greenhouse, with very little adverse impact on seed vigour. Suppression of lipolysis could therefore constitute a new method for enhancing oil yield in oilseed crops.

A hypomethylated population of Brassica rapa for forward and reverse epi-genetics

BACKGROUND

Epigenetic marks superimposed on the DNA sequence of eukaryote chromosomes provide agility and plasticity in terms of modulating gene expression, ontology, and response to the environment. Modulating the methylation status of cytosine can generate epialleles, which have been detected and characterised at specific loci in several plant systems, and have the potential to generate novel and relatively stable phenotypes. There have been no systematic attempts to explore and utilise epiallelic variation, and so extend the range of phenotypes available for selection in crop improvement. We developed an approach for generating novel epialleles by perturbation of the DNA methylation status. 5- Azacytidine (5-AzaC) provides selective targeting of 5 mCG, which in plants is associated with exonic DNA. Targeted chemical intervention using 5-AzaC has advantages over transgenic or mutant modulation of methyltransferases, allowing stochastic generation of epialleles across the genome.

RESULTS

We demonstrate the potential of stochastic chemically-induced hypomethylation to generate novel and valuable variation for crop improvement. Systematic analysis of dose-response to 5-AzaC in B. rapa guided generation of a selfed stochastically hypomethylated population, used for forward screening of several agronomic traits. Dose-response was sigmoidal for several traits, similar to that observed for chemical mutagens such as EMS. We demonstrated transgenerational inheritance of some phenotypes. BraRoAZ is a unique hypomethylated population of 1000 E2 sib lines. When compared to untreated controls, 5-Aza C-treated lines exhibited reduced immuno-staining of 5mC on pachytene chromosomes, and Methylation Sensitive Amplified Polymorphism (MSAP) profiles that were both divergent and more variable. There was coincident phenotypic variation among these lines for a range of seed yield and composition traits, including increased seed protein content and decreased oil content, as well as decreased erucic acid and corresponding increases in linoleic and/or palmitic acid. Each 5-AzaC-treated line represents a unique combination of hypomethylated epialleles.

CONCLUSIONS

The approach and populations developed are available for forward and reverse screening of epiallelic variation and subsequent functional and inheritance studies. The generation of stochastically hypomethylated populations has utility in epiallele discovery for a wide range of crop plants, and has considerable potential as an intervention strategy for crop improvement.

A hypomethylated population of Brassica rapa for forward and reverse epi-genetics

BACKGROUND

Epigenetic marks superimposed on the DNA sequence of eukaryote chromosomes provide agility and plasticity in terms of modulating gene expression, ontology, and response to the environment. Modulating the methylation status of cytosine can generate epialleles, which have been detected and characterised at specific loci in several plant systems, and have the potential to generate novel and relatively stable phenotypes. There have been no systematic attempts to explore and utilise epiallelic variation, and so extend the range of phenotypes available for selection in crop improvement. We developed an approach for generating novel epialleles by perturbation of the DNA methylation status. 5- Azacytidine (5-AzaC) provides selective targeting of 5 mCG, which in plants is associated with exonic DNA. Targeted chemical intervention using 5-AzaC has advantages over transgenic or mutant modulation of methyltransferases, allowing stochastic generation of epialleles across the genome.

RESULTS

We demonstrate the potential of stochastic chemically-induced hypomethylation to generate novel and valuable variation for crop improvement. Systematic analysis of dose-response to 5-AzaC in B. rapa guided generation of a selfed stochastically hypomethylated population, used for forward screening of several agronomic traits. Dose-response was sigmoidal for several traits, similar to that observed for chemical mutagens such as EMS. We demonstrated transgenerational inheritance of some phenotypes. BraRoAZ is a unique hypomethylated population of 1000 E2 sib lines. When compared to untreated controls, 5-Aza C-treated lines exhibited reduced immuno-staining of 5mC on pachytene chromosomes, and Methylation Sensitive Amplified Polymorphism (MSAP) profiles that were both divergent and more variable. There was coincident phenotypic variation among these lines for a range of seed yield and composition traits, including increased seed protein content and decreased oil content, as well as decreased erucic acid and corresponding increases in linoleic and/or palmitic acid. Each 5-AzaC-treated line represents a unique combination of hypomethylated epialleles.

CONCLUSIONS

The approach and populations developed are available for forward and reverse screening of epiallelic variation and subsequent functional and inheritance studies. The generation of stochastically hypomethylated populations has utility in epiallele discovery for a wide range of crop plants, and has considerable potential as an intervention strategy for crop improvement.

Identification of quantitative trait Loci for lipid metabolism in rice seeds

Plant seed oil is important for human dietary consumption and industrial application. The oil trait is controlled by quantitative trait loci (QTLs), but no QTLs for fatty acid composition are known in rice, the monocot model plant. QTL analysis was performed using F(2) and F(2:3) progeny from a cross of an indica variety and a japonica variety. Gas chromatography-mass spectrometry (GC-MS) analysis revealed significant differences between parental lines in fatty acid composition of brown rice oil, and 29 associated QTLs in F(2) and/or F(2:3) populations were identified throughout the rice genome, except chromosomes 9 and 10. Eight QTLs were repeatedly identified in both populations across different environments. Five loci pleiotropically controlled different traits, contributing to complex interactions of oil with fatty acids and between fatty acids. Nine rice orthologs of Arabidopsis genes encoding key enzymes in lipid metabolism co-localized with 11 mapped QTLs. A strong QTL for oleic (18:1) and linoleic (18:2) acid were associated with a rice ortholog of a gene encoding acyl-CoA:diacylglycerol acyltransferase (DGAT), and another for palmitic acid (16:0) mapped similarly to the acyl-ACP thioesterase (FatB) gene ortholog. Our approach rapidly and efficiently identified candidate genes for mapped QTLs controlling fatty acid composition and oil concentration, providing information for improving rice grain quality by marker assisted selection.

InterStoreDB: a generic integration resource for genetic and genomic data

Associating phenotypic traits and quantitative trait loci (QTL) to causative regions of the underlying genome is a key goal in agricultural research. InterStoreDB is a suite of integrated databases designed to assist in this process. The individual databases are species independent and generic in design, providing access to curated datasets relating to plant populations, phenotypic traits, genetic maps, marker loci and QTL, with links to functional gene annotation and genomic sequence data. Each component database provides access to associated metadata, including data provenance and parameters used in analyses, thus providing users with information to evaluate the relative worth of any associations identified. The databases include CropStoreDB, for management of population, genetic map, QTL and trait measurement data, SeqStoreDB for sequence-related data and AlignStoreDB, which stores sequence alignment information, and allows navigation between genetic and genomic datasets. Genetic maps are visualized and compared using the CMAP tool, and functional annotation from sequenced genomes is provided via an EnsEMBL-based genome browser. This framework facilitates navigation of the multiple biological domains involved in genetics and genomics research in a transparent manner within a single portal. We demonstrate the value of InterStoreDB as a tool for Brassica research. InterStoreDB is available from: http://www.interstoredb.org.

Quantitative trait loci involved in regulating seed oil composition in Arabidopsis thaliana and their evolutionary implications

Fatty acid composition is an important determinant of seed oil quality. Overall, 72 QTL for 12 fatty acid traits that control seed oil composition were identified in four recombinant inbred line (RIL) populations (Ler-0 × Sha, Ler-0 × Col-4, Ler-2 × Cvi, Ler-0 × No-0) of Arabidopsis thaliana. The identified QTL explained 3.2-79.8% of the phenotypic variance; 33 of the 59 QTL identified in the Ler-0 × Sha and the Ler-0 × Col RIL populations co-located with several a priori candidate genes for seed oil composition. QTL for fatty acids 18:1, 18:2, 22:1, and fatty acids synthesized in plastids was identified in both Ler-0 × Sha and Ler-0 × Col-4 RIL populations, and QTL for 16:0 was identified in the Ler-0 × Sha and Ler-0 × No-0 RIL populations providing strong support for the importance of these QTL in determining seed oil composition. We identified melting point QTL in three RIL populations, and fatty acid QTL collocated with two of them, suggesting that the loci could be under selection for altering the melting point of seed oils to enhance adaptation and could be useful for breeding purposes. Nuclear-cytoplasmic interactions and epistasis were rare. Analysis of the genetic correlations between these loci and other fatty acids indicated that these correlations would tend to strongly enhance selection for desirable fatty acids.

Characterization of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves of Brassica napus

• Glucosinolates are a major class of secondary metabolites found in the Brassicaceae, whose degradation products are proving to be increasingly important for human health and in crop protection. • The genetic and metabolic basis of glucosinolate accumulation was dissected through analysis of total glucosinolate concentration and its individual components in both leaves and seeds of a doubled-haploid (DH) mapping population of oilseed rape/canola (Brassica napus). • The quantitative trait loci (QTL) that had an effect on glucosinolate concentration in either or both of the organs were integrated, resulting in 105 metabolite QTL (mQTL). Pairwise correlations between individual glucosinolates and prior knowledge of the metabolic pathways involved in the biosynthesis of different glucosinolates allowed us to predict the function of genes underlying the mQTL. Moreover, this information allowed us to construct an advanced metabolic network and associated epistatic interactions responsible for the glucosinolate composition in both leaves and seeds of B. napus. • A number of previously unknown potential regulatory relationships involved in glucosinolate synthesis were identified and this study illustrates how genetic variation can affect a biochemical pathway.

Genetic variation in Mexican Jatropha curcas L. estimated with seed oil fatty acids

The genetic diversity of Mesoamerican populations of the biofuel plant Jatropha curcas, using the fatty acids of the seeds as chemical markers was studied. The oil content of the whole seed in 135 accessions from 38 sites varied between 8.020% and 54.28%. The prevalent fatty acids were oleic acid (18:1) and linoleic acid (18:2), and the proportion of unsaturated fatty acids varied between 74.5% and 83.7%. A study with cloned plants grown in common garden showed that both the content of oil as well as the proportion of fatty acids are highly inheritable, therefore these chemical markers are valid for estimating the genetic diversity of the species. An analysis of principal components showed that the fatty acids that contribute more to the variance are stearic, oleic, linoleic, methylpalmitic, gadoleic and ricinoleic. The populations were classified in ten groups when the data were analyzed for fatty acids by analysis of clusters, showing the elevated genetic variation in natural populations of this native species of Mesoamerica. A discriminant analysis separated the populations in accordance with their geographic origin, which was verified with a Mantel test. Using the Monmonier's algorithm two genetic barriers between the populations were identified. The results are discussed in light of their microevolutionary significance.

Eliminating expression of erucic acid-encoding loci allows the identification of "hidden" QTL contributing to oil quality fractions and oil content in Brassica juncea (Indian mustard)

Oil content and oil quality fractions (viz., oleic, linoleic and linolenic acid) are strongly influenced by the erucic acid pathway in oilseed Brassicas. Low levels of erucic acid in seed oil increases oleic acid content to nutritionally desirable levels, but also increases the linoleic and linolenic acid fractions and reduces oil content in Indian mustard (Brassica juncea). Analysis of phenotypic variability for oil quality fractions among a high-erucic Indian variety (Varuna), a low-erucic east-European variety (Heera) and a zero-erucic Indian variety (ZE-Varuna) developed by backcross breeding in this study indicated that lower levels of linoleic and linolenic acid in Varuna are due to substrate limitation caused by an active erucic acid pathway and not due to weaker alleles or enzyme limitation. To identify compensatory loci that could be used to increase oil content and maintain desirable levels of oil quality fractions under zero-erucic conditions, we performed Quantitative Trait Loci (QTL) mapping for the above traits on two independent F1 doubled haploid (F1DH) mapping populations developed from a cross between Varuna and Heera. One of the populations comprised plants segregating for erucic acid content (SE) and was used earlier for construction of a linkage map and QTL mapping of several yield-influencing traits in B. juncea. The second population consisted of zero-erucic acid individuals (ZE) for which, an Amplified Fragment Length Polymorphism (AFLP)-based framework linkage map was constructed in the present study. By QTL mapping for oil quality fractions and oil content in the ZE population, we detected novel loci contributing to the above traits. These loci did not co-localize with mapped locations of the fatty acid desaturase 2 (FAD2), fatty acid desaturase 3 (FAD3) or fatty acid elongase (FAE) genes unlike those of the SE population wherein major QTL were found to coincide with mapped locations of the FAE genes. Some of the new loci identified in the ZE population could be detected as 'weak' contributors (with LOD < 2.5) in the SE population in which their contribution to the traits was "masked" due to pleiotropic effects of erucic acid genes. The novel loci identified in this study could now be used to improve oil quality parameters and oil content in B. juncea under zero-erucic conditions.

Yield assessment of integument-led seed growth following targeted repair of auxin response factor 2

It is becoming increasingly vital to improve the yield of seed crops to feed an expanding population and, more recently, for biofuel production. One strategy to increase the yield is to increase the seed size, provided that there is not a concomitant decrease in seed number. In a previous study, we described a mutant in the auxin response factor 2 (ARF2) gene which produced extra cells in the seed coat and, subsequently, enlarged seeds. However, arf2 mutant plants also show severely reduced self-fertility caused, in part, by over-elongated sepals that prevent flower opening. As a low seed set increases individual seed size, a meaningful comparison of the yield in arf2 and wild-type plants could not be conducted. In this study, we show that targeted expression of wild-type ARF2 in the sepals and petals of arf2-9 mutant flowers restores flower opening and dramatically increases seed set. The restored plants retain both enlarged integuments and increased seed size, reinforcing previous evidence that arf2 mutations increase seed weight through their effect on integuments and not only via reduced fertility. We also show that the measurement of the harvest index in Arabidopsis is useful in assessing the impact of introduced traits on the yield.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F(1) hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40%). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14% to 55% of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca(2+) and Mg(2+) uptake and homeostasis at the molecular level.

A functional genomics resource for Brassica napus: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING

Two ethylmethanesulfonate (EMS) mutant populations of the semi-winter rapeseed cv. Ningyou7 were constructed with high mutant load, to provide a TILLING platform for functional genomics in Brassica napus, and for introduction of novel allelic variation in rapeseed breeding. Forward genetic screening of mutants from the M2 populations resulted in identification of a large number of novel phenotypes. Reverse genetic screening focused on the potentially multi-paralogous gene FAE1 (fatty acid elongase1), which controls seed erucic acid synthesis in rapeseed. A B. napus BAC library was screened, and loci in a reference mapping population (TNDH) were mapped to conclude that there are two paralogous copies of FAE1, one on each of the B. napus A and C genomes. A new procedure is demonstrated to identify novel mutations in situations where two or more very similar paralogous gene copies exist in a genome. The procedure involves TILLING of single plants, using existing SNPs as a positive control, and is able to distinguish novel mutations based on primer pairs designed to amplify both FAE1 paralogues simultaneously. The procedure was applied to 1344 M2 plants, with 19 mutations identified, of which three were functionally compromised with reduced seed erucic acid content.