Next generation sequencing has revealed the molecular landscape of cells in unprecedented detail. However, for the massively large-scale data produced by assays based on these technologies, informativeness is not only a function of wet-lab technology, but is critically also a function of the analytical pipelines that interpret th data. Our group has developed four statistical tools designed maximize the informativeness of these assays: 1) the Genome Structural Correction (GSC), a nonparametric model of genomic annotations used to assess the significance of relationships between features;2) the Irreproducible Discovery Rate (IDR), an analogue of the FDR that leverages information from biological replicates;3) Statmap, a comprehensive analysis pipeline for ChIP-seq and CAGE data that propagates statistical confidence from base-calling to peak-calling;and 4) Sparse Linear Isoform Discovery and abundance Estimation (SLIDE), an integrative statistical framework for the analysis of RNA-seq, cDNA, and other RNA data aimed at obtaining and quantifying de novo transcript models. These tools are designed to identify and characterize functional elements in genomes;they make minimal assumptions about the data they analyze, and therefore draw reliable conclusions and measures of statistical confidence. During the K99, we will expand and integrate our tools to extend the reach of statistical confidence throughout data interpretation. During the R00, my research will progress toward the inference and assessment of biological networks. Just as ortholog identification has become an essential step in developing animal models of human disease, multi-species network analysis promises to become a key step in interpreting the relationship between genome variation and phenotype. Many mutations, even gene deletions, do not reveal an obvious phenotype. This is due to network robustness, which often differs between closely related species. To understand these phenomena, we aim to: 1) develop standard statistical tools for network inference, and 2) develop """"""""meta models"""""""" of networks that will permit general measures of network orthology. These two aims are tightly linked: we will need critically to characterize the semantics of biological networks to model them. Currently, some models lack consistent definitions of edges and weights, resulting in untestable representations of genomics data. We will develop testable, quantitative models of biological processes, establishing a uniform semantics leveraging the rich theory of complex systems. Each of the tools above will play a key role, especially Statmap and the GSC, which will be needed to propagate statistical confidence into network analysis. Advances will have a transformative effect on our ability to map animal models of disease onto human biology. Nearly nine out of ten new drugs fail in human trials due to issues (e.g. toxicity) not present in animal models. Understanding the orthology not just of individual genes, but of entire biochemical networks will be essential to infer and correct for differences between models of disease and human biology. Solving this problem will be a major step forward in the march from """"""""base-pairs to bedside"""""""".
This proposal outlines training and mentoring plans that emphasize modern nonparametric statistical theory, developmental biology, and hands-on wet-lab techniques. The goal is to produce an independent investigator who functions as a nexus of communication between data producers and data analysts;who is able to recognize and to solve otherwise """"""""orphan"""""""" problems: important biological questions that require advances in statistical theory to be well-answered. The statistical tools that the candidate will generate during the award will lead to testable, quantitative models of biological processes, with the ultimate goal of establishing a uniform semantics for biological network analysis that leverages the increasingly rich theory of complex systems.
|Boley, Nathan; Stoiber, Marcus H; Booth, Benjamin W et al. (2014) Genome-guided transcript assembly by integrative analysis of RNA sequence data. Nat Biotechnol 32:341-6|
|Alam, Tanvir; Medvedeva, Yulia A; Jia, Hui et al. (2014) Promoter analysis reveals globally differential regulation of human long non-coding RNA and protein-coding genes. PLoS One 9:e109443|
|Gerstein, Mark B; Rozowsky, Joel; Yan, Koon-Kiu et al. (2014) Comparative analysis of the transcriptome across distant species. Nature 512:445-8|
|Brown, James B; Boley, Nathan; Eisman, Robert et al. (2014) Diversity and dynamics of the Drosophila transcriptome. Nature 512:393-9|
|Yip, Kevin Y; Cheng, Chao; Bhardwaj, Nitin et al. (2012) Classification of human genomic regions based on experimentally determined binding sites of more than 100 transcription-related factors. Genome Biol 13:R48|
|Derrien, Thomas; Johnson, Rory; Bussotti, Giovanni et al. (2012) The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res 22:1775-89|
|Banfai, Balazs; Jia, Hui; Khatun, Jainab et al. (2012) Long noncoding RNAs are rarely translated in two human cell lines. Genome Res 22:1646-57|
|ENCODE Project Consortium (2012) An integrated encyclopedia of DNA elements in the human genome. Nature 489:57-74|
|Landt, Stephen G; Marinov, Georgi K; Kundaje, Anshul et al. (2012) ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res 22:1813-31|
|Dong, Xianjun; Greven, Melissa C; Kundaje, Anshul et al. (2012) Modeling gene expression using chromatin features in various cellular contexts. Genome Biol 13:R53|