Detecting robust timedelayed regulation in Mycobacterium tuberculosis Iti
Detecting robust time-delayed regulation in Mycobacterium tuberculosis Iti Chaturvedi and Jagath C Rajapakse INCOB 2009
Gene Regulatory Networks (GRN) • GRN represents the regulatory effects (causal effects) among genes involved in a particular pathway. • Signal transduction may is transient – dynamic delayed regulations seem to exist. • Distributed nature causes intense cross-talk • M. tuberculosis is a bacteria causing TB in man and has a very slow growth rate in vitro. • The DNA repair pathway is activated when a damage to DNA occurs. • System consists of Lex. A and Rec. A and upto 40 other genes that are regulated by these two proteins. agenda Rv 2719 c lin. B dna. B lex. A ruv. C rec. A fadd 21 dna. E 2 fadd 23
Methods of Building GRN • Bayesian networks (BN): graphical models with regulations denoted by conditional probabilities (Heckerman et al 1995) • Dynamic BN (DBN): Transition network over time can model cyclic events (Friedman, N. et al 1998) • Higher-order (HDBN): Extended transition network for longer delays (Zheng et al 2006) • Skip-chain BN : Can model very long delays of arbitary length using feature functions. (Galley, M. 2006) ai are the parents of gene i
Optimization Using GA Each individual is an o-nary interaction matrix where for no interaction and for o-order interaction Crossover involves swapping several rows between two parents.
Initialize Selection Mating Initialize N individuals with 0. 7 similarity using Mutual Information GA Algorithm Rank all individuals using fitness function. Elite individual E is sent to next generation. Select two parents using roulette strategy for mating. E is optimal network Crossover Swap last few rows of the selected parents to generate two new children. Try another crossover point if population similarity is < 0. 7 Yes Terminate Invert a random cell to cause mutation Mutate No If dga < 1 for 20 generations or number of generations greater than Q
Skip-chain Model • The likelihood of a gene expression xi is given by a weighted sum of linear and skip edge scores : • Linear-chain feature functions represent local dependencies of o-order. • The skip-chain features represent long range dependencies in a GRN using a HMM agenda
Viterbi Forward Path • The skip-edge score is given as a the normalized MAP interaction: • We can use maximum likelihood to estimate state transition and emission probabilities : where denotes number of occurrences for • The most probable path is given by MAP estimate using dynamic programming :
Priors for Networks • Most higher-order Markov models are sensitive to change in pathways and associated data. – Gibbs prior is used to model target network prior. – Interaction potentials, denotes an interaction in target network and no interaction where – Here a small and large will reflect prior more and vice versa. • We use adaption to reduce over-fitting due to sparse feature specific data. • Adaption model can combine the reliable DBN with a volatile feature specific HMM for long delays.
Dirichlet Prior over Parameters • We extend the MLE to a Bayesian learning. Dirichlet is a conjugate prior for multinomial distribution. • We can maximize probability as (MAP) : • Using the linear feature as a Dirichlet conjugate prior for the skip feature of a gene we get : where • Lastly, the interpolated probability of gene based on linear and skip-edges is : where
Experiments : M. tuberculosis • Here we looked at the response of bacteria to drug-induced stress. Treatment with Mitocyin C caused DNA damage and hence led to the upregulation of associated repair genes. • Eight time points are available at NCBI Gene Expression Omnibus (GSE 1642 -GPL 1396 series) 0. 33 hr, 0. 75 hr, 1. 5 hr, 2 hr, 4 hr, 6 hr, 8 hr and 12 hr after DNA damage. • Data was discretized into 0 for down and 1 for up regulation. • The corresponding skip probabilities were calculated as described in methods. Upto seven time points of delays were allowed. • Firstly, we used 9 genes previously specified. In order to get an expanded dataset, the original dataset was subjected to ICA and the components closest to 9 genes were identified. This gave us a second dataset of 32 genes. agenda
Time delays Table 1. Predicted by DBN, HDBN, and skip-chain without priors Higher-order edges(hrs) # Genes 9 32 Model: o DBN: 1 HDBN: 3 SKIP-CHAIN: 1/5 DBN: 1 HDBN: 4 SKIP-CHAIN: 2/5 ML 1(0. 75) 2(1. 5) -14. 7 -8. 69 -6. 05 -48. 9 -39. 4 -37. 2 9 8 13 36 20 54 3(2) 2 7 (3) 6 18 14 (41) 4(4) 5(6) 20 (4) It can be seen that the ML of the underlying skip-chain prediction is much higher than the DBN or HDBN, confirming that the network fits data well.
Time delays : Priors Table 2. Predicted by skip-chain models with priors: Gibbs prior for structures and Dirichlet prior for parameters # Genes 9 32 Model: o ML Higher-order edges(hrs) 1(0. 75) 2(1. 5) 3(2) 4(4) 5(6) SKIP-CHAIN: 1 -6. 05 13 (3) SKIP-CHAIN(Gibbs): 1 -5. 8 11 (2) SKIP-CHAIN(Dirichlet): 2 -5. 2 7 13 (11) (1) SKIP-CHAIN(Gibbs and Dirichlet): 3 -3. 27 2 7 (4) (5) SKIP-CHAIN: 2 -37. 2 54 18 (41) (4) SKIP-CHAIN(Gibbs): 3 -35. 2 37 16 24 (40) (3) SKIP-CHAIN(Dirichlet): 2 -35. 05 54 16 (37) (4) SKIP-CHAIN(Gibbs and Dirichlet): 2 -34. 54 50 15 (41) (4) Using priors further increased likelihood and gave many new timedelayed interactions. The combined use of both Dirichlet and Gibbs priors is optimal.
Networks 9 genes Fig 3 : Time-delayed interactions in predicted network of 9 genes (a) DBN network, (b) HDBN network, (c) Skipchain network, (d) Skip-chain network with Gibbs prior, (e) Skip-chain network with Dirichlet prior, (f) Skip-chain network with Gibbs and Dirichlet prior A small number of transcription factors (TF) regulate the rest of the repair system. At the same time the in-degree is low, as each gene is regulated by just one TF.
Networks 32 genes Fig 4 : Time-delayed interactions in predicted network of 32 genes (a) DBN network, (b) HDBN network, (c) Skipchain network, (d) Skip-chain network with Gibbs prior, (e) Skip-chain network with Dirichlet prior, (f) Skip-chain network with Gibbs and Dirichlet prior The second dataset of 32 genes indicated that our method is good for identifying core genes. Use of priors gave better networks with fewer hubs.
Conclusion • An organism responds to changes in its environment by altering the level of expression of critical genes. • Skip-chain models address the difficulties of a HDBN by easily incorporating long time-delayed regulations. • Numerous time-delays are identified between the same pair of genes. The forward Viterbi path determines the best long-distant time delay between two genes. • Using priors gave us higher likelihood and improved the over-fitting in building the regulatory networks. • The work can be extended to skip-chain conditional random fields and fusion with protein interaction networks.
Thank You backup
- Slides: 16