Streamlining Protein Fractional Synthesis Rates Using SP3 Beads and Stable Isotope Mass Spectrometry: A Case Study on the Plant Ribosome.

Ribosomes are an archetypal ribonucleoprotein assembly. Due to ribosomal evolution and function, r-proteins share specific physicochemical similarities, making the riboproteome particularly suited for tailored proteome profiling methods. Moreover, the structural proteome of ribonucleoprotein assemblies reflects context-dependent functional features. Thus, characterizing the state of riboproteomes provides insights to uncover the context-dependent functionality of r-protein rearrangements, as they relate to what has been termed the ribosomal code, a concept that parallels that of the histone code, in which chromatin rearrangements influence gene expression. Compared to high-resolution ribosomal structures, omics methods lag when it comes to offering customized solutions to close the knowledge gap between structure and function that currently exists in riboproteomes. Purifying the riboproteome and subsequent shot-gun proteomics typically involves protein denaturation and digestion with proteases. The results are relative abundances of r-proteins at the ribosome population level. We have previously shown that, to gain insight into the stoichiometry of individual proteins, it is necessary to measure by proteomics bound r-proteins and normalize their intensities by the sum of r-protein abundances per ribosomal complex, i.e., 40S or 60S subunits. These calculations ensure that individual r-protein stoichiometries represent the fraction of each family/paralog relative to the complex, effectively revealing which r-proteins become substoichiometric in specific physiological scenarios. Here, we present an optimized method to profile the riboproteome of any organism as well as the synthesis rates of r-proteins determined by stable isotope-assisted mass spectrometry. Our method purifies the r-proteins in a reversibly denatured state, which offers the possibility for combined top-down and bottom-up proteomics. Our method offers a milder native denaturation of the r-proteome via a chaotropic GuHCl solution as compared with previous studies that use irreversible denaturation under highly acidic conditions to dissociate rRNA and r-proteins. As such, our method is better suited to conserve post-translational modifications (PTMs). Subsequently, our method carefully considers the amino acid composition of r-proteins to select an appropriate protease for digestion. We avoid non-specific protease cleavage by increasing the pH of our standardized r-proteome dilutions that enter the digestion pipeline and by using a digestion buffer that ensures an optimal pH for a reliable protease digestion process. Finally, we provide the R package ProtSynthesis to study the fractional synthesis rates of r-proteins. The package uses physiological parameters as input to determine peptide or protein fractional synthesis rates. Once the physiological parameters are measured, our equations allow a fair comparison between treatments that alter the biological equilibrium state of the system under study. Our equations correct peptide enrichment using enrichments in soluble amino acids, growth rates, and total protein accumulation. As a means of validation, our pipeline fails to find "false" enrichments in non-labeled samples while also filtering out proteins with multiple unique peptides that have different enrichment values, which are rare in our datasets. These two aspects reflect the accuracy of our tool. Our method offers the possibility of elucidating individual r-protein family/paralog abundances, PTM status, fractional synthesis rates, and dynamic assembly into ribosomal complexes if top-down and bottom-up proteomic approaches are used concomitantly, taking one step further into mapping the native and dynamic status of the r-proteome onto high-resolution ribosome structures. In addition, our method can be used to study the proteomes of all macromolecular assemblies that can be purified, although purification is the limiting step, and the efficacy and accuracy of the proteases may be limited depending on the digestion requirements. Key features • Efficient purification of the ribosomal proteome: streamlined procedure for the specific purification of the ribosomal proteome or complex Ome. • Accurate calculation of fractional synthesis rates: robust method for calculating fractional protein synthesis rates in macromolecular complexes under different physiological steady states. • Holistic ribosome methodology focused on plants: comprehensive approach that provides insights into the ribosomes and translational control of plants, demonstrated using cold acclimation [1]. • Tailored strategies for stable isotope labeling in plants: methodology focusing on materials and labeling considerations specific to free and proteinogenic amino acid analysis [2].

