Use of non-aqueous fractionation and metabolomics to study chloroplast function in Arabidopsis

"Omics" insights into plastid behavior toward improved carotenoid accumulation

Plastids are a group of diverse organelles with conserved carotenoids synthesizing and sequestering functions in plants. They optimize the carotenoid composition and content in response to developmental transitions and environmental stimuli. In this review, we describe the turbulence and reforming of transcripts, proteins, and metabolic pathways for carotenoid metabolism and storage in various plastid types upon organogenesis and external influences, which have been studied using approaches including genomics, transcriptomics, proteomics, and metabonomics. Meanwhile, the coordination of plastid signaling and carotenoid metabolism including the effects of disturbed carotenoid biosynthesis on plastid morphology and function are also discussed. The "omics" insight extends our understanding of the interaction between plastids and carotenoids and provides significant implications for designing strategies for carotenoid-biofortified crops.

Non-aqueous fractionation revealed changing subcellular metabolite distribution during apple fruit development

In developing apple fruit, metabolic compartmentation is poorly understood due to the lack of experimental data. Distinguishing subcellular compartments in fruit using non-aqueous fractionation has been technically difficult due to the excess amount of sugars present in the different subcellular compartments limiting the resolution of the technique. The work described in this study represents the first attempt to apply non-aqueous fractionation to developing apple fruit, covering the major events occurring during fruit development (cell division, cell expansion, and maturation). Here we describe the non-aqueous fractionation method to study the subcellular compartmentation of metabolites during apple fruit development considering three main cellular compartments (cytosol, plastids, and vacuole). Evidence is presented that most of the sugars and organic acids were predominantly located in the vacuole, whereas some of the amino acids were distributed between the cytosol and the vacuole. The results showed a shift in the plastid marker from the lightest fractions in the early growth stage to the dense fractions in the later fruit growth stages. This implies that the accumulation of starch content with progressing fruit development substantially influenced the distribution of plastidial fragments within the non-aqueous density gradient applied. Results from this study provide substantial baseline information on assessing the subcellular compartmentation of metabolites in apple fruit in general and during fruit growth in particular.

Overexpression of Hevea brasiliensis ethylene response factor HbERF-IXc5 enhances growth and tolerance to abiotic stress and affects laticifer differentiation

Ethylene response factor 1 (ERF1) is an essential integrator of the jasmonate and ethylene signalling pathways coordinating a large number of genes involved in plant defences. Its orthologue in Hevea brasiliensis, HbERF-IXc5, has been assumed to play a major role in laticifer metabolism and tolerance to harvesting stress for better latex production. This study sets out to establish and characterize rubber transgenic lines overexpressing HbERF-IXc5. Overexpression of HbERF-IXc5 dramatically enhanced plant growth and enabled plants to maintain some ecophysiological parameters in response to abiotic stress such as water deficit, cold and salt treatments. This study revealed that HbERF-IXc5 has rubber-specific functions compared to Arabidopsis ERF1 as transgenic plants overexpressing HbERF-IXc5 accumulated more starch and differentiated more latex cells at the histological level. The role of HbERF-IXc5 in driving the expression of some target genes involved in laticifer differentiation is discussed.

A network property necessary for concentration robustness

Maintenance of functionality of complex cellular networks and entire organisms exposed to environmental perturbations often depends on concentration robustness of the underlying components. Yet, the reasons and consequences of concentration robustness in large-scale cellular networks remain largely unknown. Here, we derive a necessary condition for concentration robustness based only on the structure of networks endowed with mass action kinetics. The structural condition can be used to design targeted experiments to study concentration robustness. We show that metabolites satisfying the necessary condition are present in metabolic networks from diverse species, suggesting prevalence of this property across kingdoms of life. We also demonstrate that our predictions about concentration robustness of energy-related metabolites are in line with experimental evidence from Escherichia coli. The necessary condition is applicable to mass action biological systems of arbitrary size, and will enable understanding the implications of concentration robustness in genetic engineering strategies and medical applications.

Absolute concentration robustness (ACR), independence of the steady-state concentration of a molecule from the environment, is difficult to predict. Here, the authors derive a network structure-based necessary condition for ACR, and suggest that metabolites satisfying the condition are prevalent.

A Benchtop Fractionation Procedure for Subcellular Analysis of the Plant Metabolome

