Proteomic profiling using mass spectrometry – does normalising by total ion current potentially mask some biological differences?

Exploratory Study on Application of MALDI-TOF-MS to Detect SARS-CoV-2 Infection in Human Saliva

SARS-CoV-2 has caused a large outbreak since its emergence in December 2019. COVID-19 diagnosis became a priority so as to isolate and treat infected individuals in order to break the contamination chain. Currently, the reference test for COVID-19 diagnosis is the molecular detection (RT-qPCR) of the virus from nasopharyngeal swab (NPS) samples. Although this sensitive and specific test remains the gold standard, it has several limitations, such as the invasive collection method, the relative high cost and the duration of the test. Moreover, the material shortage to perform tests due to the discrepancy between the high demand for tests and the production capacities puts additional constraints on RT-qPCR. Here, we propose a PCR-free method for diagnosing SARS-CoV-2 based on matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) profiling and machine learning (ML) models from salivary samples. Kinetic saliva samples were collected at enrollment and ten and thirty days later (D0, D10 and D30), to assess the classification performance of the ML models compared to the molecular tests performed on NPS specimens. Spectra were generated using an optimized protocol of saliva collection and successive quality control steps were developed to ensure the reliability of spectra. A total of 360 averaged spectra were included in the study. At D0, the comparison of MS spectra from SARS-CoV-2 positive patients (n = 105) with healthy healthcare controls (n = 51) revealed nine peaks that significantly distinguished the two groups. Among the five ML models tested, support vector machine with linear kernel (SVM-LK) provided the best performance on the training dataset (accuracy = 85.2%, sensitivity = 85.1%, specificity = 85.3%, F1-Score = 85.1%). The application of the SVM-LK model on independent datasets confirmed its performances with 88.9% and 80.8% of correct classification for samples collected at D0 and D30, respectively. Conversely, at D10, the proportion of correct classification had fallen to 64.3%. The analysis of saliva samples by MALDI-TOF MS and ML appears as an interesting supplementary tool for COVID-19 diagnosis, despite the mitigated results obtained for convalescent patients (D10).

Quantification of run order effect on chromatography - mass spectrometry profiling data

Chromatographic systems coupled with mass spectrometry detection are widely used in biological studies investigating how levels of biomolecules respond to different internal and external stimuli. Such changes are normally expected to be of low magnitude and therefore all experimental factors that can influence the analysis need to be understood and minimized. Run order effect is commonly observed and constitutes a major challenge in chromatography-mass spectrometry based profiling studies that needs to be addressed before the biological evaluation of measured data is made. So far there is no established consensus, metric or method that quickly estimates the size of this effect. In this paper we demonstrate how orthogonal projections to latent structures (OPLS®) can be used for objective quantification of the run order effect in profiling studies. The quantification metric is expressed as the amount of variation in the experimental data that is correlated to the run order. One of the primary advantages with this approach is that it provides a fast way of quantifying run-order effect for all detected features, not only internal standards. Results obtained from quantification of run order effect as provided by the OPLS can be used in the evaluation of data normalization, support the optimization of analytical protocols and identification of compounds highly influenced by instrumental drift. The application of OPLS for quantification of run order is demonstrated on experimental data from plasma profiling performed on three analytical platforms: GCMS metabolomics, LCMS metabolomics and LCMS lipidomics.

Multivariate strategy for the sample selection and integration of multi-batch data in metabolomics

INTRODUCTION

Availability of large cohorts of samples with related metadata provides scientists with extensive material for studies. At the same time, recent development of modern high-throughput 'omics' technologies, including metabolomics, has resulted in the potential for analysis of large sample sizes. Representative subset selection becomes critical for selection of samples from bigger cohorts and their division into analytical batches. This especially holds true when relative quantification of compound levels is used.

OBJECTIVES

We present a multivariate strategy for representative sample selection and integration of results from multi-batch experiments in metabolomics.

METHODS

Multivariate characterization was applied for design of experiment based sample selection and subsequent subdivision into four analytical batches which were analyzed on different days by metabolomics profiling using gas-chromatography time-of-flight mass spectrometry (GC-TOF-MS). For each batch OPLS-DA^® was used and its p(corr) vectors were averaged to obtain combined metabolic profile. Jackknifed standard errors were used to calculate confidence intervals for each metabolite in the average p(corr) profile.