The serine-glycine-one-carbon metabolic network orchestrates changes in nitrogen and sulfur metabolism and shapes plant development.

L-serine (Ser) and L-glycine (Gly) are critically important for the overall functioning of primary metabolism. We investigated the interaction of the phosphorylated pathway of Ser biosynthesis (PPSB) with the photorespiration-associated glycolate pathway of Ser biosynthesis (GPSB) using Arabidopsis thaliana PPSB-deficient lines, GPSB-deficient mutants, and crosses of PPSB with GPSB mutants. PPSB-deficient lines mainly showed retarded primary root growth. Mutation of the photorespiratory enzyme Ser-hydroxymethyltransferase 1 (SHMT1) in a PPSB-deficient background resumed primary root growth and induced a change in the plant metabolic pattern between roots and shoots. Grafting experiments demonstrated that metabolic changes in shoots were responsible for the changes in double mutant development. PPSB disruption led to a reduction in nitrogen (N) and sulfur (S) contents in shoots and a general transcriptional response to nutrient deficiency. Disruption of SHMT1 boosted the Gly flux out of the photorespiratory cycle, which increased the levels of the one-carbon (1C) metabolite 5,10-methylene-tetrahydrofolate and S-adenosylmethionine. Furthermore, disrupting SHMT1 reverted the transcriptional response to N and S deprivation and increased N and S contents in shoots of PPSB-deficient lines. Our work provides genetic evidence of the biological relevance of the Ser-Gly-1C metabolic network in N and S metabolism and in interorgan metabolic homeostasis.

Dynamics of ribosome composition and ribosomal protein phosphorylation in immune signaling in Arabidopsis thaliana.

In plants, the detection of microbe-associated molecular patterns (MAMPs) induces primary innate immunity by the activation of mitogen-activated protein kinases (MAPKs). We show here that the MAMP-activated MAPK MPK6 not only modulates defense through transcriptional regulation but also via the ribosomal protein translation machinery. To understand the effects of MPK6 on ribosomes and their constituent ribosomal proteins (RPs), polysomes, monosomes and the phosphorylation status of the RPs, MAMP-treated WT and mpk6 mutant plants were analysed. MAMP-activation induced rapid changes in RP composition of monosomes, polysomes and in the 60S ribosomal subunit in an MPK6-specific manner. Phosphoproteome analysis showed that MAMP-activation of MPK6 regulates the phosphorylation status of the P-stalk ribosomal proteins by phosphorylation of RPP0 and the concomitant dephosphorylation of RPP1 and RPP2. These events coincide with a significant decrease in the abundance of ribosome-bound RPP0s, RPP1s and RPP3s in polysomes. The P-stalk is essential in regulating protein translation by recruiting elongation factors. Accordingly, we found that RPP0C mutant plants are compromised in basal resistance to Pseudomonas syringae infection. These data suggest that MAMP-induced defense also involves MPK6-induced regulation of P-stalk proteins, highlighting a new role of ribosomal regulation in plant innate immunity.

Protein interactome of 3',5'-cAMP reveals its role in regulating the actin cytoskeleton.

