Creating yellow seed Camelina sativa with enhanced oil accumulation by CRISPR-mediated disruption of Transparent Testa 8.

Camelina (Camelina sativa L.), a hexaploid member of the Brassicaceae family, is an emerging oilseed crop being developed to meet the increasing demand for plant oils as biofuel feedstocks. In other Brassicas, high oil content can be associated with a yellow seed phenotype, which is unknown for camelina. We sought to create yellow seed camelina using CRISPR/Cas9 technology to disrupt its Transparent Testa 8 (TT8) transcription factor genes and to evaluate the resulting seed phenotype. We identified three TT8 genes, one in each of the three camelina subgenomes, and obtained independent CsTT8 lines containing frameshift edits. Disruption of TT8 caused seed coat colour to change from brown to yellow reflecting their reduced flavonoid accumulation of up to 44%, and the loss of a well-organized seed coat mucilage layer. Transcriptomic analysis of CsTT8-edited seeds revealed significantly increased expression of the lipid-related transcription factors LEC1, LEC2, FUS3, and WRI1 and their downstream fatty acid synthesis-related targets. These changes caused metabolic remodelling with increased fatty acid synthesis rates and corresponding increases in total fatty acid (TFA) accumulation from 32.4% to as high as 38.0% of seed weight, and TAG yield by more than 21% without significant changes in starch or protein levels compared to parental line. These data highlight the effectiveness of CRISPR in creating novel enhanced-oil germplasm in camelina. The resulting lines may directly contribute to future net-zero carbon energy production or be combined with other traits to produce desired lipid-derived bioproducts at high yields.

Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis.

BACKGROUND: Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions. RESULTS: To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts. CONCLUSIONS: Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

CYCLIN-DEPENDENT KINASE 8 positively regulates oil synthesis by activating WRINKLED1 transcription.

CYCLIN-DEPENDENT KINASE 8 (CDK8), a component of the kinase module of the Mediator complex in Arabidopsis, is involved in many processes, including flowering, plant defense, drought, and energy stress responses. Here, we investigated cdk8 mutants and CDK8-overexpressing lines to evaluate whether CDK8 also plays a role in regulating lipid synthesis, an energy-demanding anabolism. Quantitative lipid analysis demonstrated significant reductions in lipid synthesis rates and lipid accumulation in developing siliques and seedlings of cdk8, and conversely, elevated lipid contents in wild-type seed overexpressing CDK8. Transactivation assays show that CDK8 is necessary for maximal transactivation of the master seed oil activator WRINKLED1 (WRI1) by the seed maturation transcription factor ABSCISIC ACID INSENSITIVE3, supporting a direct regulatory role of CDK8 in oil synthesis. Thermophoretic studies show GEMINIVIRUS REP INTERACTING KINASE1, an activating kinase of KIN10 (a catalytic subunit of SUCROSE NON-FERMENTING1-RELATED KINASE1), physically interacts with CDK8, resulting in its phosphorylation and degradation in the presence of KIN10. This work defines a mechanism whereby, once activated, KIN10 downregulates WRI1 expression and suppresses lipid synthesis via promoting the degradation of CDK8. The KIN10-CDK8-dependent regulation of lipid synthesis described herein is additional to our previously reported KIN10-dependent phosphorylation and degradation of WRI1.

Engineering triacylglycerol accumulation in duckweed (Lemna japonica).

