In this proposal, we develop Bayesian methodology for high dimensional genomic data. The overarching theme in this proposal is that we develop several novel statistical methods for motif discovery in genomic sequence data. Chromatin Immunoprecipitation microarray (ChIP-chip) data allows the direct identification of transcription factor binding sites that are active in particular biological states. Jointly modeling array intensities and DNA sequence will lead to more accurate estimation of binding sites. We develop these joint models to account for multiple motifs and varied relationships between binding sites and array intensities. We also propose a novel joint model framework for direct estimation of a motif using gene expression and the DNA sequence that bypasses computationally expensive motif selection procedures. Chromatin structure, in the form of positioning of nucleosomes in DNA, has long been known to play a huge role in protein-DNA binding, however, a quantitative assessment of this role has not been available until very recently. Taking advantage of the increasing availability of accurate experimental data assessing chromatin features, we propose a novel Bayesian statistical model framework for improving motif detection through integration of nucleosome positioning and genomic sequence data. Alternative splicing of mRNA greatly expands the functional repertoire of many genes in the mammalian genome by including or excluding the exons making up the genetic coding sequence. Standard gene expression arrays fail to capture the variability of the exon composition of mRNA species, but rather give a crude measure of overall gene expression. We propose a method that detects over-representation of specific splice junctions in different biological states while adjusting for overall gene expression. The advent of high-throughput genomic technologies has ushered in a new data-driven era, allowing the ability to measure biological activity on a genome-wide scale. Chromatin Immunoprecipitation (ChIP), histone modification, and FAIRE for example are procedures that benefited from this technology, allowing one to determine relative enrichment for their isolated fragments genome wide. The recent development of Next generation sequencing (NGS) platforms offers greater dynamic range, resolution, and genomic coverage in measuring relative enrichment of DNA fragments compared to microarrays. We develop classes of statistical mixture models based on the zero-inflated negative binomial distribution to model such count data and develop an R software package called Zero-Inflated Negative Binomial Algorithm (ZINBA) to carry out the peak calling for a given dataset. 1

Public Health Relevance

We develop Bayesian methodology for high dimensional genomic data. The overarching theme in this proposal is that we develop several novel statistical methods for motif discovery in genomic sequence data. The proposed methodology has major applications in chronic diseases such as cancer, AIDS, cardiovascular disease, and environmental health. We will develop new statistical methods for ChIP-chip data, integrating chormatin structure into motif discovery, joint modeling of gene expression and sequence data, alternative mRNA splicing, and analysis of next generation sequencing (NGS) data. 1

National Institute of Health (NIH)
National Institute of General Medical Sciences (NIGMS)
Research Project (R01)
Project #
Application #
Study Section
Biostatistical Methods and Research Design Study Section (BMRD)
Program Officer
Brazhnik, Paul
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of North Carolina Chapel Hill
Biostatistics & Other Math Sci
Schools of Public Health
Chapel Hill
United States
Zip Code
Gelfond, Jonathan; Goros, Martin; Hernandez, Brian et al. (2018) A System for an Accountable Data Analysis Process in R. R J 10:6-21
Li, Wenqing; Chen, Ming-Hui; Wangy, Xiaojing et al. (2018) Bayesian Design of Non-Inferiority Clinical Trials via the Bayes Factor. Stat Biosci 10:439-459
Wu, Jing; Ibrahim, Joseph G; Chen, Ming-Hui et al. (2018) Bayesian Modeling and Inference for Nonignorably Missing Longitudinal Binary Response Data with Applications to HIV Prevention Trials. Stat Sin 28:1929-1963
Li, Tengfei; Xie, Fengchang; Feng, Xiangnan et al. (2018) Functional Linear Regression Models for Nonignorable Missing Scalar Responses. Stat Sin 28:1867-1886
Psioda, Matthew A; Ibrahim, Joseph G (2018) Bayesian design of a survival trial with a cured fraction using historical data. Stat Med 37:3814-3831
Psioda, Matthew A; Ibrahim, Joseph G (2018) Bayesian clinical trial design using historical data that inform the treatment effect. Biostatistics :
Li, Hao; Chen, Ming-Hui; Ibrahim, Joseph G et al. (2018) Bayesian inference for network meta-regression using multivariate random effects with applications to cholesterol lowering drugs. Biostatistics :
Wang, Yu-Bo; Chen, Ming-Hui; Kuo, Lynn et al. (2018) A New Monte Carlo Method for Estimating Marginal Likelihoods. Bayesian Anal 13:311-333
Lachos, Victor H; A Matos, Larissa; Castro, Luis M et al. (2018) Flexible longitudinal linear mixed models for multiple censored responses data. Stat Med :
Kong, Dehan; Ibrahim, Joseph G; Lee, Eunjee et al. (2018) FLCRM: Functional linear cox regression model. Biometrics 74:109-117

Showing the most recent 10 out of 136 publications