Statistical Methods for Integrated Gene Regulation Analyses

统计方法的快速发展主要是由于需要描述、建模和分析各种科学和工程学科产生的复杂的大规模数据集。为了充分利用现有的和即将获得的大量基因调控数据，本提案旨在(1)开发结合序列分析、基因表达数据和蛋白质结合数据的预测建模方法；(2)建立一个完整的贝叶斯模型，用于在多个相关物种中重新鉴定顺式调控模块（介导调节蛋白和DNA序列之间相互作用的多个序列基序的组合模式）。在第一个项目中，研究了许多当代统计学习方法的使用，如boost、随机森林、MARS和BART，用于检测有影响的序列信号和预测蛋白质- dna相互作用。提出了将协变量的不确定性纳入统计学习框架的多级模型，并开发了高效的推理计算算法。第二个项目的统计方面涉及通过随机变量的耦合链对多个相互作用的随机过程进行建模。有效的算法利用二维动态规划和先进的蒙特卡罗技术，如回火和等能跳被开发为具有挑战性的贝叶斯推理提出的模型。预计该研究将对分子生物学、遗传学和医学等领域产生直接影响，在这些领域中，基因调控分析发挥着关键作用。除了方法的发展，算法和软件将提供给生物学家使用他们自己的实验数据。这些项目中的许多统计组成部分，如隐马尔可夫模型的耦合和高级蒙特卡罗采样与动态规划的设计，预计将对统计学和其他计算科学做出重大贡献。

The fast development in statistical methodology is mostly driven by the necessity to describe, model and analyze complex large-scale data sets generated from various scientific and engineering disciplines. In order to make full use of available and incoming large amount of data in gene regulation, this proposal aims (1) to develop predictive modeling approaches to combine sequence analyses, gene expression data, and protein binding data; and (2) to develop a full Bayesian model for de novo identification of cis-regulatory modules (combinatorial patterns of multiple sequence motifs that mediate the interactions between regulatory proteins and DNA sequences) in multiple related species. For the first project, the use of many contemporary statistical learning methods is investigated, such as boosting, random forests, MARS and BART, for detecting influential sequence signals and predicting protein-DNA interactions. Multi-level models are proposed to incorporate the uncertainty in covariates into a statistical learning framework and efficient computational algorithms are developed for the inference. The statistical aspects of the second project involve modeling multiple interacting stochastic processes by coupling chains of random variables. Efficient algorithms that utilize two-dimensional dynamic programming and advanced Monte Carlo techniques such as tempering and equi-energy jumps are developed for the challenging Bayesian inference on the proposed model.The proposed research is expected to have direct and immediate impact on various fields in molecular biology, genetics, and medical sciences, in which gene regulation analyses play critical roles. In addition to methodological development, algorithms and software will be delivered for biologists to use on their own experimental data. Many statistical components in these projects, such as the coupling of hidden Markov models and the design of advanced Monte Carlo sampling with dynamic programming, are expected to contribute significantly to statistics and other computational sciences as well.