Although compartmentation is a key feature of eukaryotic cells, biological research is frequently limited by methods allowing for the comprehensive subcellular resolution of the metabolome. It has been widely accepted that such a resolution would be necessary in order to approximate cellular biochemistry and metabolic regulation, yet technical challenges still limit both the reproducible subcellular fractionation and the sample throughput being necessary for a statistically robust analysis. Here, we present a method and a detailed protocol which is based on the non-aqueous fractionation technique enabling the assignment of metabolites to their subcellular localization. The presented benchtop method aims at unraveling subcellular metabolome dynamics in a precise and statistically robust manner using a relatively small amount of tissue material. The method is based on the separation of cellular fractions via density gradients consisting of organic, non-aqueous solvents. By determining the relative distribution of compartment-specific marker enzymes together with metabolite profiles over the density gradient it is possible to estimate compartment-specific metabolite concentrations by correlation. To support this correlation analysis, a spreadsheet is provided executing a calculation algorithm to determine the distribution of metabolites over subcellular compartments. The calculation algorithm performs correlation of marker enzyme activity and metabolite abundance accounting for technical errors, reproducibility and the resulting error propagation. The method was developed, tested and validated in three natural accessions of Arabidopsis thaliana showing different ability to acclimate to low temperature. Particularly, amino acids were strongly shuffled between subcellular compartments in a cold-sensitive accession while a cold-tolerant accession was characterized by a stable subcellular metabolic homeostasis. Finally, we conclude that subcellular metabolome analysis is essential to unambiguously unravel regulatory strategies being involved in plant-environment interactions.

Disruption of both chloroplastic and cytosolic FBPase genes results in a dwarf phenotype and important starch and metabolite changes in Arabidopsis thaliana

In this study, evidence is provided for the role of fructose-1,6-bisphosphatases (FBPases) in plant development and carbohydrate synthesis and distribution by analysing two Arabidopsis thaliana T-DNA knockout mutant lines, cyfbp and cfbp1, and one double mutant cyfbp cfbp1 which affect each FBPase isoform, cytosolic and chloroplastic, respectively. cyFBP is involved in sucrose synthesis, whilst cFBP1 is a key enzyme in the Calvin-Benson cycle. In addition to the smaller rosette size and lower rate of photosynthesis, the lack of cFBP1 in the mutants cfbp1 and cyfbp cfbp1 leads to a lower content of soluble sugars, less starch accumulation, and a greater superoxide dismutase (SOD) activity. The mutants also had some developmental alterations, including stomatal opening defects and increased numbers of root vascular layers. Complementation also confirmed that the mutant phenotypes were caused by disruption of the cFBP1 gene. cyfbp mutant plants without cyFBP showed a higher starch content in the chloroplasts, but this did not greatly affect the phenotype. Notably, the sucrose content in cyfbp was close to that found in the wild type. The cyfbp cfbp1 double mutant displayed features of both parental lines but had the cfbp1 phenotype. All the mutants accumulated fructose-1,6-bisphosphate and triose-phosphate during the light period. These results prove that while the lack of cFBP1 induces important changes in a wide range of metabolites such as amino acids, sugars, and organic acids, the lack of cyFBP activity in Arabidopsis essentially provokes a carbon metabolism imbalance which does not compromise the viability of the double mutant cyfbp cfbp1.

Metabolic profiling of plant extracts using direct-injection electrospray ionization mass spectrometry allows for high-throughput phenotypic characterization according to genetic and environmental effects

In comparison to the exponential increase of genotyping methods, phenotyping strategies are lagging behind in agricultural sciences. Genetic improvement depends upon the abundance of quantitative phenotypic data and the statistical partitioning of variance into environmental, genetic, and random effects. A metabolic phenotyping strategy was adapted to increase sample throughput while saving reagents, reducing cost, and simplifying data analysis. The chemical profiles of stem extracts from maize plants grown under low nitrogen (LN) or control trial (CT) were analyzed using optimized protocols for direct-injection electrospray ionization mass spectrometry (DIESI-MS). Specific ions significantly decreased or increased because of environmental (LN versus CT) or genotypic effects. Biochemical profiling with DIESI-MS had a superior cost-benefit compared to other standard analytical technologies (e.g., ultraviolet, near-infrared reflectance spectroscopy, high-performance liquid chromatography, and gas chromatography with flame ionization detection) routinely used for plant breeding. The method can be successfully applied in maize, strawberry, coffee, and other crop species.

Dissecting the subcellular compartmentation of proteins and metabolites in arabidopsis leaves using non-aqueous fractionation

