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

Remodelled ribosomal populations synthesize a specific proteome in proliferating plant tissue during cold

Plant acclimation occurs through system-wide mechanisms that include proteome shifts, some of which occur at the level of protein synthesis. All proteins are synthesized by ribosomes. Rather than being monolithic, transcript-to-protein translation machines, ribosomes can be selective and cause proteome shifts. In this study, we use apical root meristems of germinating seedlings of the monocotyledonous plant barley as a model to examine changes in protein abundance and synthesis during cold acclimation. We measured metabolic and physiological parameters that allowed us to compare protein synthesis in the cold to optimal rearing temperatures. We demonstrated that the synthesis and assembly of ribosomal proteins are independent processes in root proliferative tissue. We report the synthesis and accumulation of various macromolecular complexes and propose how ribosome compositional shifts may be associated with functional proteome changes that are part of successful cold acclimation. Our study indicates that translation initiation is limiting during cold acclimation while the ribosome population is remodelled. The distribution of the triggered ribosomal protein heterogeneity suggests that altered compositions may confer 60S subunits selective association capabilities towards translation initiation complexes. To what extent selective translation depends on heterogeneous ribo-proteome compositions in barley proliferative root tissue remains a yet unresolved question.This article is part of the discussion meeting issue 'Ribosome diversity and its impact on protein synthesis, development and disease'.

Dynamic photosynthetic labeling and carbon-positional mass spectrometry monitor in vivo RUBISCO carbon assimilation rates

RIBULOSE-1,5-BISPHOSPHATE CARBOXYLASE/OXYGENASE (RUBISCO) is the most abundant enzyme and CO2 bio-sequestration system on Earth. Its in vivo activity is usually determined by 14CO2 incorporation into 3-phosphoglycerate (3PGA). However, the radiometric analysis of 3PGA does not distinguish carbon positions. Hence, RUBISCO activity that fixes carbon into the 1-C position of 3PGA and Calvin-Benson-Bassham (CBB) cycle activities that redistribute carbon into its 2-C and 3-C positions are not resolved. This study aims to develop technology that differentiates between these activities. In source fragmentation of gas chromatography-mass spectrometry (GC-MS) enables paired isotopologue distribution analyses of fragmented substructures and the complete metabolite structure. GC-MS measurements after dynamic photosynthetic 13CO2 labeling allowed quantification of the 13C fractional enrichment (E13C) and molar carbon assimilation rates (A13C) at carbon position 1-C of 3PGA by combining E13C from carbon positions 2,3-C2 and 1,2,3-C3 with quantification of 3PGA concentrations. We validated the procedure using two GC-time of flight-MS instruments, operated at nominal or high mass resolution, and tested the expected 3PGA positional labeling by in vivo glycolysis of positional labeled glucose isotopomers. Mutant analysis of the highly divergent GLYCERALDEHYDE-3-PHOSPHATE DEHYDROGENASEs (GAPDH1 and 2) from Synechocystis sp. PCC 6803 revealed full inactivation of the CBB cycle with maintained RUBISCO activity in Δgapdh2 and a CBB cycle modulating role of GAPDH1 under fluctuating CO2 supply. RUBISCO activity in the CBB-deficient Δgapdh2 can re-assimilate CO2 released by catabolic pathways. We suggest that RUBISCO activity in Synechocystis can scavenge carbon lost through the pentose phosphate pathway or other cellular decarboxylation reactions.

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.

Specific CBF transcription factors and cold-responsive genes fine-tune the early triggering response after acquisition of cold priming and memory

