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

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

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

Prologue

The work on a molecular theory of evolution started 42 years ago ......

Chemical kinetics of molecular evolution

1988 1971 1977

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

• 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

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 ± ≈ ) ( The flowreactor is a device for studying evolution in vitro and in silico

Replication and mutation in the flowreactor

Derivation of the replication-mutation equation from the flowreactor

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

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 K

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

Taq = thermus aquaticus

Accuracy of replication: Q = q1 · q2 · q3 · … · qn

The logics of DNA replication

RNA replication by Q-replicase

• C. Weissmann, The making of a phage.

FEBS Letters 40 (1974), S10-S18

1 1 2 2 2 1

and x f dt dx x f dt dx = =

2 1 2 1 2 1 2 1 2 1 2 1

, , , , f f f f x f x = − = + = = = ξ ξ η ξ ξ ζ ξ ξ

ft ft

e t e t ) ( ) ( ) ( ) ( ζ ζ η η = =

−

Complementary replication as the simplest molecular mechanism of reproduction

Christof K. Biebricher, 1941-2009

Kinetics of RNA replication

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

∑ ∑ ∑ ∑

= = = =

= = − = − =

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 K

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

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 , φ φ L

( ) ( ) ( ) ( ) ( )

) ( ) ( ; , , 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 L l l λ λ

{ } { } { }

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

L L l L = = = = = = ÷

−

{ }

1 , , 1 , ;

1

− = = Λ = ⋅ ⋅

−

n k L W L

k

L λ

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

λ λ

constant level sets of

Selection of quasispecies with f1 = 1.9, f2 = 2.0, f3 = 2.1, and p = 0.01 , parametric plot on S3

Chain length and error threshold

n p n p n p p n p Q

n

σ σ σ σ σ ln : constant ln : constant ln ) 1 ( ln 1 ) 1 (

max max

≈ ≈ − ≥ − ⋅ ⇒ ≥ ⋅ − = ⋅ K K

sequence master

• f

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

∑ ≠

= − =

m j j m n

f f σ n p p Q

Error rate p = 1-q

0.00 0.05 0.10

Quasispecies Uniform distribution

Stationary population or quasispecies as a function

• f the mutation or error

rate p

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

Quasispecies

Driving virus populations through threshold

The error threshold in replication: No mutational backflow approximation

Molecular evolution of viruses

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

W = G

• F

0 , 0 largest eigenvalue and eigenvector

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

sequence

• structure

„ complex “

mutation selection

Complexity in molecular evolution

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.

Build-up principle of binary sequence spaces

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

linear and multiplicative Smooth fitness landscapes hyperbolic

The linear fitness landscape shows no error threshold

Error threshold on the hyperbolic landscape

single peak landscape

Rugged fitness landscapes

step linear landscape

Error threshold on the single peak landscape

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

Error threshold on the step linear landscape

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

• f the quasispecies distribuiton.

3. Transition to the uniform distribution at small mutation rates. All three phenomena coincide for the quasispecies

• n the single peak fitness lanscape.

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

• f the quasispecies distribuiton.

3. Transition to the uniform distribution at small mutation rates. All three phenomena coincide for the quasispecies

• n the single peak fitness lanscape.

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

single peak landscape

Rugged fitness landscapes

• ver individual binary sequences

with n = 10

„realistic“ landscape

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

Shift of the error threshold with increasing ruggedness

• f the fitness landscape

d = 0.100

Case I: Strong Quasispecies n = 10, f0 = 1.1, fn = 1.0, s = 919

d = 0.200

d = 0.190

Case II: Dominant single transition n = 10, f0 = 1.1, fn = 1.0, s = 541

d = 0.190

d = 0.190

Case II: Dominant single transition n = 10, f0 = 1.1, fn = 1.0, s = 541

d = 0.195

d = 0.199

Case II: Dominant single transition n = 10, f0 = 1.1, fn = 1.0, s = 541

d = 0.199

d = 0.100

Case III: Multiple transitions n = 10, f0 = 1.1, fn = 1.0, s = 637

d = 0.195

d = 0.199

Case III: Multiple transitions n = 10, f0 = 1.1, fn = 1.0, s = 637

d = 0.200

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

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.

Motoo Kimura

Is the Kimura scenario correct for frequent mutations?

5 . ) ( ) ( lim

2 1

= =

→

p x p x

p

dH = 1

a p x a p x

p p

− = =

→ →

1 ) ( lim ) ( lim

2 1

dH = 2 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

Random fixation in the sense of Motoo Kimura Pairs of neutral sequences in replication networks

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

A fitness landscape including neutrality

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

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

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

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

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

Adjacency matrix

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

Coworkers

Peter Stadler, Bärbel M. Stadler, Universität Leipzig, GE Walter Fontana, Harvard Medical School, MA Martin Nowak, Harvard University, MA Christian Reidys, Nankai University, Tien Tsin, China Thomas Wiehe, Ulrike Göbel, Walter Grüner, Stefan Kopp, Jaqueline Weber, Institut für Molekulare Biotechnologie, Jena, GE Ivo L.Hofacker, Christoph Flamm, 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

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

Universität Wien

Thank you for your attention!

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