Triacylglycerol stability limits futile cycles and inhibition of carbon capture in oil-accumulating leaves.

Engineering plant vegetative tissue to accumulate triacylglycerols (TAG, e.g., oil) can increase the amount of oil harvested per acre to levels that exceed current oilseed crops. Engineered tobacco (Nicotiana tabacum) lines that accumulate 15% to 30% oil of leaf dry weight resulted in starkly different metabolic phenotypes. In-depth analysis of the leaf lipid accumulation and 14CO2 tracking describe metabolic adaptations to the leaf oil engineering. An oil-for-membrane lipid tradeoff in the 15% oil line (referred to as HO) was surprisingly not further exacerbated when lipid production was enhanced to 30% (LEC2 line). The HO line exhibited a futile cycle that limited TAG yield through exchange with starch, altered carbon flux into various metabolite pools and end products, and suggested interference of the glyoxylate cycle with photorespiration that limited CO2 assimilation by 50%. In contrast, inclusion of the LEAFY COTYLEDON 2 (LEC2) transcription factor in tobacco improved TAG stability, alleviated the TAG-to-starch futile cycle, and recovered CO2 assimilation and plant growth comparable to wild type but with much higher lipid levels in leaves. Thus, the unstable production of storage reserves and futile cycling limit vegetative oil engineering approaches. The capacity to overcome futile cycles and maintain enhanced stable TAG levels in LEC2 demonstrated the importance of considering unanticipated metabolic adaptations while engineering vegetative oil crops.

Program for Integration and Rapid Analysis of Mass Isotopomer Distributions (PIRAMID).

SUMMARY: The analysis of stable isotope labeling experiments requires accurate, efficient, and reproducible quantification of mass isotopomer distributions (MIDs), which is not a core feature of general-purpose metabolomics software tools that are optimized to quantify metabolite abundance. Here, we present PIRAMID (Program for Integration and Rapid Analysis of Mass Isotopomer Distributions), a MATLAB-based tool that addresses this need by offering a user-friendly, graphical user interface-driven program to automate the extraction of isotopic information from mass spectrometry (MS) datasets. This tool can simultaneously extract ion chromatograms for various metabolites from multiple data files in common vendor-agnostic file formats, locate chromatographic peaks based on a targeted list of characteristic ions and retention times, and integrate MIDs for each target ion. These MIDs can be corrected for natural isotopic background based on the user-defined molecular formula of each ion. PIRAMID offers support for datasets acquired from low- or high-resolution MS, and single (MS) or tandem (MS/MS) instruments. It also enables the analysis of single or dual labeling experiments using a variety of isotopes (i.e. 2H, 13C, 15N, 18O, 34S). DATA AVAILABILITY AND IMPLEMENTATION: MATLAB p-code files are freely available for non-commercial use and can be downloaded from https://mfa.vueinnovations.com/. Commercial licenses are also available. All the data presented in this publication are available under the "Help_menu" folder of the PIRAMID software.

Physcomitrium patens response to elevated CO2 is flexible and determined by an interaction between sugar and nitrogen availability.

Mosses hold a unique position in plant evolution and are crucial for protecting natural, long-term carbon storage systems such as permafrost and bogs. Due to small stature, mosses grow close to the soil surface and are exposed to high levels of CO2 , produced by soil respiration. However, the impact of elevated CO2 (eCO2 ) levels on mosses remains underexplored. We determined the growth responses of the moss Physcomitrium patens to eCO2 in combination with different nitrogen levels and characterized the underlying physiological and metabolic changes. Three distinct growth characteristics, an early transition to caulonema, the development of longer, highly pigmented rhizoids, and increased biomass, define the phenotypic responses of P. patens to eCO2 . Elevated CO2 impacts growth by enhancing the level of a sugar signaling metabolite, T6P. The quantity and form of nitrogen source influences these metabolic and phenotypic changes. Under eCO2 , P. patens exhibits a diffused growth pattern in the presence of nitrate, but ammonium supplementation results in dense growth with tall gametophores, demonstrating high phenotypic plasticity under different environments. These results provide a framework for comparing the eCO2 responses of P. patens with other plant groups and provide crucial insights into moss growth that may benefit climate change models.

