Dataset Information

Bayesian hierarchical modeling for subject-level response classification in peptide microarray immunoassays.

ABSTRACT: The peptide microarray immunoassay simultaneously screens sample serum against thousands of peptides, determining the presence of antibodies bound to array probes. Peptide microarrays tiling immunogenic regions of pathogens (e.g., envelope proteins of a virus) are an important high throughput tool for querying and mapping antibody binding. Because of the assay's many steps, from probe synthesis to incubation, peptide microarray data can be noisy with extreme outliers. In addition, subjects may produce different antibody profiles in response to an identical vaccine stimulus or infection, due to variability among subjects' immune systems. We present a robust Bayesian hierarchical model for peptide microarray experiments, pepBayes, to estimate the probability of antibody response for each subject/peptide combination. Heavy-tailed error distributions accommodate outliers and extreme responses, and tailored random effect terms automatically incorporate technical effects prevalent in the assay. We apply our model to two vaccine trial data sets to demonstrate model performance. Our approach enjoys high sensitivity and specificity when detecting vaccine induced antibody responses. A simulation study shows an adaptive thresholding classification method has appropriate false discovery rate control with high sensitivity, and receiver operating characteristics generated on vaccine trial data suggest that pepBayes clearly separates responses from non-responses.

SUBMITTER: Imholte G

PROVIDER: S-EPMC5617895 | biostudies-literature | 2016 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bayesian hierarchical modeling for subject-level response classification in peptide microarray immunoassays.

Imholte Gregory G Gottardo Raphael R

Biometrics 20160408 4

The peptide microarray immunoassay simultaneously screens sample serum against thousands of peptides, determining the presence of antibodies bound to array probes. Peptide microarrays tiling immunogenic regions of pathogens (e.g., envelope proteins of a virus) are an important high throughput tool for querying and mapping antibody binding. Because of the assay's many steps, from probe synthesis to incubation, peptide microarray data can be noisy with extreme outliers. In addition, subjects may p ...[more]

PMID: 27061097

Similar Datasets

Project description:When assessing the human risks due to exposure to environmental chemicals, traditional dose-response analyses are not straightforward when there are numerous high-quality epidemiological studies of priority cancer and non-cancer health outcomes. Given this wealth of information, selecting a single "best" study on which to base dose-response analyses is difficult and would potentially ignore much of the available data. Therefore, systematic approaches are necessary for the analysis of these rich databases. Examples are meta-analysis (and further, meta-regression), which are well established methods that consider and incorporate information from multiple studies into the estimation of risks due to exposure to environmental contaminants. In this paper, we propose a hierarchical, Bayesian meta-analysis approach for the dose-response analysis of multiple epidemiological studies. This paper is the second of two papers detailing this approach; the first covered "pre-analysis" steps necessary to prepare the data for dose-response modeling. This paper focuses on the hierarchical Bayesian approach to dose-response modeling and extrapolation of risk to populations of interest using the association between bladder cancer and oral inorganic arsenic (iAs) exposure as an illustrative case study. In particular, this paper addresses the modeling of both case-control and cohort studies with a flexible, logistic model in a hierarchical Bayesian framework that estimates study-specific slopes, as well as a pooled slope across all studies. This approach is akin to a random effects model in which no assumption is made a priori that there is a single, common slope for all included studies. Further, this paper also details extrapolation of the estimates of logistic slope to extra risk in a target population using a lifetable analysis and basic assumptions about background iAs exposure levels. In this case, the target population was the general United States population and information on all-cause mortality and incidence and mortality from bladder cancer was used to perform the lifetable analysis. The methods herein were developed for general use in investigating the association between any pollutant and observed health-effects in epidemiological studies. In order to demonstrate these methods, inorganic arsenic was chosen as a case study given the large epidemiological database that exists for this contaminant.

Project description:BACKGROUND: Identification of differentially expressed genes is a typical objective when analyzing gene expression data. Recently, Bayesian hierarchical models have become increasingly popular to solve this type of problems. These models show good performance in accommodating noise, variability and low replication of microarray data. However, the correlation between different fluorescent signals measured from a gene spot is ignored, which can diversely affect the data analysis step. In fact, the intensities of the two signals are significantly correlated across samples. The larger the log-transformed intensities are, the smaller the correlation is. RESULTS: Motivated by the complicated error relations in microarray data, we propose a multivariate hierarchical Bayesian framework for data analysis in the replicated microarray experiments. Gene expression data are modelled by a multivariate normal distribution, parameterized by the corresponding mean vectors and covariance matrixes with a conjugate prior distribution. Within the Bayesian framework, a generalized likelihood ratio test (GLRT) is also developed to infer the gene expression patterns. Simulation studies show that the proposed approach presents better operating characteristics and lower false discovery rate (FDR) than existing methods, especially when the correlation coefficient is large. The approach is illustrated with two examples of microarray analysis. The proposed method successfully detects significant genes closely related to the experimental states, which are verified by the biological information. CONCLUSIONS: The multivariate Bayesian model, compatible with the dependence between mean and variance in the univariate Bayesian model, relaxes the constant coefficient of variation assumption between measurements by adding a covariance structure. This model improves the identification of differentially expressed genes significantly since the Bayesian model fit well with the microarray data.

Dataset Information

Bayesian hierarchical modeling for subject-level response classification in peptide microarray immunoassays.

Publications

Bayesian hierarchical modeling for subject-level response classification in peptide microarray immunoassays.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets