Metabolite Profiling Reveals Predictive Biomarkers and the Absence of ß-Methyl Amino-l-alanine in Plasma from Individuals Diagnosed with Amyotrophic Lateral Sclerosis.

By employing chip-based capillary zone electrophoresis coupled to high-resolution mass spectrometry, we profiled the plasma metabolome of 134 patients diagnosed with sporadic amyotrophic lateral sclerosis (ALS) (81 males and 53 females) and 118 individuals deemed healthy (49 males and 69 females). The most significant markers (p < 0.01) were creatine, which was 49% elevated, and creatinine and methylhistidine, which were decreased by 20 and 24%, respectively, in ALS patients. The ratio of creatine versus creatinine increased 370 and 200% for male and female ALS patients, respectively. In addition, male ALS patients on an average had 5-13% lower amounts of seven essential amino acids, whereas females did not significantly differ from healthy controls. We developed two models using the metabolite abundances: (1) a classification model for the separation of ALS and healthy samples and (2) a classification model for the prediction of disease progression based on the ALS functional rating score. Utilizing a Monte Carlo cross-validation approach, a linear discriminant analysis model achieved a mean area under the receiver operating characteristic curve (AUC) of 0.85 (0.06) with a mean sensitivity of 80% (9%) and specificity of 78% (10%) for the separation of ALS and controls, respectively. A support vector machine classifier predicted progression categories with an AUC of 0.90 (0.06) with a mean sensitivity of 73% (10%) and a specificity of 86% (5%). Lastly, using a previously reported assay with a stable isotope-labeled (13C315N2) spike-in standard, we were unable to detect the exogenous neurotoxic metabolite, ß-methylamino-l-alanine, in the free or protein-bound fraction of any of the 252 plasma samples.

Peptide variability and signatures associated with disease progression in CSF collected longitudinally from ALS patients.

We employ shotgun proteomics and data-independent acquisition (DIA) mass spectrometry to analyze cerebrospinal fluid longitudinally collected from 14 amyotrophic lateral sclerosis (ALS) patients (8 males and 6 females). We perform three main analyses of these data: (1) examine the intra- and inter-patient protein variability in CSF; (2) explore the association of inflammation with rate of disease progression; and (3) develop a mixed-effects model to best explain the decrease in ALS-Functional Rating Scale (ALS-FRS) score. Overall, the CSF protein abundances are tightly regulated with the intra-individual variability contributing just 4% to the overall variance. In four patients, a moderately significant correlation (p < 0.1) was observed between inflammation and rate of disease progression. Using a least absolute shrinkage and selection operator (LASSO) variable selection, we selected 55 viable peptides for mathematical modeling via a linear mixed-effects regression. We then employed forward selection to generate a final model by minimizing Akaike's information criterion (AIC). The final model utilized changes in abundance from 28 peptides as fixed effects to model progression of the disease in these patients. These peptides were from proteins involved in stress response and innate immunity. Graphical abstract.

MSstatsQC: Longitudinal System Suitability Monitoring and Quality Control for Targeted Proteomic Experiments.

Selected Reaction Monitoring (SRM) is a powerful tool for targeted detection and quantification of peptides in complex matrices. An important objective of SRM is to obtain peptide quantifications that are (1) suitable for the investigation, and (2) reproducible across laboratories and runs. The first objective is achieved by system suitability tests (SST), which verify that mass spectrometric instrumentation performs as specified. The second objective is achieved by quality control (QC), which provides in-process quality assurance of the sample profile. A common aspect of SST and QC is the longitudinal nature of the data. Although SST and QC have received a lot of attention in the proteomic community, the currently used statistical methods are limited. This manuscript improves upon the statistical methodology for SST and QC that is currently used in proteomics. It adapts the modern methods of longitudinal statistical process control, such as simultaneous and time weighted control charts and change point analysis, to SST and QC of SRM experiments, discusses their advantages, and provides practical guidelines. Evaluations on simulated data sets, and on data sets from the Clinical Proteomics Technology Assessment for Cancer (CPTAC) consortium, demonstrated that these methods substantially improve our ability of real time monitoring, early detection and prevention of chromatographic and instrumental problems. We implemented the methods in an open-source R-based software package MSstatsQC and its web-based graphical user interface. They are available for use stand-alone, or for integration with automated pipelines. Although the examples focus on targeted proteomics, the statistical methods in this manuscript apply more generally to quantitative proteomics.