Isotopically nonstationary metabolic flux analysis of plants: recent progress and future opportunities.

Metabolic flux analysis (MFA) is a valuable tool for quantifying cellular phenotypes and to guide plant metabolic engineering. By introducing stable isotopic tracers and employing mathematical models, MFA can quantify the rates of metabolic reactions through biochemical pathways. Recent applications of isotopically nonstationary MFA (INST-MFA) to plants have elucidated nonintuitive metabolism in leaves under optimal and stress conditions, described coupled fluxes for fast-growing algae, and produced a synergistic multi-organ flux map that is a first in MFA for any biological system. These insights could not be elucidated through other approaches and show the potential of INST-MFA to correct an oversimplified understanding of plant metabolism.

Arabidopsis ACYL CARRIER PROTEIN4 and RHOMBOID LIKE10 act independently in chloroplast phosphatidate synthesis.

ACYL CARRIER PROTEIN4 (ACP4) is the most abundant ACP isoform in Arabidopsis (Arabidopsis thaliana) leaves and acts as a scaffold for de novo fatty acid biosynthesis and as a substrate for acyl-ACP-utilizing enzymes. Recently, ACP4 was found to interact with a protein-designated plastid RHOMBOID LIKE10 (RBL10) that affects chloroplast monogalactosyldiacylglycerol (MGDG) biosynthesis, but the cellular function of this interaction remains to be explored. Here, we generated and characterized acp4 rbl10 double mutants to explore whether ACP4 and RBL10 directly interact in influencing chloroplast lipid metabolism. Alterations in the content and molecular species of chloroplast lipids such as MGDG and phosphatidylglycerol were observed in the acp4 and rbl10 mutants, which are likely associated with the changes in the size and profiles of diacylglycerol (DAG), phosphatidic acid (PA), and acyl-ACP precursor pools. ACP4 contributed to the size and profile of the acyl-ACP pool and interacted with acyl-ACP-utilizing enzymes, as expected for its role in fatty acid biosynthesis and chloroplast lipid assembly. RBL10 appeared to be involved in the conversion of PA to DAG precursors for MGDG biosynthesis as evidenced by the increased 34:x PA and decreased 34:x DAG in the rbl10 mutant and the slow turnover of radiolabeled PA in isolated chloroplasts fed with [14C] acetate. Interestingly, the impaired PA turnover in rbl10 was partially reversed in the acp4 rbl10 double mutant. Collectively, this study shows that ACP4 and RBL10 affect chloroplast lipid biosynthesis by modulating substrate precursor pools and appear to act independently.

Expression of malic enzyme reveals subcellular carbon partitioning for storage reserve production in soybeans.

Central metabolism produces amino and fatty acids for protein and lipids that establish seed value. Biosynthesis of storage reserves occurs in multiple organelles that exchange central intermediates including two essential metabolites, malate, and pyruvate that are linked by malic enzyme. Malic enzyme can be active in multiple subcellular compartments, partitioning carbon and reducing equivalents for anabolic and catabolic requirements. Prior studies based on isotopic labeling and steady-state metabolic flux analyses indicated malic enzyme provides carbon for fatty acid biosynthesis in plants, though genetic evidence confirming this role is lacking. We hypothesized that increasing malic enzyme flux would alter carbon partitioning and result in increased lipid levels in soybeans. Homozygous transgenic soybean plants expressing Arabidopsis malic enzyme alleles, targeting the translational products to plastid or outside the plastid during seed development, were verified by transcript and enzyme activity analyses, organelle proteomics, and transient expression assays. Protein, oil, central metabolites, cofactors, and acyl-acyl carrier protein (ACPs) levels were quantified overdevelopment. Amino and fatty acid levels were altered resulting in an increase in lipids by 0.5-2% of seed biomass (i.e. 2-9% change in oil). Subcellular targeting of a single gene product in central metabolism impacts carbon and reducing equivalent partitioning for seed storage reserves in soybeans.

A General Method for Quantification and Discovery of Acyl Groups Attached to Acyl Carrier Proteins in Fatty Acid Metabolism Using LC-MS/MS.

