Whole cell,label free protein quantitation with data independent acquisition: Quantitation at the MS2 level

Label free quantitation by measurement of peptide fragment signal intensity (MS2 quantitation) is a technique that has seen limited use due to the stochastic nature of data dependent acquisition (DDA). However, data independent acquisition has the potential to make large scale MS2 quantitation a more viable technique. In this study we used an implementation of data independent acquisition—SWATH—to perform label free protein quantitation in a model bacterium Clostridium stercorarium. Four tryptic digests analyzed by SWATH were probed by an ion library containing information on peptide mass and retention time obtained from DDA experiments. Application of this ion library to SWATH data quantified 1030 proteins with at least two peptides quantified (～40% of predicted proteins in the C. stercorarium genome) in each replicate. Quantitative results obtained were very consistent between biological replicates (R² ～ 0.960). Protein quantitation by summation of peptide fragment signal intensities was also highly consistent between biological replicates (R² ～ 0.930), indicating that this approach may have increased viability compared to recent applications in label free protein quantitation. SWATH based quantitation was able to consistently detect differences in relative protein quantity and it provided coverage for a number of proteins that were missed in some samples by DDA analysis.     

A multi-well filtration assay for quantitation of inositol phosphates in biological samples

The quantitation of inositol phosphates (IPs), mediators of certain signal transduction processes, typically involves laborious and time consuming conventional ion-exchange chromatography procedures. We have developed a high throughput microtiter plate-based IP assay that utilizes vacuum rather than gravitational flow and has significant advantages over existing methods. The response of recombinant HEK-293 cells expressing human LHRH receptor cDNA to LHRH agonists was used as a model system to develop the assay conditions. Cell lysates containing labeled IPs were applied in 96-well plates fitted with filtration discs containing regenerated Dowex AGI-X8 resin. Specifically bound inositol phosphates were eluted with 1 M ammonium formate in 0.1 M formic acid directly into a fresh 96-well plate and an aliquot of the eluate from each well is transferred into a 96-well plate and counted. The results were comparable to those obtained with the conventional column method and the variation among replicates was significantly improved. This assay facilitates rapid quantitation of inositol phosphates from a large number of samples with relative ease and reduced generation of radioactive waste.     

Reproducibility and quantitation of amplicon sequencing-based detection

To determine the reproducibility and quantitation of the amplicon sequencing-based detection approach for analyzing microbial community structure, a total of 24 microbial communities from a long-term global change experimental site were examined. Genomic DNA obtained from each community was used to amplify 16S rRNA genes with two or three barcode tags as technical replicates in the presence of a small quantity (0.1% wt/wt) of genomic DNA from Shewanella oneidensis MR-1 as the control. The technical reproducibility of the amplicon sequencing-based detection approach is quite low, with an average operational taxonomic unit (OTU) overlap of 17.2%±2.3% between two technical replicates, and 8.2%±2.3% among three technical replicates, which is most likely due to problems associated with random sampling processes. Such variations in technical replicates could have substantial effects on estimating β-diversity but less on α-diversity. A high variation was also observed in the control across different samples (for example, 66.7-fold for the forward primer), suggesting that the amplicon sequencing-based detection approach could not be quantitative. In addition, various strategies were examined to improve the comparability of amplicon sequencing data, such as increasing biological replicates, and removing singleton sequences and less-representative OTUs across biological replicates. Finally, as expected, various statistical analyses with preprocessed experimental data revealed clear differences in the composition and structure of microbial communities between warming and non-warming, or between clipping and non-clipping. Taken together, these results suggest that amplicon sequencing-based detection is useful in analyzing microbial community structure even though it is not reproducible and quantitative. However, great caution should be taken in experimental design and data interpretation when the amplicon sequencing-based detection approach is used for quantitative analysis of the β-diversity of microbial communities.     

Highly reproducible label free quantitative proteomic analysis of RNA polymerase complexes