A novel integrated strategy for the detection and quantification of the neurotoxin ß-N-methylamino-L-alanine in environmental samples.

We describe a set of new tools for the detection and quantification of ß-N-methylamino-L-alanine (BMAA) which includes a novel stable isotope-labeled BMAA standard (13C3,15N2) and a chip-based capillary electrophoresis mass spectrometry platform for separation and detection. Baseline resolution of BMAA from its potentially confounding structural isomers N-2-aminoethylglycine (AEG) and 2,4-diaminobutyric acid (2,4-DAB) is achieved using the chip-based CE-MS system in less than 1 min. Detection and linearity of response are demonstrated across > 3.5 orders of dynamic range using parallel reaction monitoring (PRM). The lower limit of detection and quantification were calculated for BMAA detection at 40 nM (4.8 ng/mL) and 400 nM (48 ng/mL), respectively. Finally, the strategy was applied to detect BMAA in seafood samples purchased at a local market in Raleigh, NC where their harvest location was known. BMAA was detected in a sea scallop sample. Because the BMAA/stable isotope-labeled 13C3,15N2-BMAA (SIL-BMAA) ratio in the scallop sample was below the limit of quantification, a semiquantitative analysis of BMAA content was carried out, and BMAA content was estimated to be approximately 820 ng BMAA/1 g of wet scallop tissue. Identification was verified by high mass measurement accuracy of precursor (< 5 ppm) and product ions (< 10 ppm), comigration with SIL-BMAA spike-in standard, and conservation of ion abundance ratios for product ions between BMAA and SIL-BMAA. Interestingly, BMAA was not identified in the free protein fraction but only detected after protein hydrolysis which suggests that BMAA is tightly bound by and/or incorporated into proteins. Graphical abstract Utilization of novel 13C3,15N2-BMAA and chip-based CE-MS/MS for detection and quantification of BMAA in environmental samples.

Exposure to BMAA mirrors molecular processes linked to neurodegenerative disease.

The goal of this study is to investigate the molecular pathways perturbed by in vitro exposure of beta-methylamino-L-alanine (BMAA) to NSC-34 cells via contemporary proteomics. Our analysis of differentially regulated proteins reveals significant enrichment (p < 0.01) of pathways related to ER stress, protein ubiquitination, the unfolded protein response, and mitochondrial dysfunction. Upstream regulator analysis indicates that exposure to BMAA induces activation of transcription factors (X-box binding protein 1; nuclear factor 2 erythroid like 2; promyelocytic leukemia) involved in regulation of the UPR, oxidative stress, and cellular senescence. Furthermore, the authors examine the hypothesis that BMAA causes protein damage via misincorporation in place of L-Serine. The authors are unable to detect misincorporation of BMAA into protein via analysis of cellular protein, secreted protein, targeted detection of BMAA after protein hydrolysis, or through the use of in vitro protein translation kits.

An Automated Pipeline to Monitor System Performance in Liquid Chromatography-Tandem Mass Spectrometry Proteomic Experiments.

We report the development of a completely automated pipeline to monitor system suitability in bottom-up proteomic experiments. LC-MS/MS runs are automatically imported into Skyline and multiple identification-free metrics are extracted from targeted peptides. These data are then uploaded to the Panorama Skyline document repository where metrics can be viewed in a web-based interface using powerful process control techniques, including Levey-Jennings and Pareto plots. The interface is versatile and takes user input, which allows the user significant control over the visualization of the data. The pipeline is vendor and instrument-type neutral, supports multiple acquisition techniques (e.g., MS 1 filtering, data-independent acquisition, parallel reaction monitoring, and selected reaction monitoring), can track performance of multiple instruments, and requires no manual intervention aside from initial setup. Data can be viewed from any computer with Internet access and a web browser, facilitating sharing of QC data between researchers. Herein, we describe the use of this pipeline, termed Panorama AutoQC, to evaluate LC-MS/MS performance in a range of scenarios including identification of suboptimal instrument performance, evaluation of ultrahigh pressure chromatography, and identification of the major sources of variation throughout years of peptide data collection.

Improving the analytical performance and versatility of paper spray mass spectrometry via paper microfluidics.

