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

A role for fermentation in aerobic conditions as revealed by computational analysis of maize root metabolism during growth by cell elongation

The root is a well-studied example of cell specialisation, yet little is known about the metabolism that supports the transport functions and growth of different root cell types. To address this, we used computational modelling to study metabolism in the elongation zone of a maize lateral root. A functional-structural model captured the cell-anatomical features of the root and modelled how they changed as the root elongated. From these data, we derived constraints for a flux balance analysis model that predicted metabolic fluxes of the 11 concentric rings of cells in the root. We discovered a distinct metabolic flux pattern in the cortical cell rings, endodermis and pericycle (but absent in the epidermis) that involved a high rate of glycolysis and production of the fermentation end-products lactate and ethanol. This aerobic fermentation was confirmed experimentally by metabolite analysis. The use of fermentation in the model was not obligatory but was the most efficient way to meet the specific demands for energy, reducing power and carbon skeletons of expanding cells. Cytosolic acidification was avoided in the fermentative mode due to the substantial consumption of protons by lipid synthesis. These results expand our understanding of fermentative metabolism beyond that of hypoxic niches and suggest that fermentation could play an important role in the metabolism of aerobic tissues.

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.

Transcriptome-based identification and functional characterization of iridoid synthase involved in monotropein biosynthesis in blueberry

Blueberries (Vaccinium spp.) are well known for their nutritional quality, and recent work has shown that Vaccinium spp. also produce iridoids, which are specialized metabolites with potent health-promoting benefits. The iridoid glycoside monotropein, which has anti-inflammatory and antinociceptive activities, has been detected in several wild blueberry species but in only a few cultivated highbush blueberry cultivars. How monotropein is produced in blueberry and the genes involved in its biosynthesis remain to be elucidated. Using a monotropein-positive (M+) and monotropein-negative (M-) cultivar of blueberry, we employed transcriptomics and comparative genomics to identify candidate genes in the blueberry iridoid biosynthetic pathway. Orthology analysis was completed using de novo transcript assemblies for both the M+ and M- blueberry cultivars along with the known iridoid-producing plant species Catharanthus roseus to identify putative genes involved in key steps in the early iridoid biosynthetic pathway. From the identified orthologs, we functionally characterized iridoid synthase (ISY), a key enzyme involved in formation of the iridoid scaffold, from both the M+ and M- cultivars. Detection of nepetalactol suggests that ISY from both the M+ and M- cultivars produce functional enzymes that catalyze the formation of iridoids. Transcript accumulation of the putative ISY gene did not correlate with monotropein production, suggesting other genes in the monotropein biosynthetic pathway may be more directly responsible for differential accumulation of the metabolite in blueberry. Mutual rank analysis revealed that ISY is co-expressed with UDP-glucuronosyltransferase, which encodes an enzyme downstream of the ISY step. Results from this study contribute new knowledge in our understanding of iridoid biosynthesis in blueberry and could lead to development of new cultivars with increased human health benefits.

Causes and consequences of endogenous hypoxia on growth and metabolism of developing maize kernels

Maize (Zea mays) kernels are the largest cereal grains, and their endosperm is severely oxygen deficient during grain fill. The causes, dynamics, and mechanisms of acclimation to hypoxia are minimally understood. Here, we demonstrate that hypoxia develops in the small, growing endosperm, but not the nucellus, and becomes the standard state, regardless of diverse structural and genetic perturbations in modern maize (B73, popcorn, sweet corn), mutants (sweet4c, glossy6, waxy), and non-domesticated wild relatives (teosintes and Tripsacum species). We also uncovered an interconnected void space at the chalazal pericarp, providing superior oxygen supply to the placental tissues and basal endosperm transfer layer. Modeling indicated a very high diffusion resistance inside the endosperm, which, together with internal oxygen consumption, could generate steep oxygen gradients at the endosperm surface. Manipulation of oxygen supply induced reciprocal shifts in gene expression implicated in controlling mitochondrial functions (23.6 kDa Heat-Shock Protein, Voltage-Dependent Anion Channel 2) and multiple signaling pathways (core hypoxia genes, cyclic nucleotide metabolism, ethylene synthesis). Metabolite profiling revealed oxygen-dependent shifts in mitochondrial pathways, ascorbate metabolism, starch synthesis, and auxin degradation. Long-term elevated oxygen supply enhanced the rate of kernel development. Altogether, evidence here supports a mechanistic framework for the establishment of and acclimation to hypoxia in the maize endosperm.

Polyunsaturated Fatty Acids: Conversion to Lipid Mediators, Roles in Inflammatory Diseases and Dietary Sources

Polyunsaturated fatty acids (PUFAs) are important components of the diet of mammals. Their role was first established when the essential fatty acids (EFAs) linoleic acid and α-linolenic acid were discovered nearly a century ago. However, most of the biochemical and physiological actions of PUFAs rely on their conversion to 20C or 22C acids and subsequent metabolism to lipid mediators. As a generalisation, lipid mediators formed from n-6 PUFAs are pro-inflammatory while those from n-3 PUFAs are anti-inflammatory or neutral. Apart from the actions of the classic eicosanoids or docosanoids, many newly discovered compounds are described as Specialised Pro-resolving Mediators (SPMs) which have been proposed to have a role in resolving inflammatory conditions such as infections and preventing them from becoming chronic. In addition, a large group of molecules, termed isoprostanes, can be generated by free radical reactions and these too have powerful properties towards inflammation. The ultimate source of n-3 and n-6 PUFAs are photosynthetic organisms which contain Δ-12 and Δ-15 desaturases, which are almost exclusively absent from animals. Moreover, the EFAs consumed from plant food are in competition with each other for conversion to lipid mediators. Thus, the relative amounts of n-3 and n-6 PUFAs in the diet are important. Furthermore, the conversion of the EFAs to 20C and 22C PUFAs in mammals is rather poor. Thus, there has been much interest recently in the use of algae, many of which make substantial quantities of long-chain PUFAs or in manipulating oil crops to make such acids. This is especially important because fish oils, which are their main source in human diets, are becoming limited. In this review, the metabolic conversion of PUFAs into different lipid mediators is described. Then, the biological roles and molecular mechanisms of such mediators in inflammatory diseases are outlined. Finally, natural sources of PUFAs (including 20 or 22 carbon compounds) are detailed, as well as recent efforts to increase their production.

Ectopic expression of the AtCDF1 transcription factor in potato enhances tuber starch and amino acid contents and yield under open field conditions

INTRODUCTION

Cycling Dof transcription factors (CDFs) have been involved in different aspects of plant growth and development. In Arabidopsis and tomato, one member of this family (CDF1) has recently been associated with the regulation of primary metabolism and abiotic stress responses, but their roles in crop production under open field conditions remain unknown.