RESULTS

A combined, representative metabolic profile describing differences between systemic lupus erythematosus (SLE) patients and controls was obtained and used for elucidation of metabolic pathways that could be disturbed in SLE.

CONCLUSION

Design of experiment based representative sample selection ensured diversity and minimized bias that could be introduced at this step. Combined metabolic profile enabled unified analysis and interpretation.

Statistical methods for handling unwanted variation in metabolomics data

Metabolomics experiments are inevitably subject to a component of unwanted variation, due to factors such as batch effects, long runs of samples, and confounding biological variation. Although the removal of this unwanted variation is a vital step in the analysis of metabolomics data, it is considered a gray area in which there is a recognized need to develop a better understanding of the procedures and statistical methods required to achieve statistically relevant optimal biological outcomes. In this paper, we discuss the causes of unwanted variation in metabolomics experiments, review commonly used metabolomics approaches for handling this unwanted variation, and present a statistical approach for the removal of unwanted variation to obtain normalized metabolomics data. The advantages and performance of the approach relative to several widely used metabolomics normalization approaches are illustrated through two metabolomics studies, and recommendations are provided for choosing and assessing the most suitable normalization method for a given metabolomics experiment. Software for the approach is made freely available.

Current advances in systems and integrative biology

Systems biology has gained a tremendous amount of interest in the last few years. This is partly due to the realization that traditional approaches focusing only on a few molecules at a time cannot describe the impact of aberrant or modulated molecular environments across a whole system. Furthermore, a hypothesis-driven study aims to prove or disprove its postulations, whereas a hypothesis-free systems approach can yield an unbiased and novel testable hypothesis as an end-result. This latter approach foregoes assumptions which predict how a biological system should react to an altered microenvironment within a cellular context, across a tissue or impacting on distant organs. Additionally, re-use of existing data by systematic data mining and re-stratification, one of the cornerstones of integrative systems biology, is also gaining attention. While tremendous efforts using a systems methodology have already yielded excellent results, it is apparent that a lack of suitable analytic tools and purpose-built databases poses a major bottleneck in applying a systematic workflow. This review addresses the current approaches used in systems analysis and obstacles often encountered in large-scale data analysis and integration which tend to go unnoticed, but have a direct impact on the final outcome of a systems approach. Its wide applicability, ranging from basic research, disease descriptors, pharmacological studies, to personalized medicine, makes this emerging approach well suited to address biological and medical questions where conventional methods are not ideal.

Problems inherent to a meta-analysis of proteomics data: a case study on the plants' response to Cd in different cultivation conditions

This meta-analysis focuses on plant-proteome responses to cadmium (Cd) stress. Initially, some general topics related to a proteomics meta-analysis are discussed: (1) obstacles encountered during data analysis, (2) a consensus in proteomic research, (3) validation and good reporting practices for protein identification and (4) guidelines for statistical analysis of differentially abundant proteins. In a second part, the Cd responses in leaves and roots obtained from a proteomics meta-analysis are discussed in (1) a time comparison (short versus long term exposure), and (2) a culture comparison (hydroponics versus soil cultivation). Data of the meta-analysis confirmed the existence of an initial alarm phase upon Cd exposure. Whereas no metabolic equilibrium is established in hydroponically exposed plants, an equilibrium seems to be manifested in roots of plants grown in Cd-contaminated soil after long term exposure. In leaves, the carbohydrate metabolism is primarily affected independent of the exposure time and the cultivation method. In addition, a metabolic shift from CO2-fixation towards respiration is manifested, independent of the cultivation system. Finally, some ideas for the improvement of proteomics setups and for comparisons between studies are discussed.

BIOLOGICAL SIGNIFICANCE

This meta-analysis focuses on the plant responses to Cd stress in leaves and roots at the proteome level. This meta-analysis points out the encountered obstacles when performing a proteomics meta-analysis related to inherent technologies, but also related to experimental setups. Furthermore, the question is addressed whether an extrapolation of results obtained in hydroponic cultivation towards soil-grown plants is possible.

Statistical analysis of metabolomics data

Statistical matters form an integral part of a metabolomics experiment. In this chapter we describe several important aspects in the analysis of metabolomics data such as the removal of unwanted variation and the identification of differentially abundant metabolites, along with a number of other essential statistical considerations.

The age of the "ome": genome, transcriptome and proteome data set collection and analysis

