Quasispecies, virus evolution, and lethal mutagenesis on realistic - - PowerPoint PPT Presentation

quasispecies virus evolution and lethal mutagenesis on
SMART_READER_LITE
LIVE PREVIEW

Quasispecies, virus evolution, and lethal mutagenesis on realistic - - PowerPoint PPT Presentation

Quasispecies, virus evolution, and lethal mutagenesis on realistic fitness landscapes Peter Schuster Institut fr Theoretische Chemie, Universitt Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA Seminar Lecture


slide-1
SLIDE 1
slide-2
SLIDE 2

Quasispecies, virus evolution, and lethal mutagenesis

  • n realistic fitness landscapes

Peter Schuster

Institut für Theoretische Chemie, Universität Wien, Austria and The Santa Fe Institute, Santa Fe, New Mexico, USA

Seminar Lecture Universität Marburg, 14.10.2011

slide-3
SLIDE 3

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

slide-4
SLIDE 4

Chemical kinetics of molecular evolution Historical prologue

The work on a molecular theory of evolution started more than 40 years ago ......

1971

slide-5
SLIDE 5

Manfred Eigen 1927 -

  

  

   

n i i n i i i j i n i ji j

x x f Φ n j Φ x x W x

1 1 1

, , 2 , 1 ; dt d 

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
slide-6
SLIDE 6

Evolution in the test tube: G.F. Joyce, Angew.Chem.Int.Ed. 46 (2007), 6420-6436

Sol Spiegelman, 1914 - 1983

slide-7
SLIDE 7

Kinetics of RNA replication

C.K. Biebricher, M. Eigen, W.C. Gardiner, Jr. Biochemistry 22:2544-2559, 1983

Christof K. Biebricher, 1941-2009

slide-8
SLIDE 8

C.K. Biebricher, R. Luce. 1992. In vitro recombination and terminal recombination of RNA by Q replicase. The EMBO Journal 11:5129-5135. stable

does not replicate!

metastable

replicates!

slide-9
SLIDE 9

RNA replication by Q-replicase

  • C. Weissmann, The making of a phage.

FEBS Letters 40 (1974), S10-S18

Charles Weissmann 1931-

slide-10
SLIDE 10

Chemical kinetics of molecular evolution (continued)

1977 1988

slide-11
SLIDE 11

Application of quasispecies theory to the fight against viruses Esteban Domingo 1943 -

slide-12
SLIDE 12

Error threshold versus lethal mutagenesis Vol.1(6), e61, 2005, pp.450 – 460.

slide-13
SLIDE 13

1. Complexity in molecular evolution 2. The error threshold 3. Simple landscapes and error thresholds 4. ‚Realistic‘ fitness landscapes 5. Quasispecies on realistic landscapes 6. Neutrality and consensus sequences 7. Error thresholds and lethal mutagenesis

slide-14
SLIDE 14
  • 1. Complexity in molecular evolution

2. The error threshold 3. Simple landscapes and error thresholds 4. ‚Realistic‘ fitness landscapes 5. Quasispecies on realistic landscapes 6. Neutrality and consensus sequences 7. Error thresholds and lethal mutagenesis

slide-15
SLIDE 15

Chemical kinetics of replication and mutation as parallel reactions

slide-16
SLIDE 16

Factorization of the value matrix W separates mutation and fitness effects.

   

   

     

n i i n i i i j i i n i ji j i n i ji j

x x f Φ n j Φ x x f Q Φ x x W x

1 1 1 1

, , 2 , 1 ; dt d 

slide-17
SLIDE 17

Mutation-selection equation: [Ii] = xi  0, fi  0, Qij  0 solutions are obtained after integrating factor transformation by means

  • f an eigenvalue problem

f x f x n i x x f Q dt dx

n j j j n i i i j j n j ij i

     

  

   1 1 1