Duckweeds are amongst the fastest growing of higher plants, making them attractive high-biomass targets for biofuel feedstock production. Their fronds have high rates of fatty acid synthesis to meet the demand for new membranes, but triacylglycerols (TAG) only accumulate to very low levels. Here we report on the engineering of Lemna japonica for the synthesis and accumulation of TAG in its fronds. This was achieved by expression of an estradiol-inducible cyan fluorescent protein-Arabidopsis WRINKLED1 fusion protein (CFP-AtWRI1), strong constitutive expression of a mouse diacylglycerol:acyl-CoA acyltransferase2 (MmDGAT), and a sesame oleosin variant (SiOLE(*)). Individual expression of each gene increased TAG accumulation by 1- to 7-fold relative to controls, while expression of pairs of these genes increased TAG by 7- to 45-fold. In uninduced transgenics containing all three genes, TAG accumulation increased by 45-fold to 3.6% of dry weight (DW) without severely impacting growth, and by 108-fold to 8.7% of DW after incubation on medium containing 100 µm estradiol for 4 days. TAG accumulation was accompanied by an increase in total fatty acids of up to three-fold to approximately 15% of DW. Lipid droplets from fronds of all transgenic lines were visible by confocal microscopy of BODIPY-stained fronds. At a conservative 12 tonnes (dry matter) per acre and 10% (DW) TAG, duckweed could produce 350 gallons of oil/acre/year, approximately seven-fold the yield of soybean, and similar to that of oil palm. These findings provide the foundation for optimizing TAG accumulation in duckweed and present a new opportunity for producing biofuels and lipidic bioproducts.

An expanded role for the transcription factor WRINKLED1 in the biosynthesis of triacylglycerols during seed development.

The transcription factor WRINKLED1 (WRI1) is known as a master regulator of fatty acid synthesis in developing oilseeds of Arabidopsis thaliana and other species. WRI1 is known to directly stimulate the expression of many fatty acid biosynthetic enzymes and a few targets in the lower part of the glycolytic pathway. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors is controlled by WRI1. To shortlist possible gene targets for future in-planta experimental validation, here we present a strategy that combines phylogenetic foot printing of cis-regulatory elements with additional layers of evidence. Upstream regions of protein-encoding genes in A. thaliana were searched for the previously described DNA-binding consensus for WRI1, the ASML1/WRI1 (AW)-box. For about 900 genes, AW-box sites were found to be conserved across orthologous upstream regions in 11 related species of the crucifer family. For 145 select potential target genes identified this way, affinity of upstream AW-box sequences to WRI1 was assayed by Microscale Thermophoresis. This allowed definition of a refined WRI1 DNA-binding consensus. We find that known WRI1 gene targets are predictable with good confidence when upstream AW-sites are phylogenetically conserved, specifically binding WRI1 in the in vitro assay, positioned in proximity to the transcriptional start site, and if the gene is co-expressed with WRI1 during seed development. When targets predicted in this way are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in glycolysis and fatty acid biosynthesis by WRI1. Several of the newly predicted targets are in the upper glycolysis pathway and the pentose phosphate pathway. Of these, plastidic isoforms of fructokinase (FRK3) and of phosphoglucose isomerase (PGI1) are particularly corroborated by previously reported seed phenotypes of respective null mutations.

Purple acid phosphatase2 stimulates a futile cycle of lipid synthesis and degradation, and mitigates the negative growth effects of triacylglycerol accumulation in vegetative tissues.

Storage lipids (mostly triacylglycerols, TAGs) serve as an important energy and carbon reserve in plants, and hyperaccumulation of TAG in vegetative tissues can have negative effects on plant growth. Purple acid phosphatase2 (PAP2) was previously shown to affect carbon metabolism and boost plant growth. However, the effects of PAP2 on lipid metabolism remain unknown. Here, we demonstrated that PAP2 can stimulate a futile cycle of fatty acid (FA) synthesis and degradation, and mitigate negative growth effects associated with high accumulation of TAG in vegetative tissues. Constitutive expression of PAP2 in Arabidopsis thaliana enhanced both lipid synthesis and degradation in leaves and led to a substantial increase in seed oil yield. Suppressing lipid degradation in a PAP2-overexpressing line by disrupting sugar-dependent1 (SDP1), a predominant TAG lipase, significantly elevated vegetative TAG content and improved plant growth. Diverting FAs from membrane lipids to TAGs in PAP2-overexpressing plants by constitutively expressing phospholipid:diacylglycerol acyltransferase1 (PDAT1) greatly increased TAG content in vegetative tissues without compromising biomass yield. These results highlight the potential of combining PAP2 with TAG-promoting factors to enhance carbon assimilation, FA synthesis and allocation to TAGs for optimized plant growth and storage lipid accumulation in vegetative tissues.