METHODS

In this study, we compared the growth, and tuber yield and composition of plants ectopically expressing the CDF1 gene from Arabidopsis under the control of the 35S promoter with wild-type (WT) potato plants cultured in growth chamber and open field conditions.

RESULTS

In growth chambers, the 35S::AtCDF1 plants showed a greater tuber yield than the WT by increasing the biomass partition for tuber development. Under field conditions, the ectopic expression of CDF1 also promoted the sink strength of the tubers, since 35S::AtCDF1 plants exhibited significant increases in tuber size and weight resulting in higher tuber yield. A metabolomic analysis revealed that tubers of 35S::AtCDF1 plants cultured under open field conditions accumulated higher levels of glucose, starch and amino acids than WT tubers. A comparative proteomic analysis of tubers of 35S::AtCDF1 and WT plants cultured under open field conditions revealed that these changes can be accounted for changes in the expression of proteins involved in energy production and different aspects of C and N metabolism.

DISCUSSION

The results from this study advance our collective understanding of the role of CDFs and are of great interest for the purposes of improving the yield and breeding of crop plants.

Evidence that mitochondrial alternative oxidase respiration supports carbon balance in source leaves of Nicotiana tabacum

Alternative oxidase (AOX) represents a non-energy conserving pathway within the mitochondrial electron transport chain. One potential physiological role of AOX could be to manage leaf carbohydrate amounts by supporting respiratory carbon oxidation reactions. In this study, several approaches tested the hypothesis that AOX1a gene expression in Nicotiana tabacum leaf is enhanced in conditions expected to promote an increased leaf carbohydrate status. These approaches included supplying leaves with exogenous carbohydrates, comparing plants grown at different atmospheric CO₂ concentrations, comparing sink leaves with source leaves, comparing plants with different ratios of source to sink activity, and examining gene expression over the diel cycle. In each case, the pattern of AOX1a gene expression was compared with that of other genes known to respond to carbohydrates and/or other factors related to source:sink activity. These included GPT1 and GPT3 (that encode chloroplast glucose 6-phosphate/phosphate translocators), SPS (that encodes sucrose phosphate synthase), SUT1 (that encodes a sucrose/H⁺ symporter involved in phloem loading) and UCP1 (that encodes a mitochondrial uncoupling protein). The AOX1a transcript amount was higher following the leaf sink-to-source transition, and in plants with higher source relative to sink activity due to increasing plant age. Further, these effects were amplified in plants grown at elevated CO₂ to stimulate source activity, particularly at end-of-day time periods. The AOX1a transcript amount was also higher following treatment of leaves with carbohydrate, in particular sucrose. Overall, the results provide evidence that, while source leaf sucrose accumulation may signal for a down-regulation of sucrose synthesis and transport, it also signals for means to manage the excess cytosolic carbohydrate pools. This includes increased AOX respiration to support carbon oxidation pathways even if energy charge is high, in combination perhaps with some return flux of carbohydrate from cytosol to stroma through the GPT3 translocator. As discussed, these activities could contribute to maintaining plant source:sink balance, as well as photosynthetic and phloem loading capacity.

Glucose-6-P/phosphate translocator2 mediates the phosphoglucose-isomerase1-independent response to microbial volatiles

In Arabidopsis (Arabidopsis thaliana), the plastidial isoform of phosphoglucose isomerase (PGI1) mediates photosynthesis, metabolism, and development, probably due to its involvement in the synthesis of isoprenoid-derived signals in vascular tissues. Microbial volatile compounds (VCs) with molecular masses of <45 Da promote photosynthesis, growth, and starch overaccumulation in leaves through PGI1-independent mechanisms. Exposure to these compounds in leaves enhances the levels of GLUCOSE-6-PHOSPHATE/PHOSPHATE TRANSLOCATOR2 (GPT2) transcripts. We hypothesized that the PGI1-independent response to microbial volatile emissions involves GPT2 action. To test this hypothesis, we characterized the responses of wild-type (WT), GPT2-null gpt2-1, PGI1-null pgi1-2, and pgi1-2gpt2-1 plants to small fungal VCs. In addition, we characterized the responses of pgi1-2gpt2-1 plants expressing GPT2 under the control of a vascular tissue- and root tip-specific promoter to small fungal VCs. Fungal VCs promoted increases in growth, starch content, and photosynthesis in WT and gpt2-1 plants. These changes were substantially weaker in VC-exposed pgi1-2gpt2-1 plants but reverted to WT levels with vascular and root tip-specific GPT2 expression. Proteomic analyses did not detect enhanced levels of GPT2 protein in VC-exposed leaves and showed that knocking out GPT2 reduced the expression of photosynthesis-related proteins in pgi1-2 plants. Histochemical analyses of GUS activity in plants expressing GPT2-GUS under the control of the GPT2 promoter showed that GPT2 is mainly expressed in root tips and vascular tissues around hydathodes. Overall, the data indicated that the PGI1-independent response to microbial VCs involves resetting of the photosynthesis-related proteome in leaves through long-distance GPT2 action.

Integrated transcriptome and proteome analysis of developing embryo reveals the mechanisms underlying the high levels of oil accumulation in Carya cathayensis Sarg

Hickory (Carya cathayensis Sarg.) is an extraordinary nut-bearing deciduous arbor with high content of oil in its embryo. However, the molecular mechanism underlying high oil accumulation is mostly unknown. Here, we reported that the lipid droplets and oil accumulation gradually increased with the embryo development and the oil content was up to ~76% at maturity. Furthermore, transcriptome and proteome analysis of developing hickory embryo identified 32,907 genes and 9857 proteins. Time-series analysis of gene expressions showed that these genes were divided into 12 clusters and lipid metabolism-related genes were enriched in Cluster 3, with the highest expression levels at 95 days after pollination (S2). Differentially expressed genes and proteins indicated high correlation, and both were enriched in the lipid metabolism. Notably, the genes involved in biosynthesis, transport of fatty acid/lipid and lipid droplets formation had high expression levels at S2, while the expression levels of other genes required for suberin/wax/cutin biosynthesis and lipid degradation were very low at all the sampling time points, ultimately promoting the accumulation of oil. Quantitative reverse-transcription PCR analysis also verified the results of RNA-seq. The co-regulatory networks of lipid metabolism were further constructed and WRINKLED1 (WRI1) was a core transcriptional factor located in the nucleus. Of note, CcWRI1A/B could directly activate the expression of some genes (CcBCCP2A, CcBCCP2B, CcFATA and CcFAD3) required for fatty acid synthesis. These results provided in-depth evidence for revealing the molecular mechanism of high oil accumulation in hickory embryo.

Transcriptome Analysis and Identification of Lipid Genes in Physaria lindheimeri, a Genetic Resource for Hydroxy Fatty Acids in Seed Oil