Acyl carrier proteins (ACPs) are the scaffolds for fatty acid biosynthesis in living systems, rendering them essential to a comprehensive understanding of lipid metabolism. However, accurate quantitative methods to assess individual acyl-ACPs do not exist. We developed a robust method to quantify acyl-ACPs to the picogram level. We successfully identified acyl-ACP elongation intermediates (3-hydroxyacyl-ACPs and 2,3-trans-enoyl-ACPs) and unexpected medium-chain (C10:1, C14:1) and polyunsaturated long-chain (C16:3) acyl-ACPs, indicating both the sensitivity of the method and how current descriptions of lipid metabolism and ACP function are incomplete. Such ACPs are likely important to medium-chain lipid production for fuels and highlight poorly understood lipid remodeling events in the chloroplast. The approach is broadly applicable to type II fatty acid synthase systems found in plants and bacteria as well as mitochondria from mammals and fungi because it capitalizes on a highly conserved Asp-Ser-Leu-Asp amino acid sequence in ACPs to which acyl groups attach. Our method allows for sensitive quantification using liquid chromatography-tandem mass spectrometry with de novo-generated standards and an isotopic dilution strategy and will fill a gap in our understanding, providing insights through quantitative exploration of fatty acid biosynthesis processes for optimal biofuels, renewable feedstocks, and medical studies in health and disease.

Sterile Spikelets Contribute to Yield in Sorghum and Related Grasses.

Sorghum (Sorghum bicolor) and its relatives in the grass tribe Andropogoneae bear their flowers in pairs of spikelets in which one spikelet (seed-bearing or sessile spikelet [SS]) of the pair produces a seed and the other is sterile or male (staminate). This division of function does not occur in other major cereals such as wheat (Triticum aestivum) or rice (Oryza sativa). Additionally, one bract of the SS spikelet often produces a long extension, the awn, that is in the same position as, but independently derived from, that of wheat and rice. The function of the sterile spikelet is unknown and that of the awn has not been tested in Andropogoneae. We used radioactive and stable isotopes of carbon, RNA sequencing of metabolically important enzymes, and immunolocalization of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) to show that the sterile spikelet assimilates carbon, which is translocated to the largely heterotrophic SS. The awn shows no evidence of photosynthesis. These results apply to distantly related species of Andropogoneae. Removal of sterile spikelets in sorghum significantly decreases seed weight (yield) by ∼9%. Thus, the sterile spikelet, but not the awn, affects yield in the cultivated species and fitness in the wild species.

Metabolic flux analysis of the non-transitory starch tradeoff for lipid production in mature tobacco leaves.

The metabolic plasticity of tobacco leaves has been demonstrated via the generation of transgenic plants that can accumulate over 30% dry weight as triacylglycerols. In investigating the changes in carbon partitioning in these high lipid-producing (HLP) leaves, foliar lipids accumulated stepwise over development. Interestingly, non-transient starch was observed to accumulate with plant age in WT but not HLP leaves, with a drop in foliar starch concurrent with an increase in lipid content. The metabolic carbon tradeoff between starch and lipid was studied using 13CO2-labeling experiments and isotopically nonstationary metabolic flux analysis, not previously applied to the mature leaves of a crop. Fatty acid synthesis was investigated through assessment of acyl-acyl carrier proteins using a recently derived quantification method that was extended to accommodate isotopic labeling. Analysis of labeling patterns and flux modeling indicated the continued production of unlabeled starch, sucrose cycling, and a significant contribution of NADP-malic enzyme to plastidic pyruvate production for the production of lipids in HLP leaves, with the latter verified by enzyme activity assays. The results suggest an inherent capacity for a developmentally regulated carbon sink in tobacco leaves and may in part explain the uniquely successful leaf lipid engineering efforts in this crop.

Temporal changes in metabolism late in seed development affect biomass composition.