Non-aqueous fractionation is a technique for the enrichment of different subcellular compartments derived from lyophilized material. It was developed to study the subcellular distribution of metabolites. Here we analyzed the distribution of about 1,000 proteins and 70 metabolites, including 22 phosphorylated intermediates in wild-type Arabidopsis rosette leaves, using non-aqueous gradients divided into 12 fractions. Good separation of plastidial, cytosolic, and vacuolar metabolites and proteins was achieved, but cytosolic, mitochondrial, and peroxisomal proteins clustered together. There was considerable heterogeneity in the fractional distribution of transcription factors, ribosomal proteins, and subunits of the vacuolar-ATPase, indicating diverse compartmental location. Within the plastid, sub-organellar separation of thylakoids and stromal proteins was observed. Metabolites from the Calvin-Benson cycle, photorespiration, starch and sucrose synthesis, glycolysis, and the tricarboxylic acid cycle grouped with their associated proteins of the respective compartment. Non-aqueous fractionation thus proved to be a powerful method for the study of the organellar, and in some cases sub-organellar, distribution of proteins and their association with metabolites. It remains the technique of choice for the assignment of subcellular location to metabolites in intact plant tissues, and thus the technique of choice for doing combined metabolite-protein analysis on a single tissue sample.

Analysis of isotopic labeling in peptide fragments by tandem mass spectrometry

Phenotype in multicellular organisms is the consequence of dynamic metabolic events that occur in a spatially dependent fashion. This spatial and temporal complexity presents challenges for investigating metabolism; creating a need for improved methods that effectively probe biochemical events such as amino acid biosynthesis. Isotopic labeling can provide a temporal-spatial recording of metabolic events through, for example, the description of enriched amino acids in the protein pool. Proteins are therefore an important readout of metabolism and can be assessed with modern mass spectrometers. We compared the measurement of isotopic labeling in MS2 spectra obtained from tandem mass spectrometry under either higher energy collision dissociation (HCD) or collision induced dissociation (CID) at varied energy levels. Developing soybean embryos cultured with or without 13C-labeled substrates, and Escherichia coli MG1655 enriched by feeding 7% uniformly labeled glucose served as a source of biological material for protein evaluation. CID with low energies resulted in a disproportionate amount of heavier isotopologues remaining in the precursor isotopic distribution. HCD resulted in fewer quantifiable products; however deviation from predicted distributions were small relative to the CID-based comparisons. Fragment ions have the potential to provide information on the labeling of amino acids in peptides, but our results indicate that without further development the use of this readout in quantitative methods such as metabolic flux analysis is limited.

Isotopically nonstationary MFA (INST-MFA) of autotrophic metabolism

Metabolic flux analysis (MFA) is a powerful approach for quantifying plant central carbon metabolism based upon a combination of extracellular flux measurements and intracellular isotope labeling measurements. In this chapter, we present the method of isotopically nonstationary (13)C MFA (INST-MFA), which is applicable to autotrophic systems that are at metabolic steady state but are sampled during the transient period prior to achieving isotopic steady state following the introduction of (13)CO2. We describe protocols for performing the necessary isotope labeling experiments, sample collection and quenching, nonaqueous fractionation and extraction of intracellular metabolites, and mass spectrometry (MS) analysis of metabolite labeling. We also outline the steps required to perform computational flux estimation using INST-MFA. By combining several recently developed experimental and computational techniques, INST-MFA provides an important new platform for mapping carbon fluxes that is especially applicable to autotrophic organisms, which are not amenable to steady-state (13)C MFA experiments.

Arabidopsis tic62 trol mutant lacking thylakoid-bound ferredoxin-NADP+ oxidoreductase shows distinct metabolic phenotype

Ferredoxin-NADP+ oxidoreductase (FNR), functioning in the last step of the photosynthetic electron transfer chain, exists both as a soluble protein in the chloroplast stroma and tightly attached to chloroplast membranes. Surface plasmon resonance assays showed that the two FNR isoforms, LFNR1 and LFNR2, are bound to the thylakoid membrane via the C-terminal domains of Tic62 and TROL proteins in a pH-dependent manner. The tic62 trol double mutants contained a reduced level of FNR, exclusively found in the soluble stroma. Although the mutant plants showed no visual phenotype or defects in the function of photosystems under any conditions studied, a low ratio of NADPH/NADP+ was detected. Since the CO₂ fixation capacity did not differ between the tic62 trol plants and wild-type, it seems that the plants are able to funnel reducing power to most crucial reactions to ensure survival and fitness of the plants. However, the activity of malate dehydrogenase was down-regulated in the mutant plants. Apparently, the plastid metabolism is able to cope with substantial changes in directing the electrons from the light reactions to stromal metabolism and thus only few differences are visible in steady-state metabolite pool sizes of the tic62 trol plants.

