Multi-omics analysis of green lineage osmotic stress pathways unveils crucial roles of different cellular compartments.

Maintenance of water homeostasis is a fundamental cellular process required by all living organisms. Here, we use the single-celled green alga Chlamydomonas reinhardtii to establish a foundational understanding of osmotic-stress signaling pathways through transcriptomics, phosphoproteomics, and functional genomics approaches. Comparison of pathways identified through these analyses with yeast and Arabidopsis allows us to infer their evolutionary conservation and divergence across these lineages. 76 genes, acting across diverse cellular compartments, were found to be important for osmotic-stress tolerance in Chlamydomonas through their functions in cytoskeletal organization, potassium transport, vesicle trafficking, mitogen-activated protein kinase and chloroplast signaling. We show that homologs for five of these genes have conserved functions in stress tolerance in Arabidopsis and reveal a novel PROFILIN-dependent stage of acclimation affecting the actin cytoskeleton that ensures tissue integrity upon osmotic stress. This study highlights the conservation of the stress response in algae and land plants, and establishes Chlamydomonas as a unicellular plant model system to dissect the osmotic stress signaling pathway.

Radical-Mediated Covalent Azidylation of Hydrophobic Microdomains in Water-Soluble Proteins.

Hydrophobic microdomains, also known as hydrophobic patches, are essential for many important biological functions of water-soluble proteins. These include ligand or substrate binding, protein-protein interactions, proper folding after translation, and aggregation during denaturation. Unlike transmembrane domains, which are easily recognized from stretches of contiguous hydrophobic sidechains in amino acids via primary protein sequence, these three-dimensional hydrophobic patches cannot be easily predicted. The lack of experimental strategies for directly determining their locations hinders further understanding of their structure and function. Here, we posit that the small triatomic anion N3- (azide) is attracted to these patches and, in the presence of an oxidant, forms a radical that covalently modifies C-H bonds of nearby amino acids. Using two model proteins (BSA and lysozyme) and a cell-free lysate from the model higher plant Arabidopsis thaliana, we find that radical-mediated covalent azidylation occurs within buried catalytic active sites and ligand binding sites and exhibits similar behavior to established hydrophobic probes. The results herein suggest a model in which the azido radical is acting as an "affinity reagent" for nonaqueous three-dimensional protein microenvironments and is consistent with both the nonlocalized electron density of the azide moiety and the known high reactivity of azido radicals widely used in organic chemistry syntheses. We propose that the azido radical is a facile means of identifying hydrophobic microenvironments in soluble proteins and, in addition, provides a simple new method for attaching chemical handles to proteins without the need for genetic manipulation or specialized reagents.

Mass spectrometric based analysis of whole eggs dissolved in formic acid.

We have developed a method for complete dissolution of whole eggs in formic acid that provides a new approach to analyzing egg biomolecules. As expected from prior work with extracted lipids, phosphatidylcholine represents the most abundant 31P NMR signal. A simplified methanol/chloroform partitioning method for separating the dissolved egg solution into metabolites, lipids and protein was performed and after ultra-high mass resolution and tandem MS fragmentation analyses several phosphatidylcholine molecules containing different fatty acid chain lengths as well as number and position of double bonds was detected. The MS based proteomic analysis further revealed 6 Gallus sequences annotated as 'uncharacterized' because they show no sequence homology with any other protein found in nature and thus, may represent proteins uniquely evolved to perform functions specific to chickens. Overall, this procedure is a rapid and facile means of characterizing in a high throughput and comprehensive manner, the molecular components of whole eggs.

Mass spectrometry-based technologies for probing the 3D world of plant proteins.