Hydroxy fatty acids (HFAs) have numerous industrial applications but are absent in most vegetable oils. Physaria lindheimeri accumulating 85% HFA in its seed oil makes it a valuable resource for engineering oilseed crops for HFA production. To discover lipid genes involved in HFA synthesis in P. lindheimeri, transcripts from developing seeds at various stages, as well as leaf and flower buds, were sequenced. Ninety-seven percent clean reads from 552,614,582 raw reads were assembled to 129,633 contigs (or transcripts) which represented 85,948 unique genes. Gene Ontology analysis indicated that 60% of the contigs matched proteins involved in biological process, cellular component or molecular function, while the remaining matched unknown proteins. We identified 42 P. lindheimeri genes involved in fatty acid and seed oil biosynthesis, and 39 of them shared 78-100% nucleotide identity with Arabidopsis orthologs. We manually annotated 16 key genes and 14 of them contained full-length protein sequences, indicating high coverage of clean reads to the assembled contigs. A detailed profiling of the 16 genes revealed various spatial and temporal expression patterns. The further comparison of their protein sequences uncovered amino acids conserved among HFA-producing species, but these varied among non-HFA-producing species. Our findings provide essential information for basic and applied research on HFA biosynthesis.

Toxicant-Induced Metabolic Alterations in Lipid and Amino Acid Pathways Are Predictive of Acute Liver Toxicity in Rats

Liver disease and disorders associated with aberrant hepatocyte metabolism can be initiated via drug and environmental toxicant exposures. In this study, we tested the hypothesis that gene and metabolic profiling can reveal commonalities in liver response to different toxicants and provide the capability to identify early signatures of acute liver toxicity. We used Sprague Dawley rats and three classical hepatotoxicants: acetaminophen (2 g/kg), bromobenzene (0.4 g/kg), and carbon tetrachloride (0.3 g/kg), to identify early perturbations in liver metabolism after a single acute exposure dose. We measured changes in liver genes and plasma metabolites at two time points (5 and 10 h) and used genome-scale metabolic models to identify commonalities in liver responses across the three toxicants. We found strong correlations for gene and metabolic profiles between the toxicants, indicative of similarities in the liver response to toxicity. We identified several injury-specific pathways in lipid and amino acid metabolism that changed similarly across the three toxicants. Our findings suggest that several plasma metabolites in lipid and amino acid metabolism are strongly associated with the progression of liver toxicity, and as such, could be targeted and clinically assessed for their potential as early predictors of acute liver toxicity.

Reorganization of Acyl Flux through the Lipid Metabolic Network in Oil-Accumulating Tobacco Leaves

The triacylglycerols (TAGs; i.e. oils) that accumulate in plants represent the most energy-dense form of biological carbon storage, and are used for food, fuels, and chemicals. The increasing human population and decreasing amount of arable land have amplified the need to produce plant oil more efficiently. Engineering plants to accumulate oils in vegetative tissues is a novel strategy, because most plants only accumulate large amounts of lipids in the seeds. Recently, tobacco (Nicotiana tabacum) leaves were engineered to accumulate oil at 15% of dry weight due to a push (increased fatty acid synthesis)-and-pull (increased final step of TAG biosynthesis) engineering strategy. However, to accumulate both TAG and essential membrane lipids, fatty acid flux through nonengineered reactions of the endogenous metabolic network must also adapt, which is not evident from total oil analysis. To increase our understanding of endogenous leaf lipid metabolism and its ability to adapt to metabolic engineering, we utilized a series of in vitro and in vivo experiments to characterize the path of acyl flux in wild-type and transgenic oil-accumulating tobacco leaves. Acyl flux around the phosphatidylcholine acyl editing cycle was the largest acyl flux reaction in wild-type and engineered tobacco leaves. In oil-accumulating leaves, acyl flux into the eukaryotic pathway of glycerolipid assembly was enhanced at the expense of the prokaryotic pathway. However, a direct Kennedy pathway of TAG biosynthesis was not detected, as acyl flux through phosphatidylcholine preceded the incorporation into TAG. These results provide insight into the plasticity and control of acyl lipid metabolism in leaves.

Cannabis glandular trichomes alter morphology and metabolite content during flower maturation

The cannabis leaf is iconic, but it is the flowers of cannabis that are consumed for the psychoactive and medicinal effects of their specialized metabolites. Cannabinoid metabolites, together with terpenes, are produced in glandular trichomes. Superficially, stalked and sessile trichomes in cannabis only differ in size and whether they have a stalk. The objectives of this study were: to define each trichome type using patterns of autofluorescence and secretory cell numbers, to test the hypothesis that stalked trichomes develop from sessile-like precursors, and to test whether metabolic specialization occurs in cannabis glandular trichomes. A two-photon microscopy technique using glandular trichome intrinsic autofluorescence was developed which demonstrated that stalked glandular trichomes possessed blue autofluorescence correlated with high cannabinoid levels. These stalked trichomes had 12-16 secretory disc cells and strongly monoterpene-dominant terpene profiles. In contrast, sessile trichomes on mature flowers and vegetative leaves possessed red-shifted autofluorescence, eight secretory disc cells and less monoterpene-dominant terpene profiles. Moreover, intrinsic autofluorescence patterns and disc cell numbers supported a developmental model where stalked trichomes develop from apparently sessile trichomes. Transcriptomes of isolated floral trichomes revealed strong expression of cannabinoid and terpene biosynthetic genes, as well as uncharacterized genes highly co-expressed with CBDA synthase. Identification and characterization of two previously unknown and highly expressed monoterpene synthases highlighted the metabolic specialization of stalked trichomes for monoterpene production. These unique properties and highly expressed genes of cannabis trichomes determine the medicinal, psychoactive and sensory properties of cannabis products.

Common metabolic networks contribute to carbon sink strength of sorghum internodes: implications for bioenergy improvement

BACKGROUND

Sorghum bicolor (L.) is an important bioenergy source. The stems of sweet sorghum function as carbon sinks and accumulate large amounts of sugars and lignocellulosic biomass and considerable amounts of starch, therefore providing a model of carbon allocation and accumulation for other bioenergy crops. While omics data sets for sugar accumulation have been reported in different genotypes, the common features of primary metabolism in sweet genotypes remain unclear. To obtain a cohesive and comparative picture of carbohydrate metabolism between sorghum genotypes, we compared the phenotypes and transcriptome dynamics of sugar-accumulating internodes among three different sweet genotypes (Della, Rio, and SIL-05) and two non-sweet genotypes (BTx406 and R9188).

RESULTS