The current state of human genetic studies is both a marvel and a morass. A marvel in that with the completion of the human genome sequence, projects that used to take years now take months or weeks; however, this creates a wealth of data concomitant to a black hole of meaning. In terms of the well used analogy: the human genome sequence is a library in an ancient language with no Rosetta stone. Researchers have readily exploited the human genome map and thousands of candidate gene studies for a multitude of diseases have been performed. However, many of those studies have found that the variants associated with disease risk are not obvious coding changes. The question now becomes: what do these associations mean? One approach to the downstream mapping of associations is to use additional information to map which variant might truly be causative of risk and what that risk variant is doing. This review will summarize the current state of both data set collection and analysis for the understanding of DNA variants and their downstream effects on transcripts and proteins. This article is part of a Special Issue entitled 'Transcriptome'.

Improved signal processing and normalization for biomarker protein detection in broad-mass-range TOF mass spectra from clinical samples

PURPOSE

To demonstrate robust detection of biomarkers in broad-mass-range TOF-MS data.

EXPERIMENTAL DESIGN

Spectra were obtained for two serum protein profiling studies: (i) 2-200 kDa for 132 patients, 67 healthy and 65 diagnosed as having adult T-cell leukemia and (ii) 2-100 kDa for 140 patients, 70 pairs, each with matched prostate-specific antigen (PSA) levels and biopsy-confirmed diagnoses of one benign and one prostate cancer. Signal processing was performed on raw spectra and peak data were normalized using four methods. Feature selection was performed using Bayesian Network Analysis and a classifier was tested on withheld data. Identification of candidate biomarkers was pursued.

RESULTS

Integrated peak intensities were resolved over full spectra. Normalization using local noise values was superior to global methods in reducing peak correlations, reducing replicate variability and improving feature selection stability. For the leukemia data set, potential disease biomarkers were detected and were found to be predictive for withheld data. Preliminary assignments of protein IDs were consistent with published results and LC-MS/MS identification. No prostate-specific-antigen-independent biomarkers were detected in the prostate cancer data set.

CONCLUSIONS AND CLINICAL RELEVANCE

Signal processing, local signal-to-noise (SNR) normalization and Bayesian Network Analysis feature selection facilitate robust detection and identification of biomarker proteins in broad-mass-range clinical TOF-MS data.

Recent advances in computational analysis of mass spectrometry for proteomic profiling

The proteome, defined as an organism's proteins and their actions, is a highly complex end-effector of molecular and cellular events. Differing amounts of proteins in a sample can be indicators of an individual's health status; thus, it is valuable to identify key proteins that serve as 'biomarkers' for diseases. Since the proteome cannot be simply inferred from the genome due to pre- and posttranslational modifications, a direct approach toward mapping the proteome must be taken. The difficulty in evaluating a large number of individual proteins has been eased with the development of high-throughput methods based on mass spectrometry (MS) of peptide or protein mixtures, bypassing the time-consuming, laborious process of protein purification. However, proteomic profiling by MS requires extensive computational analysis. This article describes key issues and recent advances in computational analysis of mass spectra for biomarker identification.

Protein mass spectra data analysis for clinical biomarker discovery: a global review

The identification of new diagnostic or prognostic biomarkers is one of the main aims of clinical cancer research. In recent years there has been a growing interest in using high throughput technologies for the detection of such biomarkers. In particular, mass spectrometry appears as an exciting tool with great potential. However, to extract any benefit from the massive potential of clinical proteomic studies, appropriate methods, improvement and validation are required. To better understand the key statistical points involved with such studies, this review presents the main data analysis steps of protein mass spectra data analysis, from the pre-processing of the data to the identification and validation of biomarkers.

The promise of bone cancer proteomics

Mass spectrometric analysis of the low-molecular-weight (LMW) range of the serum/plasma proteome is revealing the existence of large numbers of previously unknown peptides and protein fragments, predicted to be derived from circulating low-abundance proteins. While genomics and proteomics are the primary discovery research tool, recent innovations in high-throughput proteomics are now standard practice for biomarker and target discovery. Surface-enhanced laser desorption/ionization time-of-flight (SELDI-TOF) mass spectrometry (MS) is the current mainstay for serum or plasma analysis, although other methods are emerging as alternative high-throughput approaches. From a proteomics perspective, the bone cancers, such as myeloma, breast and prostate cancer bony metastases, and osteosarcoma, are likely among the least studied. As recent advances in proteomic technology have thrust the bone cancer field into the era of proteomics, a review of the current status of the proteome as it relates to the skeletal consequences of malignancy seems reasonable.

