error thresholds on realistic fitness landscapes
play

Error thresholds on realistic fitness landscapes Peter Schuster - PowerPoint PPT Presentation

Error thresholds on realistic fitness landscapes Peter Schuster Institut fr Theoretische Chemie, Universitt Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Evolutionary Dynamics: Mutation, Selection, and the Origin of


  1. Error thresholds on realistic fitness landscapes Peter Schuster Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Evolutionary Dynamics: Mutation, Selection, and the Origin of Information University of Utrecht, 07.04.2010

  2. Web-Page for further information: http://www.tbi.univie.ac.at/~pks

  3. Prologue The work on a molecular theory of evolution started 42 years ago ...... 1971 1977 1988 Chemical kinetics of molecular evolution

  4. 1. Open system and constant population size 2. Chemical kinetics of replication and mutation 3. Complexity of fitness landscapes 4. Quasispecies on realistic landscapes 5. Neutrality and replication 6. Lethal variants and mutagenesis

  5. 1. Open system and constant population size 2. Chemical kinetics of replication and mutation 3. Complexity of fitness landscapes 4. Quasispecies on realistic landscapes 5. Neutrality and replication 6. Lethal variants and mutagenesis

  6. Stock solution : [A] = a = a 0 Flow rate : r = � R -1 The population size N , the number of polynucleotide molecules, is controlled by the flow r ≈ ± N t N N ( ) The flowreactor is a device for studying evolution in vitro and in silico

  7. Replication and mutation in the flowreactor

  8. Derivation of the replication-mutation equation from the flowreactor

  9. 1. Open system and constant population size 2. Chemical kinetics of replication and mutation 3. Complexity of fitness landscapes 4. Quasispecies on realistic landscapes 5. Neutrality and replication 6. Lethal variants and mutagenesis

  10. x d ∑ n j = − = W x x Φ j n ; 1 , 2 , , K ji i j = i 1 dt ∑ ∑ n n = Φ f x x i i i = = i i 1 1 Manfred Eigen 1927 - Mutation and (correct) replication as parallel chemical reactions M. Eigen. 1971. Naturwissenschaften 58:465, M. Eigen & P. Schuster.1977. Naturwissenschaften 64:541, 65:7 und 65:341

  11. Taq = thermus aquaticus Accuracy of replication: Q = q 1 · q 2 · q 3 · … · q n The logics of DNA replication

  12. RNA replication by Q � -replicase C. Weissmann, The making of a phage . FEBS Letters 40 (1974), S10-S18

  13. dx dx = = 1 f x 2 f x and 2 2 1 1 dt dt = ξ = ξ ζ = ξ + ξ η = ξ − ξ = x f x f f f f , , , , 1 2 1 2 1 2 1 2 1 2 1 2 − η = η ft t e ( ) ( 0 ) ζ = ζ ft t e ( ) ( 0 ) Complementary replication as the simplest molecular mechanism of reproduction

  14. Christof K. Biebricher, 1941-2009 Kinetics of RNA replication C.K. Biebricher, M. Eigen, W.C. Gardiner, Jr. Biochemistry 22 :2544-2559, 1983

  15. x d ∑ ∑ n n j = − = − = W x x Φ Q f x x Φ j n ; 1 , 2 , , K ji i j ji i i j = = i i dt 1 1 ∑ ∑ n n = Φ f x x i i i = = i 1 i 1 Factorization of the value matrix W separates mutation and fitness effects.

  16. Mutation-selection equation : [I i ] = x i � 0, f i � 0, Q ij � 0 dx ∑ ∑ ∑ n n n = − φ = = φ = = i Q f x x i n x f x f , 1 , 2 , , ; 1 ; L ij j j i i j j = = = dt j 1 i 1 j 1 solutions are obtained after integrating factor transformation by means of an eigenvalue problem ( ) ( ) ∑ − n 1 ⋅ ⋅ λ c t 0 exp l ( ) ∑ n ik k k = = = = k x t 0 i n c h x ; 1 , 2 , , ; ( 0 ) ( 0 ) L ( ) ( ) ∑ ∑ i − k ki i n n 1 = i ⋅ ⋅ λ 1 c t 0 exp l jk k k = = j k 1 0 { } { } { } − ÷ = = = = = = 1 W f Q i j n L i j n L H h i j n ; , 1 , 2 , , ; ; , 1 , 2 , , ; ; , 1 , 2 , , L l L L i ij ij ij { } − ⋅ ⋅ = Λ = λ = − 1 L W L k n ; 0 , 1 , , 1 L k

  17. Perron-Frobenius theorem applied to the value matrix W λ W is primitive: (i) is real and strictly positive 0 (ii) λ > λ ≠ k for all 0 k 0 λ (iii) is associated with strictly positive eigenvectors 0 (iv) is a simple root of the characteristic equation of W λ 0 (v-vi) etc. W is irreducible: (i), (iii), (iv), etc. as above λ ≥ λ ≠ k (ii) for all 0 k 0

  18. constant level sets of � Selection of quasispecies with f 1 = 1.9, f 2 = 2.0, f 3 = 2.1, and p = 0.01 , parametric plot on S 3

  19. Chain length and error threshold ⋅ σ = − ⋅ σ ≥ ⇒ ⋅ − ≥ − σ n Q p n p ( 1 ) 1 ln ( 1 ) ln σ ln ≈ p n constant : K max p σ ln ≈ n p K constant : max n = − n Q p ( 1 ) replicatio n accuracy K p error rate K n chain length K f = m σ superiorit y of master sequence K ∑ ≠ f j j m

  20. Quasispecies Uniform distribution Stationary population or quasispecies as a function of the mutation or error rate p 0.00 0.05 0.10 Error rate p = 1-q

  21. Eigenvalues of the matrix W as a function of the error rate p

  22. Quasispecies Driving virus populations through threshold The error threshold in replication: No mutational backflow approximation

  23. Molecular evolution of viruses

  24. 1. Open system and constant population size 2. Chemical kinetics of replication and mutation 3. Complexity of fitness landscapes 4. Quasispecies on realistic landscapes 5. Neutrality and replication 6. Lethal variants and mutagenesis

  25. � 0 , � 0 � largest eigenvalue and eigenvector diagonalization of matrix W „ complicated but not complex “ � = G W F mutation matrix fitness landscape „ complex “ ( complex ) sequence structure � „ complex “ mutation selection Complexity in molecular evolution

  26. Sewall Wright. 1931. Evolution in Mendelian populations. Genetics 16:97-159. -- --. 1932. The roles of mutation, inbreeding, crossbreeding, and selection in evolution. In: D.F.Jones, ed. Proceedings of the Sixth International Congress on Genetics, Vol.I. Brooklyn Botanical Garden. Ithaca, NY, pp. 356-366. -- --. 1988. Surfaces of selective value revisited. The American Naturalist 131:115-131.

  27. Build-up principle of binary sequence spaces

  28. Build-up principle of four letter (AUGC) sequence spaces

  29. linear and multiplicative hyperbolic Smooth fitness landscapes

  30. The linear fitness landscape shows no error threshold

  31. Error threshold on the hyperbolic landscape

  32. step linear landscape single peak landscape Rugged fitness landscapes

  33. Error threshold on the single peak landscape

  34. Error threshold on a single peak fitness landscape with n = 50 and � = 10

  35. Error threshold on the step linear landscape

  36. The error threshold can be separated into three phenomena: 1. Decrease in the concentration of the master sequence to very small values. 2. Sharp change in the stationary concentration of the quasispecies distribuiton. 3. Transition to the uniform distribution at small mutation rates. All three phenomena coincide for the quasispecies on the single peak fitness lanscape.

  37. The error threshold can be separated into three phenomena: 1. Decrease in the concentration of the master sequence to very small values. 2. Sharp change in the stationary concentration of the quasispecies distribuiton. 3. Transition to the uniform distribution at small mutation rates. All three phenomena coincide for the quasispecies on the single peak fitness lanscape.

  38. 1. Open system and constant population size 2. Chemical kinetics of replication and mutation 3. Complexity of fitness landscapes 4. Quasispecies on realistic landscapes 5. Neutrality and replication 6. Lethal variants and mutagenesis

  39. single peak landscape „realistic“ landscape Rugged fitness landscapes over individual binary sequences with n = 10

  40. Error threshold: Individual sequences n = 10, � = 2, s = 491 and d = 0, 1.0, 1.875

  41. Shift of the error threshold with increasing ruggedness of the fitness landscape

  42. d = 0.100 d = 0.200 n = 10, f 0 = 1.1, f n = 1.0, s = 919 Case I : Strong Quasispecies

  43. Case II : Dominant single transition d = 0.190 d = 0.190 n = 10, f 0 = 1.1, f n = 1.0, s = 541

  44. Case II : Dominant single transition d = 0.190 d = 0.195 n = 10, f 0 = 1.1, f n = 1.0, s = 541

  45. Case II : Dominant single transition d = 0.199 d = 0.199 n = 10, f 0 = 1.1, f n = 1.0, s = 541

  46. d = 0.100 d = 0.195 n = 10, f 0 = 1.1, f n = 1.0, s = 637 Case III : Multiple transitions

  47. d = 0.199 d = 0.200 n = 10, f 0 = 1.1, f n = 1.0, s = 637 Case III : Multiple transitions

  48. 1. Open system and constant population size 2. Chemical kinetics of replication and mutation 3. Complexity of fitness landscapes 4. Quasispecies on realistic landscapes 5. Neutrality and replication 6. Lethal variants and mutagenesis

  49. Motoo Kimuras population genetics of neutral evolution. Evolutionary rate at the molecular level. Nature 217 : 624-626, 1955. The Neutral Theory of Molecular Evolution . Cambridge University Press. Cambridge, UK, 1983.

  50. Motoo Kimura Is the Kimura scenario correct for frequent mutations?

Recommend


More recommend