Critical assessment of alignment procedures for LC-MS proteomics and metabolomics measurements

Analytical approaches to metabolomics and applications to systems biology

Phenotypic expression of renal diseases encompasses a complex interaction between genetic, environmental, and local tissue factors. The level of complexity requires integrated understanding of perturbations in the network of genes, proteins, and metabolites. Metabolomics attempts to systematically identify and quantitate metabolites from biological samples. The small molecules represent the end result of complexity of biological processes in a given cell, tissue, or organ, and thus form attractive candidates to understand disease phenotypes. Metabolites represent a diverse group of low-molecular-weight structures including lipids, amino acids, peptides, nucleic acids, and organic acids, which makes comprehensive analysis a difficult analytical challenge. The recent rapid development of a variety of analytical platforms based on mass spectrometry and nuclear magnetic resonance have enabled separation, characterization, detection, and quantification of such chemically diverse structures. Continued development of bioinformatics and analytical strategies will accelerate widespread use and integration of metabolomics into systems biology. Here, we will discuss analytical and bioinformatic techniques and highlight recent studies that use metabolomics in understanding pathophysiology of disease processes.

Generic workflow for quality assessment of quantitative label-free LC-MS analysis

As high-resolution instruments are becoming standard in proteomics laboratories, label-free quantification using precursor measurements is becoming a viable option, and is consequently rapidly gaining popularity. Several software solutions have been presented for label-free analysis, but to our knowledge no conclusive studies regarding the sensitivity and reliability of each step of the analysis procedure has been described. Here, we use real complex samples to assess the reliability of label-free quantification using four different software solutions. A generic approach to quality test quantitative label-free LC-MS is introduced. Measures for evaluation are defined for feature detection, alignment and quantification. All steps of the analysis could be considered adequately performed by the utilized software solutions, although differences and possibilities for improvement could be identified. The described method provides an effective testing procedure, which can help the user to quickly pinpoint where in the workflow changes are needed.

SIMA: simultaneous multiple alignment of LC/MS peak lists

MOTIVATION

Alignment of multiple liquid chromatography/mass spectrometry (LC/MS) experiments is a necessity today, which arises from the need for biological and technical repeats. Due to limits in sampling frequency and poor reproducibility of retention times, current LC systems suffer from missing observations and non-linear distortions of the retention times across runs. Existing approaches for peak correspondence estimation focus almost exclusively on solving the pairwise alignment problem, yielding straightforward but suboptimal results for multiple alignment problems.

RESULTS

We propose SIMA, a novel automated procedure for alignment of peak lists from multiple LC/MS runs. SIMA combines hierarchical pairwise correspondence estimation with simultaneous alignment and global retention time correction. It employs a tailored multidimensional kernel function and a procedure based on maximum likelihood estimation to find the retention time distortion function that best fits the observed data. SIMA does not require a dedicated reference spectrum, is robust with regard to outliers, needs only two intuitive parameters and naturally incorporates incomplete correspondence information. In a comparison with seven alternative methods on four different datasets, we show that SIMA yields competitive and superior performance on real-world data.

AVAILABILITY

A C++ implementation of the SIMA algorithm is available from http://hci.iwr.uni-heidelberg.de/MIP/Software.

Label-free peptide profiling of Orbitrap™ full mass spectra

BACKGROUND

We developed a new version of the open source software package Peptrix that can yet compare large numbers of Orbitrap™ LC-MS data. The peptide profiling results for Peptrix on MS1 spectra were compared with those obtained from a small selection of open source and commercial software packages: msInspect, Sieve™ and Progenesis™. The properties compared in these packages were speed, total number of detected masses, redundancy of masses, reproducibility in numbers and CV of intensity, overlap of masses, and differences in peptide peak intensities. Reproducibility measurements were taken for the different MS1 software applications by measuring in triplicate a complex peptide mixture of immunoglobulin on the Orbitrap™ mass spectrometer. Values of peptide masses detected from the high intensity peaks of the MS1 spectra by peptide profiling were verified with values of the MS2 fragmented and sequenced masses that resulted in protein identifications with a significant score.

FINDINGS

Peptrix finds about the same number of peptide features as the other packages, but peptide masses are in some cases approximately 5 to 10 times less redundant present in the peptide profile matrix. The Peptrix profile matrix displays the largest overlap when comparing the number of masses in a pair between two software applications. The overlap of peptide masses between software packages of low intensity peaks in the spectra is remarkably low with about 50% of the detected masses in the individual packages. Peptrix does not differ from the other packages in detecting 96% of the masses that relate to highly abundant sequenced proteins. MS1 peak intensities vary between the applications in a non linear way as they are not processed using the same method.

CONCLUSIONS

