Peak quantification in surface-enhanced laser desorption/ionization by using mixture models

Enhancing Element Identification by Expectation-Maximization Method in Atom Probe Tomography

This paper describes an alternative way to assign elemental identity to atoms collected by atom probe tomography (APT). This method is based on Bayesian assignation of label through the expectation-maximization method (well known in data analysis). Assuming the correct shape of mass over charge peaks in mass spectra, the probability of each atom to be labeled as a given element is determined, and is used to enhance data visualization and composition mapping in APT analyses. The method is particularly efficient for small count experiments with a low signal to noise ratio, and can be used on small subsets of analyzed volumes, and is complementary to single-ion decomposition methods. Based on the selected model and experimental examples, it is shown that the method enhances our ability to observe and extract information from the raw dataset. The experimental case of the superimposition of the Si peak and N peak in a steel is presented.

A study of the properties of Gaussian mixture model for stable isotope standard quantification in MALDI-TOF MS

The quantification of peptides in Matrix assisted laser desorption/ionization time-of-flight mass spectrum analysis coupled with stable isotope standards has been used to quantify native peptides under many experimental conditions. This approach has difficulties quantifying samples containing peptides with ion currents in overlapping (convolved) spectra. In a previous article we proposed a reparametrized Gaussian mixture model based on the known characteristics of the peptides that could also accommodate overlapping spectra. We demonstrated the application of our model in a series of single and overlapping peptides quantification experiments. Here, we focus solely on studying the properties of our approach and examine the characteristics of the GMM approach in convolved peptides using simulated spectra and provide a method for simulating these spectra.

Signal Partitioning Algorithm for Highly Efficient Gaussian Mixture Modeling in Mass Spectrometry

Mixture - modeling of mass spectra is an approach with many potential applications including peak detection and quantification, smoothing, de-noising, feature extraction and spectral signal compression. However, existing algorithms do not allow for automated analyses of whole spectra. Therefore, despite highlighting potential advantages of mixture modeling of mass spectra of peptide/protein mixtures and some preliminary results presented in several papers, the mixture modeling approach was so far not developed to the stage enabling systematic comparisons with existing software packages for proteomic mass spectra analyses. In this paper we present an efficient algorithm for Gaussian mixture modeling of proteomic mass spectra of different types (e.g., MALDI-ToF profiling, MALDI-IMS). The main idea is automated partitioning of protein mass spectral signal into fragments. The obtained fragments are separately decomposed into Gaussian mixture models. The parameters of the mixture models of fragments are then aggregated to form the mixture model of the whole spectrum. We compare the elaborated algorithm to existing algorithms for peak detection and we demonstrate improvements of peak detection efficiency obtained by using Gaussian mixture modeling. We also show applications of the elaborated algorithm to real proteomic datasets of low and high resolution.

Partitioning the genetic diversity of a virus family: approach and evaluation through a case study of picornaviruses

The recent advent of genome sequences as the only source available to classify many newly discovered viruses challenges the development of virus taxonomy by expert virologists who traditionally rely on extensive virus characterization. In this proof-of-principle study, we address this issue by presenting a computational approach (DEmARC) to classify viruses of a family into groups at hierarchical levels using a sole criterion-intervirus genetic divergence. To quantify genetic divergence, we used pairwise evolutionary distances (PEDs) estimated by maximum likelihood inference on a multiple alignment of family-wide conserved proteins. PEDs were calculated for all virus pairs, and the resulting distribution was modeled via a mixture of probability density functions. The model enables the quantitative inference of regions of distance discontinuity in the family-wide PED distribution, which define the levels of hierarchy. For each level, a limit on genetic divergence, below which two viruses join the same group, was objectively selected among a set of candidates by minimizing violations of intragroup PEDs to the limit. In a case study, we applied the procedure to hundreds of genome sequences of picornaviruses and extensively evaluated it by modulating four key parameters. It was found that the genetics-based classification largely tolerates variations in virus sampling and multiple alignment construction but is affected by the choice of protein and the measure of genetic divergence. In an accompanying paper (C. Lauber and A. E. Gorbalenya, J. Virol. 86:3905-3915, 2012), we analyze the substantial insight gained with the genetics-based classification approach by comparing it with the expert-based picornavirus taxonomy.

Image analysis tools and emerging algorithms for expression proteomics