Arabidopsis and primary photosynthetic metabolism - more than the icing on the cake

Historically speaking, Arabidopsis was not the plant of choice for investigating photosynthesis, with physiologists and biochemists favouring other species such as Chlorella, spinach and pea. However, its inherent advantages for forward genetics rapidly led to its adoption for photosynthesis research. In the last ten years, the availability of the Arabidopsis genome sequence - still the gold-standard for plant genomes - and the rapid expansion of genetic and genomic resources have further increased its importance. Research in Arabidopsis has not only provided comprehensive information about the enzymes and other proteins involved in photosynthesis, but has also allowed transcriptional responses, protein levels and compartmentation to be analysed at a global level for the first time. Emerging technical and theoretical advances offer another leap forward in our understanding of post-translational regulation and the control of metabolism. To illustrate the impact of Arabidopsis, we provide a historical review of research in primary photosynthetic metabolism, highlighting the role of Arabidopsis in elucidation of the pathway of photorespiration and the regulation of RubisCO, as well as elucidation of the pathways of starch turnover and studies of the significance of starch for plant growth.

Challenges for biomarker discovery in body fluids using SELDI-TOF-MS

Protein profiling using SELDI-TOF-MS has gained over the past few years an increasing interest in the field of biomarker discovery. The technology presents great potential if some parameters, such as sample handling, SELDI settings, and data analysis, are strictly controlled. Practical considerations to set up a robust and sensitive strategy for biomarker discovery are presented. This paper also reviews biological fluids generally available including a description of their peculiar properties and the preanalytical challenges inherent to sample collection and storage. Finally, some new insights for biomarker identification and validation challenges are provided.

Optimal preparation methods for automated matrix-assisted laser desorption/ionization time-of-flight mass spectrometry profiling of low molecular weight proteins and peptides

Mass spectrometry (MS) profiling of the proteome and peptidome for disease-associated patterns is a new concept in clinical diagnostics. The technique, however, is highly sensitive to external sources of variation leading to potentially unacceptable numbers of false positive and false negative results. Before MS profiling can be confidently implemented in a medical setting, standard experimental methods must be developed that minimize technical variance. Past studies of variance have focused largely on pre-analytical variation (i.e., sample collection, handling, etc.). Here, we examined how factors at the analytical stage including the matrix and solid-phase extraction influence MS profiling. Firstly, a standard peptide/protein sample was measured automatically by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) MS across five consecutive days using two different preparation methods, dried droplet and sample/matrix, of four types of matrix: alpha-cyano-4-hydroxycinnamic acid (HCCA), sinapinic acid (SA), 2,5-dihydroxybenzoic acid (DHB) and 2,5-dihydroxyacetophenone (DHAP). The results indicated that the matrix preparation greatly influenced a number of key parameters of the spectra including repeatability (within-day variability), reproducibility (inter-day variability), resolution, signal strength, background intensity and detectability. Secondly, an investigation into the variance associated with C8 magnetic bead extraction of the standard sample prior to automated MS profiling demonstrated that the process did not adversely affect these same parameters. In fact, the spectra were generally more robust following extraction. Thirdly, the best performing matrix preparations were evaluated using C8 magnetic bead extracted human plasma. We conclude that the DHAP prepared according to the dried-droplet method is the most appropriate matrix to use when performing automated MS profiling.

Normalization and statistical analysis of quantitative proteomics data generated by metabolic labeling