Over the past two decades, mass spectrometric (MS)-based proteomics technologies have facilitated the study of signaling pathways throughout biology. Nowhere is this needed more than in plants, where an evolutionary history of genome duplications has resulted in large gene families involved in posttranslational modifications and regulatory pathways. For example, at least 5% of the Arabidopsis thaliana genome (ca. 1,200 genes) encodes protein kinases and protein phosphatases that regulate nearly all aspects of plant growth and development. MS-based technologies that quantify covalent changes in the side-chain of amino acids are critically important, but they only address one piece of the puzzle. A more crucially important mechanistic question is how noncovalent interactions-which are more difficult to study-dynamically regulate the proteome's 3D structure. The advent of improvements in protein 3D technologies such as cryo-electron microscopy, nuclear magnetic resonance, and X-ray crystallography has allowed considerable progress to be made at this level, but these methods are typically limited to analyzing proteins, which can be expressed and purified in milligram quantities. Newly emerging MS-based technologies have recently been developed for studying the 3D structure of proteins. Importantly, these methods do not require protein samples to be purified and require smaller amounts of sample, opening the wider proteome for structural analysis in complex mixtures, crude lysates, and even in intact cells. These MS-based methods include covalent labeling, crosslinking, thermal proteome profiling, and limited proteolysis, all of which can be leveraged by established MS workflows, as well as newly emerging methods capable of analyzing intact macromolecules and the complexes they form. In this review, we discuss these recent innovations in MS-based "structural" proteomics to provide readers with an understanding of the opportunities they offer and the remaining challenges for understanding the molecular underpinnings of plant structure and function.

Proteome Damage Inflicted by Ionizing Radiation: Advancing a Theme in the Research of Miroslav Radman.

Oxidative proteome damage has been implicated as a major contributor to cell death and aging. Protein damage and aging has been a particular theme of the recent research of Miroslav Radman. However, the study of how cellular proteins are damaged by oxidative processes is still in its infancy. Here we examine oxidative changes in the proteomes of four bacterial populations-wild type E. coli, two isolates from E. coli populations evolved for high levels of ionizing radiation (IR) resistance, and D. radiodurans-immediately following exposure to 3000 Gy of ionizing radiation. By a substantial margin, the most prominent intracellular oxidation events involve hydroxylation of methionine residues. Significant but much less frequent are carbonylation events on tyrosine and dioxidation events on tryptophan. A few proteins are exquisitely sensitive to targeted oxidation events, notably the active site of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) in E. coli. Extensive experimental evolution of E. coli for IR resistance has decreased overall proteome sensitivity to oxidation but not to the level seen in D. radiodurans. Many observed oxidation events may reflect aspects of protein structure and/or exposure of protein surfaces to water. Proteins such as GAPDH and possibly Ef-Tu may have an evolved sensitivity to oxidation by H2O2.

Corrigendum: Physiology of Highly Radioresistant Escherichia coli After Experimental Evolution for 100 Cycles of Selection.

[This corrects the article DOI: 10.3389/fmicb.2020.582590.].

Physiology of Highly Radioresistant Escherichia coli After Experimental Evolution for 100 Cycles of Selection.

Ionizing radiation (IR) is lethal to most organisms at high doses, damaging every cellular macromolecule via induction of reactive oxygen species (ROS). Utilizing experimental evolution and continuing previous work, we have generated the most IR-resistant Escherichia coli populations developed to date. After 100 cycles of selection, the dose required to kill 99% the four replicate populations (IR9-100, IR10-100, IR11-100, and IR12-100) has increased from 750 Gy to approximately 3,000 Gy. Fitness trade-offs, specialization, and clonal interference are evident. Long-lived competing sub-populations are present in three of the four lineages. In IR9, one lineage accumulates the heme precursor, porphyrin, leading to generation of yellow-brown colonies. Major genomic alterations are present. IR9 and IR10 exhibit major deletions and/or duplications proximal to the chromosome replication terminus. Contributions to IR resistance have expanded beyond the alterations in DNA repair systems documented previously. Variants of proteins involved in ATP synthesis (AtpA), iron-sulfur cluster biogenesis (SufD) and cadaverine synthesis (CadA) each contribute to IR resistance in IR9-100. Major genomic and physiological changes are emerging. An isolate from IR10 exhibits protein protection from ROS similar to the extremely radiation resistant bacterium Deinococcus radiodurans, without evident changes in cellular metal homeostasis. Selection is continuing with no limit to IR resistance in evidence as our E. coli populations approach levels of IR resistance typical of D. radiodurans.