Field experiments showed that Della and Rio had similar dynamics and internode patterns of sugar concentration, albeit distinct other phenotypes. Interestingly, cellulose synthases for primary cell wall and key genes in starch synthesis and degradation were coordinately upregulated in sweet genotypes. Sweet sorghums maintained active monolignol biosynthesis compared to the non-sweet genotypes. Comparative RNA-seq results support the role of candidate Tonoplast Sugar Transporter gene (TST), but not the Sugars Will Eventually be Exported Transporter genes (SWEETs) in the different sugar accumulations between sweet and non-sweet genotypes.

CONCLUSIONS

Comparisons of the expression dynamics of carbon metabolic genes across the RNA-seq data sets identify several candidate genes with contrasting expression patterns between sweet and non-sweet sorghum lines, including genes required for cellulose and monolignol synthesis (CesA, PTAL, and CCR), starch metabolism (AGPase, SS, SBE, and G6P-translocator SbGPT2), and sucrose metabolism and transport (TPP and TST2). The common transcriptome features of primary metabolism identified here suggest the metabolic networks contributing to carbon sink strength in sorghum internodes, prioritize the candidate genes for manipulating carbon allocation with bioenergy purposes, and provide a comparative and cohesive picture of the complexity of carbon sink strength in sorghum stem.

Gene Networks Underlying Cannabinoid and Terpenoid Accumulation in Cannabis

Glandular trichomes are specialized anatomical structures that accumulate secretions with important biological roles in plant-environment interactions. These secretions also have commercial uses in the flavor, fragrance, and pharmaceutical industries. The capitate-stalked glandular trichomes of Cannabis sativa (cannabis), situated on the surfaces of the bracts of the female flowers, are the primary site for the biosynthesis and storage of resins rich in cannabinoids and terpenoids. In this study, we profiled nine commercial cannabis strains with purportedly different attributes, such as taste, color, smell, and genetic origin. Glandular trichomes were isolated from each of these strains, and cell type-specific transcriptome data sets were acquired. Cannabinoids and terpenoids were quantified in flower buds. Statistical analyses indicated that these data sets enable the high-resolution differentiation of strains by providing complementary information. Integrative analyses revealed a coexpression network of genes involved in the biosynthesis of both cannabinoids and terpenoids from imported precursors. Terpene synthase genes involved in the biosynthesis of the major monoterpenes and sesquiterpenes routinely assayed by cannabis testing laboratories were identified and functionally evaluated. In addition to cloning variants of previously characterized genes, specifically CsTPS14CT [(-)-limonene synthase] and CsTPS15CT (β-myrcene synthase), we functionally evaluated genes that encode enzymes with activities not previously described in cannabis, namely CsTPS18VF and CsTPS19BL (nerolidol/linalool synthases), CsTPS16CC (germacrene B synthase), and CsTPS20CT (hedycaryol synthase). This study lays the groundwork for developing a better understanding of the complex chemistry and biochemistry underlying resin accumulation across commercial cannabis strains.

Induced, Imprinted, and Primed Responses to Changing Environments: Does Metabolism Store and Process Information?

Metabolism is the system layer that determines growth by the rate of matter uptake and conversion into biomass. The scaffold of enzymatic reaction rates drives the metabolic network in a given physico-chemical environment. In response to the diverse environmental stresses, plants have evolved the capability of integrating macro- and micro-environmental events to be prepared, i.e., to be primed for upcoming environmental challenges. The hierarchical view on stress signaling, where metabolites are seen as final downstream products, has recently been complemented by findings that metabolites themselves function as stress signals. We present a systematic concept of metabolic responses that are induced by environmental stresses and persist in the plant system. Such metabolic imprints may prime metabolic responses of plants for subsequent environmental stresses. We describe response types with examples of biotic and abiotic environmental stresses and suggest that plants use metabolic imprints, the metabolic changes that last beyond recovery from stress events, and priming, the imprints that function to prepare for upcoming stresses, to integrate diverse environmental stress histories. As a consequence, even genetically identical plants should be studied and understood as phenotypically plastic organisms that continuously adjust their metabolic state in response to their individually experienced local environment. To explore the occurrence and to unravel functions of metabolic imprints, we encourage researchers to extend stress studies by including detailed metabolic and stress response monitoring into extended recovery phases.

Transcriptome and metabolome reveal distinct carbon allocation patterns during internode sugar accumulation in different sorghum genotypes

Sweet sorghum accumulates large amounts of soluble sugar in its stem. However, a system-based understanding of this carbohydrate allocation process is lacking. Here, we compared the dynamic transcriptome and metabolome between the conversion line R9188 and its two parents, sweet sorghum RIO and grain sorghum BTx406 that have contrasting sugar-accumulating phenotypes. We identified two features of sucrose metabolism, stable concentrations of sugar phosphates in RIO and opposite trend of trehalose-6-phosphate (T6P) between RIO vs R9188/BTx406. Integration of transcriptome and metabolome revealed R9188 is partially active in starch metabolism together with medium sucrose level, whereas sweet sorghum had the highest sucrose concentration and remained highly active in sucrose, starch, and cell wall metabolism post-anthesis. Similar expression pattern of genes involved in sucrose degradation decreased the pool of sugar phosphates for precursors of starch and cell wall synthesis in R9188 and BTx406. Differential T6P signal between RIO vs R9188/BTx406 is associated with introgression of T6P regulators from BTx406 into R9188, including C-group bZIP and trehalose 6-phosphate phosphatase (TPP). The inverted T6P signalling in R9188 appears to down-regulate sucrose and starch metabolism partly through transcriptome reprogramming, whereas introgressed metabolic genes could be related to reduced cell wall metabolism. Our results show that coordinated primary metabolic pathways lead to high sucrose demand and accumulation in sweet sorghum, providing us with targets for genetic improvements of carbohydrate allocation in bioenergy crops.

Identification of genes associated with ricinoleic acid accumulation in Hiptage benghalensis via transcriptome analysis

BACKGROUND

Ricinoleic acid is a high-value hydroxy fatty acid with broad industrial applications. Hiptage benghalensis seed oil contains a high amount of ricinoleic acid (~ 80%) and represents an emerging source of this unusual fatty acid. However, the mechanism of ricinoleic acid accumulation in H. benghalensis is yet to be explored at the molecular level, which hampers the exploration of its potential in ricinoleic acid production.

RESULTS