Identification of protein interactors is ideally suited for the functional characterization of small molecules. 3',5'-cAMP is an evolutionary ancient signaling metabolite largely uncharacterized in plants. To tap into the physiological roles of 3',5'-cAMP, we used a chemo-proteomics approach, thermal proteome profiling (TPP), for the unbiased identification of 3',5'-cAMP protein targets. TPP measures shifts in the protein thermal stability upon ligand binding. Comprehensive proteomics analysis yielded a list of 51 proteins significantly altered in their thermal stability upon incubation with 3',5'-cAMP. The list contained metabolic enzymes, ribosomal subunits, translation initiation factors, and proteins associated with the regulation of plant growth such as CELL DIVISION CYCLE 48. To functionally validate obtained results, we focused on the role of 3',5'-cAMP in regulating the actin cytoskeleton suggested by the presence of actin among the 51 identified proteins. 3',5'-cAMP supplementation affected actin organization by inducing actin-bundling. Consistent with these results, the increase in 3',5'-cAMP levels, obtained either by feeding or by chemical modulation of 3',5'-cAMP metabolism, was sufficient to partially rescue the short hypocotyl phenotype of the actin2 actin7 mutant, severely compromised in actin level. The observed rescue was specific to 3',5'-cAMP, as demonstrated using a positional isomer 2',3'-cAMP, and true for the nanomolar 3',5'-cAMP concentrations reported for plant cells. In vitro characterization of the 3',5'-cAMP-actin pairing argues against a direct interaction between actin and 3',5'-cAMP. Alternative mechanisms by which 3',5'-cAMP would affect actin dynamics, such as by interfering with calcium signaling, are discussed. In summary, our work provides a specific resource, 3',5'-cAMP interactome, as well as functional insight into 3',5'-cAMP-mediated regulation in plants.

Global transcriptome profiling reveals differential regulatory, metabolic and hormonal networks during somatic embryogenesis in Coffea arabica.

BACKGROUND: Somatic embryogenesis (SE) is one of the most promising processes for large-scale dissemination of elite varieties. However, for many plant species, optimizing SE protocols still relies on a trial and error approach. We report the first global scale transcriptome profiling performed at all developmental stages of SE in coffee to unravel the mechanisms that regulate cell fate and totipotency. RESULTS: RNA-seq of 48 samples (12 developmental stages × 4 biological replicates) generated 90 million high quality reads per sample, approximately 74% of which were uniquely mapped to the Arabica genome. First, the statistical analysis of transcript data clearly grouped SE developmental stages into seven important phases (Leaf, Dedifferentiation, Primary callus, Embryogenic callus, Embryogenic cell clusters, Redifferentiation and Embryo) enabling the identification of six key developmental phase switches, which are strategic for the overall biological efficiency of embryo regeneration. Differential gene expression and functional analysis showed that genes encoding transcription factors, stress-related genes, metabolism-related genes and hormone signaling-related genes were significantly enriched. Second, the standard environmental drivers used to control SE, i.e. light, growth regulators and cell density, were clearly perceived at the molecular level at different developmental stages. Third, expression profiles of auxin-related genes, transcription factor-related genes and secondary metabolism-related genes were analyzed during SE. Gene co-expression networks were also inferred. Auxin-related genes were upregulated during dedifferentiation and redifferentiation while transcription factor-related genes were switched on from the embryogenic callus and onward. Secondary metabolism-related genes were switched off during dedifferentiation and switched back on at the onset of redifferentiation. Secondary metabolites and endogenous IAA content were tightly linked with their respective gene expression. Lastly, comparing Arabica embryogenic and non-embryogenic cell transcriptomes enabled the identification of biological processes involved in the acquisition of embryogenic capacity. CONCLUSIONS: The present analysis showed that transcript fingerprints are discriminating signatures of cell fate and are under the direct influence of environmental drivers. A total of 23 molecular candidates were successfully identified overall the 12 developmental stages and can be tested in many plant species to optimize SE protocols in a rational way.

The Effect of Cold Stress on the Root-Specific Lipidome of Two Wheat Varieties with Contrasting Cold Tolerance.

Complex glycerolipidome analysis of wheat upon low temperature stress has been reported for above-ground tissues only. There are no reports on the effects of cold stress on the root lipidome nor on tissue-specific responses of cold stress wheat roots. This study aims to investigate the changes of lipid profiles in the different developmental zones of the seedling roots of two wheat varieties with contrasting cold tolerance exposed to chilling and freezing temperatures. We analyzed 273 lipid species derived from 21 lipid classes using a targeted profiling approach based on MS/MS data acquired from schedule parallel reaction monitoring assays. For both the tolerant Young and sensitive Wyalkatchem species, cold stress increased the phosphatidylcholine and phosphatidylethanolamine compositions, but decreased the monohexosyl ceramide compositions in the root zones. We show that the difference between the two varieties with contrasting cold tolerance could be attributed to the change in the individual lipid species, rather than the fluctuation of the whole lipid classes. The outcomes gained from this study may advance our understanding of the mechanisms of wheat adaptation to cold and contribute to wheat breeding for the improvement of cold-tolerance.