Ionizing Radiation-induced Proteomic Oxidation in Escherichia coli.

Recent work has begun to investigate the role of protein damage in cell death because of ionizing radiation (IR) exposure, but none have been performed on a proteome-wide basis, nor have they utilized MS (MS) to determine chemical identity of the amino acid side chain alteration. Here, we use Escherichia coli to perform the first MS analysis of IR-treated intact cells on a proteome scale. From quintuplicate IR-treated (1000 Gy) and untreated replicates, we successfully quantified 13,262 peptides mapping to 1938 unique proteins. Statistically significant, but low in magnitude (<2-fold), IR-induced changes in peptide abundance were observed in 12% of all peptides detected, although oxidative alterations were rare. Hydroxylation (+15.99 Da) was the most prevalent covalent adduct detected. In parallel with these studies on E. coli, identical experiments with the IR-resistant bacterium, Deinococcus radiodurans, revealed orders of magnitude less effect of IR on the proteome. In E. coli, the most significant target of IR by a wide margin was glyceraldehyde 3'-phosphate dehydrogenase (GAPDH), in which the thiol side chain of the catalytic Cys residue was oxidized to sulfonic acid. The same modification was detected in IR-treated human breast carcinoma cells. Sensitivity of GAPDH to reactive oxygen species (ROS) has been described previously in microbes and here, we present GAPDH as an immediate, primary target of IR-induced oxidation across all domains of life.

Proteomic Characterization of Lignocellulolytic Enzymes Secreted by the Insect-Associated Fungus Daldinia decipiens oita, Isolated from a Forest in Northern Japan.

Wood-devastating insects utilize their symbiotic microbes with lignocellulose-degrading abilities to extract energy from recalcitrant woods. It is well known that free-living lignocellulose-degrading fungi secrete various carbohydrate-active enzymes (CAZymes) to degrade plant cell wall components, mainly cellulose, hemicellulose, and lignin. However, CAZymes from insect-symbiotic fungi have not been well documented except for a few examples. In this study, an insect-associated fungus, Daldinia decipiens oita, was isolated as a potential symbiotic fungus of female Xiphydria albopicta captured from Hokkaido forest. This fungus was grown in seven different media containing a single carbon source, glucose, cellulose, xylan, mannan, pectin, poplar, or larch, and the secreted proteins were identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS). A total of 128 CAZymes, including domains of 92 glycoside hydrolases, 15 carbohydrate esterases, 5 polysaccharide lyases, 17 auxiliary activities, and 11 carbohydrate-binding modules, were identified, and these are involved in degradation of cellulose and hemicellulose but not lignin. Together with the results of polysaccharide-degrading activity measurements, we concluded that D. decipiens oita tightly regulates the expression of these CAZymes in response to the tested plant cell wall materials. Overall, this study described the detailed proteomic approach of a woodwasp-associated fungus and revealed that the new isolate, D. decipiens oita, secretes diverse CAZymes to efficiently degrade lignocellulose in the symbiotic environment.IMPORTANCE Recent studies show the potential impacts of insect symbiont microbes on biofuel application with regard to their degradation capability of a recalcitrant plant cell wall. In this study, we describe a novel fungal isolate, D. decipiens oita, as a single symbiotic fungus from the Xiphydria woodwasp found in the northern forests of Japan. Our detailed secretome analyses of D. decipiens oita, together with activity measurements, reveal that this insect-associated fungus exhibits high and broad activities for plant cell wall material degradation, suggesting potential applications within the biomass conversion industry for plant mass degradation.