Peptrix is capable of peptide profiling using Orbitrap™ files and finding differential expressed peptides in body fluid and tissue samples. The number of peptide masses detected in Orbitrap™ files can be increased by using more MS1 peptide profiling applications, including Peptrix, since it appears from the comparison of Peptrix with the other applications that all software packages have likely a high false negative rate of low intensity peptide peaks (missing peptides).

Bioinformatics for mass spectrometry-based metabolomics

The broad view of the state of biological systems cannot be complete without the added value of integrating proteomic and genomic data with metabolite measurement. By definition, metabolomics aims at quantifying not less than the totality of small molecules present in a biofluid, tissue, organism, or any material beyond living systems. To cope with the complexity of the task, mass spectrometry (MS) is the most promising analytical environment to fulfill increasing appetite for more accurate and larger view of the metabolome while providing sufficient data generation throughput. Bioinformatics and associated disciplines naturally play a central role in bridging the gap between fast evolving technology and domain experts. Here, we describe the strategies to translate crude MS information into features characteristics of metabolites, and resources available to guide scientists along the metabolomics pipeline. A particular emphasis is put on pragmatic solutions to interpret the outcome of metabolomics experiments at the level of signal processing, statistical treatment, and biochemical understanding.

Bioinformatics for qualitative and quantitative proteomics

Mass spectrometry is today a key analytical technique to elucidate the amount and content of proteins expressed in a certain cellular context. The degree of automation in proteomics has yet to reach that of genomic techniques, but even current technologies make a manual inspection of the data infeasible. This article addresses the key algorithmic problems bioinformaticians face when handling modern proteomic samples and shows common solutions to them. We provide examples on how algorithms can be combined to build relatively complex analysis pipelines, point out certain pitfalls and aspects worth considering and give a list of current state-of-the-art tools.

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.

Comparative evaluation of software for retention time alignment of gas chromatography/time-of-flight mass spectrometry-based metabonomic data

In chromatography-based metabonomic research, retention time (RT) alignment of chromatographic peaks poses a challenge for the accurate profiling of biomarkers. Although a number of RT alignment software has been reported, the performance of these software packages have not been comprehensively evaluated. This study aimed to evaluate the RT alignment accuracy of publicly available and commercial RT alignment software. Two gas chromatography/mass spectrometry (GC/MS) datasets acquired from a mixture of standard metabolites and human bladder cancer urine samples, were used to assess three publicly available software packages, MetAlign, MZmine and TagFinder, and two commercial applications comprising the Calibration feature and Statistical Compare of ChromaTOF software. The overall RT alignment accuracies in aligning standard compounds mixture were 93, 92, 74, 73 and 42% for Calibration feature, MZmine, MetAlign, Statistical Compare and TagFinder, respectively. Additionally, unique trends were observed for the individual software with regards to the different experimental conditions related to extent and direction of RT shifts. Conflicting performance was observed for human urine samples suggesting that RT misalignments still occurred despite the use of RT alignment software. While RT alignment remains an inevitable step in data preprocessing, metabonomic researchers are recommended to perform manual check on the RT alignment of important biomarkers as part of their validation process.

Computational mass spectrometry for metabolomics: identification of metabolites and small molecules

The identification of compounds from mass spectrometry (MS) data is still seen as a major bottleneck in the interpretation of MS data. This is particularly the case for the identification of small compounds such as metabolites, where until recently little progress has been made. Here we review the available approaches to annotation and identification of chemical compounds based on electrospray ionization (ESI-MS) data. The methods are not limited to metabolomics applications, but are applicable to any small compounds amenable to MS analysis. Starting with the definition of identification, we focus on the analysis of tandem mass and MS(n) spectra, which can provide a wealth of structural information. Searching in libraries of reference spectra provides the most reliable source of identification, especially if measured on comparable instruments. We review several choices for the distance functions. The identification without reference spectra is even more challenging, because it requires approaches to interpret tandem mass spectra with regard to the molecular structure. Both commercial and free tools are capable of mining general-purpose compound libraries, and identifying candidate compounds. The holy grail of computational mass spectrometry is the de novo deduction of structure hypotheses for compounds, where method development has only started thus far. In a case study, we apply several of the available methods to the three compounds, kaempferol, reserpine, and verapamil, and investigate whether this results in reliable identifications.

MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data

Background

Mass spectrometry (MS) coupled with online separation methods is commonly applied for differential and quantitative profiling of biological samples in metabolomic as well as proteomic research. Such approaches are used for systems biology, functional genomics, and biomarker discovery, among others. An ongoing challenge of these molecular profiling approaches, however, is the development of better data processing methods. Here we introduce a new generation of a popular open-source data processing toolbox, MZmine 2.

Results