Since their origins in academic endeavours in the 1970s, computational analysis tools have matured into a number of established commercial packages that underpin research in expression proteomics. In this paper we describe the image analysis pipeline for the established 2-DE technique of protein separation, and by first covering signal analysis for MS, we also explain the current image analysis workflow for the emerging high-throughput 'shotgun' proteomics platform of LC coupled to MS (LC/MS). The bioinformatics challenges for both methods are illustrated and compared, whereas existing commercial and academic packages and their workflows are described from both a user's and a technical perspective. Attention is given to the importance of sound statistical treatment of the resultant quantifications in the search for differential expression. Despite wide availability of proteomics software, a number of challenges have yet to be overcome regarding algorithm accuracy, objectivity and automation, generally due to deterministic spot-centric approaches that discard information early in the pipeline, propagating errors. We review recent advances in signal and image analysis algorithms in 2-DE, MS, LC/MS and Imaging MS. Particular attention is given to wavelet techniques, automated image-based alignment and differential analysis in 2-DE, Bayesian peak mixture models, and functional mixed modelling in MS, and group-wise consensus alignment methods for LC/MS.

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.

An approach to automated frequency-domain feature extraction in nuclear magnetic resonance spectroscopy

For the analysis of metabolite systems, nuclear magnetic resonance (NMR) spectroscopy has become an important quantitative monitoring technology. Automated quantitation methods are highly desired and mainly characterized by the tasks of model selection and parameter approximation. This paper proposes a promising automated two stage approach in the frequency-domain, in which signaling peaks are first identified and filtered from noise based on curvature properties of the spectrum, and then proportionally approximated based on the analytical solution of a Lorentz-function. Remarkably, in opposition to common least-squares approaches, the proposed approximation scheme does not rely on partial derivatives, and furthermore, the runtime is independent to the number of spectral datapoints. Simulations provide promising empirical evidence for successful peak selection and parameter approximation, with the results for the latter highly outperforming the Levenberg-Marquardt algorithm in terms of error minimization and robustness.

On-chip detection of low-molecular-weight recombinant proteins in plant crude extracts by SELDI-TOF MS

This chapter presents a general procedure for the on-chip detection and quantitation of low-molecular-weight recombinant proteins in transgenic plant crude extracts by surface-enhanced laser desorption ionization time-of-flight mass spectrometry (SELDI-TOF MS). A protocol is first described to detect the protein of interest in crude protein extracts of transgenic plant lines, by differential protein mapping against similar extracts from a control, nontransgenic line. A complementary protocol is then presented to generate a standard curve with the SELDI system, allowing the protein to be quantified in different transgenic lines. Overall, this procedure may be carried out within a few hours, without the need for prior purification or enrichment of the recombinant protein.

Reversible jump MCMC approach for peak identification for stroke SELDI mass spectrometry using mixture model

Mass spectrometry (MS) has shown great potential in detecting disease-related biomarkers for early diagnosis of stroke. To discover potential biomarkers from large volume of noisy MS data, peak detection must be performed first. This article proposes a novel automatic peak detection method for the stroke MS data. In this method, a mixture model is proposed to model the spectrum. Bayesian approach is used to estimate parameters of the mixture model, and Markov chain Monte Carlo method is employed to perform Bayesian inference. By introducing a reversible jump method, we can automatically estimate the number of peaks in the model. Instead of separating peak detection into substeps, the proposed peak detection method can do baseline correction, denoising and peak identification simultaneously. Therefore, it minimizes the risk of introducing irrecoverable bias and errors from each substep. In addition, this peak detection method does not require a manually selected denoising threshold. Experimental results on both simulated dataset and stroke MS dataset show that the proposed peak detection method not only has the ability to detect small signal-to-noise ratio peaks, but also greatly reduces false detection rate while maintaining the same sensitivity.

SELDI-TOF mass spectra: A view on sources of variation

Adequate interpretation of mass spectrometry data can yield valuable biomarkers. However, spectrum interpretation is a complicated task. This paper reviews the various factors that determine a sample's spectrum and demonstrates the role of these factors in the interpretation process. We derive a simulation model that adequately predicts the expected spectrum based on known sample content and, in the reverse mode, obtain an analysis model that adequately fits an observed spectrum based on the hypothesized sources of variation.