Elucidation of Triacylglycerol Overproduction in the C4 Bioenergy Crop Sorghum bicolor by Constraint-Based Analysis.

Upregulation of triacylglycerols (TAGs) in vegetative plant tissues such as leaves has the potential to drastically increase the energy density and biomass yield of bioenergy crops. In this context, constraint-based analysis has the promise to improve metabolic engineering strategies. Here we present a core metabolism model for the C4 biomass crop Sorghum bicolor (iTJC1414) along with a minimal model for photosynthetic CO2 assimilation, sucrose and TAG biosynthesis in C3 plants. Extending iTJC1414 to a four-cell diel model we simulate C4 photosynthesis in mature leaves with the principal photo-assimilatory product being replaced by TAG produced at different levels. Independent of specific pathways and per unit carbon assimilated, energy content and biosynthetic demands in reducing equivalents are about 1.3 to 1.4 times higher for TAG than for sucrose. For plant generic pathways, ATP- and NADPH-demands per CO2 assimilated are higher by 1.3- and 1.5-fold, respectively. If the photosynthetic supply in ATP and NADPH in iTJC1414 is adjusted to be balanced for sucrose as the sole photo-assimilatory product, overproduction of TAG is predicted to cause a substantial surplus in photosynthetic ATP. This means that if TAG synthesis was the sole photo-assimilatory process, there could be an energy imbalance that might impede the process. Adjusting iTJC1414 to a photo-assimilatory rate that approximates field conditions, we predict possible daily rates of TAG accumulation, dependent on varying ratios of carbon partitioning between exported assimilates and accumulated oil droplets (TAG, oleosin) and in dependence of activation of futile cycles of TAG synthesis and degradation. We find that, based on the capacity of leaves for photosynthetic synthesis of exported assimilates, mature leaves should be able to reach a 20% level of TAG per dry weight within one month if only 5% of the photosynthetic net assimilation can be allocated into oil droplets. From this we conclude that high TAG levels should be achievable if TAG synthesis is induced only during a final phase of the plant life cycle.

Cellular Plasticity in Response to Suppression of Storage Proteins in the Brassica napus Embryo.

The tradeoff between protein and oil storage in oilseed crops has been tested here in oilseed rape (Brassica napus) by analyzing the effect of suppressing key genes encoding protein storage products (napin and cruciferin). The phenotypic outcomes were assessed using NMR and mass spectrometry imaging, microscopy, transcriptomics, proteomics, metabolomics, lipidomics, immunological assays, and flux balance analysis. Surprisingly, the profile of storage products was only moderately changed in RNA interference transgenics. However, embryonic cells had undergone remarkable architectural rearrangements. The suppression of storage proteins led to the elaboration of membrane stacks enriched with oleosin (sixfold higher protein abundance) and novel endoplasmic reticulum morphology. Protein rebalancing and amino acid metabolism were focal points of the metabolic adjustments to maintain embryonic carbon/nitrogen homeostasis. Flux balance analysis indicated a rather minor additional demand for cofactors (ATP and NADPH). Thus, cellular plasticity in seeds protects against perturbations to its storage capabilities and, hence, contributes materially to homeostasis. This study provides mechanistic insights into the intriguing link between lipid and protein storage, which have implications for biotechnological strategies directed at improving oilseed crops.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models.

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For 13C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

WRINKLED1 Regulates BIOTIN ATTACHMENT DOMAIN-CONTAINING Proteins that Inhibit Fatty Acid Synthesis.