Covalent Modification of Amino Acids and Peptides Induced by Ionizing Radiation from an Electron Beam Linear Accelerator Used in Radiotherapy.

To identify modifications to amino acids that are directly induced by ionizing radiation, free amino acids and 3-residue peptides were irradiated using a linear accelerator (Linac) radiotherapy device. Mass spectrometry was performed to detail the relative sensitivity to radiation as well as identify covalent, radiation-dependent adducts. The order of reactivity of the 20 common amino acids was generally in agreement with published literature except for His (most reactive of the 20) and Cys (less reactive). Novel and previously identified modifications on the free amino acids were detected. Amino acids were far less reactive when flanked by glycine residues in a tripeptide. Order of reactivity, with GVG most and GEG least, was substantially altered, as were patterns of modification. Radiation reactivity of amino acids is clearly and strongly affected by conversion of the α-amino and α-carboxyl groups to peptide bonds, and the presence of neighboring amino acid residues.

Probing a Plant Plasma Membrane Receptor Kinase's Three-Dimensional Structure Using Mass Spectrometry-Based Protein Footprinting.

FERONIA (FER), one of the 17 malectin-like receptor-like kinases encoded in the Arabidopsis genome, acts as a receptor for a 5 kDa growth-inhibiting secreted protein hormone, rapid alkalinization factor 1 (RALF1). Upon binding the peptide ligand, FER is involved in a variety of signaling pathways eliciting ovule fertilization and vegetative root cell expansion. Here, we report the use of mass spectrometry-based, carbodiimide-mediated protein carboxyl group (aspartic and glutamic acid) footprinting to map solvent accessible amino acids of the ectodomain of FER (ectoFER), including those involved in RALF1 binding and/or allosteric changes. Aspartate and glutamate residues labeled in this procedure were located in various regions, including the N-terminus, malectin-like domains, and juxtamembrane region, and these correlated well with a three-dimensional structural model of ectoFER predicted from the crystal structure of a related receptor. Covalent cross-linking experiments also revealed the N-terminus of ectoFER linked to the highly conserved C-terminus of RALF1. RALF1 binding assays performed with truncation mutants of ectoFER further implicated the receptor N-terminal and juxtamembrane regions in the binding of RALF1. In conclusion, our results of mass spectrometry-based footprinting methods provide a framework for understanding ligand-induced changes in solvent accessibility and their positions within the three-dimensional structure of a plant receptor kinase.

Plasma-Generated OH Radical Production for Analyzing Three-Dimensional Structure in Protein Therapeutics.

Protein three-dimensional structure dynamically changes in solution depending on the presence of ligands and interacting proteins. Methods for detecting these changes in protein conformation include 'protein footprinting,' using mass spectrometry. We describe herein a new technique, PLIMB (Plasma Induced Modification of Biomolecules), that generates µs bursts of hydroxyl radicals from water, to measure changes in protein structure via altered solvent accessibility of amino acid side chains. PLIMB was first benchmarked with model compounds, and then applied to a biological problem, i.e., ligand (EGF) induced changes in the conformation of the external (ecto) domain of Epidermal Growth Factor Receptor (EGFR). Regions in which oxidation decreased upon adding EGF fall along the dimerization interface, consistent with models derived from crystal structures. These results demonstrate that plasma-generated hydroxyl radicals from water can be used to map protein conformational changes, and provide a readily accessible means of studying protein structure in solution.

A cell-free method for expressing and reconstituting membrane proteins enables functional characterization of the plant receptor-like protein kinase FERONIA.