Plants need to adapt to fluctuating temperatures throughout their lifetime. Previous research showed that Arabidopsis memorizes a first cold stress (priming) and improves its primed freezing tolerance further when subjected to a second similar stress after a lag phase. This study investigates primary metabolomic and transcriptomic changes during early cold priming or triggering after 3 days at 4°C interrupted by a memory phase. DREB1 family transcription factors DREB1C/CBF2, DREB1D/CBF4, DREB1E/DDF2, and DREB1F/DDF1 were strongly significantly induced throughout the entire triggering. During triggering, genes encoding Late Embryogenesis Abundant (LEA), antifreeze proteins or detoxifiers of reactive oxygen species (ROS) were higher expressed compared with priming. Examples of early triggering responders were xyloglucan endotransglucosylase/hydrolase genes encoding proteins involved in cell wall remodeling, while late responders were identified to act in fine-tuning the stress response and developmental regulation. Induction of non-typical members of the DREB subfamily of ERF/AP2 transcription factors, the relatively small number of induced CBF regulon genes and a slower accumulation of selected cold stress associated metabolites indicate that a cold triggering stimulus might be sensed as milder stress in plants compared with priming. Further, strong induction of CBF4 throughout triggering suggests a unique function of this gene for the response to alternating temperatures.

Establishment of a GC-MS-based 13 C-positional isotopomer approach suitable for investigating metabolic fluxes in plant primary metabolism

¹³ C-Metabolic flux analysis (¹³ C-MFA) has greatly contributed to our understanding of plant metabolic regulation. However, the generation of detailed in vivo flux maps remains a major challenge. Flux investigations based on nuclear magnetic resonance have resolved small networks with high accuracy. Mass spectrometry (MS) approaches have broader potential, but have hitherto been limited in their power to deduce flux information due to lack of atomic level position information. Herein we established a gas chromatography (GC) coupled to MS-based approach that provides ¹³ C-positional labelling information in glucose, malate and glutamate (Glu). A map of electron impact (EI)-mediated MS fragmentation was created and validated by ¹³ C-positionally labelled references via GC-EI-MS and GC-atmospheric pressure chemical ionization-MS technologies. The power of the approach was revealed by analysing previous ¹³ C-MFA data from leaves and guard cells, and ¹³ C-HCO₃ labelling of guard cells harvested in the dark and after the dark-to-light transition. We demonstrated that the approach is applicable to established GC-EI-MS-based ¹³ C-MFA without the need for experimental adjustment, but will benefit in the future from paired analyses by the two GC-MS platforms. We identified specific glucose carbon atoms that are preferentially labelled by photosynthesis and gluconeogenesis, and provide an approach to investigate the phosphoenolpyruvate carboxylase (PEPc)-derived ¹³ C-incorporation into malate and Glu. Our results suggest that gluconeogenesis and the PEPc-mediated CO₂ assimilation into malate are activated in a light-independent manner in guard cells. We further highlight that the fluxes from glycolysis and PEPc toward Glu are restricted by the mitochondrial thioredoxin system in illuminated leaves.

Metabolic Profiling and Metabolite Correlation Network Analysis Reveal That Fusarium solani Induces Differential Metabolic Responses in Lotus japonicus and Lotus tenuis against Severe Phosphate Starvation

Root fungal endophytes are essential mediators of plant nutrition under mild stress conditions. However, variations in the rhizosphere environment, such as nutrient depletion, could result in a stressful situation for both partners, shifting mutualistic to nonconvenient interactions. Mycorrhizal fungi and dark septate endophytes (DSEs) have demonstrated their ability to facilitate phosphate (Pi) acquisition. However, few studies have investigated other plant-fungal interactions that take place in the root environment with regard to phosphate nutrition. In the present research work, we aimed to analyze the effect of extreme Pi starvation and the fungal endophyte Fusarium solani on the model Lotus japonicus and the crop L. tenuis. We conducted metabolomics analysis based on gas chromatography-mass spectrometry (GC-MS) on plant tissues under optimal conditions, severe Pi starvation and F.solani presence. By combining statistical and correlation network analysis strategies, we demonstrated the differential outcomes of the two plant species against the combination of treatments. The combination of nutritional stress and Fusarium presence activated significant modifications in the metabolism of L. japonicus affecting the levels of sugars, polyols and some amino acids. Our results display potential markers for further inspection of the factors related to plant nutrition and plant-fungal interactions.