WRINKLED1 (WRI1) is a transcriptional activator that binds to a conserved sequence (designated as AW box) boxes in the promoters of many genes from central metabolism and fatty acid (FA) synthesis, resulting in their transcription. BIOTIN ATTACHMENT DOMAIN-CONTAINING (BADC) proteins lack a biotin-attachment domain and are therefore inactive, but in the presence of excess FA, BADC1 and BADC3 are primarily responsible for the observed long-term irreversible inhibition of ACETYL-COA CARBOXYLASE, and consequently FA synthesis. Here, we tested the interaction of WRI1 with BADC genes in Arabidopsis (Arabidopsis thaliana) and found purified WRI1 bound with high affinity to canonical AW boxes from the promoters of all three BADC genes. Consistent with this observation, both expression of BADC1, BADC2, and BADC3 genes and BADC1 protein levels were reduced in wri1-1 relative to the wild type, and elevated upon WRI1 overexpression. The double mutant badc1 badc2 phenocopied wri1-1 with respect to both reduction in root length and elevation of indole-3-acetic acid-Asp levels relative to the wild type. Overexpression of BADC1 in wri1-1 decreased indole-3-acetic acid-Asp content and partially rescued its short-root phenotype, demonstrating a role for BADCs in seedling establishment. That WRI1 positively regulates genes encoding both FA synthesis and BADC proteins (i.e. conditional inhibitors of FA synthesis), represents a coordinated mechanism to achieve lipid homeostasis in which plants couple the transcription of their FA synthetic capacity with their capacity to biochemically downregulate it.

Expression of a Lychee PHOSPHATIDYLCHOLINE:DIACYLGLYCEROL CHOLINEPHOSPHOTRANSFERASE with an Escherichia coli CYCLOPROPANE SYNTHASE Enhances Cyclopropane Fatty Acid Accumulation in Camelina Seeds.

Cyclopropane fatty acids (CPAs) are useful feedstocks for biofuels and bioproducts such as lubricants and biodiesel. Our goal is to identify factors that can facilitate the accumulation of CPA in seed triacylglycerol (TAG) storage oil. We hypothesized that the poor metabolism of CPA through the TAG biosynthetic network could be overcome by the addition of enzymes from species that naturally accumulate CPA in their seed oil, such as lychee (Litchi chinensis), which contains approximately 40% CPA in TAG. Our previous work on engineering CPA accumulation in crop and model plants identified a metabolic bottleneck between phosphatidylcholine (PC), the site of CPA biosynthesis, diacylglycerol (DAG), and TAG. Here, we report the cloning and heterologous expression in camelina (Camelina sativa) of a lychee PHOSPHATIDYLCHOLINE:DIACYLGLYCEROL CHOLINEPHOSPHOTRANSFERASE (PDCT), which encodes the enzyme that catalyzes the transfer of the phosphocholine headgroup from PC to DAG. Camelina lines coexpressing LcPDCT and Escherichia coli CYCLOPROPANE SYNTHASE (EcCPS) showed up to a 50% increase of CPA in mature seed, relative to the EcCPS background. Stereospecific lipid compositional analysis showed that the expression of LcPDCT strongly reduced the level of C18:1 substrate at PC-sn-1 and PC-sn-2 (i.e. the sites of CPA synthesis), while the levels of CPA increased in PC-sn-2, DAG-sn-1 and DAG-sn-2, and both sn-1/3 and sn-2 positions in TAG. Taken together, these data suggest that the addition of PDCT facilitates more efficient movement of CPA from PC to DAG and establishes LcPDCT as a useful factor to combine with others to enhance CPA accumulation in plant seed oil.

Identification of bottlenecks in the accumulation of cyclic fatty acids in camelina seed oil.