There are more than 600 receptor-like kinases (RLKs) in Arabidopsis, but due to challenges associated with the characterization of membrane proteins, only a few have known biological functions. The plant RLK FERONIA is a peptide receptor and has been implicated in plant growth regulation, but little is known about its molecular mechanism of action. To investigate the properties of this enzyme, we used a cell-free wheat germ-based expression system in which mRNA encoding FERONIA was co-expressed with mRNA encoding the membrane scaffold protein variant MSP1D1. With the addition of the lipid cardiolipin, assembly of these proteins into nanodiscs was initiated. FERONIA protein kinase activity in nanodiscs was higher than that of soluble protein and comparable with other heterologously expressed protein kinases. Truncation experiments revealed that the cytoplasmic juxtamembrane domain is necessary for maximal FERONIA activity, whereas the transmembrane domain is inhibitory. An ATP analogue that reacts with lysine residues inhibited catalytic activity and labeled four lysines; mutagenesis demonstrated that two of these, Lys-565 and Lys-663, coordinate ATP in the active site. Mass spectrometric phosphoproteomic measurements further identified phosphorylation sites that were examined using phosphomimetic mutagenesis. The results of these experiments are consistent with a model in which kinase-mediated phosphorylation within the C-terminal region is inhibitory and regulates catalytic activity. These data represent a step further toward understanding the molecular basis for the protein kinase catalytic activity of FERONIA and show promise for future characterization of eukaryotic membrane proteins.

Rapid Oligo-Galacturonide Induced Changes in Protein Phosphorylation in Arabidopsis.

The wall-associated kinases (WAKs)(1)are receptor protein kinases that bind to long polymers of cross-linked pectin in the cell wall. These plasma-membrane-associated protein kinases also bind soluble pectin fragments called oligo-galacturonides (OGs) released from the wall after pathogen attack and damage. WAKs are required for cell expansion during development but bind water soluble OGs generated from walls with a higher affinity than the wall-associated polysaccharides. OGs activate a WAK-dependent, distinct stress-like response pathway to help plants resist pathogen attack. In this report, a quantitative mass-spectrometric-based phosphoproteomic analysis was used to identify Arabidopsis cellular events rapidly induced by OGsin planta Using N(14/)N(15)isotopicin vivometabolic labeling, we screened 1,000 phosphoproteins for rapid OG-induced changes and found 50 proteins with increased phosphorylation, while there were none that decreased significantly. Seven of the phosphosites within these proteins overlap with those altered by another signaling molecule plants use to indicate the presence of pathogens (the bacterial "elicitor" peptide Flg22), indicating distinct but overlapping pathways activated by these two types of chemicals. Genetic analysis of genes encoding 10 OG-specific and two Flg22/OG-induced phosphoproteins reveals that null mutations in eight proteins compromise the OG response. These phosphorylated proteins with genetic evidence supporting their role in the OG response include two cytoplasmic kinases, two membrane-associated scaffold proteins, a phospholipase C, a CDPK, an unknown cadmium response protein, and a motor protein. Null mutants in two proteins, the putative scaffold protein REM1.3, and a cytoplasmic receptor like kinase ROG2, enhance and suppress, respectively, a dominantWAKallele. Altogether, the results of these chemical and genetic experiments reveal the identity of several phosphorylated proteins involved in the kinase/phosphatase-mediated signaling pathway initiated by cell wall changes.

Functional characterization of PCRK1, a putative protein kinase with a role in immunity.

In Arabidopsis, defense signaling is triggered by the perception of conserved molecular patterns by pattern recognition receptors (PRRs). Signal transduction from the PRRs requires members of a family of Receptor-Like Cytoplasmic Kinases (RLCKs). Previously, we described one such RLCK, PTI Compromised Receptor-Like Cytoplasmic Kinase 1 (PCRK1) that is important for immunity induced by Microbe Associated Molecular Patterns (MAMPs) as well as Damage Associated Molecular Patterns (DAMPs). In this study, we measured the growth of Pma ES4326 in double mutants carrying pcrk1 together with the salicylic acid (SA) biosynthesis mutation sid2-2 or the jasmonic acid (JA) receptor mutation coi1-1, showing that the function of PCRK1 is SA independent but may be partially dependent on JA. Mutation of phosphorylated serine residues S232, S233 and S237 compromised the immune signaling function of PCRK1.