The use of quantitative proteomics methods to study protein complexes has the potential to provide in-depth information on the abundance of different protein components as well as their modification state in various cellular conditions. To interrogate protein complex quantitation using shotgun proteomic methods, we have focused on the analysis of protein complexes using label-free multidimensional protein identification technology and studied the reproducibility of biological replicates. For these studies, we focused on three highly related and essential multi-protein enzymes, RNA polymerase I, II, and III from Saccharomyces cerevisiae. We found that label-free quantitation using spectral counting is highly reproducible at the protein and peptide level when analyzing RNA polymerase I, II, and III. In addition, we show that peptide sampling does not follow a random sampling model, and we show the need for advanced computational models to predict peptide detection probabilities. In order to address these issues, we used the APEX protocol to model the expected peptide detectability based on whole cell lysate acquired using the same multidimensional protein identification technology analysis used for the protein complexes. Neither method was able to predict the peptide sampling levels that we observed using replicate multidimensional protein identification technology analyses. In addition to the analysis of the RNA polymerase complexes, our analysis provides quantitative information about several RNAP associated proteins including the RNAPII elongation factor complexes DSIF and TFIIF. Our data shows that DSIF and TFIIF are the most highly enriched RNAP accessory factors in Rpb3-TAP purifications and demonstrate our ability to measure low level associated protein abundance across biological replicates. In addition, our quantitative data supports a model in which DSIF and TFIIF interact with RNAPII in a dynamic fashion in agreement with previously published reports.     

Comparison of endotoxin exposure assessment by bioaerosol impinger and filter-sampling methods

Environmental assessment data collected in two prior occupational hygiene studies of swine barns and sawmills allowed the comparison of concurrent, triplicate, side-by-side endotoxin measurements using air sampling filters and bioaerosol impingers. Endotoxin concentrations in impinger solutions and filter eluates were assayed using the Limulus amebocyte lysate assay. In sawmills, impinger sampling yielded significantly higher endotoxin concentration measurements and lower variances than filter sampling with IOM inhalable dust samplers. Analysis of variance for repeated measures showed that this association remained after controlling for other factors such as replicate, sawmill, sawmill operation, wood type, and interaction terms. Endotoxin concentrations in the swine barns were 10-fold higher on average than in sawmills. These samples demonstrated comparable endotoxin concentration estimates for impinger and filter methods although the variability was lower using the impinger method. In both occupational settings, side-by-side replicates were more uniform for the impinger samples than for the filter samples. This study demonstrates that impinger sampling is an acceptable method for quantitation of area endotoxin concentrations. Further, when sampling is performed with impingers for airborne microorganism quantitation, these same impinger solutions can yield valid endotoxin exposure estimates, negating the need for additional filter sampling.     

Direct determination of anabolic steroid conjugates in human urine by combined high-performance liquid chromatography and tandem mass spectrometry

A novel screening procedure for the sulfate and glucuronide conjugates of testosterone (T) and epitestosterone (E) in human urine was developed based on liquid-solid extraction and microbore high-performance liquid chromatography combined on-line with ion-spray tandem mass spectrometry. Confirmation of the sulfate and glucuronide conjugates of testosterone and epitestosterone isolated frrm normal human urine was acheived by selected reaction monitoring of characteristic product ions of the parent compounds. Endogenous levels of the steroid conjugates are detected in normal male urine and an increase is observed when the sample is fortified with authentic analytical standards of the conjugates. Calibration curves of all steroid conjugates in urine are linear over a range of twenty. Deuterated internal standards of testosterone glucuronide and epitestosterone sulfate were used for quantitation of the endogenous conjugates. T/E ratios were determined based on the glucuronide fractions of six replicates from a normal male and were shown to be statistically reproducible and below the accepted T/E threshold of 6:1. Sulfate conjugates were shown to be present at significantly lower levels in the urine. The method has potential as an alternative for monitoring anabolic steroid conjugates in human urin.     

Comparison of Endotoxin Exposure Assessment by Bioaerosol Impinger and Filter-Sampling Methods

