A halotyrosine antibody that detects increased protein modifications in asthma patients.

Airway inflammation has a pathophysiological role in asthma. Eosinophils, which are commonly increased in asthmatic airways, express eosinophil peroxidase and thereby produce hypobromite and bromotyrosine. Bromotyrosine is believed to be a specific marker for eosinophil activity, but developing an antibody against monobromotyrosine, the predominant brominated tyrosine residue found in vivo has proven difficult. We evaluated whether a 3-bromobenozoic acid hapten antigen produced antibodies that recognized halogenated tyrosine residues. Studies with small-molecule inhibitors or brominated or chlorinated protein suggested that a mouse monoclonal antibody (BTK-94C) selectively bound free and protein mono- and dibromotyrosine and, to a lesser degree, chlorotyrosine, and thus was designated a general halotyrosine antibody. We evaluated if this antibody had potential for characterizing human asthma using an enzyme-linked immunosorbent assay (ELISA) microarray platform to examine the halogenation of 23 proteins in three independent sets of sputum samples (52 samples total). In 15 healthy control or asthmatic subjects, ICAM, PDGF and RANTES had greater proportional amounts of halogenation in asthmatic subjects and the halogenation signal was associated with the severity of exercise-induced airway hyperresponsiveness. In 17 severe asthma patients treated with placebo or mepolizumab to suppress eosinophils, drug-related decreases in halogenation were observed with p values ranging from 0.006 to 0.11 for these 3 proteins. Analysis of 20 subjects that either had neutrophilic asthma or were healthy controls demonstrated a broad increase in halotyrosine (possibly chlorotyrosine) in neutrophilic asthmatics. Overall, these results suggest that an ELISA utilizing BTK-94C could prove useful for assessing airway inflammation in asthma patients.

Oxidatively modified proteins as plasma biomarkers in breast cancer.

BACKGROUND: Post-translational protein modifications (PTMs) are increased in breast tumors. OBJECTIVE: We explored whether PTMs on proteins secreted by the breast could be detected in plasma and potentially used for the early detection of breast cancer. METHODS: We used a custom ELISA microarray platform to measure 4-hydroxynonenal (HNE), glutathione (GSH), nitrotyrosine and halotyrosine adducts in 27 secreted proteins, for a total of 108 candidate biomarkers. Two independent sets of human plasma samples were measured, for a total of 160 samples. The results were analyzed for consistent cancer-associated changes across the two sample sets. Plasma samples for both cases and benign controls were collected at the time of tissue diagnosis after referral from a positive screen (such as mammography). The results from both studies were evaluated using ANOVA and t-tests or receiver operator curves (ROC). RESULTS: Levels of GSH-modified ceruloplasmin and HNE-modified PDGF were significantly altered in plasma samples from cancer patients relative to benign controls. Healthy controls, which were only included in the first set of samples, were similar to the benign controls for both of these markers. A combination of three glutathionylated proteins produced the best area under the ROC curve, with a value of 76%. CONCLUSIONS: Specific PTMs in individual proteins may be useful for distinguishing between women with breast cancer and those with benign breast disease. These oxidative changes in plasma proteins may reflect redox changes in breast cancer. Additional studies on oxidative modifications in individual proteins are warranted.

Plasma biomarker profiles differ depending on breast cancer subtype but RANTES is consistently increased.

BACKGROUND: Current biomarkers for breast cancer have little potential for detection. We determined whether breast cancer subtypes influence circulating protein biomarkers. METHODS: A sandwich ELISA microarray platform was used to evaluate 23 candidate biomarkers in plasma samples that were obtained from subjects with either benign breast disease or invasive breast cancer. All plasma samples were collected at the time of biopsy, after a referral due to a suspicious screen (e.g., mammography). Cancer samples were evaluated on the basis of breast cancer subtypes, as defined by the HER2 and estrogen receptor statuses. RESULTS: Ten proteins were statistically altered in at least one breast cancer subtype, including four epidermal growth factor receptor ligands, two matrix metalloproteases, two cytokines, and two angiogenic factors. Only one cytokine, RANTES, was significantly increased (P < 0.01 for each analysis) in all four subtypes, with areas under the curve (AUC) for receiver operating characteristic values that ranged from 0.76 to 0.82, depending on cancer subtype. The best AUC values were observed for analyses that combined data from multiple biomarkers, with values ranging from 0.70 to 0.99, depending on the cancer subtype. Although the results for RANTES are consistent with previous publications, the multi-assay results need to be validated in independent sample sets. CONCLUSIONS: Different breast cancer subtypes produce distinct biomarker profiles, and circulating protein biomarkers have potential to differentiate between true- and false-positive screens for breast cancer. IMPACT: Subtype-specific biomarker panels may be useful for detecting breast cancer or as an adjunct assay to improve the accuracy of current screening methods.