To explore the molecular mechanism of ricinoleic acid biosynthesis and regulation, H. benghalensis seeds were harvested at five developing stages (13, 16, 19, 22, and 25 days after pollination) for lipid analysis. The results revealed that the rapid accumulation of ricinoleic acid occurred at the early-mid-seed development stages (16-22 days after pollination). Subsequently, the gene transcription profiles of the developing seeds were characterized via a comprehensive transcriptome analysis with second-generation sequencing and single-molecule real-time sequencing. Differential expression patterns were identified in 12,555 transcripts, including 71 enzymes in lipid metabolic pathways, 246 putative transcription factors (TFs) and 124 long noncoding RNAs (lncRNAs). Twelve genes involved in diverse lipid metabolism pathways, including fatty acid biosynthesis and modification (hydroxylation), lipid traffic, triacylglycerol assembly, acyl editing and oil-body formation, displayed high expression levels and consistent expression patterns with ricinoleic acid accumulation in the developing seeds, suggesting their primary roles in ricinoleic acid production. Subsequent co-expression network analysis identified 57 TFs and 35 lncRNAs, which are putatively involved in the regulation of ricinoleic acid biosynthesis. The transcriptome data were further validated by analyzing the expression profiles of key enzyme-encoding genes, TFs and lncRNAs with quantitative real-time PCR. Finally, a network of genes associated with ricinoleic acid accumulation in H. benghalensis was established.

CONCLUSIONS

This study was the first step toward the understating of the molecular mechanisms of ricinoleic acid biosynthesis and oil accumulation in H. benghalensis seeds and identified a pool of novel genes regulating ricinoleic acid accumulation. The results set a foundation for developing H. benghalensis into a novel ricinoleic acid feedstock at the transcriptomic level and provided valuable candidate genes for improving ricinoleic acid production in other plants.

Identification of genes associated with ricinoleic acid accumulation in Hiptage benghalensis via transcriptome analysis

BACKGROUND

Ricinoleic acid is a high-value hydroxy fatty acid with broad industrial applications. Hiptage benghalensis seed oil contains a high amount of ricinoleic acid (~ 80%) and represents an emerging source of this unusual fatty acid. However, the mechanism of ricinoleic acid accumulation in H. benghalensis is yet to be explored at the molecular level, which hampers the exploration of its potential in ricinoleic acid production.

RESULTS

To explore the molecular mechanism of ricinoleic acid biosynthesis and regulation, H. benghalensis seeds were harvested at five developing stages (13, 16, 19, 22, and 25 days after pollination) for lipid analysis. The results revealed that the rapid accumulation of ricinoleic acid occurred at the early-mid-seed development stages (16-22 days after pollination). Subsequently, the gene transcription profiles of the developing seeds were characterized via a comprehensive transcriptome analysis with second-generation sequencing and single-molecule real-time sequencing. Differential expression patterns were identified in 12,555 transcripts, including 71 enzymes in lipid metabolic pathways, 246 putative transcription factors (TFs) and 124 long noncoding RNAs (lncRNAs). Twelve genes involved in diverse lipid metabolism pathways, including fatty acid biosynthesis and modification (hydroxylation), lipid traffic, triacylglycerol assembly, acyl editing and oil-body formation, displayed high expression levels and consistent expression patterns with ricinoleic acid accumulation in the developing seeds, suggesting their primary roles in ricinoleic acid production. Subsequent co-expression network analysis identified 57 TFs and 35 lncRNAs, which are putatively involved in the regulation of ricinoleic acid biosynthesis. The transcriptome data were further validated by analyzing the expression profiles of key enzyme-encoding genes, TFs and lncRNAs with quantitative real-time PCR. Finally, a network of genes associated with ricinoleic acid accumulation in H. benghalensis was established.

CONCLUSIONS

This study was the first step toward the understating of the molecular mechanisms of ricinoleic acid biosynthesis and oil accumulation in H. benghalensis seeds and identified a pool of novel genes regulating ricinoleic acid accumulation. The results set a foundation for developing H. benghalensis into a novel ricinoleic acid feedstock at the transcriptomic level and provided valuable candidate genes for improving ricinoleic acid production in other plants.

Carbon Orientation in the Diatom Phaeodactylum tricornutum: The Effects of Carbon Limitation and Photon Flux Density

Diatoms adapt to changing environmental conditions in very efficient ways. Among the mechanisms that can be activated, the reorientation of carbon metabolism is crucial because it allows the storage of energy into energy-dense molecules, typically lipids. Beside their roles in physiology, lipids are commercially interesting compounds. Therefore studies dealing with this topic are relevant for both basic and applied science. Although the molecular mechanisms involved in the reorientation of carbon metabolism as a response to a deficiency in nutrients such as nitrogen or phosphorus has been partially elucidated, the impacts of carbon availability on the implementation of the reorientation mechanisms remain unclear. Indeed, it has not been determined if the same types of mechanisms are activated under carbon and other nutrient deficiencies or limitations. The first aim of this work was to get insights into the physiological, biological and molecular processes triggered by progressive carbon starvation in the model diatom Phaeodactylum tricornutum. The second aim was to investigate the effects of the growth light intensity on these processes. For such a purpose three different photon flux densities 30, 300, and 1000 μmol photons m^-2 s^-1 were used. The results presented here demonstrate that under carbon limitation, diatom cells still reorient carbon metabolism toward either phosphoenolpyruvate or pyruvate, which serves as a hub for the production of more complex molecules. The distribution of carbon atoms between the different pathways was partially affected by the growth photon flux density because low light (LL) provides conditions for the accumulation of chrysolaminarin, while medium light mostly stimulated lipid synthesis. A significant increase in the amount of proteins was observed under high light (HL).

Integrated analysis of transcriptomic and proteomic data from tree peony (P. ostii) seeds reveals key developmental stages and candidate genes related to oil biosynthesis and fatty acid metabolism

Tree peony (Paeonia section Moutan DC.) seeds are an excellent source of beneficial natural compounds that promote health, and they contain high levels of alpha-linolenic acid (ALA). In recent years, tree peony has been emerging as an oil crop. Therefore, combined analysis of the transcriptome and proteome of tree peony (P. ostii) seeds at 25, 32, 39, 53, 67, 81, 88, 95, and 109 days after pollination (DAP) was conducted to better understand the transcriptional and translational regulation of seed development and oil biosynthesis. A total of 38,482 unigenes and 2841 proteins were identified. A total of 26,912 differentially expressed genes (DEGs) and 592 differentially expressed proteins (DEPs) were clustered into three groups corresponding to the rapid growth, seed inclusion enrichment and conversion, and late dehydration and mature stages of seed development. Fifteen lipid metabolism pathways were identified at both the transcriptome and proteome levels. Pathway enrichment analysis revealed that a period of rapid fatty acid biosynthesis occurred at 53-88 DAP. Furthermore, 211 genes and 35 proteins associated with the fatty acid metabolism pathway, 63 genes and 11 proteins associated with the biosynthesis of unsaturated fatty acids (UFAs), and 115 genes and 24 proteins associated with ALA metabolism were identified. Phylogenetic analysis revealed that 16 putative fatty acid desaturase (FAD)-encoding genes clustered into four FAD groups, eight of which exhibited the highest expression at 53 DAP, suggesting that they play an important role in ALA accumulation. RT-qPCR analysis indicated that the temporal expression patterns of oil biosynthesis genes were largely similar to the RNA-seq results. The expression patterns of fatty acid metabolism- and seed development-related proteins determined by MRM were also highly consistent with the results obtained in the proteomic analysis. Correlation analysis indicated significant differences in the number and abundance of DEGs and DEPs but a high level of consistency in expression patterns and metabolic pathways. The results of the present study represent the first combined transcriptomic and proteomic analysis of tree peony seeds and provide insight into tree peony seed development and oil accumulation.

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.