Environmental assessment data collected in two prior occupational hygiene studies of swine barns and sawmills allowed the comparison of concurrent, triplicate, side-by-side endotoxin measurements using air sampling filters and bioaerosol impingers. Endotoxin concentrations in impinger solutions and filter eluates were assayed using the Limulus amebocyte lysate assay. In sawmills, impinger sampling yielded significantly higher endotoxin concentration measurements and lower variances than filter sampling with IOM inhalable dust samplers. Analysis of variance for repeated measures showed that this association remained after controlling for other factors such as replicate, sawmill, sawmill operation, wood type, and interaction terms. Endotoxin concentrations in the swine barns were 10-fold higher on average than in sawmills. These samples demonstrated comparable endotoxin concentration estimates for impinger and filter methods although the variability was lower using the impinger method. In both occupational settings, side-by-side replicates were more uniform for the impinger samples than for the filter samples. This study demonstrates that impinger sampling is an acceptable method for quantitation of area endotoxin concentrations. Further, when sampling is performed with impingers for airborne microorganism quantitation, these same impinger solutions can yield valid endotoxin exposure estimates, negating the need for additional filter sampling.     

Quantitative plant proteomics

Quantitation is an inherent requirement in comparative proteomics and there is no exception to this for plant proteomics. Quantitative proteomics has high demands on the experimental workflow, requiring a thorough design and often a complex multi-step structure. It has to include sufficient numbers of biological and technical replicates and methods that are able to facilitate a quantitative signal read-out. Quantitative plant proteomics in particular poses many additional challenges but because of the nature of plants it also offers some potential advantages. In general, analysis of plants has been less prominent in proteomics. Low protein concentration, difficulties in protein extraction, genome multiploidy, high Rubisco abundance in green tissue, and an absence of well-annotated and completed genome sequences are some of the main challenges in plant proteomics. However, the latter is now changing with several genomes emerging for model plants and crops such as potato, tomato, soybean, rice, maize and barley. This review discusses the current status in quantitative plant proteomics (MS-based and non-MS-based) and its challenges and potentials. Both relative and absolute quantitation methods in plant proteomics from DIGE to MS-based analysis after isotope labeling and label-free quantitation are described and illustrated by published studies. In particular, we describe plant-specific quantitative methods such as metabolic labeling methods that can take full advantage of plant metabolism and culture practices, and discuss other potential advantages and challenges that may arise from the unique properties of plants.     

A Software Toolkit and Interface for Performing Stable Isotope Labeling and Top3 Quantification Using Progenesis LC-MS

Abstract Numerous software packages exist to provide support for quantifying peptides and proteins from mass spectrometry (MS) data. However, many support only a subset of experimental methods or instrument types, meaning that laboratories often have to use multiple software packages. The Progenesis LC-MS software package from Nonlinear Dynamics is a software solution for label-free quantitation. However, many laboratories using Progenesis also wish to employ stable isotope-based methods that are not natively supported in Progenesis. We have developed a Java programming interface that can use the output files produced by Progenesis, allowing the basic MS features quantified across replicates to be used in a range of different experimental methods. We have developed post-processing software (the Progenesis Post-Processor) to embed Progenesis in the analysis of stable isotope labeling data and top3 pseudo-absolute quantitation. We have also created export ability to the new data standard, mzQuantML, produced by the Proteomics Standards Initiative to facilitate the development and standardization process. The software is provided to users with a simple graphical user interface for accessing the different features. The underlying programming interface may also be used by Java developers to develop other routines for analyzing data produced by Progenesis.     

Systematic discovery of ectopic pregnancy serum biomarkers using 3-D protein profiling coupled with label-free quantitation

Ectopic pregnancy (EP) and normal intrauterine pregnancy (IUP) serum proteomes were quantitatively compared to systematically identify candidate biomarkers. A 3-D biomarker discovery strategy consisting of abundant protein immunodepletion, SDS gels, LC-MS/MS, and label-free quantitation of MS signal intensities identified 70 candidate biomarkers with differences between groups greater than 2.5-fold. Further statistical analyses of peptide quantities were used to select the most promising 12 biomarkers for further study, which included known EP biomarkers, novel EP biomarkers (ADAM12 and ISM2), and five specific isoforms of the pregnancy specific beta-1-glycoprotein family. Technical replicates showed good reproducibility and protein intensities from the label-free discovery analysis compared favorably with reported abundance levels of several known reference serum proteins over at least 3 orders of magnitude. Similarly, relative abundances of candidate biomarkers from the label-free discovery analysis were consistent with relative abundances from pilot validation assays performed for five of the 12 most promising biomarkers using label-free multiple reaction monitoring of both the patient serum pools used for discovery and the individual samples that constituted these pools. These results demonstrate robust, reproducible, in-depth 3-D serum proteome discovery, and subsequent pilot-scale validation studies can be achieved readily using label-free quantitation strategies.     