Rapid Phosphoproteomic Effects of Abscisic Acid (ABA) on Wild-Type and ABA Receptor-Deficient A. thaliana Mutants.

Abscisic acid (ABA)¹ is a plant hormone that controls many aspects of plant growth, including seed germination, stomatal aperture size, and cellular drought response. ABA interacts with a unique family of 14 receptor proteins. This interaction leads to the activation of a family of protein kinases, SnRK2s, which in turn phosphorylate substrates involved in many cellular processes. The family of receptors appears functionally redundant. To observe a measurable phenotype, four of the fourteen receptors have to be mutated to create a multilocus loss-of-function quadruple receptor (QR) mutant, which is much less sensitive to ABA than wild-type (WT) plants. Given these phenotypes, we asked whether or not a difference in ABA response between the WT and QR backgrounds would manifest on a phosphorylation level as well. We tested WT and QR mutant ABA response using isotope-assisted quantitative phosphoproteomics to determine what ABA-induced phosphorylation changes occur in WT plants within 5 min of ABA treatment and how that phosphorylation pattern is altered in the QR mutant. We found multiple ABA-induced phosphorylation changes that occur within 5 min of treatment, including three SnRK2 autophosphorylation events and phosphorylation on SnRK2 substrates. The majority of robust ABA-dependent phosphorylation changes observed were partially diminished in the QR mutant, whereas many smaller ABA-dependent phosphorylation changes observed in the WT were not responsive to ABA in the mutant. A single phosphorylation event was increased in response to ABA treatment in both the WT and QR mutant. A portion of the discovery data was validated using selected reaction monitoring-based targeted measurements on a triple quadrupole mass spectrometer. These data suggest that different subsets of phosphorylation events depend upon different subsets of the ABA receptor family to occur. Altogether, these data expand our understanding of the model by which the family of ABA receptors directs rapid phosphoproteomic changes.

Phosphoproteomic Analyses Reveal Early Signaling Events in the Osmotic Stress Response.

Elucidating how plants sense and respond to water loss is important for identifying genetic and chemical interventions that may help sustain crop yields in water-limiting environments. Currently, the molecular mechanisms involved in the initial perception and response to dehydration are not well understood. Modern mass spectrometric methods for quantifying changes in the phosphoproteome provide an opportunity to identify key phosphorylation events involved in this process. Here, we have used both untargeted and targeted isotope-assisted mass spectrometric methods of phosphopeptide quantitation to characterize proteins in Arabidopsis (Arabidopsis thaliana) whose degree of phosphorylation is rapidly altered by hyperosmotic treatment. Thus, protein phosphorylation events responsive to 5 min of 0.3 m mannitol treatment were first identified using 15N metabolic labeling and untargeted mass spectrometry with a high-resolution ion-trap instrument. The results from these discovery experiments were then validated using targeted Selected Reaction Monitoring mass spectrometry with a triple quadrupole. Targeted Selected Reaction Monitoring experiments were conducted with plants treated under nine different environmental perturbations to determine whether the phosphorylation changes were specific for osmosignaling or involved cross talk with other signaling pathways. The results indicate that regulatory proteins such as members of the mitogen-activated protein kinase family are specifically phosphorylated in response to osmotic stress. Proteins involved in 5' messenger RNA decapping and phosphatidylinositol 3,5-bisphosphate synthesis were also identified as targets of dehydration-induced phosphoregulation. The results of these experiments demonstrate the utility of targeted phosphoproteomic analysis in understanding protein regulation networks and provide new insight into cellular processes involved in the osmotic stress response.

Expression of a translationally fused TAP-tagged plasma membrane proton pump in Arabidopsis thaliana.

