Mass Spectrometry utilizing labeling allows multiple specimens to be subjected to mass spectrometry simultaneously. As a result, between-experiment variability is reduced. Here we describe use of fundamental concepts of statistical experimental design in the labeling framework in order to minimize variability and avoid biases. We demonstrate how to export data in the format that is most efficient for statistical analysis. We demonstrate how to assess the need for normalization, perform normalization, and check whether it worked. We describe how to build a model explaining the observed values and test for differential protein abundance along with descriptive statistics and measures of reliability of the findings. Concepts are illustrated through the use of three case studies utilizing the iTRAQ 4-plex labeling protocol.

MOTIVATION: Using stable isotopes in global proteome scans, labeled molecules from one sample are pooled with unlabeled molecules from another sample and subsequently subjected to mass-spectral analysis. Stable-isotope methodologies make use of the fact that identical molecules of different stable-isotope compositions are differentiated in a mass spectrometer and are represented in a mass spectrum as distinct isotopic clusters with a known mass shift. We describe two multivariable linear regression models for (16)O/(18)O stable-isotope labeled data that jointly model pairs of resolved isotopic clusters from the same peptide and quantify the abundance present in each of the two biological samples while concurrently accounting for peptide-specific incorporation rates of the heavy isotope. The abundance measure for each peptide from the two biological samples is then used in down-stream statistical analyses, e.g. differential expression analysis. Because the multivariable regression models are able to correct for the abundance of the labeled peptide that appear as an unlabeled peptide due to the inability to exchange the natural C-terminal oxygen for the heavy isotope, they are particularly advantageous for a two-step digestion/labeling procedure. We discuss how estimates from the regression model are used to quantify the variability of the estimated abundance measures for the paired samples. Although discussed in the context of (16)O/(18)O stable-isotope labeled data, the multivariable regression models are generalizable to other stable-isotope labeled technologies.

MOTIVATION: The need for normalization in microarray experiments has been well documented in the literature. Currently, many analysis methods treat normalization and analysis as a series of steps, with summarized data carried forward to the next step. RESULTS: We present a unified algorithm which incorporates normalization and class comparison in one analysis using probe level perfect match and mismatch data. The algorithm is based on calibration models common to most biological assays, and the resulting chip-specific parameters have a natural interpretation. We show that the algorithm fits into the statistical generalized linear models framework, describe a practical fitting strategy and present results of the algorithm applied to an example dataset as well as based on metrics used in affycomp. The algorithm ranks amongst the top third of the affycomp competitors, performing best in measures of bias. AVAILABILITY: R functions are available on request from the authors.

MOTIVATION: Our goal was to develop a normalization technique that yields results similar to cyclic loess normalization and with speed comparable to quantile normalization. RESULTS: Fastlo yields normalized values similar to cyclic loess and quantile normalization and is fast; it is at least an order of magnitude faster than cyclic loess and approaches the speed of quantile normalization. Furthermore, fastlo is more versatile than both cyclic loess and quantile normalization because it is model-based. AVAILABILITY: The Splus/R function for fastlo normalization is available from the authors.