[Formula: see text]: ComplexOme-Structural Network Interpreter used to study spatial enrichment in metazoan ribosomes.

BACKGROUND: Upon environmental stimuli, ribosomes are surmised to undergo compositional rearrangements due to abundance changes among proteins assembled into the complex, leading to modulated structural and functional characteristics. Here, we present the ComplexOme-Structural Network Interpreter ([Formula: see text]), a computational method to allow testing whether ribosomal proteins (rProteins) that exhibit abundance changes under specific conditions are spatially confined to particular regions within the large ribosomal complex. RESULTS: [Formula: see text] translates experimentally determined structures into graphs, with nodes representing proteins and edges the spatial proximity between them. In its first implementation, [Formula: see text] considers rProteins and ignores rRNA and other objects. Spatial regions are defined using a random walk with restart methodology, followed by a procedure to obtain a minimum set of regions that cover all proteins in the complex. Structural coherence is achieved by applying weights to the edges reflecting the physical proximity between purportedly contacting proteins. The weighting probabilistically guides the random-walk path trajectory. Parameter tuning during region selection provides the option to tailor the method to specific biological questions by yielding regions of different sizes with minimum overlaps. In addition, other graph community detection algorithms may be used for the [Formula: see text] workflow, considering that they yield different sized, non-overlapping regions. All tested algorithms result in the same node kernels under equivalent regions. Based on the defined regions, available abundance change information of proteins is mapped onto the graph and subsequently tested for enrichment in any of the defined spatial regions. We applied [Formula: see text] to the cytosolic ribosome structures of Saccharomyces cerevisiae, Oryctolagus cuniculus, and Triticum aestivum using datasets with available quantitative protein abundance change information. We found that in yeast, substoichiometric rProteins depleted from translating polysomes are significantly constrained to a ribosomal region close to the tRNA entry and exit sites. CONCLUSIONS: [Formula: see text] offers a computational method to partition multi-protein complexes into structural regions and a statistical approach to test for spatial enrichments of any given subsets of proteins. [Formula: see text] is applicable to any multi-protein complex given appropriate structural and abundance-change data. [Formula: see text] is publicly available as a GitHub repository https://github.com/MSeidelFed/COSNet_i and can be installed using the python installer pip.

The Metabolic Response of Brachypodium Roots to the Interaction with Beneficial Bacteria Is Affected by the Plant Nutritional Status.

The potential of plant growth promoting (PGP) bacteria in improving the performance of plants in suboptimal environments is increasingly acknowledged, but little information is available on the mechanisms underlying this interaction, particularly when plants are subjected to a combination of stresses. In this study, we investigated the effects of the inoculation with the PGP bacteria Azospirillum brasilense (Azospirillum) on the metabolism of the model cereal Brachypodium distachyon (Brachypodium) grown at low temperatures and supplied with insufficient phosphorus. Investigating polar metabolite and lipid fluctuations during early plant development, we found that the bacteria initially elicited a defense response in Brachypodium roots, while at later stages Azospirillum reduced the stress caused by phosphorus deficiency and improved root development of inoculated plants, particularly by stimulating the growth of branch roots. We propose that the interaction of the plant with Azospirillum was influenced by its nutritional status: bacteria were sensed as pathogens while plants were still phosphorus sufficient, but the interaction became increasingly beneficial for the plants as their phosphorus levels decreased. Our results provide new insights on the dynamics of the cereal-PGP bacteria interaction, and contribute to our understanding of the role of beneficial microorganisms in the growth of cereal crops in suboptimal environments.