The negative association between protein and oil production in soybean (Glycine max) seed is well-documented. However, this inverse relationship is based primarily on the composition of mature seed, which reflects the cumulative result of events over the course of soybean seed development and therefore does not convey information specific to metabolic fluctuations during developmental growth regimes. In this study, we assessed maternal nutrient supply via measurement of seed coat exudates and metabolite levels within the cotyledon throughout development to identify trends in the accumulation of central carbon and nitrogen metabolic intermediates. Active metabolic activity during late seed development was probed through transient labeling with 13C substrates. The results indicated: (1) a drop in lipid contents during seed maturation with a concomitant increase in carbohydrates, (2) a transition from seed filling to maturation phases characterized by quantitatively balanced changes in carbon use and CO2 release, (3) changes in measured carbon and nitrogen resources supplied maternally throughout development, (4) 13C metabolite production through gluconeogenic steps for sustained carbohydrate accumulation as the maternal nutrient supply diminishes, and (5) oligosaccharide biosynthesis within the seed coat during the maturation phase. These results highlight temporal engineering targets for altering final biomass composition to increase the value of soybeans and a path to breaking the inverse correlation between seed protein and oil content.

An efficient LC-MS method for isomer separation and detection of sugars, phosphorylated sugars, and organic acids.

Assessing central carbon metabolism in plants can be challenging due to the dynamic range in pool sizes, with low levels of important phosphorylated sugars relative to more abundant sugars and organic acids. Here, we report a sensitive liquid chromatography-mass spectrometry method for analysing central metabolites on a hybrid column, where both anion-exchange and hydrophilic interaction chromatography (HILIC) ligands are embedded in the stationary phase. The liquid chromatography method was developed for enhanced selectivity of 27 central metabolites in a single run with sensitivity at femtomole levels observed for most phosphorylated sugars. The method resolved phosphorylated hexose, pentose, and triose isomers that are otherwise challenging. Compared with a standard HILIC approach, these metabolites had improved peak areas using our approach due to ion enhancement or low ion suppression in the biological sample matrix. The approach was applied to investigate metabolism in high lipid-producing tobacco leaves that exhibited increased levels of acetyl-CoA, a precursor for oil biosynthesis. The application of the method to isotopologue detection and quantification was considered through evaluating 13C-labeled seeds from Camelina sativa. The method provides a means to analyse intermediates more comprehensively in central metabolism of plant tissues.

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.

Interactions of C4 Subtype Metabolic Activities and Transport in Maize Are Revealed through the Characterization of DCT2 Mutants.

C4 photosynthesis in grasses requires the coordinated movement of metabolites through two specialized leaf cell types, mesophyll (M) and bundle sheath (BS), to concentrate CO2 around Rubisco. Despite the importance of transporters in this process, few have been identified or rigorously characterized. In maize (Zea mays), DCT2 has been proposed to function as a plastid-localized malate transporter and is preferentially expressed in BS cells. Here, we characterized the role of DCT2 in maize leaves using Activator-tagged mutant alleles. Our results indicate that DCT2 enables the transport of malate into the BS chloroplast. Isotopic labeling experiments show that the loss of DCT2 results in markedly different metabolic network operation and dramatically reduced biomass production. In the absence of a functioning malate shuttle, dct2 lines survive through the enhanced use of the phosphoenolpyruvate carboxykinase carbon shuttle pathway that in wild-type maize accounts for ∼ 25% of the photosynthetic activity. The results emphasize the importance of malate transport during C4 photosynthesis, define the role of a primary malate transporter in BS cells, and support a model for carbon exchange between BS and M cells in maize.

Synergism between Inositol Polyphosphates and TOR Kinase Signaling in Nutrient Sensing, Growth Control, and Lipid Metabolism in Chlamydomonas.

The networks that govern carbon metabolism and control intracellular carbon partitioning in photosynthetic cells are poorly understood. Target of Rapamycin (TOR) kinase is a conserved growth regulator that integrates nutrient signals and modulates cell growth in eukaryotes, though the TOR signaling pathway in plants and algae has yet to be completely elucidated. We screened the unicellular green alga Chlamydomonas reinhardtii using insertional mutagenesis to find mutants that conferred hypersensitivity to the TOR inhibitor rapamycin. We characterized one mutant, vip1-1, that is predicted to encode a conserved inositol hexakisphosphate kinase from the VIP family that pyrophosphorylates phytic acid (InsP6) to produce the low abundance signaling molecules InsP7 and InsP8 Unexpectedly, the rapamycin hypersensitive growth arrest of vip1-1 cells was dependent on the presence of external acetate, which normally has a growth-stimulatory effect on Chlamydomonas. vip1-1 mutants also constitutively overaccumulated triacylglycerols (TAGs) in a manner that was synergistic with other TAG inducing stimuli such as starvation. vip1-1 cells had reduced InsP7 and InsP8, both of which are dynamically modulated in wild-type cells by TOR kinase activity and the presence of acetate. Our data uncover an interaction between the TOR kinase and inositol polyphosphate signaling systems that we propose governs carbon metabolism and intracellular pathways that lead to storage lipid accumulation.