Modified fatty acids (mFA) have diverse uses; for example, cyclopropane fatty acids (CPA) are feedstocks for producing coatings, lubricants, plastics and cosmetics. The expression of mFA-producing enzymes in crop and model plants generally results in lower levels of mFA accumulation than in their natural-occurring source plants. Thus, to further our understanding of metabolic bottlenecks that limit mFA accumulation, we generated transgenic Camelina sativa lines co-expressing Escherichia coli cyclopropane synthase (EcCPS) and Sterculia foetida lysophosphatidic acid acyltransferase (SfLPAT). In contrast to transgenic CPA-accumulating Arabidopsis, CPA accumulation in camelina caused only minor changes in seed weight, germination rate, oil accumulation and seedling development. CPA accumulated to much higher levels in membrane than storage lipids, comprising more than 60% of total fatty acid in both phosphatidylcholine (PC) and phosphatidylethanolamine (PE) versus 26% in diacylglycerol (DAG) and 12% in triacylglycerol (TAG) indicating bottlenecks in the transfer of CPA from PC to DAG and from DAG to TAG. Upon co-expression of SfLPAT with EcCPS, di-CPA-PC increased by ~50% relative to lines expressing EcCPS alone with the di-CPA-PC primarily observed in the embryonic axis and mono-CPA-PC primarily in cotyledon tissue. EcCPS-SfLPAT lines revealed a redistribution of CPA from the sn-1 to sn-2 positions within PC and PE that was associated with a doubling of CPA accumulation in both DAG and TAG. The identification of metabolic bottlenecks in acyl transfer between site of synthesis (phospholipids) and deposition in storage oils (TAGs) lays the foundation for the optimizing CPA accumulation through directed engineering of oil synthesis in target crops.

Mathematical models of plant metabolism.

Among various modeling approaches in plant metabolic research, applications of Constraint-Based modeling are fast increasing in recent years, apparently driven by current advances in genomics and genome sequencing. Constraint-Based modeling, the functional analysis of metabolic networks at the whole cell or genome scale, is more difficult to apply to plants than to microbes. Here we discuss recent developments in Constraint-Based modeling in plants with focus on issues of model reconstruction and flux prediction. Another topic is the emerging application of integration of Constraint-Based modeling with omics data to increase predictive power. Furthermore, advances in experimental measurements of cellular fluxes by (13)C-Metabolic Flux Analysis are highlighted, including instationary (13)C-MFA used to probe autotrophic metabolism in photosynthetic tissue in the light.

Quantitative Multilevel Analysis of Central Metabolism in Developing Oilseeds of Oilseed Rape during in Vitro Culture.

Seeds provide the basis for many food, feed, and fuel products. Continued increases in seed yield, composition, and quality require an improved understanding of how the developing seed converts carbon and nitrogen supplies into storage. Current knowledge of this process is often based on the premise that transcriptional regulation directly translates via enzyme concentration into flux. In an attempt to highlight metabolic control, we explore genotypic differences in carbon partitioning for in vitro cultured developing embryos of oilseed rape (Brassica napus). We determined biomass composition as well as 79 net fluxes, the levels of 77 metabolites, and 26 enzyme activities with specific focus on central metabolism in nine selected germplasm accessions. Overall, we observed a tradeoff between the biomass component fractions of lipid and starch. With increasing lipid content over the spectrum of genotypes, plastidic fatty acid synthesis and glycolytic flux increased concomitantly, while glycolytic intermediates decreased. The lipid/starch tradeoff was not reflected at the proteome level, pointing to the significance of (posttranslational) metabolic control. Enzyme activity/flux and metabolite/flux correlations suggest that plastidic pyruvate kinase exerts flux control and that the lipid/starch tradeoff is most likely mediated by allosteric feedback regulation of phosphofructokinase and ADP-glucose pyrophosphorylase. Quantitative data were also used to calculate in vivo mass action ratios, reaction equilibria, and metabolite turnover times. Compounds like cyclic 3',5'-AMP and sucrose-6-phosphate were identified to potentially be involved in so far unknown mechanisms of metabolic control. This study provides a rich source of quantitative data for those studying central metabolism.