Characterization of the Heat-Stable Proteome during Seed Germination in Arabidopsis with Special Focus on LEA Proteins

During seed germination, desiccation tolerance is lost in the radicle with progressing radicle protrusion and seedling establishment. This process is accompanied by comprehensive changes in the metabolome and proteome. Germination of Arabidopsis seeds was investigated over 72 h with special focus on the heat-stable proteome including late embryogenesis abundant (LEA) proteins together with changes in primary metabolites. Six metabolites in dry seeds known to be important for seed longevity decreased during germination and seedling establishment, while all other metabolites increased simultaneously with activation of growth and development. Thermo-stable proteins were associated with a multitude of biological processes. In the heat-stable proteome, a relatively similar proportion of fully ordered and fully intrinsically disordered proteins (IDP) was discovered. Highly disordered proteins were found to be associated with functional categories development, protein, RNA and stress. As expected, the majority of LEA proteins decreased during germination and seedling establishment. However, four germination-specific dehydrins were identified, not present in dry seeds. A network analysis of proteins, metabolites and amino acids generated during the course of germination revealed a highly connected LEA protein network.

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.

The Impact of Metabolic Scion-Rootstock Interactions in Different Grapevine Tissues and Phloem Exudates

In viticulture, grafting is used to propagate Phylloxera-susceptible European grapevines, thereby using resistant American rootstocks. Although scion–rootstock reciprocal signaling is essential for the formation of a proper vascular union and for coordinated growth, our knowledge of graft partner interactions is very limited. In order to elucidate the scale and the content of scion–rootstock metabolic interactions, we profiled the metabolome of eleven graft combination in leaves, stems, and phloem exudate from both above and below the graft union 5–6 months after grafting. We compared the metabolome of scions vs. rootstocks of homografts vs. heterografts and investigated the reciprocal effect of the rootstock on the scion metabolome. This approach revealed that (1) grafting has a minor impact on the metabolome of grafted grapevines when tissues and genotypes were compared, (2) heterografting affects rootstocks more than scions, (3) the presence of a heterologous grafting partner increases defense-related compounds in both scion and rootstocks in shorter and longer distances from the graft, and (4) leaves were revealed as the best tissue to search for grafting-related metabolic markers. These results will provide a valuable metabolomics resource for scion–rootstock interaction studies and will facilitate future efforts on the identification of metabolic markers for important agronomic traits in grafted grapevines.

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.

Cognitive analysis of metabolomics data for systems biology

Cognitive computing is revolutionizing the way big data are processed and integrated, with artificial intelligence (AI) natural language processing (NLP) platforms helping researchers to efficiently search and digest the vast scientific literature. Most available platforms have been developed for biomedical researchers, but new NLP tools are emerging for biologists in other fields and an important example is metabolomics. NLP provides literature-based contextualization of metabolic features that decreases the time and expert-level subject knowledge required during the prioritization, identification and interpretation steps in the metabolomics data analysis pipeline. Here, we describe and demonstrate four workflows that combine metabolomics data with NLP-based literature searches of scientific databases to aid in the analysis of metabolomics data and their biological interpretation. The four procedures can be used in isolation or consecutively, depending on the research questions. The first, used for initial metabolite annotation and prioritization, creates a list of metabolites that would be interesting for follow-up. The second workflow finds literature evidence of the activity of metabolites and metabolic pathways in governing the biological condition on a systems biology level. The third is used to identify candidate biomarkers, and the fourth looks for metabolic conditions or drug-repurposing targets that the two diseases have in common. The protocol can take 1-4 h or more to complete, depending on the processing time of the various software used.

Acclimatisation of guard cell metabolism to long-term salinity