Phospholipase Dζ Enhances Diacylglycerol Flux into Triacylglycerol.

Plant seeds are the primary source of triacylglycerols (TAG) for food, feed, fuel, and industrial applications. As TAG is produced from diacylglycerol (DAG), successful engineering strategies to enhance TAG levels have focused on the conversion of DAG to TAG. However, the production of TAG can be limited by flux through the enzymatic reactions that supply DAG. In this study, two Arabidopsis phospholipase Dζ genes (AtPLDζ1 and AtPLDζ2 ) were coexpressed in Camelina sativa to test whether the conversion of phosphatidylcholine to DAG impacts TAG levels in seeds. The resulting transgenic plants produced 2% to 3% more TAG as a component of total seed biomass and had increased 18:3 and 20:1 fatty acid levels relative to wild type. Increased DAG and decreased PC levels were examined through the kinetics of lipid assembly by [14C]acetate and [14C]glycerol incorporation into glycerolipids. [14C]acetate was rapidly incorporated into TAG in both wild-type and overexpression lines, indicating a significant flux of nascent and elongated acyl-CoAs into the sn-3 position of TAG. Stereochemical analysis revealed that newly synthesized fatty acids were preferentially incorporated into the sn-2 position of PC, but the sn-1 position of de novo DAG and indicated similar rates of nascent acyl groups into the Kennedy pathway and acyl editing. [14C]glycerol studies demonstrated PC-derived DAG is the major source of DAG for TAG synthesis in both tissues. The results emphasize that the interconversions of DAG and PC pools can impact oil production and composition.

Assessing compartmentalized flux in lipid metabolism with isotopes.

Metabolism in plants takes place across multiple cell types and within distinct organelles. The distributions equate to spatial heterogeneity; though the limited means to experimentally assess metabolism frequently involve homogenizing tissues and mixing metabolites from different locations. Most current isotope investigations of metabolism therefore lack the ability to resolve spatially distinct events. Recognition of this limitation has resulted in inspired efforts to advance metabolic flux analysis and isotopic labeling techniques. Though a number of these efforts have been applied to studies in central metabolism; recent advances in instrumentation and techniques present an untapped opportunity to make similar progress in lipid metabolism where the use of stable isotopes has been more limited. These efforts will benefit from sophisticated radiolabeling reports that continue to enrich our knowledge on lipid biosynthetic pathways and provide some direction for stable isotope experimental design and extension of MFA. Evidence for this assertion is presented through the review of several elegant stable isotope studies and by taking stock of what has been learned from radioisotope investigations when spatial aspects of metabolism were considered. The studies emphasize that glycerolipid production occurs across several locations with assembly of lipids in the ER or plastid, fatty acid biosynthesis occurring in the plastid, and the generation of acetyl-CoA and glycerol-3-phosphate taking place at multiple sites. Considering metabolism in this context underscores the cellular and subcellular organization that is important to enhanced production of glycerolipids in plants. An attempt is made to unify salient features from a number of reports into a diagrammatic model of lipid metabolism and propose where stable isotope labeling experiments and further flux analysis may help address questions in the field. This article is part of a Special Issue entitled: Plant Lipid Biology edited by Kent D. Chapman and Ivo Feussner.

Cytosolic phosphorylating glyceraldehyde-3-phosphate dehydrogenases affect Arabidopsis cellular metabolism and promote seed oil accumulation.