The Arabidopsis thaliana plasma membrane proton ATPase genes, AHA1 and AHA2, are the two most highly expressed isoforms of an 11 gene family and are collectively essential for embryo development. We report the translational fusion of a tandem affinity-purification tag to the 5' end of the AHA1 open reading frame in a genomic clone. Stable expression of TAP-tagged AHA1 in Arabidopsis rescues the embryonic lethal phenotype of endogenous double aha1/aha2 knockdowns. Western blots of SDS-PAGE and Blue Native gels show enrichment of AHA1 in plasma membrane fractions and indicate a hexameric quaternary structure. TAP-tagged AHA1 rescue lines exhibited reduced vertical root growth. Analysis of the plasma membrane and soluble proteomes identified several plasma membrane-localized proteins with alterred abundance in TAP-tagged AHA1 rescue lines compared to wild type. Using affinity-purification mass spectrometry, we uniquely identified two additional AHA isoforms, AHA9 and AHA11, which copurified with TAP-tagged AHA1. In conclusion, we have generated transgenic Arabidopsis lines in which a TAP-tagged AHA1 transgene has complemented all essential endogenous AHA1 and AHA2 functions and have shown that these plants can be used to purify AHA1 protein and to identify in planta interacting proteins by mass spectrometry.

A peptide hormone and its receptor protein kinase regulate plant cell expansion.

Plant cells are immobile; thus, plant growth and development depend on cell expansion rather than cell migration. The molecular mechanism by which the plasma membrane initiates changes in the cell expansion rate remains elusive. We found that a secreted peptide, RALF (rapid alkalinization factor), suppresses cell elongation of the primary root by activating the cell surface receptor FERONIA in Arabidopsis thaliana. A direct peptide-receptor interaction is supported by specific binding of RALF to FERONIA and reduced binding and insensitivity to RALF-induced growth inhibition in feronia mutants. Phosphoproteome measurements demonstrate that the RALF-FERONIA interaction causes phosphorylation of plasma membrane H(+)-adenosine triphosphatase 2 at Ser(899), mediating the inhibition of proton transport. The results reveal a molecular mechanism for RALF-induced extracellular alkalinization and a signaling pathway that regulates cell expansion.

A pipeline for 15N metabolic labeling and phosphoproteome analysis in Arabidopsis thaliana.

Within the past two decades, the biological application of mass spectrometric technology has seen great advances in terms of innovations in hardware, software, and reagents. Concurrently, the burgeoning field of proteomics has followed closely (Yates et al., Annu Rev Biomed Eng 11:49-79, 2009)-and with it, importantly, the ability to globally assay altered levels of posttranslational modifications in response to a variety of stimuli. Though many posttranslational modifications have been described, a major focus of these efforts has been protein-level phosphorylation of serine, threonine, and tyrosine residues (Schreiber et al., Proteomics 8:4416-4432, 2008). The desire to examine changes across signal transduction cascades and networks in their entirety using a single mass spectrometric analysis accounts for this push-namely, preservation and enrichment of the transient yet informative phosphoryl side group. Analyzing global changes in phosphorylation allows inferences surrounding cascades/networks as a whole to be made. Towards this same end, much work has explored ways to permit quantitation and combine experimental samples such that more than one replicate or experimental condition can be identically processed and analyzed, cutting down on experimental and instrument variability, in addition to instrument run time. One such technique that has emerged is metabolic labeling (Gouw et al., Mol Cell Proteomics 9:11-24, 2010), wherein biological samples are labeled in living cells with nonradioactive heavy isotopes such as (15)N or (13)C. Since metabolic labeling in living organisms allows one to combine the material to be processed at the earliest possible step, before the tissue is homogenized, it provides a unique and excellent method for comparing experimental samples in a high-throughput, reproducible fashion with minimal technical variability. This chapter describes a pipeline used for labeling living Arabidopsis thaliana plants with nitrogen-15 ((15)N) and how this can be used, in conjunction with a technique for enrichment of phosphorylated peptides (phosphopeptides), to determine changes in A. thaliana's phosphoproteome on an untargeted, global scale.