Transcript abundance on its own cannot be used to infer fluxes in central metabolism.

An attempt has been made to define the extent to which metabolic flux in central plant metabolism is reflected by changes in the transcriptome and metabolome, based on an analysis of in vitro cultured immature embryos of two oilseed rape (Brassica napus) accessions which contrast for seed lipid accumulation. Metabolic flux analysis (MFA) was used to constrain a flux balance metabolic model which included 671 biochemical and transport reactions within the central metabolism. This highly confident flux information was eventually used for comparative analysis of flux vs. transcript (metabolite). Metabolite profiling succeeded in identifying 79 intermediates within the central metabolism, some of which differed quantitatively between the two accessions and displayed a significant shift corresponding to flux. An RNA-Seq based transcriptome analysis revealed a large number of genes which were differentially transcribed in the two accessions, including some enzymes/proteins active in major metabolic pathways. With a few exceptions, differential activity in the major pathways (glycolysis, TCA cycle, amino acid, and fatty acid synthesis) was not reflected in contrasting abundances of the relevant transcripts. The conclusion was that transcript abundance on its own cannot be used to infer metabolic activity/fluxes in central plant metabolism. This limitation needs to be borne in mind in evaluating transcriptome data and designing metabolic engineering experiments.

FAD2 and FAD3 desaturases form heterodimers that facilitate metabolic channeling in vivo.

Plant desaturases comprise two independently evolved classes, a structurally well characterized soluble class responsible for the production of monoenes in the plastids of higher plants and the poorly structurally characterized integral membrane class that has members in the plastid and endoplasmic reticulum that are responsible for producing mono- and polyunsaturated fatty acids. Both require iron and oxygen for activity and are inhibited by azide and cyanide underscoring their common chemical imperatives. We previously showed that the Δ(9) acyl-CoA integral membrane desaturase Ole1p from Saccharomyces cerevisiae exhibits dimeric organization, like the soluble plastidial acyl-ACP desaturases. Here we use two independent bimolecular complementation assays, i.e. yeast two-hybrid analysis and Arabidopsis leaf protoplast split luciferase assay, to demonstrate that members of the plant integral membrane fatty acid desaturase (FAD) family, FAD2, FAD3, FAD6, FAD7, and FAD8, self-associate. Further, the endoplasmic reticulum-localized desaturase FAD2 can associate with FAD3, as can the plastid-localized FAD6 desaturase with either FAD7 or FAD8. These pairings appear to be specific because pairs such as FAD3 and FAD7 (or FAD8) and FAD2 and FAD6 do not interact despite their high amino acid similarity. These results are consistent also with their known endoplasmic reticulum and plastid subcellular localizations. Chemical cross-linking experiments confirm that FAD2 and FAD3 can form dimers like the yeast Ole1p and, when coexpressed, can form FAD2-FAD3 heterodimers. Metabolic flux analysis of yeast coexpressing FAD2 and FAD3 indicates that heterodimers can form a metabolic channel in which 18:1-PC is converted to 18:3-PC without releasing a free 18:2-PC intermediate.

Flux variability analysis: application to developing oilseed rape embryos using toolboxes for constraint-based modeling.

Flux variability analysis enables comprehensive exploration of alternate optimal routes in a metabolic network. This method is especially useful with models such as bna572 for the developing oilseed rape embryo which is highly compartmentalized. Here, we describe a protocol for carrying out flux variability analysis on reactions and network projections of bna572 using well-established software (CellNetAnalyzer and COBRA) for constraint-based analysis of stoichiometric network reconstructions.

Integration of a constraint-based metabolic model of Brassica napus developing seeds with (13)C-metabolic flux analysis.