Analysis of high-throughput ELISA microarray data.

Our research group develops analytical methods and software for the high-throughput analysis of quantitative enzyme-linked immunosorbent assay (ELISA) microarrays. ELISA microarrays differ from DNA microarrays in several fundamental aspects and most algorithms for analysis of DNA microarray data are not applicable to ELISA microarrays. In this review, we provide an overview of the steps involved in ELISA microarray data analysis and how the statistically sound algorithms we have developed provide an integrated software suite to address the needs of each data-processing step. The algorithms discussed are available in a set of open-source software tools (http://www.pnl.gov/statistics/ProMAT).

ELISA-BASE: an integrated bioinformatics tool for analyzing and tracking ELISA microarray data.

SUMMARY: ELISA-BASE is an open source database for capturing, organizing and analyzing enzyme-linked immunosorbent assay (ELISA) microarray data. ELISA-BASE is an extension of the BioArray Software Environment (BASE) database system. AVAILABILITY: http://www.pnl.gov/statistics/ProMAT/ELISA-BASE.stm.

An analysis pipeline for the inference of protein-protein interaction networks.

We present a platform for the reconstruction of protein-protein interaction networks inferred from Mass Spectrometry (MS) bait-prey data. The Software Environment for Biological Network Inference (SEBINI), an environment for the deployment of network inference algorithms that use high-throughput data, forms the platform core. Among the many algorithms available in SEBINI is the Bayesian Estimator of Probabilities of Protein-Protein Associations (BEPro3) algorithm, which is used to infer interaction networks from such MS affinity isolation data. Also, the pipeline incorporates the Collective Analysis of Biological Interaction Networks (CABIN) software. We have thus created a structured workflow for protein-protein network inference and supplemental analysis.

A general system for studying protein-protein interactions in Gram-negative bacteria.

One of the most promising methods for large-scale studies of protein interactions is isolation of an affinity-tagged protein with its in vivo interaction partners, followed by mass spectrometric identification of the copurified proteins. Previous studies have generated affinity-tagged proteins using genetic tools or cloning systems that are specific to a particular organism. To enable protein-protein interaction studies across a wider range of Gram-negative bacteria, we have developed a methodology based on expression of affinity-tagged "bait" proteins from a medium copy-number plasmid. This construct is based on a broad-host-range vector backbone (pBBR1MCS5). The vector has been modified to incorporate the Gateway DEST vector recombination region, to facilitate cloning and expression of fusion proteins bearing a variety of affinity, fluorescent, or other tags. We demonstrate this methodology by characterizing interactions among subunits of the DNA-dependent RNA polymerase complex in two metabolically versatile Gram-negative microbial species of environmental interest, Rhodopseudomonas palustris CGA010 and Shewanella oneidensis MR-1. Results compared favorably with those for both plasmid and chromosomally encoded affinity-tagged fusion proteins expressed in a model organism, Escherichia coli.

A Bayesian estimator of protein-protein association probabilities.

UNLABELLED: The Bayesian Estimator of Protein-Protein Association Probabilities (BEPro aff3) is a software tool for estimating probabilities of protein-protein association between bait and prey protein pairs using data from multiple-bait, multiple-replicate, protein liquid chromatography tandem mass spectrometry LC-MS/MS affinity isolation experiments. AVAILABILITY: BEPro (3) is public domain software, has been tested on WIndows XP, Linux and Mac OS, and is freely available from http://www.pnl.gov/statistics/BEPro3. SUPPLEMENTARY INFORMATION: A user guide, example dataset with analysis and additional documentation are included with the BEPro (3) download.

MASIC: a software program for fast quantitation and flexible visualization of chromatographic profiles from detected LC-MS(/MS) features.

Quantitative analysis of liquid chromatography (LC)-mass spectrometry (MS) and tandem mass spectrometry (MS/MS) data is essential to many proteomics studies. We have developed MASIC(2) to accurately measure peptide abundances and LC elution times in LC-MS/MS analyses. This software program uses an efficient processing algorithm to quickly generate mass specific selected ion chromatograms from a dataset and provides an interactive browser that allows users to examine individual chromatograms with a variety of options.

Statistically inferring protein-protein associations with affinity isolation LC-MS/MS assays.

Affinity isolation of protein complexes followed by protein identification by LC-MS/MS is an increasingly popular approach for mapping protein interactions. However, systematic and random assay errors from multiple sources must be considered to confidently infer authentic protein-protein interactions. To address this issue, we developed a general, robust statistical method for inferring authentic interactions from protein prey-by-bait frequency tables using a binomial-based likelihood ratio test (LRT) coupled with Bayes' Odds estimation. We then applied our LRT-Bayes' algorithm experimentally using data from protein complexes isolated from Rhodopseudomonas palustris. Our algorithm, in conjunction with the experimental protocol, inferred with high confidence authentic interacting proteins from abundant, stable complexes, but few or no authentic interactions for lower-abundance complexes. The algorithm can discriminate against a background of prey proteins that are detected in association with a large number of baits as an artifact of the measurement. We conclude that the experimental protocol including the LRT-Bayes' algorithm produces results with high confidence but moderate sensitivity. We also found that Monte Carlo simulation is a feasible tool for checking modeling assumptions, estimating parameters, and evaluating the significance of results in protein association studies.

An estimate of biofilm properties using an acoustic microscope.

Noninvasive measurements over a biofilm, a three-dimensional (3-D) community of microorganisms immobilized at a substratum, were made using an acoustic microscope operating at frequencies up to 70 MHz. The microscope scanned a 2.5-mm by 2.5-mm region of a living biofilm having a nominal thickness of 100 microm. Spatial variation of surface heterogeneity, thickness, interior structure, and biomass were estimated. Thickness was estimated as the product of the speed of sound of the medium and the interim between the highest signal peak and that of the substratum plane without biofilm. The thickest portions of biofilm were 145 microm; however, slender structures attributed as streamers extended above, with one obtaining a 274-microm height above the substratum. Three-dimensional iso-contours of amplitude were used to estimate the internal structure of the biofilm. Backscatter amplitude was examined at five zones of increasing height from the substratum to examine biomass distribution. Ultrasound-based estimates of thickness were corroborated with optical microscopy. The experimental acoustic and optical systems, methods used to estimate biofilm properties, and potential applications for the resulting data are discussed.

ProMAT: protein microarray analysis tool.

SUMMARY: ProMAT is a software tool for statistically analyzing data from enzyme-linked immunosorbent assay microarray experiments. The software estimates standard curves, sample protein concentrations and their uncertainties for multiple assays. ProMAT generates a set of comprehensive figures for assessing results and diagnosing process quality. The tool is available for Windows or Mac, and is distributed as open-source Java and R code. AVAILABILITY: ProMAT is available at http://www.pnl.gov/statistics/ProMAT. ProMAT requires Java version 1.5.0 and R version 1.9.1 (or more recent versions). ProMAT requires either Windows XP or Mac OS 10.4 or newer versions.

Estimating probabilities of peptide database identifications to LC-FTICR-MS observations.

BACKGROUND: The field of proteomics involves the characterization of the peptides and proteins expressed in a cell under specific conditions. Proteomics has made rapid advances in recent years following the sequencing of the genomes of an increasing number of organisms. A prominent technology for high throughput proteomics analysis is the use of liquid chromatography coupled to Fourier transform ion cyclotron resonance mass spectrometry (LC-FTICR-MS). Meaningful biological conclusions can best be made when the peptide identities returned by this technique are accompanied by measures of accuracy and confidence. METHODS: After a tryptically digested protein mixture is analyzed by LC-FTICR-MS, the observed masses and normalized elution times of the detected features are statistically matched to the theoretical masses and elution times of known peptides listed in a large database. The probability of matching is estimated for each peptide in the reference database using statistical classification methods assuming bivariate Gaussian probability distributions on the uncertainties in the masses and the normalized elution times. RESULTS: A database of 69,220 features from 32 LC-FTICR-MS analyses of a tryptically digested bovine serum albumin (BSA) sample was matched to a database populated with 97% false positive peptides. The percentage of high confidence identifications was found to be consistent with other database search procedures. BSA database peptides were identified with high confidence on average in 14.1 of the 32 analyses. False positives were identified on average in just 2.7 analyses. CONCLUSION: Using a priori probabilities that contrast peptides from expected and unexpected proteins was shown to perform better in identifying target peptides than using equally likely a priori probabilities. This is because a large percentage of the target peptides were similar to unexpected peptides which were included to be false positives. The use of triplicate analyses with a "2 out of 3" reporting rule was shown to have excellent rejection of false positives.

ELISA microarray technology as a high-throughput system for cancer biomarker validation.

A large gap currently exists between the ability to discover potential biomarkers and the ability to assess the real value of these proteins for cancer screening. One major challenge in biomarker validation is the inherent variability in biomarker levels. This variability stems from the diversity across the human population and the considerable molecular heterogeneity between individual tumors, even those that originate from a single tissue. An additional challenge with cancer screening is that most cancers are rare in the general population, meaning that assay specificity must be very high. Otherwise, the number of false positives will be much greater than the number of true positives. Due to these challenges associated with biomarker validation, it is necessary to analyze thousands of samples in order to obtain a clear idea of the utility of a screening assay. Enzyme-linked immunosorbent assay (ELISA) microarray technology can simultaneously quantify levels of multiple proteins and, thus, has the potential to accelerate validation of protein biomarkers for clinical use. This review will discuss current ELISA microarray technology and potential advances that could help to achieve the reproducibility and throughput that are required to evaluate cancer biomarkers.

Diagnostic oligonucleotide microarray fingerprinting of Bacillus isolates.

A genome-independent microarray and new statistical techniques were used to genotype Bacillus strains and quantitatively compare DNA fingerprints with the known taxonomy of the genus. A synthetic DNA standard was used to understand process level variability and lead to recommended standard operating procedures for microbial forensics and clinical diagnostics.

Automated Microarray Image Analysis Toolbox for MATLAB.

UNLABELLED: The Automated Microarray Image Analysis (AMIA) Toolbox for MATLAB is a flexible, open-source, microarray image analysis tool that allows the user to customize analyses of microarray image sets. This tool provides several methods to identify and quantify spot statistics, as well as extensive diagnostic statistics and images to evaluate data quality and array processing. The open, modular nature of AMIA provides access to implementation details and encourages modification and extension of AMIA's capabilities. AVAILABILITY: The AMIA Toolbox is freely available at http://www.pnl.gov/statistics/amia. The AMIA Toolbox requires MATLAB 6.5 (R13) (MathWorks, Inc. Natick, MA), as well as the Statistics Toolbox 4.1 and Image Processing Toolbox 4.1 for MATLAB or more recent versions. CONTACT: amanda.white@pnl.gov

The utility of accurate mass and LC elution time information in the analysis of complex proteomes.

The combination of mass and normalized elution time (NET) of a peptide identified by liquid chromatography-mass spectrometry (LC-MS) measurements can serve as a unique signature for that peptide. However, the specificity of an LC-MS measurement depends upon the complexity of the proteome (i.e., the number of possible peptides) and the accuracy of the LC-MS measurements. In this work, theoretical tryptic digests of all predicted proteins from the genomes of three organisms of varying complexity were evaluated for specificity. Accuracy of the LC-MS measurement of mass-NET pairs (on a 0 to 1.0 NET scale) was described by bivariate normal sampling distributions centered on the peptide signatures. Measurement accuracy (i.e., mass and NET standard deviations of +/-0.1, 1, 5, and 10 ppm, and +/-0.01 and 0.05, respectively) was varied to evaluate improvements in process quality. The spatially localized confidence score, a conditional probability of peptide uniqueness, formed the basis for the peptide identification. Application of this approach to organisms with comparatively small proteomes, such as Deinococcus radiodurans, shows that modest mass and elution time accuracies are generally adequate for confidently identifying most peptides. For more complex proteomes, more accurate measurements are required. However, the study suggests that the majority of proteins for even the human proteome should be identifiable with reasonable confidence by using LC-MS measurements with mass accuracies within +/-1 ppm and high efficiency separations having elution time measurements within +/-0.01 NET.

Comparing bacterial DNA microarray fingerprints.

Epidemiologic and forensic investigations often require assays to detect subtle genetic differences between closely related microorganisms. Typically, gel electrophoresis is used to compare randomly amplified DNA fragments between microbial samples, where the patterns of DNA fragment sizes are viewed as genotype 'fingerprints'. The limited genomic sample captured on a gel, however, is not always sufficient to discriminate closely related strains. This paper examines the application of microarray technology to DNA fingerprinting as a high-resolution alternative to gel-based methods. The so-called universal microarray, which uses short oligonucleotide probes that do not target specific genes or species, is intended to be applicable to all microorganisms because it does not require prior knowledge of genomic sequence. In principle, closely related strains can be distinguished if enough independent oligonucleotide probes are used on the microarray, i.e., if the genome is sufficiently sampled. In practice, we confront noisy data, imperfectly matched hybridizations, and a high-dimensional inference problem. We describe the statistical problems of microarray fingerprinting, outline similarities with and differences from more conventional microarray applications, and illustrate a statistical measurement error model to fingerprint 10 closely related strains from three Bacillus species, and 3 strains from non-Bacillus species.

Quantitative oligonucleotide microarray fingerprinting of Salmonella enterica isolates.

We report on a genome-independent microbial fingerprinting method using nucleic acid microarrays for microbial forensics and epidemiology applications and demonstrate that the microarray method provides high resolution differentiation between closely related microorganisms, using Salmonella enterica strains as the test case. In replicate trials we used a simple 192 probe nonamer array to construct a fingerprint library of 25 closely related Salmonella isolates. Controlling false discovery rate for multiple testing at alpha = 0.05, at least 295 of 300 pairs of S.enterica isolate fingerprints were found to be statistically distinct using a modified Hotelling T2 test. Although most pairs of Salmonella fingerprints are found to be distinct, forensic applications will also require a protocol for library construction and reliable microbial classification against a fingerprint library. We outline additional steps required to produce such a protocol.

Sequence versus structure for the direct detection of 16S rRNA on planar oligonucleotide microarrays.

A two-probe proximal chaperone detection system consisting of a species-specific capture probe for the microarray and a labeled, proximal chaperone probe for detection was recently described for direct detection of intact rRNAs from environmental samples on oligonucleotide arrays. In this study, we investigated the physical spacing and nucleotide mismatch tolerance between capture and proximal chaperone detector probes that are required to achieve species-specific 16S rRNA detection for the dissimilatory metal and sulfate reducer 16S rRNAs. Microarray specificity was deduced by analyzing signal intensities across replicate microarrays with a statistical analysis-of-variance model that accommodates well-to-well and slide-to-slide variations in microarray signal intensity. Chaperone detector probes located in immediate proximity to the capture probe resulted in detectable, nonspecific binding of nontarget rRNA, presumably due to base-stacking effects. Species-specific rRNA detection was achieved by using a 22-nt capture probe and a 15-nt detector probe separated by 10 to 14 nt along the primary sequence. Chaperone detector probes with up to three mismatched nucleotides still resulted in species-specific capture of 16S rRNAs. There was no obvious relationship between position or number of mismatches and within- or between-genus hybridization specificity. From these results, we conclude that relieving secondary structure is of principal concern for the successful capture and detection of 16S rRNAs on planar surfaces but that the sequence of the capture probe is more important than relieving secondary structure for achieving specific hybridization.