A key concept of the MZmine 2 software design is the strict separation of core functionality and data processing modules, with emphasis on easy usability and support for high-resolution spectra processing. Data processing modules take advantage of embedded visualization tools, allowing for immediate previews of parameter settings. Newly introduced functionality includes the identification of peaks using online databases, MSⁿ data support, improved isotope pattern support, scatter plot visualization, and a new method for peak list alignment based on the random sample consensus (RANSAC) algorithm. The performance of the RANSAC alignment was evaluated using synthetic datasets as well as actual experimental data, and the results were compared to those obtained using other alignment algorithms.

Conclusions

MZmine 2 is freely available under a GNU GPL license and can be obtained from the project website at: http://mzmine.sourceforge.net/. The current version of MZmine 2 is suitable for processing large batches of data and has been applied to both targeted and non-targeted metabolomic analyses.

Mass spectrometry-based metabolomics of yeast

Driven by the advent of metabolomics, recent years have seen renewed interest in the investigation of yeast metabolism. Here we provide a practical guide to metabolomic analysis of yeast using liquid chromatography-mass spectrometry (LC-MS). We begin with background on LC-MS and its utility in studying yeast metabolism. We then describe key issues involved at each step of a typical yeast metabolomics experiment: in experimental design, cell culture, metabolite extraction, LC-MS, and data processing and analysis. Throughout, we highlight interdependencies between the steps that are relevant to developing an integrated workflow which effectively leverages LC-MS to reveal yeast biology.

Characterization of 1H NMR spectroscopic data and the generation of synthetic validation sets

MOTIVATION

Common contemporary practice within the nuclear magnetic resonance (NMR) metabolomics community is to evaluate and validate novel algorithms on empirical data or simplified simulated data. Empirical data captures the complex characteristics of experimental data, but the optimal or most correct analysis is unknown a priori; therefore, researchers are forced to rely on indirect performance metrics, which are of limited value. In order to achieve fair and complete analysis of competing techniques more exacting metrics are required. Thus, metabolomics researchers often evaluate their algorithms on simplified simulated data with a known answer. Unfortunately, the conclusions obtained on simulated data are only of value if the data sets are complex enough for results to generalize to true experimental data. Ideally, synthetic data should be indistinguishable from empirical data, yet retain a known best analysis.

RESULTS

We have developed a technique for creating realistic synthetic metabolomics validation sets based on NMR spectroscopic data. The validation sets are developed by characterizing the salient distributions in sets of empirical spectroscopic data. Using this technique, several validation sets are constructed with a variety of characteristics present in 'real' data. A case study is then presented to compare the relative accuracy of several alignment algorithms using the increased precision afforded by these synthetic data sets.

AVAILABILITY

These data sets are available for download at http://birg.cs.wright.edu/nmr_synthetic_data_sets.

Analytical and statistical approaches to metabolomics research

Metabolomics, the global profiling of metabolites in different living systems, has experienced a rekindling of interest partially due to the improved detection capabilities of the instrumental techniques currently being used in this area of biomedical research. The analytical methods of choice for the analysis of metabolites in search of disease biomarkers in biological specimens, and for the study of various low molecular weight metabolic pathways include NMR spectroscopy, GC/MS, CE/MS, and HPLC/MS. Global metabolite analysis and profiling of two different sets of data results in a plethora of data that is difficult to manage or interpret manually because of their subtle differences. Multivariate statistical methods and pattern-recognition programs were developed to handle the acquired data and to search for the discriminating features between data acquired from two sample sets, healthy and diseased. Metabolomics have been used in toxicology, plant physiology, and biomedical research. In this paper, we discuss various aspects of metabolomic research including sample collection, handling, storage, requirements for sample analysis, peak alignment, data interpretation using statistical approaches, metabolite identification, and finally recommendations for successful analysis.

Metabolite signal identification in accurate mass metabolomics data with MZedDB, an interactive m/z annotation tool utilising predicted ionisation behaviour 'rules'

Background

Metabolomics experiments using Mass Spectrometry (MS) technology measure the mass to charge ratio (m/z) and intensity of ionised molecules in crude extracts of complex biological samples to generate high dimensional metabolite 'fingerprint' or metabolite 'profile' data. High resolution MS instruments perform routinely with a mass accuracy of < 5 ppm (parts per million) thus providing potentially a direct method for signal putative annotation using databases containing metabolite mass information. Most database interfaces support only simple queries with the default assumption that molecules either gain or lose a single proton when ionised. In reality the annotation process is confounded by the fact that many ionisation products will be not only molecular isotopes but also salt/solvent adducts and neutral loss fragments of original metabolites. This report describes an annotation strategy that will allow searching based on all potential ionisation products predicted to form during electrospray ionisation (ESI).

Results

