1 / 75

341: Introduction to Bioinformatics

341: Introduction to Bioinformatics. Dr. Nataša Pržulj Department of Comput ing Imperial College London natasha@imperial.ac.uk Winter 2011. Topics. Introduction to biology (cell, DNA, RNA, genes, proteins) Sequencing and genomics (sequencing technology, sequence alignment algorithms)

analiese
Download Presentation

341: Introduction to Bioinformatics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 341: Introduction to Bioinformatics Dr. Nataša Pržulj Department of Computing Imperial College London natasha@imperial.ac.uk Winter 2011

  2. Topics • Introduction to biology (cell, DNA, RNA, genes, proteins) • Sequencing and genomics (sequencing technology, sequence alignment algorithms) • Functional genomics and microarray analysis (array technology, statistics, clustering and classification) • Introduction to biological networks • Introduction to graph theory • Network properties • Global: network/node centralities • Local: network motifs and graphlets • Network models • Network/node clustering • Network comparison/alignment • Software tools for network analysis • Interplay between topology and biology 2 2

  3. Topics • Introduction to biology (cell, DNA, RNA, genes, proteins) • Sequencing and genomics (sequencing technology, sequence alignment algorithms) • Functional genomics and microarray analysis (array technology, statistics, clustering and classification) • Introduction to biological networks • Introduction to graph theory • Network properties • Global: network/node centralities • Local: network motifs and graphlets • Network models • Network/node clustering • Network comparison/alignment • Software tools for network analysis • Interplay between topology and biology 3 3

  4. Network properties: summary of last class Network Comparisons: • Large network comparison is computationally hard due to NP-completeness of the underlying subgraph isomorphism problem: • Given 2 graphs G and H as input, determine whether G contains a subgraph that is isomorphic to H. • Thus, network comparisons rely on easily computable heuristics (approximate solutions), called “network properties” • Network properties can roughly & historically be divided in two categories: • Global network properties: give an overall view of the network, but might not be detailed enough to capture complex topological characteristics of large networks. • Local network properties: more detailed network descriptors which usually encompass larger number of constraints, thus reducing degrees of freedom in which the networks being compared can vary. 4

  5. Network properties: summary of last class 1. Global Network Properties • Readings: Chapter 3 of “Analysis of biological networks” by Junker and Schreiber. • Global Network Properties: • Degree distribution • Average clustering coefficient • Clustering spectrum • Average Diameter • Spectrum of shortest path lengths • Centralities

  6. Network properties: summary of last class • 2. Local Network Properties • Readings: Chapter 5 of “Analysis of Biological Networks” by Junker and Schreiber. • Network motifs • Graphlets Two network comparison measures based on graphlets: • 2.1) Relative Graphlet Frequency Distance between two networks • 2.2) Graphlet Degree Distribution Agreement between two networks

  7. 1) Network motifs (Uri Alon’s group, ’02-’04) http://www.weizmann.ac.il/mcb/UriAlon/ Also, see Pajek, MAVisto, and FANMOD

  8. 2) Graphlets 2.1) Reltive graphlet frequency distance between two networks N. Przulj, D. G. Corneil, and I. Jurisica, “Modeling Interactome: Scale Free or Geometric?,” Bioinformatics, vol. 20, num. 18, pg. 3508-3515, 2004.

  9. 2) Graphlets 2.1) Graphlet degree distribution agreement between two networks N. Przulj, “Biological Network Comparison Using Graphlet Degree Distribution,” ECCB, Bioinformatics, vol. 23, pg. e177-e183, 2007.

  10. 2) Graphlets 2.1) Graphlet degree distribution agreement between two networks Graphlet Degree (GD) vectors, or “node signatures” T. Milenkovic and N. Przulj, “Uncovering Biological Network Function via Graphlet Degree Signatures”, Cancer Informatics, vol. 4, pg. 257-273, 2008.

  11. 2) Graphlets 2.1) Graphlet degree distribution agreement between two networks Signature Similarity Measure between nodes u and v T. Milenkovic and N. Przulj, “Uncovering Biological Network Function via Graphlet Degree Signatures”, Cancer Informatics, vol. 4, pg. 257-273, 2008.

  12. Software that implements many of these network properties and compares networks with respect to them: GraphCrunch http://bio-nets.doc.ic.ac.uk/graphcrunch/

  13. Software that implements many of these network properties and compares networks with respect to them: GraphCrunch http://bio-nets.doc.ic.ac.uk/graphcrunch2/

  14. Software that implements many of these network properties and compares networks with respect to them: GraphCrunch http://bio-nets.doc.ic.ac.uk/graphcrunch2/

  15. Another Software: Cytoscape http://www.cytoscape.org/

  16. Examples of signatures and signature similarities: T. Milenković and N. Pržulj, “Uncovering Biological Network Function via Graphlet Degree Signatures,” Cancer Informatics, 2008:6 257-273, 2008 (Highly Visible).

  17. Examples of signatures and signature similarities: SMD1 YBR095C 40% PMA1 T. Milenković and N. Pržulj, “Uncovering Biological Network Function via Graphlet Degree Signatures,” Cancer Informatics, 2008:6 257-273, 2008 (Highly Visible).

  18. Examples of signatures and signature similarities: T. Milenković and N. Pržulj, “Uncovering Biological Network Function via Graphlet Degree Signatures,” Cancer Informatics, 2008:6 257-273, 2008 (Highly Visible).

  19. Examples of signatures and signature similarities: 90%* SMD1 RPO26 SMB1 *Statistically significant threshold at ~85% T. Milenković and N. Pržulj, “Uncovering Biological Network Function via Graphlet Degree Signatures,” Cancer Informatics, 2008:6 257-273, 2008 (Highly Visible).

  20. Later we will see how to use this and other techniques to link network structure with biological function

  21. Generalize Degree Distribution of a network • The degree distribution measures: • the number of nodes “touching” k edges for each value of k N. Przulj, “Biological Network Comparison Using Graphlet Degree Distribution,” Bioinformatics, vol. 23, pg. e177-e183, 2007.

  22. N. Przulj, “Biological Network Comparison Using Graphlet Degree Distribution,” Bioinformatics, vol. 23, pg. e177-e183, 2007.

  23. N. Przulj, “Biological Network Comparison Using Graphlet Degree Distribution,” Bioinformatics, vol. 23, pg. e177-e183, 2007.

  24. / sqrt(2) ( to make it between 0 and 1) This is called Graphlet Degree Distribution (GDD) Agreement between networks G and H.

  25. Software that implements many of these network properties and compares networks with respect to them: GraphCrunch http://bio-nets.doc.ic.ac.uk/graphcrunch/

  26. Software that implements many of these network properties and compares networks with respect to them: GraphCrunch http://bio-nets.doc.ic.ac.uk/graphcrunch2/

  27. Topics • Introduction to biology (cell, DNA, RNA, genes, proteins) • Sequencing and genomics (sequencing technology, sequence alignment algorithms) • Functional genomics and microarray analysis (array technology, statistics, clustering and classification) • Introduction to biological networks • Introduction to graph theory • Network properties • Network/node centralities • Network motifs • Network models • Network/node clustering • Network comparison/alignment • Software tools for network analysis • Interplay between topology and biology 27 27

  28. What is a network (graph) model?

  29. Does the model network fit the data? • Use network properties: • Local • Global • Why? • “Hardness” of graph theoretic problems • E.g. NP-completeness of subgraph isomorphism • Cannot exactly compare/align networks • Use heuristics (approximate solutions) • Exact comparison inappropriate in biology • Due to biological variation • Noise revise models as data sets evolve

  30. Why model networks? • Understand laws  reproduction/predictions • Network models have already been used in biological applications: • Network motifs (Shen-Orr et al., Nature Genetics 2002, Milo et al., Science 2002) • De-noising of PPI network data (Kuchaiev et al., PLoS Comp. Biology, 2009) • Guiding biological experiments (Lappe and Holm, Nature Biotechnology, 2004) • Development of computationally easy algorithms for PPI nets that are computationally intensive on graphs in general(Przulj et al., Bioinformatics, 2006)

  31. Network models We will cover the following network models: • Erdos–Renyi random graphs • Generalized random graphs (with the same degree distribution as the data networks) • Small-world networks • Scale-free networks • Hierarchical model • Geometric random graphs • Stickiness index-based network model

  32. Erdos–Renyi random graphs (ER) • Model a data network G(V,E) with |V|=n and |E|=m • An ER graph that models G is constructed as follows: • It has n nodes • Edges are added between pairs of nodes uniformly at random with the same probability p • Two (equivalent) methods for constructing ER graphs: • Gn,p: pick p so that the resulting model network has m edges • Gn,m: pick randomly m pairs of nodes and add edges between them with probability 1

  33. Erdos–Renyi random graphs (ER) • Number of edges, |E|=m, in Gn,pis: • Average degree is:

  34. Erdos–Renyi random graphs (ER) • Many properties of ER can be proven theoretically (See: Bollobas, "Random Graphs," 2002) • Example: • When m=n/2,suddenly the giant component emerges, i.e.: • One connected component of the network has O(n) nodes • The next largest connected component has O(log(n)) nodes

  35. Erdos–Renyi random graphs (ER) • The degree distribution is binomial: • For large n, this can be approximated with Poisson distribution: where z is the average degree • However, currently available biological networkshave power-law degree distribution

  36. Erdos–Renyi random graphs (ER) • Clustering coefficient, C, of ER is low (for low p) • C=p, since probability p of connecting any two nodes in an ER graph is the same, regardless of whether the nodes are neighbors • However, biological networkshave high clustering coefficients

  37. Erdos–Renyi random graphs (ER) • Average diameter of ER graphs is small • It is equal to • Biological networks alsohave small average diameters • Summary

  38. Generalized random graphs (ER-DD) • Preserve the degree distribution of data (“ER-DD”) • Constructed as follows: • An ER-DD network has n nodes (so does the data) • Edges are added between pairs of nodes using the “stubs method”

  39. Generalized random graphs (ER-DD) • The “stubs method” for constructing ER-DD graphs: • The number of “stubs” (to be filled by edges) is assigned to each node in the model network according to the degree distribution of the real network to be modeled • Edges are created between pairs of nodes with “available” stubs picked at random • After an edge is created, the number of stubs left available at the corresponding “end nodes” of the edges is decreased by one • Multiple edges between the same pair of nodes are not allowed

  40. Generalized random graphs (ER-DD) • Summary • 2 global network properties are matched by ER-DD • How about local network properties (graphlet frequencies)? • Low-density graphlets are over-represented in ER and ER-DD • However, data have lots of dense graphlets, since they have high clustering coefficients

  41. Small-world networks (SW) • Watts and Strogatz, 1998 • Created from regular ring lattices by random rewiring of a small percentage of their edges • E.g.

  42. Small-world networks (SW) • SW networks have: • High clustering coefficients – introduced by “ring regularity” • Large average diameters of regular lattices – fixed by randomly re-wiring a small percentage of edges • Summary

  43. Scale-free networks (SF) • Power-law degree distributions: P(k) = k−γ • γ > 0; 2 < γ < 3

  44. Scale-free networks (SF) • Power-law degree distributions: P(k) = k−γ • γ > 0; 2 < γ < 3

  45. Scale-free networks (SF) • Different models exist, e.g.: • Preferential Attachment Model (SF-BA) (Barabasi-Albert, 1999) • Gene Duplication and Mutation Model (SF-GD) (Vazquez et al., 2003)

  46. Scale-free networks (SF) • Preferential Attachment Model (SF-BA) • “Growth” model: nodes are added to an existing network • New nodes preferentially attach to existing nodes with probability proportional to the degrees of the existing nodes; e.g.: • This is repeated until the size of SF network matches the size of the data • “Rich getting richer” • The starting network strongly influences the properties of the resulting network (F. Hormozdiari, et al., PLoS Computational Biology, 3(7):e118, July 2007. ) • SF-BA: particularly effective at describing Internet

  47. Scale-free networks (SF) • Gene Duplication and Mutation Model (SF-GD) • Biologically motivated • Attempts to mimic gene duplication and mutation processes

  48. Scale-free networks (SF) • Gene Duplication and Mutation Model (SF-GD) • At each time step, a node is added to the network as follows:

  49. Scale-free networks (SF) • Summary

More Related