Two paper-based microfluidic techniques, photolithography and wax patterning, were investigated for their potential to improve upon the sensitivity, reproducibility, and versatility of paper spray mass spectrometry. The main limitation of photolithography was the significant signal (approximately three orders of magnitude) above background which was attributed to the chemicals used in the photoresist process. Hydrophobic barriers created via wax patterning were discovered to have approximately 2 orders of magnitude less background signal compared to analogous barriers created using photolithography. A minimum printed wax barrier thickness of approximately 0.3 mm was necessary to consistently retain commonly used paper spray solvents (1 : 1 water : acetonitrile/methanol) and avoid leakage. Constricting capillary flow via wax-printed channels yielded both a significant increase in signal and detection time for detection of model analytes. This signal increase, which was attributed to restricting the radial flow of analyte/solvent on paper (i.e., a concentrating effect), afforded a significant increase in sensitivity (p ≪ 0.05) for the detection of pesticides spiked into residential tap water using a five-point calibration curve. Finally, unique mixing designs using wax patterning can be envisioned to perform on-paper analyte derivatization.

Tools for monitoring system suitability in LC MS/MS centric proteomic experiments.

With advances in liquid chromatography coupled to tandem mass spectrometry technologies combined with the continued goals of biomarker discovery, clinical applications of established biomarkers, and integrating large multiomic datasets (i.e. "big data"), there remains an urgent need for robust tools to assess instrument performance (i.e. system suitability) in proteomic workflows. To this end, several freely available tools have been introduced that monitor a number of peptide identification (ID) and/or peptide ID free metrics. Peptide ID metrics include numbers of proteins, peptides, or peptide spectral matches identified from a complex mixture. Peptide ID free metrics include retention time reproducibility, full width half maximum, ion injection times, and integrated peptide intensities. The main driving force in the development of these tools is to monitor both intra- and interexperiment performance variability and to identify sources of variation. The purpose of this review is to summarize and evaluate these tools based on versatility, automation, vendor neutrality, metrics monitored, and visualization capabilities. In addition, the implementation of a robust system suitability workflow is discussed in terms of metrics, type of standard, and frequency of evaluation along with the obstacles to overcome prior to incorporating a more proactive approach to overall quality control in liquid chromatography coupled to tandem mass spectrometry based proteomic workflows.

Reagent for Evaluating Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) Performance in Bottom-Up Proteomic Experiments.

We present a novel proteomic standard for assessing liquid chromatography-tandem mass spectrometry (LC-MS/MS) instrument performance, in terms of chromatographic reproducibility and dynamic range within a single LC-MS/MS injection. The peptide mixture standard consists of six peptides that were specifically synthesized to cover a wide range of hydrophobicities (grand average hydropathy (GRAVY) scores of -0.6 to 1.9). A combination of stable isotope labeled amino acids ((13)C and (15)N) were inserted to create five isotopologues. By combining these isotopologues at different ratios, they span four orders of magnitude within each distinct peptide sequence. Each peptide, from lightest to heaviest, increases in abundance by a factor of 10. We evaluate several metrics on our quadrupole orbitrap instrument using the 6 × 5 LC-MS/MS reference mixture spiked into a complex lysate background as a function of dynamic range, including mass measurement accuracy (MMA) and the linear range of quantitation of MS1 and parallel reaction monitoring experiments. Detection and linearity of the instrument routinely spanned three orders of magnitude across the gradient (500 fmol to 0.5 fmol on column) and no systematic trend was observed for MMA of targeted peptides as a function of abundance by analysis of variance analysis (p = 0.17). Detection and linearity of the fifth isotopologue (i.e., 0.05 fmol on column) was dependent on the peptide and instrument scan type (MS1 vs PRM). We foresee that this standard will serve as a powerful method to conduct both intra-instrument performance monitoring/evaluation, technology development, and inter-instrument comparisons.

Data-independent-acquisition mass spectrometry for identification of targeted-peptide site-specific modifications.

