Olga Vitek

Targeted proteomics is a method of choice for accurate and high-throughput quantification of predefined sets of proteins. Many workflows use isotope-labeled reference peptides for every target protein, which is time consuming and costly. We report a statistical approach for quantifying full protein panels with a reduced set of reference peptides. This label-sparse approach achieves accurate quantification while reducing experimental cost and time. It is implemented in the software tool SparseQuant.

Motivation: RNA-seq experiments produce digital counts of reads that are affected by both biological and technical variation. To distinguish the systematic changes in expression between conditions from noise, the counts are frequently modeled by the Negative Binomial distribution. However, in experiments with small sample size, the per-gene estimates of the dispersion parameter are unreliable. Method: We propose a simple and effective approach for estimating the dispersions. First, we obtain the initial estimates for each gene using the method of moments. Second, the estimates are regularized, i.e. shrunk towards a common value that minimizes the average squared difference between the initial estimates and the shrinkage estimates. The approach does not require extra modeling assumptions, is easy to compute and is compatible with the exact test of differential expression. Results: We evaluated the proposed approach using 10 simulated and experimental datasets and compared its performance with that of currently popular packages edgeR, DESeq, baySeq, BBSeq and SAMseq. For these datasets, sSeq performed favorably for experiments with small sample size in sensitivity, specificity and computational time. Availability:http://www.stat.purdue.edu/∼ovitek/Software.html and Bioconductor. Contact:ovitek@purdue.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Metabolic profiling of urine presents challenges because of the extensive random variation of metabolite concentrations and the dilution resulting from changes in the overall urine volume. Thus statistical analysis methods play a particularly important role; however, appropriate choices of these methods are not straightforward. Here we investigate constant and variance-stabilization normalization of raw and peak picked spectra, for use with exploratory analysis (principal component analysis) and confirmatory analysis (ordinary and Empirical Bayes t-test) in ¹H NMR-based metabolic profiling of urine. We compare the performance of these methods using urine samples spiked with known metabolites according to a Latin square design. We find that analysis of peak picked and logarithm-transformed spectra is preferred, and that signal processing and statistical analysis steps are interdependent. While variance-stabilizing transformation is preferred in conjunction with principal component analysis, constant normalization is more appropriate for use with a t-test. Empirical Bayes t-test provides more reliable conclusions when the number of samples in each group is relatively small. Performance of these methods is illustrated using a clinical metabolomics experiment on patients with type 1 diabetes to evaluate the effect of insulin deprivation.