Spatially Enriched Paralog Rearrangements Argue Functionally Diverse Ribosomes Arise during Cold Acclimation in Arabidopsis.

Ribosome biogenesis is essential for plants to successfully acclimate to low temperature. Without dedicated steps supervising the 60S large subunits (LSUs) maturation in the cytosol, e.g., Rei-like (REIL) factors, plants fail to accumulate dry weight and fail to grow at suboptimal low temperatures. Around REIL, the final 60S cytosolic maturation steps include proofreading and assembly of functional ribosomal centers such as the polypeptide exit tunnel and the P-Stalk, respectively. In consequence, these ribosomal substructures and their assembly, especially during low temperatures, might be changed and provoke the need for dedicated quality controls. To test this, we blocked ribosome maturation during cold acclimation using two independent reil double mutant genotypes and tested changes in their ribosomal proteomes. Additionally, we normalized our mutant datasets using as a blank the cold responsiveness of a wild-type Arabidopsis genotype. This allowed us to neglect any reil-specific effects that may happen due to the presence or absence of the factor during LSU cytosolic maturation, thus allowing us to test for cold-induced changes that happen in the early nucleolar biogenesis. As a result, we report that cold acclimation triggers a reprogramming in the structural ribosomal proteome. The reprogramming alters the abundance of specific RP families and/or paralogs in non-translational LSU and translational polysome fractions, a phenomenon known as substoichiometry. Next, we tested whether the cold-substoichiometry was spatially confined to specific regions of the complex. In terms of RP proteoforms, we report that remodeling of ribosomes after a cold stimulus is significantly constrained to the polypeptide exit tunnel (PET), i.e., REIL factor binding and functional site. In terms of RP transcripts, cold acclimation induces changes in RP families or paralogs that are significantly constrained to the P-Stalk and the ribosomal head. The three modulated substructures represent possible targets of mechanisms that may constrain translation by controlled ribosome heterogeneity. We propose that non-random ribosome heterogeneity controlled by specialized biogenesis mechanisms may contribute to a preferential or ultimately even rigorous selection of transcripts needed for rapid proteome shifts and successful acclimation.

Membrane-Enriched Proteomics Link Ribosome Accumulation and Proteome Reprogramming With Cold Acclimation in Barley Root Meristems.

Due to their sessile nature, plants rely on root systems to mediate many biotic and abiotic cues. To overcome these challenges, the root proteome is shaped to specific responses. Proteome-wide reprogramming events are magnified in meristems due to their active protein production. Using meristems as a test system, here, we study the major rewiring that plants undergo during cold acclimation. We performed tandem mass tag-based bottom-up quantitative proteomics of two consecutive segments of barley seminal root apexes subjected to suboptimal temperatures. After comparing changes in total and ribosomal protein (RP) fraction-enriched contents with shifts in individual protein abundances, we report ribosome accumulation accompanied by an intricate translational reprogramming in the distal apex zone. Reprogramming ranges from increases in ribosome biogenesis to protein folding factors and suggests roles for cold-specific RP paralogs. Ribosome biogenesis is the largest cellular investment; thus, the vast accumulation of ribosomes and specific translation-related proteins during cold acclimation could imply a divergent ribosomal population that would lead to a proteome shift across the root. Consequently, beyond the translational reprogramming, we report a proteome rewiring. First, triggered protein accumulation includes spliceosome activity in the root tip and a ubiquitous upregulation of glutathione production and S-glutathionylation (S-GSH) assemblage machineries in both root zones. Second, triggered protein depletion includes intrinsically enriched proteins in the tip-adjacent zone, which comprise the plant immune system. In summary, ribosome and translation-related protein accumulation happens concomitantly to a proteome reprogramming in barley root meristems during cold acclimation. The cold-accumulated proteome is functionally implicated in feedbacking transcript to protein translation at both ends and could guide cold acclimation.

Arabidopsis REI-LIKE proteins activate ribosome biogenesis during cold acclimation.