Assessing Flux Distribution Associated with Metabolic Specialization of Glandular Trichomes

Many aromatic plants accumulate mixtures of secondary (or specialized) metabolites in anatomical structures called glandular trichomes (GTs). Different GT types may also synthesize different mixtures of secreted metabolites, and this contributes to the enormous chemical diversity reported to occur across species. Over the past two decades, significant progress has been made in characterizing the genes and enzymes that are responsible for the unique metabolic capabilities of GTs in different lineages of flowering plants. Less is known about the processes that regulate flux distribution through precursor pathways toward metabolic end-products. We discuss here the results from a meta-analysis of genome-scale models that were developed to capture the unique metabolic capabilities of different GT types.

Divergent N Deficiency-Dependent Senescence and Transcriptome Response in Developmentally Old and Young Brassica napus Leaves

In the spring oilseed rape (OSR) cultivar 'Mozart' grown under optimal N supply (NO) or mild N deficiency (NL) the transcriptome changes associated with progressing age until early senescence in developmentally old lower canopy leaves (leaf #4) and younger higher canopy leaves (leaf #8) were investigated. Twelve weeks old NO and NL plants appeared phenotypically and transcriptomically identical, but thereafter distinct nutrition-dependent differences in gene expression patterns in lower and upper canopy leaves emerged. In NO leaves #4 of 14-week-old compared to 13-week-old plants, ∼600 genes were up- or downregulated, whereas in NL leaves #4 ∼3000 genes were up- or downregulated. In contrast, in 15-week-old compared to 13-week-old upper canopy leaves #8 more genes were up- or downregulated in optimally N-supplied plants (∼2000 genes) than in N-depleted plants (∼750 genes). This opposing effect of N depletion on gene regulation was even more prominent among photosynthesis-related genes (PSGs). Between week 13 and 14 in leaves #4, 99 of 110 PSGs were downregulated in NL plants, but none in NO plants. In contrast, from weeks 13 to 16 in leaves #8 of NL plants only 11 PSGs were downregulated in comparison to 66 PSGs in NO plants. Different effects of N depletion in lower versus upper canopy leaves were also apparent in upregulation of autophagy genes and NAC transcription factors. More than half of the regulated NAC and WRKY transcription factor, autophagy and protease genes were specifically regulated in NL leaves #4 or NO leaves #8 and thus may contribute to differences in senescence and nutrient mobilization in these leaves. We suggest that in N-deficient plants the upper leaves retain their N resources longer than in amply fertilized plants and remobilize them only after shedding of the lower leaves.

Comparative transcriptome and metabolome analysis suggests bottlenecks that limit seed and oil yields in transgenic Camelina sativa expressing diacylglycerol acyltransferase 1 and glycerol-3-phosphate dehydrogenase

BACKGROUND

Camelina sativa has attracted much interest as alternative renewable resources for biodiesel, other oil-based industrial products and a source for edible oils. Its unique oil attributes attract research to engineering new varieties of improved oil quantity and quality. The overexpression of enzymes catalyzing the synthesis of the glycerol backbone and the sequential conjugation of fatty acids into this backbone is a promising approach for increasing the levels of triacylglycerol (TAG). In a previous study, we co-expressed the diacylglycerol acyltransferase (DGAT1) and glycerol-3-phosphate dehydrogenase (GPD1), involved in TAG metabolism, in Camelina seeds. Transgenic plants exhibited a higher-percentage seed oil content, a greater seed mass, and overall improved seed and oil yields relative to wild-type plants. To further increase seed oil content in Camelina, we utilized metabolite profiling, in conjunction with transcriptome profiling during seed development to examine potential rate-limiting step(s) in the production of building blocks for TAG biosynthesis.

RESULTS

Transcriptomic analysis revealed approximately 2518 and 3136 transcripts differentially regulated at significant levels in DGAT1 and GPD1 transgenics, respectively. These transcripts were found to be involved in various functional categories, including alternative metabolic routes in fatty acid synthesis, TAG assembly, and TAG degradation. We quantified the relative contents of over 240 metabolites. Our results indicate major metabolic switches in transgenic seeds associated with significant changes in the levels of glycerolipids, amino acids, sugars, and organic acids, especially the TCA cycle and glycolysis intermediates.

CONCLUSIONS

From the transcriptomic and metabolomic analysis of DGAT1, GPD1 and DGAT1 + GPD1 expressing lines of C. sativa, we conclude that TAG production is limited by (1) utilization of fixed carbon from the source tissues supported by the increase in glycolysis pathway metabolites and decreased transcripts levels of transcription factors controlling fatty acids synthesis; (2) TAG accumulation is limited by the activity of lipases/hydrolases that hydrolyze TAG pool supported by the increase in free fatty acids and monoacylglycerols. This comparative transcriptomics and metabolomics approach is useful in understanding the regulation of TAG biosynthesis, identifying bottlenecks, and the corresponding genes controlling these pathways identified as limitations, for generating Camelina varieties with improved seed and oil yields.

Fluxomics links cellular functional analyses to whole-plant phenotyping

Fluxes through metabolic pathways reflect the integration of genetic and metabolic regulations. While it is attractive to measure all the mRNAs (transcriptome), all the proteins (proteome), and a large number of the metabolites (metabolome) in a given cellular system, linking and integrating this information remains difficult. Measurement of metabolome-wide fluxes (termed the fluxome) provides an integrated functional output of the cell machinery and a better tool to link functional analyses to plant phenotyping. This review presents and discusses sets of methodologies that have been developed to measure the fluxome. First, the principles of metabolic flux analysis (MFA), its 'short time interval' version Inst-MFA, and of constraints-based methods, such as flux balance analysis and kinetic analysis, are briefly described. The use of these powerful methods for flux characterization at the cellular scale up to the organ (fruits, seeds) and whole-plant level is illustrated. The added value given by fluxomics methods for unravelling how the abiotic environment affects flux, the process, and key metabolic steps are also described. Challenges associated with the development of fluxomics and its integration with 'omics' for thorough plant and organ functional phenotyping are discussed. Taken together, these will ultimately provide crucial clues for identifying appropriate target plant phenotypes for breeding.

Label-free deep shotgun proteomics reveals protein dynamics during tomato fruit tissues development

Current innovations in mass-spectrometry-based technologies allow deep coverage of protein expression. Despite its immense value and in contrast to transcriptomics, only a handful of studies in crop plants engaged with global proteome assays. Here, we present large-scale shotgun proteomics profiling of tomato fruit across two key tissues and five developmental stages. A total of 7738 individual protein groups were identified and reliably measured at least in one of the analyzed tissues or stages. The depth of our assay enabled identification of 61 differentially expressed transcription factors, including renowned ripening-related regulators and elements of ethylene signaling. Significantly, we measured proteins involved in 83% of all predicted enzymatic reactions in the tomato metabolic network. Hence, proteins representing almost the complete set of reactions in major metabolic pathways were identified, including the cytosolic and plastidic isoprenoid and the phenylpropanoid pathways. Furthermore, the data allowed us to discern between protein isoforms according to expression patterns, which is most significant in light of the weak transcript-protein expression correspondence. Finally, visualization of changes in protein abundance associated with a particular process provided us with a unique view of skin and flesh tissues in developing fruit. This study adds a new dimension to the existing genomic, transcriptomic and metabolomic resources. It is therefore likely to promote translational and post-translational research in tomato and additional species, which is presently focused on transcription.

Transcriptome Analysis of Stem and Globally Comparison with Other Tissues in Brassica napus

Brassica napus is one of the most important oilseed crops in the world. However, there is currently no enough stem transcriptome information and comparative transcriptome analysis of different tissues, which impedes further functional genomics research on B. napus. In this study, the stem transcriptome of B. napus was characterized by RNA-seq technology. Approximately 13.4 Gb high-quality clean reads with an average length of 100 bp were generated and used for comparative transcriptome analysis with the existing transcriptome sequencing data of roots, leaves, flower buds, and immature embryos of B. napus. All the transcripts were annotated against GO and KEGG databases. The common genes in five tissues, differentially expressed genes (DEGs) of the common genes between stems and other tissues, and tissue-specific genes were detected, and the main biochemical activities and pathways implying the common genes, DEGs and tissue-specific genes were investigated. Accordingly, the common transcription factors (TFs) in the five tissues and tissue-specific TFs were identified, and a TFs-based regulation network between TFs and the target genes involved in 'Phenylpropanoid biosynthesis' pathway were constructed to show several important TFs and key nodes in the regulation process. Collectively, this study not only provided an available stem transcriptome resource in B. napus, but also revealed valuable comparative transcriptome information of five tissues of B. napus for future investigation on specific processes, functions and pathways.

Identification of multiple lipid genes with modifications in expression and sequence associated with the evolution of hydroxy fatty acid accumulation in Physaria fendleri

Two Brassicaceae species, Physaria fendleri and Camelina sativa, are genetically very closely related to each other and to Arabidopsis thaliana. Physaria fendleri seeds contain over 50% hydroxy fatty acids (HFAs), while Camelina sativa and Arabidopsis do not accumulate HFAs. To better understand how plants evolved new biochemical pathways with the capacity to accumulate high levels of unusual fatty acids, transcript expression and protein sequences of developing seeds of Physaria fendleri, wild-type Camelina sativa, and Camelina sativa expressing a castor bean (Ricinus communis) hydroxylase were analyzed. A number of potential evolutionary adaptations within lipid metabolism that probably enhance HFA production and accumulation in Physaria fendleri, and, in their absence, limit accumulation in transgenic tissues were revealed. These adaptations occurred in at least 20 genes within several lipid pathways from the onset of fatty acid synthesis and its regulation to the assembly of triacylglycerols. Lipid genes of Physaria fendleri appear to have co-evolved through modulation of transcriptional abundances and alterations within protein sequences. Only a handful of genes showed evidence for sequence adaptation through gene duplication. Collectively, these evolutionary changes probably occurred to minimize deleterious effects of high HFA amounts and/or to enhance accumulation for physiological advantage. These results shed light on the evolution of pathways for novel fatty acid production in seeds, help explain some of the current limitations to accumulation of HFAs in transgenic plants, and may provide improved strategies for future engineering of their production.

Multiscale Metabolic Modeling of C4 Plants: Connecting Nonlinear Genome-Scale Models to Leaf-Scale Metabolism in Developing Maize Leaves

C4 plants, such as maize, concentrate carbon dioxide in a specialized compartment surrounding the veins of their leaves to improve the efficiency of carbon dioxide assimilation. Nonlinear relationships between carbon dioxide and oxygen levels and reaction rates are key to their physiology but cannot be handled with standard techniques of constraint-based metabolic modeling. We demonstrate that incorporating these relationships as constraints on reaction rates and solving the resulting nonlinear optimization problem yields realistic predictions of the response of C4 systems to environmental and biochemical perturbations. Using a new genome-scale reconstruction of maize metabolism, we build an 18000-reaction, nonlinearly constrained model describing mesophyll and bundle sheath cells in 15 segments of the developing maize leaf, interacting via metabolite exchange, and use RNA-seq and enzyme activity measurements to predict spatial variation in metabolic state by a novel method that optimizes correlation between fluxes and expression data. Though such correlations are known to be weak in general, we suggest that developmental gradients may be particularly suited to the inference of metabolic fluxes from expression data, and we demonstrate that our method predicts fluxes that achieve high correlation with the data, successfully capture the experimentally observed base-to-tip transition between carbon-importing tissue and carbon-exporting tissue, and include a nonzero growth rate, in contrast to prior results from similar methods in other systems.

Modeling Rice Metabolism: From Elucidating Environmental Effects on Cellular Phenotype to Guiding Crop Improvement

Crop productivity is severely limited by various biotic and abiotic stresses. Thus, it is highly needed to understand the underlying mechanisms of environmental stress response and tolerance in plants, which could be addressed by systems biology approach. To this end, high-throughput omics profiling and in silico modeling can be considered to explore the environmental effects on phenotypic states and metabolic behaviors of rice crops at the systems level. Especially, the advent of constraint-based metabolic reconstruction and analysis paves a way to characterize the plant cellular physiology under various stresses by combining the mathematical network models with multi-omics data. Rice metabolic networks have been reconstructed since 2013 and currently six such networks are available, where five are at genome-scale. Since their publication, these models have been utilized to systematically elucidate the rice abiotic stress responses and identify agronomic traits for crop improvement. In this review, we summarize the current status of the existing rice metabolic networks and models with their applications. Furthermore, we also highlight future directions of rice modeling studies, particularly stressing how these models can be used to contextualize the affluent multi-omics data that are readily available in the public domain. Overall, we envisage a number of studies in the future, exploiting the available metabolic models to enhance the yield and quality of rice and other food crops.

Genome wide analysis of protein production load in Trichoderma reesei

BACKGROUND

The filamentous fungus Trichoderma reesei (teleomorph Hypocrea jecorina) is a widely used industrial host organism for protein production. In industrial cultivations, it can produce over 100 g/l of extracellular protein, mostly constituting of cellulases and hemicellulases. In order to improve protein production of T. reesei the transcriptional regulation of cellulases and secretory pathway factors have been extensively studied. However, the metabolism of T. reesei under protein production conditions has not received much attention.

RESULTS

To understand the physiology and metabolism of T. reesei under protein production conditions we carried out a well-controlled bioreactor experiment with extensive analysis. We used minimal media to make the data amenable for modelling and three strain pairs to cover different protein production levels. With RNA-sequencing transcriptomics we detected the concentration of the carbon source as the most important determinant of the transcriptome. As the major transcriptional response concomitant to protein production we detected the induction of selected genes that were putatively regulated by xyr1 and were related to protein transport, amino acid metabolism and transcriptional regulation. We found novel metabolic responses such as production of glycerol and a cellotriose-like compound. We then used this cultivation data for flux balance analysis of T. reesei metabolism and demonstrate for the first time the use of genome wide stoichiometric metabolic modelling for T. reesei. We show that our model can predict protein production rate and provides novel insight into the metabolism of protein production. We also provide this unprecedented cultivation and transcriptomics data set for future modelling efforts.

CONCLUSIONS

The use of stoichiometric modelling can open a novel path for the improvement of protein production in T. reesei. Based on this we propose sulphur assimilation as a major limiting factor of protein production. As an organism with exceptional protein production capabilities modelling of T. reesei can provide novel insight also to other less productive organisms.

Time-Series Analyses of Transcriptomes and Proteomes Reveal Molecular Networks Underlying Oil Accumulation in Canola

Understanding the regulation of lipid metabolism is vital for genetic engineering of canola (Brassica napus L.) to increase oil yield or modify oil composition. We conducted time-series analyses of transcriptomes and proteomes to uncover the molecular networks associated with oil accumulation and dynamic changes in these networks in canola. The expression levels of genes and proteins were measured at 2, 4, 6, and 8 weeks after pollination (WAP). Our results show that the biosynthesis of fatty acids is a dominant cellular process from 2 to 6 WAP, while the degradation mainly happens after 6 WAP. We found that genes in almost every node of fatty acid synthesis pathway were significantly up-regulated during oil accumulation. Moreover, significant expression changes of two genes, acetyl-CoA carboxylase and acyl-ACP desaturase, were detected on both transcriptomic and proteomic levels. We confirmed the temporal expression patterns revealed by the transcriptomic analyses using quantitative real-time PCR experiments. The gene set association analysis show that the biosynthesis of fatty acids and unsaturated fatty acids are the most significant biological processes from 2-4 WAP and 4-6 WAP, respectively, which is consistent with the results of time-series analyses. These results not only provide insight into the mechanisms underlying lipid metabolism, but also reveal novel candidate genes that are worth further investigation for their values in the genetic engineering of canola.

Inference and Prediction of Metabolic Network Fluxes

In this Update, we cover the basic principles of the estimation and prediction of the rates of the many interconnected biochemical reactions that constitute plant metabolic networks. This includes metabolic flux analysis approaches that utilize the rates or patterns of redistribution of stable isotopes of carbon and other atoms to estimate fluxes, as well as constraints-based optimization approaches such as flux balance analysis. Some of the major insights that have been gained from analysis of fluxes in plants are discussed, including the functioning of metabolic pathways in a network context, the robustness of the metabolic phenotype, the importance of cell maintenance costs, and the mechanisms that enable energy and redox balancing at steady state. We also discuss methodologies to exploit 'omic data sets for the construction of tissue-specific metabolic network models and to constrain the range of permissible fluxes in such models. Finally, we consider the future directions and challenges faced by the field of metabolic network flux phenotyping.

Metabolic Architecture of the Cereal Grain and Its Relevance to Maximize Carbon Use Efficiency

Here, we have characterized the spatial heterogeneity of the cereal grain's metabolism and demonstrated how, by integrating a distinct set of metabolic strategies, the grain has evolved to become an almost perfect entity for carbon storage. In vivo imaging revealed light-induced cycles in assimilate supply toward the ear/grain of barley (Hordeum vulgare) and wheat (Triticum aestivum). In silico modeling predicted that, in the two grain storage organs (the endosperm and embryo), the light-induced shift in solute influx does cause adjustment in metabolic flux without changes in pathway utilization patterns. The enveloping, leaf-like pericarp, in contrast, shows major shifts in flux distribution (starch metabolism, photosynthesis, remobilization, and tricarboxylic acid cycle activity) allow to refix 79% of the CO2 released by the endosperm and embryo, allowing the grain to achieve an extraordinary high carbon conversion efficiency of 95%. Shading experiments demonstrated that ears are autonomously able to raise the influx of solutes in response to light, but with little effect on the steady-state levels of metabolites or transcripts or on the pattern of sugar distribution within the grain. The finding suggests the presence of a mechanism(s) able to ensure metabolic homeostasis in the face of short-term environmental fluctuation. The proposed multicomponent modeling approach is informative for predicting the metabolic effects of either an altered level of incident light or a momentary change in the supply of sucrose. It is therefore of potential value for assessing the impact of either breeding and/or biotechnological interventions aimed at increasing grain yield.

Metabolic Plasticity and Inter-Compartmental Interactions in Rice Metabolism: An Analysis from Reaction Deletion Study

More than 20% of the total caloric intake of human population comes from rice. The expression of rice genes and hence, the concentration of enzymatic proteins might vary due to several biotic and abiotic stresses. It in turn, can influence the overall metabolism and survivability of rice plant. Thus, understanding the rice cellular metabolism, its plasticity and potential readjustments under different perturbations can help rice biotechnologists to design efficient rice cultivars. Here, using the flux balance analysis (FBA) method, with the help of in-silico reaction deletion strategy, we study the metabolic plasticity of genome-scale metabolic model of rice leaf. A set of 131 reactions, essential for the production of primary biomass precursors is identified; deletion of any of them can inhibit the overall biomass production. Usability Index (I_U) for the rest of the reactions are estimated and based on this parameter, they are classified into three categories—maximally-favourable, quasi-favourable and unfavourable for the primary biomass production. The lower value of 1 − I_U of a reaction suggests that the cell cannot easily bypass it for biomass production. While some of the alternative paths are energetically equally efficient, others demand for higher photon. The variations in (i) ATP/NADPH ratio, (ii) exchange of metabolites through chloroplastic transporters and (iii) total biomass production are also presented here. Mutual metabolic dependencies of different cellular compartments are also demonstrated.

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.

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.