; 1 ; , , 2 , 1 ,   

         

) ( ) ( ; , , 2 , 1 ; exp exp

1 1 1 1

   

     

      

n i i ki k n j k k n k jk k k n k ik i

x h c n i t c t c t x     

     

n j i h H L n j i L n j i Q f W

ij ij ij i

, , 2 , 1 , ; ; , , 2 , 1 , ; ; , , 2 , 1 , ;

1

          

 

1 , , 1 , ;

1

      

n k L W L

k

 

slide-18
SLIDE 18
slide-19
SLIDE 19

Perron-Frobenius theorem applied to the value matrix W

W is primitive: (i) is real and strictly positive (ii) (iii) is associated with strictly positive eigenvectors (iv) is a simple root of the characteristic equation of W (v-vi) etc. W is irreducible: (i), (iii), (iv), etc. as above (ii)

all for   k

k

 

  

all for   k

k

 

slide-20
SLIDE 20

Complexity in molecular evolution

W = G  F 0 , 0  largest eigenvalue and eigenvector

diagonalization of matrix W „ complicated but not complex “ fitness landscape mutation matrix „ complex “

genotype  phenotype mutation selection

slide-21
SLIDE 21

1. Complexity in molecular evolution

  • 2. The error threshold

3. Simple landscapes and error thresholds 4. ‚Realistic‘ fitness landscapes 5. Quasispecies on realistic landscapes 6. Neutrality and consensus sequences 7. Error thresholds and lethal mutagenesis

slide-22
SLIDE 22

The no-mutational backflow or zeroth order approximation

slide-23
SLIDE 23

The no-mutational backflow or zeroth order approximation

slide-24
SLIDE 24

Uniform error rate model

n jj X X d X X d n ji

p Q p p Q

j i H j i H

) 1 ( and ) 1 (

) , ( ) , (

   

p ... mutation rate per site and replication n ... chain length of the polynucleotide dH(Xi,Xj) ... Hamming distance

slide-25
SLIDE 25

 

 

       

                   

N m i i i i m m m m m n m m n m n m m m m mm m m mm m mm m m

f x x f f f p p x p Q x f Q t t f Q x x

, 1 / 1 cr 1 ) ( 1 1 ) ( ) ( ) (

) 1 ( 1 and ) ( 1 and ) 1 ( 1 ) 1 ( 1 1 1 ) ( and ) ( dt d         

The ‚no-mutational-backflow‘ or zeroth order approximation

slide-26
SLIDE 26

Chain length and error threshold

n p n p n p p n p Q

m m m m n m mm

     ln : constant ln : constant ln ) 1 ( ln 1 ) 1 (

max max

             

sequence master

  • f

y superiorit ) 1 ( length chain rate error accuracy n replicatio ) 1 (    

 

   

m j m j j m m n mm

x f x f σ n p p Q

slide-27
SLIDE 27

quasispecies

The error threshold in replication and mutation

driving virus populations through threshold

slide-28
SLIDE 28

1. Complexity in molecular evolution 2. The error threshold

  • 3. Simple landscapes and error thresholds

4. ‚Realistic‘ fitness landscapes 5. Quasispecies on realistic landscapes 6. Neutrality and consensus sequences 7. Error thresholds and lethal mutagenesis

slide-29
SLIDE 29

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.

slide-30
SLIDE 30

Sewall Wrights fitness landscape as metaphor for Darwinian evolution

Sewall Wright. 1932. The roles of mutation, inbreeding, crossbreeding and selection in evolution. In: D.F.Jones, ed. Int. Proceedings of the Sixth International Congress on Genetics. Vol.1, 356-366. Ithaca, NY.

slide-31
SLIDE 31

The simplified model

slide-32
SLIDE 32
slide-33
SLIDE 33

Model fitness landscapes I

single peak landscape step linear landscape

slide-34
SLIDE 34

Error threshold on the single peak landscape

slide-35
SLIDE 35

Error threshold on the step linear landscape

slide-36
SLIDE 36

Model fitness landscapes II linear and multiplicative hyperbolic

both are often used in population genetics

slide-37
SLIDE 37

The linear fitness landscape shows no error threshold

slide-38
SLIDE 38

Error threshold on the hyperbolic landscape

slide-39
SLIDE 39

Make things as simple as possible, but not simpler !

Albert Einstein

Albert Einstein‘s razor, precise refence is unknown.

slide-40
SLIDE 40

1. Complexity in molecular evolution 2. The error threshold 3. Simple landscapes and error thresholds

  • 4. ‚Realistic‘ fitness landscapes

5. Quasispecies on realistic landscapes 6. Neutrality and consensus sequences 7. Error thresholds and lethal mutagenesis

slide-41
SLIDE 41

Realistic fitness landscapes 1.Ruggedness: nearby lying genotypes may develop into very different phenotypes 2.Neutrality: many different genotypes give rise to phenotypes with identical selection behavior 3.Combinatorial explosion: the number of possible genomes is prohibitive for systematic searches

Facit: Any successful and applicable theory of molecular evolution must be able to predict evolutionary dynamics from a small or at least in practice measurable number of fitness values.

slide-42
SLIDE 42

Rugged fitness landscapes

  • ver individual binary sequences

with n = 10

single peak landscape „realistic“ landscape

slide-43
SLIDE 43

Random distribution of fitness values: d = 0.5 and s = 919

slide-44
SLIDE 44

Random distribution of fitness values: d = 1.0 and s = 919

slide-45
SLIDE 45

Random distribution of fitness values: d = 1.0 and s = 637

slide-46
SLIDE 46

1. Complexity in molecular evolution 2. The error threshold 3. Simple landscapes and error thresholds 4. ‚Realistic‘ fitness landscapes

  • 5. Quasispecies on realistic landscapes

6. Neutrality and consensus sequences 7. Error thresholds and lethal mutagenesis

slide-47
SLIDE 47

Error threshold: Individual sequences n = 10,  = 2, s = 491 and d = 0, 0.5, 0.9375

slide-48
SLIDE 48

Error threshold on ‚realistic‘ landscapes n = 10, f0 = 1.1, fn = 1.0, d = 0.5

s = 541 s = 919 s = 637

slide-49
SLIDE 49

Error threshold on ‚realistic‘ landscapes n = 10, f0 = 1.1, fn = 1.0, d = 0.5

s = 541 s = 637 s = 919

slide-50
SLIDE 50

s = 541 s = 919 s = 637

Error threshold on ‚realistic‘ landscapes n = 10, f0 = 1.1, fn = 1.0, d = 0.995

slide-51
SLIDE 51

s = 919 s = 541 s = 637

Error threshold on ‚realistic‘ landscapes n = 10, f0 = 1.1, fn = 1.0, d = 1.0

slide-52
SLIDE 52

Two questions: 1. Can we predict mutational behavior of quasispecies from fitness landscapes?

  • 2. What is the evolutionary consequence
  • f the occurrence of mutationally stable

and unstable quasispecies?

slide-53
SLIDE 53

Landscape analysis through the evaluation of single point mutation neighborhoods

slide-54
SLIDE 54

Landscape analysis through the evaluation of single point mutation neighborhoods

slide-55
SLIDE 55

Landscape analysis through the evaluation of single point mutation neighborhoods

slide-56
SLIDE 56

Landscape analysis through the evaluation of single point mutation neighborhoods

slide-57
SLIDE 57

Landscape analysis through the evaluation of single point mutation neighborhoods

slide-58
SLIDE 58

Landscape analysis through the evaluation of single point mutation neighborhoods

slide-59
SLIDE 59

Determination of the dominant mutation flow: d = 1 , s = 637

slide-60
SLIDE 60

Determination of the dominant mutation flow: d = 1 , s = 919

slide-61
SLIDE 61

1. Complexity in molecular evolution 2. The error threshold 3. Simple landscapes and error thresholds 4. ‚Realistic‘ fitness landscapes 5. Quasispecies on realistic landscapes

  • 6. Neutrality and consensus sequences

7. Error thresholds and lethal mutagenesis

slide-62
SLIDE 62

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.

slide-63
SLIDE 63

Motoo Kimura

Is the Kimura scenario correct for frequent mutations?

slide-64
SLIDE 64

Pairs of neutral sequences in replication networks

  • P. Schuster, J. Swetina. 1988. Bull. Math. Biol. 50:635-650

5 . ) ( ) ( lim

2 1

 

p x p x

p

dH = 1

) 1 ( 1 ) ( lim ) 1 ( ) ( lim

2 1

      

 

p x p x

p p

dH = 2

Random fixation in the sense of Motoo Kimura

dH  3

1 ) ( lim , ) ( lim

  • r

) ( lim , 1 ) ( lim

2 1 2 1

   

   

p x p x p x p x

p p p p

slide-65
SLIDE 65

A fitness landscape including neutrality

slide-66
SLIDE 66

Neutral network: Individual sequences n = 10,  = 1.1, d = 1.0

slide-67
SLIDE 67

Consensus sequence of a quasispecies of two strongly coupled sequences of Hamming distance dH(Xi,,Xj) = 1.

slide-68
SLIDE 68

Neutral network: Individual sequences n = 10,  = 1.1, d = 1.0

slide-69
SLIDE 69

Consensus sequence of a quasispecies of two strongly coupled sequences of Hamming distance dH(Xi,,Xj) = 2.

slide-70
SLIDE 70

N = 7 Neutral networks with increasing :  = 0.10, s = 229

Adjacency matrix

slide-71
SLIDE 71

All ‚realistic‘ fitness landscapes – rugged with and without neutrality – show error thresholds for replication.

slide-72
SLIDE 72

1. Complexity in molecular evolution 2. The error threshold 3. Simple landscapes and error thresholds 4. ‚Realistic‘ fitness landscapes 5. Quasispecies on realistic landscapes 6. Neutrality and consensus sequences

  • 7. Error thresholds and lethal mutagenesis
slide-73
SLIDE 73

The flowreactor is a device for studying evolution in vitro and in silico Stock solution: [A] = a = a0 Flow rate: r = R

  • 1

The population size N , the number of polynucleotide molecules, is controlled by the flow r N N t N   ) ( Enzyme immobilized

slide-74
SLIDE 74
slide-75
SLIDE 75
slide-76
SLIDE 76
slide-77
SLIDE 77
slide-78
SLIDE 78
slide-79
SLIDE 79

Extiction threshold pext

capacity carrying with ln 1 ln 1

ext

 r a f n r a f n p     

slide-80
SLIDE 80
slide-81
SLIDE 81
  • H. Tejero, A. Marín, and F. Montero. 2010. J. theoret. Biol. 262:733-741.

Three species model: Xm ... master, Xk ... viable, Xj ... lethal

slide-82
SLIDE 82

Error threshold and extinction threshold

  • H. Tejero, A. Marín, and
  • F. Montero. 2010.

J.Theoret.Biol.262:733-741

slide-83
SLIDE 83

Theory cannot remove complexity, but it shows what kind of „regular“ behavior can be expected and what experiments have to be done to get a grasp on the irregularities. Manfred Eigen,

Preface to E. Domingo, C.R. Parrish, J.J.Holland, eds. Origin and Evolution of Viruses. Academic Press 2008

slide-84
SLIDE 84

Coworkers

Peter Stadler, Bärbel M. Stadler, Universität Leipzig, GE Paul E. Phillipson, University of Colorado at Boulder, CO Heinz Engl, Philipp Kügler, James Lu, Stefan Müller, RICAM Linz, AT Jord Nagel, Kees Pleij, Universiteit Leiden, NL Walter Fontana, Harvard Medical School, MA Martin Nowak, Harvard University, MA Christian Reidys, Nankai University, Tien Tsin, China Christian Forst, Los Alamos National Laboratory, NM Thomas Wiehe, Ulrike Göbel, Walter Grüner, Stefan Kopp, Jaqueline Weber, Institut für Molekulare Biotechnologie, Jena, GE Ivo L.Hofacker, Christoph Flamm, Andreas Svrček-Seiler, Universität Wien, AT Kurt Grünberger, Michael Kospach , Andreas Wernitznig, Stefanie Widder, Stefan Wuchty, Jan Cupal, Stefan Bernhart, Lukas Endler, Ulrike Langhammer, Rainer Machne, Ulrike Mückstein, Erich Bornberg-Bauer, Universität Wien, AT

Universität Wien

slide-85
SLIDE 85

Universität Wien

Acknowledgement of support

Fonds zur Förderung der wissenschaftlichen Forschung (FWF) Projects No. 09942, 10578, 11065, 13093 13887, and 14898 Wiener Wissenschafts-, Forschungs- und Technologiefonds (WWTF) Project No. Mat05 Jubiläumsfonds der Österreichischen Nationalbank Project No. Nat-7813 European Commission: Contracts No. 98-0189, 12835 (NEST) Austrian Genome Research Program – GEN-AU: Bioinformatics Network (BIN) Österreichische Akademie der Wissenschaften Siemens AG, Austria Universität Wien and the Santa Fe Institute

slide-86
SLIDE 86

Thank you for your attention!

slide-87
SLIDE 87

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

slide-88
SLIDE 88