We present a novel strategy based on data-independent acquisition coupled to targeted data extraction for the detection and identification of site-specific modifications of targeted peptides in a completely unbiased manner. This method requires prior knowledge of the site of the modification along the peptide backbone from the protein of interest, but not the mass of the modification. The procedure, named multiplex adduct peptide profiling (MAPP), consists of three steps: 1) A fragment-ion tag is extracted from the data, consisting of the b-type and y-type ion series from the N and C-terminus, respectively, up to the amino-acid position that is believed to be modified; 2) MS1 features are matched to the fragment-ion tag in retention-time space, using the isolation window as a pre-filter to enable calculation of the mass of the modification; and 3) modified fragment ions are overlaid with the unmodified fragment ions to verify the mass calculated in step 2. We discuss the development, applications, and limitations of this new method for detection of unknown peptide modifications. We present an application of the method in profiling adducted peptides derived from abundant proteins in biological fluids with the ultimate objective of detecting biomarkers of exposure to reactive species.

Development and characterization of a novel plug and play liquid chromatography-mass spectrometry (LC-MS) source that automates connections between the capillary trap, column, and emitter.

We report the development and characterization of a novel, vendor-neutral ultra-high pressure-compatible (~10,000 p.s.i.) LC-MS source. This device is the first to make automated connections with user-packed capillary traps, columns, and capillary emitters. The source uses plastic rectangular inserts (referred to here as cartridges) where individual components (i.e. trap, column, or emitter) can be exchanged independent of one another in a plug and play manner. Automated robotic connections are made between the three cartridges using linear translation powered by stepper motors to axially compress each cartridge by applying a well controlled constant compression force to each commercial LC fitting. The user has the versatility to tailor the separation (e.g. the length of the column, type of stationary phase, and mode of separation) to the experimental design of interest in a cost-effective manner. The source is described in detail, and several experiments are performed to evaluate the robustness of both the system and the exchange of the individual trap and emitter cartridges. The standard deviation in the retention time of four targeted peptides from a standard digest interlaced with a soluble Caenorhabditis elegans lysate ranged between 3.1 and 5.3 s over 3 days of analyses. Exchange of the emitter cartridge was found to have an insignificant effect on the abundance of various peptides. In addition, the trap cartridge can be replaced with minimal effects on retention time (<20 s).

Evidence of the biochemical basis of host virulence in the greenbug aphid, Schizaphis graminum (Homoptera: Aphididae).

Biotypes of aphids and many other insect pests are defined based on the phenotypic response of host plants to the insect pest without considering their intrinsic characteristics and genotypes. Plant breeders have spent considerable effort developing aphid-resistant, small-grain varieties to limit insecticide control of the greenbug, Schizaphis graminum. However, new S. graminum biotypes frequently emerge that break resistance. Mechanisms of virulence on the aphid side of the plant-insect interaction are not well understood. S. graminum biotype H is highly virulent on most small grain varieties. This characteristic makes biotype H ideal for comparative proteomics to investigate the basis of biotype virulence in aphids. In this study, we used comparative proteomics to identify protein expression differences associated with virulence. Aphid proteins involved in the tricarboxylic acid cycle, immune system, cell division, and antiapoptosis pathways were found to be up-regulated in biotype H relative to other biotypes. Proteins from the bacterial endosymbiont of aphids were also differentially expressed in biotype H. Guided by the proteome results, we tested whether biotype H had a fitness advantage compared with other S. graminum biotypes and found that biotype H had a higher reproductive fitness as compared with two other biotypes on a range of different wheat germplasms. Finally, we tested whether aphid genetics can be used to further dissect the genetic mechanisms of biotype virulence in aphids. The genetic data showed that sexual reproduction is a source of biotypic variation observed in S. graminum.

Platform-independent and label-free quantitation of proteomic data using MS1 extracted ion chromatograms in skyline: application to protein acetylation and phosphorylation.