Evaluating Strategies to Normalise Biological Replicates of Western Blot Data

Western blot data are widely used in quantitative applications such as statistical testing and mathematical modelling. To ensure accurate quantitation and comparability between experiments, Western blot replicates must be normalised, but it is unclear how the available methods affect statistical properties of the data. Here we evaluate three commonly used normalisation strategies: (i) by fixed normalisation point or control; (ii) by sum of all data points in a replicate; and (iii) by optimal alignment of the replicates. We consider how these different strategies affect the coefficient of variation (CV) and the results of hypothesis testing with the normalised data. Normalisation by fixed point tends to increase the mean CV of normalised data in a manner that naturally depends on the choice of the normalisation point. Thus, in the context of hypothesis testing, normalisation by fixed point reduces false positives and increases false negatives. Analysis of published experimental data shows that choosing normalisation points with low quantified intensities results in a high normalised data CV and should thus be avoided. Normalisation by sum or by optimal alignment redistributes the raw data uncertainty in a mean-dependent manner, reducing the CV of high intensity points and increasing the CV of low intensity points. This causes the effect of normalisations by sum or optimal alignment on hypothesis testing to depend on the mean of the data tested; for high intensity points, false positives are increased and false negatives are decreased, while for low intensity points, false positives are decreased and false negatives are increased. These results will aid users of Western blotting to choose a suitable normalisation strategy and also understand the implications of this normalisation for subsequent hypothesis testing.     

A simple peak detection and label-free quantitation algorithm for chromatography-mass spectrometry

Background

Label-free quantitation of mass spectrometric data is one of the simplest and least expensive methods for differential expression profiling of proteins and metabolites. The need for high accuracy and performance computational label-free quantitation methods is still high in the biomarker and drug discovery research field. However, recent most advanced types of LC-MS generate huge amounts of analytical data with high scan speed, high accuracy and resolution, which is often impossible to interpret manually. Moreover, there are still issues to be improved for recent label-free methods, such as how to reduce false positive/negatives of the candidate peaks, how to expand scalability and how to enhance and automate data processing. AB3D (A simple label-free quantitation algorithm for Biomarker Discovery in Diagnostics and Drug discovery using LC-MS) has addressed these issues and has the capability to perform label-free quantitation using MS1 for proteomics study.

Results

We developed an algorithm called AB3D, a label free peak detection and quantitative algorithm using MS1 spectral data. To test our algorithm, practical applications of AB3D for LC-MS data sets were evaluated using 3 datasets. Comparisons were then carried out between widely used software tools such as MZmine 2, MSight, SuperHirn, OpenMS and our algorithm AB3D, using the same LC-MS datasets. All quantitative results were confirmed manually, and we found that AB3D could properly identify and quantify known peptides with fewer false positives and false negatives compared to four other existing software tools using either the standard peptide mixture or the real complex biological samples of Bartonella quintana (strain JK31). Moreover, AB3D showed the best reliability by comparing the variability between two technical replicates using a complex peptide mixture of HeLa and BSA samples. For performance, the AB3D algorithm is about 1.2 - 15 times faster than the four other existing software tools.

Conclusions

AB3D is a simple and fast algorithm for label-free quantitation using MS1 mass spectrometry data for large scale LC-MS data analysis with higher true positive and reasonable false positive rates. Furthermore, AB3D demonstrated the best reproducibility and is about 1.2- 15 times faster than those of existing 4 software tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0376-0) contains supplementary material, which is available to authorized users.     

Technical, experimental, and biological variations in isobaric tags for relative and absolute quantitation (iTRAQ)

