1. Bayesian negative binomial mixture regression models for the analysis of sequence count and methylation data
- Author
-
Marina Vannucci, Alberto Cassese, Michele Guindani, Qiwei Li, FPN Methodologie & Statistiek, and RS: FPN M&S I
- Subjects
Statistics and Probability ,Multivariate statistics ,Computer science ,Bayesian probability ,Breast Neoplasms ,Feature selection ,Computational biology ,01 natural sciences ,General Biochemistry, Genetics and Molecular Biology ,010104 statistics & probability ,03 medical and health sciences ,Prior probability ,Biomarkers, Tumor ,Humans ,Gene Regulatory Networks ,RNA-SEQ ,0101 mathematics ,030304 developmental biology ,Count data ,GENE-EXPRESSION ,0303 health sciences ,Models, Statistical ,Markov random field ,Base Sequence ,General Immunology and Microbiology ,Applied Mathematics ,Integrative analysis ,Linear model ,Bayes Theorem ,General Medicine ,DNA Methylation ,Mixture model ,Binomial Distribution ,VARIABLE SELECTION ,Negative binomial ,Data Interpretation, Statistical ,DIFFERENTIAL EXPRESSION ANALYSIS ,DNA methylation ,Regression Analysis ,Female ,General Agricultural and Biological Sciences ,Mixture models - Abstract
In this article, we develop a Bayesian hierarchical mixture regression model for studying the association between a multivariate response, measured as counts on a set of features, and a set of covariates. We have available RNA-Seq and DNA methylation data measured on breast cancer patients at different stages of the disease. We account for the heterogeneity and over-dispersion of count data (here, RNA-Seq data) by considering a mixture of negative binomial distributions and incorporate the covariates (here, methylation data) into the model via a linear modeling construction on the mean components. Our modeling construction includes several innovative characteristics. First, it employs selection techniques that allow the identification of a small subset of features that best discriminate the samples while simultaneously selecting a set of covariates associated to each feature. Second, it incorporates known dependencies into the feature selection process via the use of Markov random field (MRF) priors. On simulated data, we show how incorporating existing information via the prior model can improve the accuracy of feature selection. In the analysis of RNA-Seq and DNA methylation data on breast cancer, we incorporate knowledge on relationships among genes via a gene-gene network, which we extract from the KEGG database. Our data analysis identifies genes which are discriminatory of cancer stages and simultaneously selects significant associations between those genes and DNA methylation sites. A biological interpretation of our findings reveals several biomarkers that can help understanding the effect of DNA methylation on gene expression transcription across cancer stages.
- Published
- 2019