Comparative proteomics is a powerful analytical method for learning about the responses of biological systems to changes in growth parameters. To make confident inferences about biological responses, proteomics approaches must incorporate appropriate statistical measures of quantitative data. In the present work we applied microarray-based normalization and statistical analysis (significance testing) methods to analyze quantitative proteomics data generated from the metabolic labeling of a marine bacterium (Sphingopyxis alaskensis). Quantitative data were generated for 1,172 proteins, representing 1,736 high confidence protein identifications (54% genome coverage). To test approaches for normalization, cells were grown at a single temperature, metabolically labeled with (14)N or (15)N, and combined in different ratios to give an artificially skewed data set. Inspection of ratio versus average (MA) plots determined that a fixed value median normalization was most suitable for the data. To determine an appropriate statistical method for assessing differential abundance, a -fold change approach, Student's t test, unmoderated t test, and empirical Bayes moderated t test were applied to proteomics data from cells grown at two temperatures. Inverse metabolic labeling was used with multiple technical and biological replicates, and proteomics was performed on cells that were combined based on equal optical density of cultures (providing skewed data) or on cell extracts that were combined to give equal amounts of protein (no skew). To account for arbitrarily complex experiment-specific parameters, a linear modeling approach was used to analyze the data using the limma package in R/Bioconductor. A high quality list of statistically significant differentially abundant proteins was obtained by using lowess normalization (after inspection of MA plots) and applying the empirical Bayes moderated t test. The approach also effectively controlled for the number of false discoveries and corrected for the multiple testing problem using the Storey-Tibshirani false discovery rate (Storey, J. D., and Tibshirani, R. (2003) Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. U.S.A. 100, 9440-9445). The approach we have developed is generally applicable to quantitative proteomics analyses of diverse biological systems.

Sample size determination in clinical proteomic profiling experiments using mass spectrometry for class comparison

Mass spectrometric profiling approaches such as MALDI-TOF and SELDI-TOF are increasingly being used in disease marker discovery, particularly in the lower molecular weight proteome. However, little consideration has been given to the issue of sample size in experimental design. The aim of this study was to develop a protocol for the use of sample size calculations in proteomic profiling studies using MS. These sample size calculations can be based on a simple linear mixed model which allows the inclusion of estimates of biological and technical variation inherent in the experiment. The use of a pilot experiment to estimate these components of variance is investigated and is shown to work well when compared with larger studies. Examination of data from a number of studies using different sample types and different chromatographic surfaces shows the need for sample- and preparation-specific sample size calculations.

Integrated multi-level quality control for proteomic profiling studies using mass spectrometry

Background

Proteomic profiling using mass spectrometry (MS) is one of the most promising methods for the analysis of complex biological samples such as urine, serum and tissue for biomarker discovery. Such experiments are often conducted using MALDI-TOF (matrix-assisted laser desorption/ionisation time-of-flight) and SELDI-TOF (surface-enhanced laser desorption/ionisation time-of-flight) MS. Using such profiling methods it is possible to identify changes in protein expression that differentiate disease states and individual proteins or patterns that may be useful as potential biomarkers. However, the incorporation of quality control (QC) processes that allow the identification of low quality spectra reliably and hence allow the removal of such data before further analysis is often overlooked. In this paper we describe rigorous methods for the assessment of quality of spectral data. These procedures are presented in a user-friendly, web-based program. The data obtained post-QC is then examined using variance components analysis to quantify the amount of variance due to some of the factors in the experimental design.

Results

Using data from a SELDI profiling study of serum from patients with different levels of renal function, we show how the algorithms described in this paper may be used to detect systematic variability within and between sample replicates, pooled samples and SELDI chips and spots. Manual inspection of those spectral data that were identified as being of poor quality confirmed the efficacy of the algorithms. Variance components analysis demonstrated the relatively small amount of technical variance attributable to day of profile generation and experimental array.

Conclusion

Using the techniques described in this paper it is possible to reliably detect poor quality data within proteomic profiling experiments undertaken by MS. The removal of these spectra at the initial stages of the analysis substantially improves the confidence of putative biomarker identification and allows inter-experimental comparisons to be carried out with greater confidence.

The application of SELDI-TOF mass spectrometry to mammalian cell culture

Surface Enhanced Laser Desorption/Ionisation Time-of-Fight Mass Spectrometry (SELDI-TOF MS) is a technique by which protein profiles can be rapidly produced from a wide variety of biological samples. By employing chromatographic surfaces combined with the specificity and reproducibility of mass spectrometry it has allowed for profiles from complex biological samples to be analysed. Profiling and biomarker identification have been employed widely throughout the biological sciences. To date, however, the benefits of SELDI-TOF MS have not been realised in the area of mammalian cell culture. The advantages in identifying markers for cell stresses, apoptosis and other culture parameters mean that these tools could help greatly to enhance monitoring and control of bioreaction process and improve the production of therapeutics. Better characterisation of culture systems through proteome analysis will allow for improved productivity and better yields.