Metabolite 'structures' harvested from publicly accessible databases were converted into a common format to generate a comprehensive archive in MZedDB. 'Rules' were derived from chemical information that allowed MZedDB to generate a list of adducts and neutral loss fragments putatively able to form for each structure and calculate, on the fly, the exact molecular weight of every potential ionisation product to provide targets for annotation searches based on accurate mass. We demonstrate that data matrices representing populations of ionisation products generated from different biological matrices contain a large proportion (sometimes > 50%) of molecular isotopes, salt adducts and neutral loss fragments. Correlation analysis of ESI-MS data features confirmed the predicted relationships of m/z signals. An integrated isotope enumerator in MZedDB allowed verification of exact isotopic pattern distributions to corroborate experimental data.

Conclusion

We conclude that although ultra-high accurate mass instruments provide major insight into the chemical diversity of biological extracts, the facile annotation of a large proportion of signals is not possible by simple, automated query of current databases using computed molecular formulae. Parameterising MZedDB to take into account predicted ionisation behaviour and the biological source of any sample improves greatly both the frequency and accuracy of potential annotation 'hits' in ESI-MS data.

Challenges in applying chemometrics to LC–MS-based global metabolite profile data

Metabolite profiling can provide insights into the metabolic status of complex living systems through the non-targeted analysis of metabolites in any biological sample. Metabolite profiling is complementary to genomics, transcriptomics and proteomics, and its applications span epidemiology, disease diagnosis, nutrition, pharmaceutical research, and toxicology. Metabolic phenotypes are a reflection of an organism’s environment, lifestyle, diet, gut microfloral composition and are also influenced by genetic factors, with important implications in genome-wide-association studies. Specialized analytical platforms, such as NMR spectroscopy and MS, are required to interrogate such metabolic complexity. The increased sophistication of such techniques has lead to a demand for improved data analysis approaches, including preprocessing and advanced chemometric techniques. This article discusses data generation, preprocessing, multivariate analysis and data interpretation for LC-MS-based metabolite profiling, focusing on challenges encountered and potential solutions.

Reference map for liquid chromatography-mass spectrometry-based quantitative proteomics

The accurate mass and time (AMT) tag strategy has been recognized as a powerful tool for high-throughput analysis in liquid chromatography-mass spectrometry (LC-MS)-based proteomics. Due to the complexity of the human proteome, this strategy requires highly accurate mass measurements for confident identifications. We have developed a method of building a reference map that allows relaxed criteria for mass errors yet delivers high confidence for peptide identifications. The samples used for generating the peptide database were produced by collecting cysteine-containing peptides from T47D cells and then fractionating the peptides using strong cationic exchange chromatography (SCX). LC-tandem mass spectrometry (MS/MS) data from the SCX fractions were combined to create a comprehensive reference map. After the reference map was built, it was possible to skip the SCX step in further proteomic analyses. We found that the reference-driven identification increases the overall throughput and proteomic coverage by identifying peptides with low intensity or complex interference. The use of the reference map also facilitates the quantitation process by allowing extraction of peptide intensities of interest and incorporating models of theoretical isotope distribution.

MetAlign: interface-driven, versatile metabolomics tool for hyphenated full-scan mass spectrometry data preprocessing

Hyphenated full-scan MS technology creates large amounts of data. A versatile easy to handle automation tool aiding in the data analysis is very important in handling such a data stream. MetAlign softwareas described in this manuscripthandles a broad range of accurate mass and nominal mass GC/MS and LC/MS data. It is capable of automatic format conversions, accurate mass calculations, baseline corrections, peak-picking, saturation and mass-peak artifact filtering, as well as alignment of up to 1000 data sets. A 100 to 1000-fold data reduction is achieved. MetAlign software output is compatible with most multivariate statistics programs.

Retention time alignment algorithms for LC/MS data must consider non-linear shifts

MOTIVATION

Proteomics has particularly evolved to become of high interest for the field of biomarker discovery and drug development. Especially the combination of liquid chromatography and mass spectrometry (LC/MS) has proven to be a powerful technique for analyzing protein mixtures. Clinically orientated proteomic studies will have to compare hundreds of LC/MS runs at a time. In order to compare different runs, sophisticated preprocessing steps have to be performed. An important step is the retention time (rt) alignment of LC/MS runs. Especially non-linear shifts in the rt between pairs of LC/MS runs make this a crucial and non-trivial problem.

RESULTS

For the purpose of demonstrating the particular importance of correcting non-linear rt shifts, we evaluate and compare different alignment algorithms. We present and analyze two versions of a new algorithm that is based on regression techniques, once assuming and estimating only linear shifts and once also allowing for the estimation of non-linear shifts. As an example for another type of alignment method we use an established alignment algorithm based on shifting vectors that we adapted to allow for correcting non-linear shifts also. In a simulation study, we show that rt alignment procedures that can estimate non-linear shifts yield clearly better alignments. This is even true under mild non-linear deviations.

AVAILABILITY

R code for the regression-based alignment methods and simulated datasets are available at http://www.statistik.tu-dortmund.de/genetik-publikationen-alignment.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.