The use of large-scale or genome-scale metabolic reconstructions for modeling and simulation of plant metabolism and integration of those models with large-scale omics and experimental flux data is becoming increasingly important in plant metabolic research. Here we report an updated version of bna572, a bottom-up reconstruction of oilseed rape (Brassica napus L.; Brassicaceae) developing seeds with emphasis on representation of biomass-component biosynthesis. New features include additional seed-relevant pathways for isoprenoid, sterol, phenylpropanoid, flavonoid, and choline biosynthesis. Being now based on standardized data formats and procedures for model reconstruction, bna572+ is available as a COBRA-compliant Systems Biology Markup Language (SBML) model and conforms to the Minimum Information Requested in the Annotation of Biochemical Models (MIRIAM) standards for annotation of external data resources. Bna572+ contains 966 genes, 671 reactions, and 666 metabolites distributed among 11 subcellular compartments. It is referenced to the Arabidopsis thaliana genome, with gene-protein-reaction (GPR) associations resolving subcellular localization. Detailed mass and charge balancing and confidence scoring were applied to all reactions. Using B. napus seed specific transcriptome data, expression was verified for 78% of bna572+ genes and 97% of reactions. Alongside bna572+ we also present a revised carbon centric model for (13)C-Metabolic Flux Analysis ((13)C-MFA) with all its reactions being referenced to bna572+ based on linear projections. By integration of flux ratio constraints obtained from (13)C-MFA and by elimination of infinite flux bounds around thermodynamically infeasible loops based on COBRA loopless methods, we demonstrate improvements in predictive power of Flux Variability Analysis (FVA). Using this combined approach we characterize the difference in metabolic flux of developing seeds of two B. napus genotypes contrasting in starch and oil content.

Survey of the total fatty acid and triacylglycerol composition and content of 30 duckweed species and cloning of a Δ6-desaturase responsible for the production of γ-linolenic and stearidonic acids in Lemna gibba.

BACKGROUND: Duckweeds, i.e., members of the Lemnoideae family, are amongst the smallest aquatic flowering plants. Their high growth rate, aquatic habit and suitability for bio-remediation make them strong candidates for biomass production. Duckweeds have been studied for their potential as feedstocks for bioethanol production; however, less is known about their ability to accumulate reduced carbon as fatty acids (FA) and oil. RESULTS: Total FA profiles of thirty duckweed species were analysed to assess the natural diversity within the Lemnoideae. Total FA content varied between 4.6% and 14.2% of dry weight whereas triacylglycerol (TAG) levels varied between 0.02% and 0.15% of dry weight. Three FA, 16:0 (palmitic), 18:2Δ9,12 (Linoleic acid, or LN) and 18:3Δ9,12,15 (α-linolenic acid, or ALA) comprise more than 80% of total duckweed FA. Seven Lemna and two Wolffiela species also accumulate polyunsaturated FA containing Δ6-double bonds, i.e., GLA and SDA. Relative to total FA, TAG is enriched in saturated FA and deficient in polyunsaturated FA, and only five Lemna species accumulate Δ6-FA in their TAG. A putative Δ6-desaturase designated LgDes, with homology to a family of front-end Δ6-FA and Δ8-spingolipid desaturases, was identified in the assembled DNA sequence of Lemna gibba. Expression of a synthetic LgDes gene in Nicotiana benthamiana resulted in the accumulation of GLA and SDA, confirming it specifies a Δ6-desaturase. CONCLUSIONS: Total accumulation of FA varies three-fold across the 30 species of Lemnoideae surveyed. Nine species contain GLA and SDA which are synthesized by a Δ6 front-end desaturase, but FA composition is otherwise similar. TAG accumulates up to 0.15% of total dry weight, comparable to levels found in the leaves of terrestrial plants. Polyunsaturated FA is underrepresented in TAG, and the Δ6-FA GLA and SDA are found in the TAG of only five of the nine Lemna species that produce them. When present, GLA is enriched and SDA diminished relative to their abundance in the total FA pool.