The cytosolic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPC) catalyzes a key reaction in glycolysis, but its contribution to plant metabolism and growth are not well defined. Here, we show that two cytosolic GAPCs play important roles in cellular metabolism and seed oil accumulation. Knockout or overexpression of GAPCs caused significant changes in the level of intermediates in the glycolytic pathway and the ratios of ATP/ADP and NAD(P)H/NAD(P). Two double knockout seeds had ∼3% of dry weight decrease in oil content compared with that of the wild type. In transgenic seeds under the constitutive 35S promoter, oil content was increased up to 42% of dry weight compared with 36% in the wild type and the fatty acid composition was altered; however, these transgenic lines exhibited decreased fertility. Seed-specific overexpression lines had >3% increase in seed oil without compromised seed yield or fecundity. The results demonstrate that GAPC levels play important roles in the overall cellular production of reductants, energy, and carbohydrate metabolites and that GAPC levels are directly correlated with seed oil accumulation. Changes in cellular metabolites and cofactor levels highlight the complexity and tolerance of Arabidopsis thaliana cells to the metabolic perturbation. Further implications for metabolic engineering of seed oil production are discussed.

Rerouting of carbon flux in a glycogen mutant of cyanobacteria assessed via isotopically non-stationary 13 C metabolic flux analysis.

Cyanobacteria, which constitute a quantitatively dominant phylum, have attracted attention in biofuel applications due to favorable physiological characteristics, high photosynthetic efficiency and amenability to genetic manipulations. However, quantitative aspects of cyanobacterial metabolism have received limited attention. In the present study, we have performed isotopically non-stationary 13 C metabolic flux analysis (INST-13 C-MFA) to analyze rerouting of carbon in a glycogen synthase deficient mutant strain (glgA-I glgA-II) of the model cyanobacterium Synechococcus sp. PCC 7002. During balanced photoautotrophic growth, 10-20% of the fixed carbon is stored in the form of glycogen via a pathway that is conserved across the cyanobacterial phylum. Our results show that deletion of glycogen synthase gene orchestrates cascading effects on carbon distribution in various parts of the metabolic network. Carbon that was originally destined to be incorporated into glycogen gets partially diverted toward alternate storage molecules such as glucosylglycerol and sucrose. The rest is partitioned within the metabolic network, primarily via glycolysis and tricarboxylic acid cycle. A lowered flux toward carbohydrate synthesis and an altered distribution at the glucose-1-phosphate node indicate flexibility in the network. Further, reversibility of glycogen biosynthesis reactions points toward the presence of futile cycles. Similar redistribution of carbon was also predicted by Flux Balance Analysis. The results are significant to metabolic engineering efforts with cyanobacteria where fixed carbon needs to be re-routed to products of interest. Biotechnol. Bioeng. 2017;114: 2298-2308. © 2017 Wiley Periodicals, Inc.

Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation.

Improving plant productivity is an important aim for metabolic engineering. There are few comprehensive methods that quantitatively describe leaf metabolism, although such information would be valuable for increasing photosynthetic capacity, enhancing biomass production, and rerouting carbon flux toward desirable end products. Isotopically nonstationary metabolic flux analysis (INST-MFA) has been previously applied to map carbon fluxes in photoautotrophic bacteria, which involves model-based regression of transient (13)C-labeling patterns of intracellular metabolites. However, experimental and computational difficulties have hindered its application to terrestrial plant systems. We performed in vivo isotopic labeling of Arabidopsis thaliana rosettes with (13)CO2 and estimated fluxes throughout leaf photosynthetic metabolism by INST-MFA. Plants grown at 200 µmol m(-2)s(-1) light were compared with plants acclimated for 9 d at an irradiance of 500 µmol⋅m(-2)⋅s(-1). Approximately 1,400 independent mass isotopomer measurements obtained from analysis of 37 metabolite fragment ions were regressed to estimate 136 total fluxes (54 free fluxes) under each condition. The results provide a comprehensive description of changes in carbon partitioning and overall photosynthetic flux after long-term developmental acclimation of leaves to high light. Despite a doubling in the carboxylation rate, the photorespiratory flux increased from 17 to 28% of net CO2 assimilation with high-light acclimation (Vc/Vo: 3.5:1 vs. 2.3:1, respectively). This study highlights the potential of (13)C INST-MFA to describe emergent flux phenotypes that respond to environmental conditions or plant physiology and cannot be obtained by other complementary approaches.