Reconstruction of Gene Regulatory Networks from RNA-Seq Data

1. Reconstruction of Gene Regulatory Networks from RNA-Seq Data Jianlin Jack Cheng Computer Science Department University of Missouri, Columbia ACM-BCB, 2014

2. Big Data Challenge in Genomic Era Omics Data Biological Experiments DNA/RNA Sequencing Genomics Transcriptomics Proteomics Metabolomics … Biological System Mass Spectrometry Analysis Knowledge

3. Expression Profiles of Genes under Multiple Conditions / Time Points

4. Gene Regulatory Networks (GRN) Transcription factor (TF) regulates a gene GRN of yeast in rich medium Gene regulatory module TF1 TF2 TF3 Bar-Joseph et al., 2003

5. Bayesian Probabilistic Modeling • Assign genes into co-regulated modules • Construct regulatory relations of each module Posterior Likelihood Prior

6. Gene Regulatory Network Modeling Join Zhu et al., 2013

7. Gene Regulatory Logic of a Gene Module as a Decision Tree High Expression Transcription factors and binary regulatory tree gene 1 gene 2 gene 3 …. …. gene n One Gene Module Low Expression Biological Conditions (Treatments) in Columns

8. Regulatory Tree Construction • Pick a TF • Divide conditions into two subsets based expression states • Calculate probability Gaussian Mixture g1 g2 . gi . . gn μ2, σ2 μ1, σ1 Zhu et al., 2013

9. Regulatory Tree Construction • Repeat at next level Gaussian Mixture g1 g2 . gi . . gn Zhu et al., 2013

10. Regulatory Tree Construction Algorithm • Pick a TF • Divide conditions based on TF states • Calculate likelihood • Select TF maximizing likelihood • Repeat g1 g2 . gi . . gn Gaussian Mixture Zhu et al., 2013

11. Gene Re-Assignment Regulatory Tree of a Module gi

12. RNA-Seq Data of Soybean Nodulation • An important source of protein and oil • Nitrogen fixation enabled by soybean-rhizobia symbiotic interactions Nodule

13. Gene Regulatory Modules of Differentially Expressed Genes One out of 10 modules A TF functioning in nodulation according to literature. NSP, whose homologous protein is a nodulation signaling in rice. Zhu et al., 2013

14. Application to Other Species Helix-loop-helix transcription factor 2 • Arabidopsis • Drosophila • Mouse • Human • … Soybean proteins affect TWIST2 – a novel protein related to Kidney disease?

15. Acknowledgements Students • Deb Bhattacharya • Renzhi Cao • Jie Hou • Jilong Li • Matt Spencer • Trieu Tuan • Mingzhu Zhu Collaborators Jim Birchler, Bill Folk,Kevin Fritsche, Michael Greenlief, ZezongGu, Mark Hannink, Trupti Joshi, Dennis Lubahn, ValeriMossine, Alan Parrish, Frank Schmidt, Gary Stacey, Grace Sun, John Walker, Dong Xu

16. Binding Site Analysis • MEME + TomTom to identify two binding sites: BetabetaAlphazinc, finger and Leucine Zipper • TFs in GRAS family contain proteins binding to the motifs.

17. Function Enrichment Validation Function predicted by MULTICOM-PDCN P-value calculated by hypergeometric distribution. Some functions are related to formation of nodule organ. Zhu et al., 2013

18. Protein Interaction and Literature Validation I: TF-TF interactions by STRING, L: Literature Function Support Zhu et al., 2013

19. Computational Model Evaluation

20. GRN of Human Prostate Cancer Under Botanical Treatments Lu et al., submitted

21. Li et al., submitted.