Stomatal movements are enabled by changes in guard cell turgor facilitated via transient accumulation of inorganic and organic ions imported from the apoplast or biosynthesized within guard cells. Under salinity, excess salt ions accumulate within plant tissues resulting in osmotic and ionic stress. To elucidate whether (a) Na⁺ and Cl^- concentrations increase in guard cells in response to long-term NaCl exposure and how (b) guard cell metabolism acclimates to the anticipated stress, we profiled the ions and primary metabolites of leaves, the apoplast and isolated guard cells at darkness and during light, that is, closed and fully opened stomata. In contrast to leaves, the primary metabolism of guard cell preparations remained predominantly unaffected by increased salt ion concentrations. Orchestrated reductions of stomatal aperture and guard cell osmolyte synthesis were found, but unlike in leaves, no increases of stress responsive metabolites or compatible solutes occurred. Diverging regulation of guard cell metabolism might be a prerequisite to facilitate the constant adjustment of turgor that affects aperture. Moreover, the photoperiod-dependent sucrose accumulation in the apoplast and guard cells changed to a permanently replete condition under NaCl, indicating that stress-related photosynthate accumulation in leaves contributes to the permanent closing response of stomata under stress.

Homologs of Circadian Clock Proteins Impact the Metabolic Switch Between Light and Dark Growth in the Cyanobacterium Synechocystis sp. PCC 6803

The putative circadian clock system of the facultative heterotrophic cyanobacterial strain Synechocystis sp. PCC 6803 comprises the following three Kai-based systems: a KaiABC-based potential oscillator that is linked to the SasA-RpaA two-component output pathway and two additional KaiBC systems without a cognate KaiA component. Mutants lacking the genes encoding the KaiAB1C1 components or the response regulator RpaA show reduced growth in light/dark cycles and do not show heterotrophic growth in the dark. In the present study, the effect of these mutations on central metabolism was analyzed by targeted and non-targeted metabolite profiling. The strongest metabolic changes were observed in the dark in ΔrpaA and, to a lesser extent, in the ΔkaiAB1C1 mutant. These observations included the overaccumulation of 2-phosphoglycolate, which correlated with the overaccumulation of the RbcL subunit in the mutants, and taken together, these data suggest enhanced RubisCO activity in the dark. The imbalanced carbon metabolism in the ΔrpaA mutant extended to the pyruvate family of amino acids, which showed increased accumulation in the dark. Hence, the deletion of the response regulator rpaA had a more pronounced effect on metabolism than the deletion of the kai genes. The larger impact of the rpaA mutation is in agreement with previous transcriptomic analyses and likely relates to a KaiAB1C1-independent function as a transcription factor. Collectively, our data demonstrate an important role of homologs of clock proteins in Synechocystis for balanced carbon and nitrogen metabolism during light-to-dark transitions.

Assessing Dynamic Changes of Taste-Related Primary Metabolism During Ripening of Durian Pulp Using Metabolomic and Transcriptomic Analyses

Durian is an economically important fruit of Southeast Asia. There is, however, a lack of in-depth information on the alteration of its metabolic networks during ripening. Here, we annotated 94 ripening-associated metabolites from the pulp of durian cv. Monthong fruit at unripe and ripe stages, using capillary electrophoresis- and gas chromatography- time-of-flight mass spectrometry, specifically focusing on taste-related metabolites. During ripening, sucrose content increased. Change in raffinose-family oligosaccharides are reported herein for the first time. The malate and succinate contents increased, while those of citrate, an abundant organic acid, were unchanged. Notably, most amino acids increased, including isoleucine, leucine, and valine, whereas aspartate decreased, and glutamate was unchanged. Furthermore, transcriptomic analysis was performed to analyze the dynamic changes in sugar metabolism, glycolysis, TCA cycle, and amino acid pathways to identify key candidate genes. Taken together, our results elucidate the fundamental taste-related metabolism of durian, which can be exploited to develop durian metabolic and genetic markers in the future.