Despite advances in metabolic and postmetabolic labeling methods for quantitative proteomics, there remains a need for improved label-free approaches. This need is particularly pressing for workflows that incorporate affinity enrichment at the peptide level, where isobaric chemical labels such as isobaric tags for relative and absolute quantitation and tandem mass tags may prove problematic or where stable isotope labeling with amino acids in cell culture labeling cannot be readily applied. Skyline is a freely available, open source software tool for quantitative data processing and proteomic analysis. We expanded the capabilities of Skyline to process ion intensity chromatograms of peptide analytes from full scan mass spectral data (MS1) acquired during HPLC MS/MS proteomic experiments. Moreover, unlike existing programs, Skyline MS1 filtering can be used with mass spectrometers from four major vendors, which allows results to be compared directly across laboratories. The new quantitative and graphical tools now available in Skyline specifically support interrogation of multiple acquisitions for MS1 filtering, including visual inspection of peak picking and both automated and manual integration, key features often lacking in existing software. In addition, Skyline MS1 filtering displays retention time indicators from underlying MS/MS data contained within the spectral library to ensure proper peak selection. The modular structure of Skyline also provides well defined, customizable data reports and thus allows users to directly connect to existing statistical programs for post hoc data analysis. To demonstrate the utility of the MS1 filtering approach, we have carried out experiments on several MS platforms and have specifically examined the performance of this method to quantify two important post-translational modifications: acetylation and phosphorylation, in peptide-centric affinity workflows of increasing complexity using mouse and human models.

Garrod's fourth inborn error of metabolism solved by the identification of mutations causing pentosuria.

Pentosuria is one of four conditions hypothesized by Archibald Garrod in 1908 to be inborn errors of metabolism. Mutations responsible for the other three conditions (albinism, alkaptonuria, and cystinuria) have been identified, but the mutations responsible for pentosuria remained unknown. Pentosuria, which affects almost exclusively individuals of Ashkenazi Jewish ancestry, is characterized by high levels of the pentose sugar L-xylulose in blood and urine and deficiency of the enzyme L-xylulose reductase. The condition is autosomal-recessive and completely clinically benign, but in the early and mid-20th century attracted attention because it was often confused with diabetes mellitus and inappropriately treated with insulin. Persons with pentosuria were identified from records of Margaret Lasker, who studied the condition in the 1930s to 1960s. In the DCXR gene encoding L-xylulose reductase, we identified two mutations, DCXR c.583ΔC and DCXR c.52(+1)G > A, each predicted to lead to loss of enzyme activity. Of nine unrelated living pentosuric subjects, six were homozygous for DCXR c.583ΔC, one was homozygous for DCXR c.52(+1)G > A, and two were compound heterozygous for the two mutant alleles. L-xylulose reductase was not detectable in protein lysates from subjects' cells and high levels of xylulose were detected in their sera, confirming the relationship between the DCXR genotypes and the pentosuric phenotype. The combined frequency of the two mutant DCXR alleles in 1,067 Ashkenazi Jewish controls was 0.0173, suggesting a pentosuria frequency of approximately one in 3,300 in this population. Haplotype analysis indicated that the DCXR c.52(+1)G > A mutation arose more recently than the DCXR c.583ΔC mutation.

A framework for quality control in quantitative proteomics.

A thorough evaluation of the quality, reproducibility, and variability of bottom-up proteomics data is necessary at every stage of a workflow from planning to analysis. We share real-world case studies applying adaptable quality control (QC) measures to assess sample preparation, system function, and quantitative analysis. System suitability samples are repeatedly measured longitudinally with targeted methods, and we share examples where they are used on three instrument platforms to identify severe system failures and track function over months to years. Internal QCs incorporated at protein and peptide-level allow our team to assess sample preparation issues and to differentiate system failures from sample-specific issues. External QC samples prepared alongside our experimental samples are used to verify the consistency and quantitative potential of our results during batch correction and normalization before assessing biological phenotypes. We combine these controls with rapid analysis using Skyline, longitudinal QC metrics using AutoQC, and server-based data deposition using PanoramaWeb. We propose that this integrated approach to QC be used as a starting point for groups to facilitate rapid quality control assessment to ensure that valuable instrument time is used to collect the best quality data possible.

The development of selected reaction monitoring methods for targeted proteomics via empirical refinement.

Software advancements in the last several years have had a significant impact on proteomics from method development to data analysis. Herein, we detail a method, which uses our in-house developed software tool termed Skyline, for empirical refinement of candidate peptides from targeted proteins. The method consists of four main steps from generation of a testable hypothesis, method development, peptide refinement, to peptide validation. The ultimate goal is to identify the best performing peptide in terms of ionization efficiency, reproducibility, specificity, and chromatographic characteristics to monitor as a proxy for protein abundance. It is important to emphasize that this method allows the user to perform this refinement procedure in the sample matrix and organism of interest with the instrumentation available. Finally, the method is demonstrated in a case study to determine the best peptide to monitor the abundance of surfactant protein B in lung aspirates.