Arabidopsis REIL proteins are cytosolic ribosomal 60S-biogenesis factors. After shift to 10 °C, reil mutants deplete and slowly replenish non-translating eukaryotic ribosome complexes of root tissue, while controlling the balance of non-translating 40S- and 60S-subunits. Reil mutations respond by hyper-accumulation of non-translating subunits at steady-state temperature; after cold-shift, a KCl-sensitive 80S sub-fraction remains depleted. We infer that Arabidopsis may buffer fluctuating translation by pre-existing non-translating ribosomes before de novo synthesis meets temperature-induced demands. Reil1 reil2 double mutants accumulate 43S-preinitiation and pre-60S-maturation complexes and alter paralog composition of ribosomal proteins in non-translating complexes. With few exceptions, e.g. RPL3B and RPL24C, these changes are not under transcriptional control. Our study suggests requirement of de novo synthesis of eukaryotic ribosomes for long-term cold acclimation, feedback control of NUC2 and eIF3C2 transcription and links new proteins, AT1G03250, AT5G60530, to plant ribosome biogenesis. We propose that Arabidopsis requires biosynthesis of specialized ribosomes for cold acclimation.

Separation and Paired Proteome Profiling of Plant Chloroplast and Cytoplasmic Ribosomes.

Conventional preparation methods of plant ribosomes fail to resolve non-translating chloroplast or cytoplasmic ribosome subunits from translating fractions. We established preparation of these ribosome complexes from Arabidopsis thaliana leaf, root, and seed tissues by optimized sucrose density gradient centrifugation of protease protected plant extracts. The method co-purified non-translating 30S and 40S ribosome subunits separated non-translating 50S from 60S subunits, and resolved assembled monosomes from low oligomeric polysomes. Combining ribosome fractionation with microfluidic rRNA analysis and proteomics, we characterized the rRNA and ribosomal protein (RP) composition. The identity of cytoplasmic and chloroplast ribosome complexes and the presence of ribosome biogenesis factors in the 60S-80S sedimentation interval were verified. In vivo cross-linking of leaf tissue stabilized ribosome biogenesis complexes, but induced polysome run-off. Omitting cross-linking, the established paired fractionation and proteome analysis monitored relative abundances of plant chloroplast and cytoplasmic ribosome fractions and enabled analysis of RP composition and ribosome associated proteins including transiently associated biogenesis factors.

Multiplexed Profiling and Data Processing Methods to Identify Temperature-Regulated Primary Metabolites Using Gas Chromatography Coupled to Mass Spectrometry.

This book chapter describes the analytical procedures required for the profiling of a metabolite fraction enriched for primary metabolites. The profiling is based on routine gas chromatography coupled to mass spectrometry (GC-MS). The generic profiling method is adapted to plant material, specifically to the analysis of plant material that was exposed to temperature stress. The method can be combined with stable isotope labeling and tracing experiments and is equally applicable to preparations of plant material and microbial photosynthetic organisms. The described methods are modular and can be multiplexed, that is, the same sample or a paired identical backup sample can be analyzed sequentially by more than one of the described procedures. The modules include rapid sampling and metabolic inactivation protocols for samples in a wide weight range, sample extraction procedures, chemical derivatization steps that are required to make the metabolite fraction amenable to gas chromatographic analysis, routine GC-MS methods, and procedures of data processing and data mining. A basic and extendable set of standardizations for metabolite recovery and retention index alignment of the resulting GC-MS chromatograms is included. The methods have two applications: (1) The rapid screening for changes of relative metabolite pools sizes under temperature stress and (2) the verification by exact quantification using GC-MS protocols that are extended by internal and external standardization.

Systematic Review of Plant Ribosome Heterogeneity and Specialization.