Isotope labelling of Rubisco subunits provides in vivo information on subcellular biosynthesis and exchange of amino acids between compartments

The architecture of plant metabolism includes substantial duplication of metabolite pools and enzyme catalyzed reactions in different subcellular compartments. This poses challenges for understanding the regulation of metabolism particularly in primary metabolism and amino acid biosynthesis. To explore the extent to which amino acids are made in single compartments and to gain insight into the metabolic precursors from which they derive, we used steady state (13) C labelling and analysed labelling in protein amino acids from plastid and cytosol. Ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) is a major component of green tissues and its large and small subunits are synthesized from different pools of amino acids in the plastid and cytosol, respectively. Developing Brassica napus embryos were cultured in the presence of [U-(13) C]-sucrose, [U-(13) C]-glucose, [U-(13) C]-glutamine or [U-(13) C]-alanine to generate proteins. The large subunits (LSU) and small subunits (SSU) of Rubisco were isolated and the labelling in their constituent amino acids was analysed by gas chromatography-mass spectrometry. Amino acids including alanine, glycine and serine exhibited different (13) C enrichment in the LSU and SSU, demonstrating that these pools have different metabolic origins and are not isotopically equilibrated between the plastid and cytosol on the time scale of cellular growth. Potential extensions of this novel approach to other macromolecules, organelles and cell types of eukaryotes are discussed.

Subcellular analysis of starch metabolism in developing barley seeds using a non-aqueous fractionation method

Compartmentation of metabolism in developing seeds is poorly understood due to the lack of data on metabolite distributions at the subcellular level. In this report, a non-aqueous fractionation method is described that allows subcellular concentrations of metabolites in developing barley endosperm to be calculated. (i) Analysis of subcellular volumes in developing endosperm using micrographs shows that plastids and cytosol occupy 50.5% and 49.9% of the total cell volume, respectively, while vacuoles and mitochondria can be neglected. (ii) By using non-aqueous fractionation, subcellular distribution between the cytosol and plastid of the levels of metabolites involved in sucrose degradation, starch synthesis, and respiration were determined. With the exception of ADP and AMP which were mainly located in the plastid, most other metabolites of carbon and energy metabolism were mainly located outside the plastid in the cytosolic compartment. (iii) In developing barley endosperm, the ultimate precursor of starch, ADPglucose (ADPGlc), was mainly located in the cytosol (80-90%), which was opposite to the situation in growing potato tubers where ADPGlc was almost exclusively located in the plastid (98%). This reflects the different subcellular distribution of ADPGlc pyrophosphorylase (AGPase) in these tissues. (iv) Cytosolic concentrations of ADPGlc were found to be close to the published K(m) values of AGPase and the ADPGlc/ADP transporter at the plastid envelope. Also the concentrations of the reaction partners glucose-1-phosphate, ATP, and inorganic pyrophosphate were close to the respective K(m) values of AGPase. (v) Knock-out of cytosolic AGPase in Riso16 mutants led to a strong decrease in ADPGlc level, in both the cytosol and plastid, whereas knock-down of the ADPGlc/ADP transporter led to a large shift in the intracellular distribution of ADPGlc. (v) The thermodynamic structure of the pathway of sucrose to starch was determined by calculating the mass-action ratios of all the steps in the pathway. The data show that AGPase is close to equilibrium, in both the cytosol and plastid, whereas the ADPGlc/ADP transporter is strongly displaced from equilibrium in vivo. This is in contrast to most other tissues, including leaves and potato tubers. (vi) Results indicate transport rather than synthesis of ADPGlc to be the major regulatory site of starch synthesis in barley endosperm. The reversibility of AGPase in the plastid has important implications for the regulation of carbon partitioning between different biosynthetic pathways.

Conducting molecular biomarker discovery studies in plants

Molecular biomarkers are molecules whose concentrations in a biological system inform about the current phenotypical state and, more importantly, may also be predictive of future phenotypic trait endpoints. The identification of biomarkers has gained much attention in targeted plant breeding since technologies have become available that measure many molecules across different levels of molecular organization and at decreasing costs. In this chapter, we outline the general strategy and workflow of conducting biomarker discovery studies. Critical aspects of study design as well as the statistical data analysis and model building will be highlighted.