We assess the reliability of isobaric-tags for relative and absolute quantitation (iTRAQ), based on different types of replicate analyses taking into account technical, experimental, and biological variations. In total, 10 iTRAQ experiments were analyzed across three domains of life involving Saccharomyces cerevisiae KAY446, Sulfolobus solfataricus P2, and Synechocystis sp. PCC 6803. The coverage of protein expression of iTRAQ analysis increases as the variation tolerance increases. In brief, a cutoff point at +/-50% variation (+/-0.50) would yield 88% coverage in quantification based on an analysis of biological replicates. Technical replicate analysis produces a higher coverage level of 95% at a lower cutoff point of +/-30% variation. Experimental or iTRAQ variations exhibit similar behavior as biological variations, which suggest that most of the measurable deviations come from biological variations. These findings underline the importance of replicate analysis as a validation tool and benchmarking technique in protein expression analysis.     

Cat Interferon inhibits Feline Leukaemia Virus Infection in Cell Culture

TRANSMISSION of feline leukaemia can be accomplished with tissue extracts from cases which occur naturally¹. Virus particles which are morphologically indistinguishable from the murine and avian C-type viruses are present in cats with the transmitted disease². Feline leukaemia virus (FeLV) replicates in cat cell cultures³ and infected cells are demonstrable by the indirect immunofiuorescent antibody test which detects FeLV group-specific antigen as granular punctate fluorescence in the cytoplasm of acetone fixed cells⁴; this method allows easy quantitation of the antiviral effect of interferon. We report the production and assay of feline interferon using the fluorescent antibody test with FeLV infected cat cell cultures.     

Use of DNA ladders for reproducible protein fractionation by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) for quantitative proteomics

In proteomics, one-dimensional (1D) sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) is widely used for protein fractionation prior to mass spectrometric analysis to enhance the dynamic range of analysis and to improve the identification of low-abundance proteins. Such protein prefractionation works well for quantitation strategies if the proteins are labeled prior to separation. However, because of the poor reproducibility of cutting gel slices, especially when small amounts of samples are analyzed, its application in label-free and peptide-labeling quantitative proteomics methods has been greatly limited. To overcome this limitation, we developed a new strategy in which a DNA ladder is mixed with the protein sample before PAGE separation. After PAGE separation, the DNA ladder is stained to allow for easy, precise, and reproducible gel cutting. To this end, a novel visible DNA-staining method was developed. This staining method is fast, sensitive, and compatible with mass spectrometry. To evaluate the reproducibility of DNA-ladder-assisted gel cutting for quantitative protein fractionation, we used stable isotope labeling with amino acids in cell culture (SILAC). Our results show that the quantitative error associated with fractionation can be minimized using the DNA-assisted fractionation and multiple replicates of gel cutting. In conclusion, 1D PAGE fractionation in combination with DNA ladders can be used for label-free comparative proteomics without compromising quantitation.     

Evaluation of DNA pooling for the estimation of microsatellite allele frequencies: a case study using striped bass (Morone saxatilis)

Using striped bass (Morone saxatilis) and six multiplexed microsatellite markers, we evaluated procedures for estimating allele frequencies by pooling DNA from multiple individuals, a method suggested as cost-effective relative to individual genotyping. Using moment-based estimators, we estimated allele frequencies in experimental DNA pools and found that the three primary laboratory steps, DNA quantitation and pooling, PCR amplification, and electrophoresis, accounted for 23, 48, and 29%, respectively, of the technical variance of estimates in pools containing DNA from 2-24 individuals. Exact allele-frequency estimates could be made for pools of sizes 2-8, depending on the locus, by using an integer-valued estimator. Larger pools of size 12 and 24 tended to yield biased estimates; however, replicates of these estimates detected allele frequency differences among pools with different allelic compositions. We also derive an unbiased estimator of Hardy-Weinberg disequilibrium coefficients that uses multiple DNA pools and analyze the cost-efficiency of DNA pooling. DNA pooling yields the most potential cost savings when a large number of loci are employed using a large number of individuals, a situation becoming increasingly common as microsatellite loci are developed in increasing numbers of taxa.     

Model-based variance-stabilizing transformation for Illumina microarray data