Estimating relative abundances of proteins from shotgun proteomics data.

BACKGROUND: Spectral counting methods provide an easy means of identifying proteins with differing abundances between complex mixtures using shotgun proteomics data. The crux spectral-counts command, implemented as part of the Crux software toolkit, implements four previously reported spectral counting methods, the spectral index (SI(N)), the exponentially modified protein abundance index (emPAI), the normalized spectral abundance factor (NSAF), and the distributed normalized spectral abundance factor (dNSAF). RESULTS: We compared the reproducibility and the linearity relative to each protein's abundance of the four spectral counting metrics. Our analysis suggests that NSAF yields the most reproducible counts across technical and biological replicates, and both SI(N) and NSAF achieve the best linearity. CONCLUSIONS: With the crux spectral-counts command, Crux provides open-source modular methods to analyze mass spectrometry data for identifying and now quantifying peptides and proteins. The C++ source code, compiled binaries, spectra and sequence databases are available at http://noble.gs.washington.edu/proj/crux-spectral-counts.

Learning score function parameters for improved spectrum identification in tandem mass spectrometry experiments.

The identification of proteins from spectra derived from a tandem mass spectrometry experiment involves several challenges: matching each observed spectrum to a peptide sequence, ranking the resulting collection of peptide-spectrum matches, assigning statistical confidence estimates to the matches, and identifying the proteins. The present work addresses algorithms to rank peptide-spectrum matches. Many of these algorithms, such as PeptideProphet, IDPicker, or Q-ranker, follow a similar methodology that includes representing peptide-spectrum matches as feature vectors and using optimization techniques to rank them. We propose a richer and more flexible feature set representation that is based on the parametrization of the SEQUEST XCorr score and that can be used by all of these algorithms. This extended feature set allows a more effective ranking of the peptide-spectrum matches based on the target-decoy strategy, in comparison to a baseline feature set devoid of these XCorr-based features. Ranking using the extended feature set gives 10-40% improvement in the number of distinct peptide identifications relative to a range of q-value thresholds. While this work is inspired by the model of the theoretical spectrum and the similarity measure between spectra used specifically by SEQUEST, the method itself can be applied to the output of any database search. Further, our approach can be trivially extended beyond XCorr to any linear operator that can serve as similarity score between experimental spectra and peptide sequences.

Evaluation of front-end higher energy collision-induced dissociation on a benchtop dual-pressure linear ion trap mass spectrometer for shotgun proteomics.

We report the implementation of front-end higher energy collision-induced dissociation (fHCD) on a benchtop dual-pressure linear ion trap. Software and hardware modifications were employed, described in detail vide-infra, to allow isolated ions to undergo collisions with ambient gas molecules in an intermediate multipole (q00) of the instrument. Results comparing the performance of fHCD and resonance excitation collision-induced dissociation (RE-CID) in terms of injection time, total number of scans, efficiency, mass measurement accuracy (MMA), unique peptide identifications, and spectral quality of labile modified peptides are presented. fHCD is approximately 23% as efficient as RE-CID, and depending on the search algorithm, it identifies 6.6% more or 15% less peptides (q < 0.01) from a soluble whole-cell lysate ( Caenorhabditis elegans ) than RE-CID using Mascot or Sequest search algorithms, respectively. fHCD offers a clear advantage for the analysis of phosphorylated and glycosylated (O-GlcNAc) peptides as the average cross-correlation score (XCorr) for spectra using fHCD was statistically greater (p < 0.05) than for spectra collected using RE-CID.

A method to determine the kinetics of multiple proteins in human infants with respiratory distress syndrome.

We report a method to measure in vivo turnover of four proteins from sequential tracheal aspirates obtained from human newborn infants with respiratory distress syndrome using targeted proteomics. We detected enrichment for all targeted proteins approximately 3 h from the start of infusion of [5,5,5-(2)H(3)] leucine, secretion times that varied from 1.2 to 2.5 h, and half lives that ranged between 10 and 21 h. Complement factor B, a component of the alternative pathway of complement activation, had an approximately twofold-longer half-life than the other three proteins. In addition, the kinetics of mature and carboxy-terminal tryptic peptides from the same protein (surfactant protein B) were not statistically different (p = 0.49).