Plants dedicate a high amount of energy and resources to the production of ribosomes. Historically, these multi-protein ribosome complexes have been considered static protein synthesis machines that are not subject to extensive regulation but only read mRNA and produce polypeptides accordingly. New and increasing evidence across various model organisms demonstrated the heterogeneous nature of ribosomes. This heterogeneity can constitute specialized ribosomes that regulate mRNA translation and control protein synthesis. A prominent example of ribosome heterogeneity is seen in the model plant, Arabidopsis thaliana, which, due to genome duplications, has multiple paralogs of each ribosomal protein (RP) gene. We support the notion of plant evolution directing high RP paralog divergence toward functional heterogeneity, underpinned in part by a vast resource of ribosome mutants that suggest specialization extends beyond the pleiotropic effects of single structural RPs or RP paralogs. Thus, Arabidopsis is a highly suitable model to study this phenomenon. Arabidopsis enables reverse genetics approaches that could provide evidence of ribosome specialization. In this review, we critically assess evidence of plant ribosome specialization and highlight steps along ribosome biogenesis in which heterogeneity may arise, filling the knowledge gaps in plant science by providing advanced insights from the human or yeast fields. We propose a data analysis pipeline that infers the heterogeneity of ribosome complexes and deviations from canonical structural compositions linked to stress events. This analysis pipeline can be extrapolated and enhanced by combination with other high-throughput methodologies, such as proteomics. Technologies, such as kinetic mass spectrometry and ribosome profiling, will be necessary to resolve the temporal and spatial aspects of translational regulation while the functional features of ribosomal subpopulations will become clear with the combination of reverse genetics and systems biology approaches.

Unravelling the Metabolic and Hormonal Machinery During Key Steps of Somatic Embryogenesis: A Case Study in Coffee.

Somatic embryogenesis (SE) is one of the most promising processes for large-scale dissemination of elite varieties. However, for many plant species, optimizing SE protocols still relies on a trial-and-error approach. Using coffee as a model plant, we report here the first global analysis of metabolome and hormone dynamics aiming to unravel mechanisms regulating cell fate and totipotency. Sampling from leaf explant dedifferentiation until embryo development covered 15 key stages. An in-depth statistical analysis performed on 104 metabolites revealed that massive re-configuration of metabolic pathways induced SE. During initial dedifferentiation, a sharp decrease in phenolic compounds and caffeine levels was also observed while auxins, cytokinins and ethylene levels were at their highest. Totipotency reached its highest expression during the callus stages when a shut-off in hormonal and metabolic pathways related to sugar and energetic substance hydrolysis was evidenced. Abscisic acid, leucine, maltotriose, myo-inositol, proline, tricarboxylic acid cycle metabolites and zeatin appeared as key metabolic markers of the embryogenic capacity. Combining metabolomics with multiphoton microscopy led to the identification of chlorogenic acids as markers of embryo redifferentiation. The present analysis shows that metabolite fingerprints are signatures of cell fate and represent a starting point for optimizing SE protocols in a rational way.

Discovery of food identity markers by metabolomics and machine learning technology.

Verification of food authenticity establishes consumer trust in food ingredients and components of processed food. Next to genetic or protein markers, chemicals are unique identifiers of food components. Non-targeted metabolomics is ideally suited to screen food markers when coupled to efficient data analysis. This study explored feasibility of random forest (RF) machine learning, specifically its inherent feature extraction for non-targeted metabolic marker discovery. The distinction of chia, linseed, and sesame that have gained attention as "superfoods" served as test case. Chemical fractions of non-processed seeds and of wheat cookies with seed ingredients were profiled. RF technology classified original seeds unambiguously but appeared overdesigned for material with unique secondary metabolites, like sesamol or rosmarinic acid in the Lamiaceae, chia. Most unique metabolites were diluted or lost during cookie production but RF technology classified the presence of the seed ingredients in cookies with 6.7% overall error and revealed food processing markers, like 4-hydroxybenzaldehyde for chia and succinic acid monomethylester for linseed additions. RF based feature extraction was adequate for difficult classifications but marker selection should not be without human supervision. Combination with alternative data analysis technologies is advised and further testing of a wide range of seeds and food processing methods.