Variance stabilization is a step in the preprocessing of microarray data that can greatly benefit the performance of subsequent statistical modeling and inference. Due to the often limited number of technical replicates for Affymetrix and cDNA arrays, achieving variance stabilization can be difficult. Although the Illumina microarray platform provides a larger number of technical replicates on each array (usually over 30 randomly distributed beads per probe), these replicates have not been leveraged in the current log2 data transformation process. We devised a variance-stabilizing transformation (VST) method that takes advantage of the technical replicates available on an Illumina microarray. We have compared VST with log2 and Variance-stabilizing normalization (VSN) by using the Kruglyak bead-level data (2006) and Barnes titration data (2005). The results of the Kruglyak data suggest that VST stabilizes variances of bead-replicates within an array. The results of the Barnes data show that VST can improve the detection of differentially expressed genes and reduce false-positive identifications. We conclude that although both VST and VSN are built upon the same model of measurement noise, VST stabilizes the variance better and more efficiently for the Illumina platform by leveraging the availability of a larger number of within-array replicates. The algorithms and Supplementary Data are included in the lumi package of Bioconductor, available at: www.bioconductor.org.     

A mixture model approach to the tests of concordance and discordance between two large-scale experiments with two-sample groups

MOTIVATION: Due to advances in experimental technologies, such as microarray, mass spectrometry and nuclear magnetic resonance, it is feasible to obtain large-scale data sets, in which measurements for a large number of features can be simultaneously collected. However, the sample sizes of these data sets are usually small due to their relatively high costs, which leads to the issue of concordance among different data sets collected for the same study: features should have consistent behavior in different data sets. There is a lack of rigorous statistical methods for evaluating this concordance or discordance. METHODS: Based on a three-component normal-mixture model, we propose two likelihood ratio tests for evaluating the concordance and discordance between two large-scale data sets with two sample groups. The parameter estimation is achieved through the expectation-maximization (E-M) algorithm. A normal-distribution-quantile-based method is used for data transformation. RESULTS: To evaluate the proposed tests, we conducted some simulation studies, which suggested their satisfactory performances. As applications, the proposed tests were applied to three SELDI-MS data sets with replicates. One data set has replicates from different platforms and the other two have replicates from the same platform. We found that data generated by SELDI-MS showed satisfactory concordance between replicates from the same platform but unsatisfactory concordance between replicates from different platforms. AVAILABILITY: The R codes are freely available at http://home.gwu.edu/~ylai/research/Concordance.     

Single quadrupole mass spectrometry for pre-clinical pharmacokinetic analysis: quantitation of carvedilol in dog plasma

The pharmaceutical industry standard for bioanalysis is LC/MS/MS. There are, however, many instances where a single quadrupole detector could successfully be used to provide adequate sensitivity and selectivity for quantitation of drug substances in biological matrices. This paper presents one example of how a single quadrupole detector can be employed in a sensitive and selective analytical method for quantitation of carvedilol. A Synergi Hydro-RP (50 mm x 2 mm i.d.; 4 microm) column was used with acetonitile:water:formic acid mobile phase (32:68:0.01, v/v) at a flow rate of 200 microL/min into a single quadrupole mass spectrometer with an electrospray interface in the positive SIM mode. Using a 300 microL plasma aliquot and a liquid-liquid extraction procedure the limit of quantitation for the assay was 1 ng/mL. The assay utility was demonstrated in the analysis of carvedilol pharmacokinetic profiles in beagle dogs following oral carvedilol administration.     

A cell-cycle-dependent DNA polymerase activity that replicates intact DNA in chromatin

An insoluble DNA polymerase activity that replicates the intact chromatin template at 85% of the rate found in vivo has been partially characterized. HeLa cells, encapsulated in agarose microbeads, are lysed using an isotonic salt concentration: the resulting encapsulated nuclei contain polymerase associated with a nucleoskeleton and the unbroken template. This preparation can be manipulated freely without aggregation or breaking the DNA and yet is accessible to enzymes and other probes. The major activity, which is sensitive to aphidicolin, is found only in S-phase nuclei and replicates DNA semi-conservatively, forming intermediates that are ligated efficiently into larger products.