Skip to main content
Sandbox Physics
Evolution as statistical physics

L3 · Chapter 5

Does a big sector mean selection? Opening angles and luck in colonies

A colony seeded with twelve colours ends with one of them holding forty per cent of the edge. Is that colour fitter? It might be. A faster-growing mutant opens a wedge whose angle depends on how much faster it grows. But the previous chapter showed that drift at the edge also produces big sectors with no selection at all. This chapter separates the two with the numbers from the Lab.

Selection gives the boundaries a drift

For a neutral sector, each boundary wanders with no preferred side, so the expected change in sector width XX is zero. Hallatschek and Nelson add selection as a mean drift: each boundary of a beneficial sector moves outwards by m⊥m_{\perp} per unit of front advance, while the random part keeps the diffusion constant DXD_X of the neutral case.

Once the sector is wide enough that drift wins, its boundaries become straight lines and the sector becomes a wedge. The drift then equals the tangent of half the opening angle Φ\Phi.

⟨ΔX⟩=2m⊥ Δr,⟨ΔX2⟩=4DX Δr,m⊥=v⊥v=tan⁡Φ2\langle\Delta X\rangle=2m_{\perp}\,\Delta r,\qquad \langle\Delta X^2\rangle=4D_X\,\Delta r,\qquad m_{\perp}=\frac{v_{\perp}}{v}=\tan\frac{\Phi}{2}
Mean and variance of the change in sector width over a front advance Δr\Delta r, and the drift written as an angle.

Source: O. Hallatschek and D. R. Nelson

A geometric rule for the angle

Suppose the only effect of the mutation is that mutant regions grow into empty space at speed v⋆v^{\star} instead of vv. Where the wedge meets the wild-type front, both populations must reach the kink at the same time. That gives the cosine rule on the right.

How v⋆/vv^{\star}/v depends on the growth advantage ss is a property of the growth process. In the Lab every gold cell divides at rate 1+s1+s instead of 11. A colony made only of gold cells is the same random process run 1+s1+s times faster, so v⋆/v=1+sv^{\star}/v=1+s exactly. Hallatschek and Nelson give this linear rule for fronts where drift dominates, and a square-root rule for smooth Fisher waves.

cos⁡Φ2=vv⋆=11+s\cos\frac{\Phi}{2}=\frac{v}{v^{\star}}=\frac{1}{1+s}
Opening angle of a beneficial sector when the mutant front is 1+s1+s times faster.

Source: O. Hallatschek and D. R. Nelson

How often a fitter founder survives

A new mutant starts as a sector one cell wide, and drift can close it before the drift of its boundaries opens it. With the biased walk above, the chance that a sector of width x0x_0 on a front of width LL eventually takes over has the same form as Kimura’s formula from chapter 2. Hallatschek and Nelson point this out in a footnote: the product of population size and selection coefficient is replaced by m⊥L/(2DX)m_{\perp}L/(2D_X).

Sectors much narrower than the establishment length DX/m⊥D_X/m_{\perp} survive with probability close to m⊥x0/DXm_{\perp}x_0/D_X; wider ones almost always survive. In the Lab strip, one gold founder reached the far end of the strip in 3.3±0.93.3\pm0.9 per cent of 400400 seeded runs when neutral, 13.5±1.713.5\pm1.7 per cent at s=0.1s=0.1 and 29.5±2.329.5\pm2.3 per cent at s=0.2s=0.2.

u(x0)=1−e−m⊥x0/DX1−e−m⊥L/DXu(x_0)=\frac{1-e^{-m_{\perp}x_0/D_X}}{1-e^{-m_{\perp}L/D_X}}
Long-time survival probability of a beneficial sector of initial width x0x_0 at a straight front of width LL.

Source: O. Hallatschek and D. R. Nelson

Worked example: reading an advantage from an angle

Problem. A gold block 1010 cells wide, with s=0.1s=0.1, starts on a straight front 700700 cells wide. What half-angle does the geometric rule predict? If you only had the measured angle, what advantage would you infer?

Prediction. cos⁡(Φ/2)=1/1.1\cos(\Phi/2)=1/1.1, so tan⁡(Φ/2)=1.12−1=0.458\tan(\Phi/2)=\sqrt{1.1^2-1}=0.458 and Φ=49∘\Phi=49^{\circ}.

Measurement. In 2020 of 2424 seeded runs the wedge survived. Its half-angle tangent, measured from the change in wedge width between bands of rows, rose from 0.296±0.0200.296\pm0.020 between rows 5050 and 150150 to 0.392±0.0140.392\pm0.014 between rows 450450 and 600600. At s=0.2s=0.2 it rose from 0.5100.510 to 0.6030.603 against 0.6630.663 predicted.

Inference. Reading 0.3920.392 as the final angle gives 1+s=1+0.3922=1.0741+s=\sqrt{1+0.392^2}=1.074, so s=0.074s=0.074, a quarter below the true value. The wedge is still opening at the end of the strip. The rule assumes a steady angle after transients, and in this lattice model the transient is long.

Answer. The geometric rule predicts tan⁡(Φ/2)=0.458\tan(\Phi/2)=0.458, and angles measured over a few hundred rows give only a lower bound on the advantage.

Source: O. Hallatschek and D. R. Nelson

Luck alone makes big sectors

Now remove selection. On the Lab strip of width 200200 and height 240240, twelve neutral colours start in equal numbers, a fair share of 1/12=0.0831/12=0.083 each. Over 200200 seeded colonies the largest colour at the final front held a median share of 0.410.41, between 0.310.31 and 0.610.61 in four colonies out of five. In 97.597.5 per cent of them it held at least a quarter of the front.

So a single colour holding a large share of the edge is weak evidence of selection. Two kinds of evidence are stronger. One is the same opening angle in many independent sectors of the same mutant. The other is a survival rate of single founders well above the neutral rate. Hallatschek and Nelson propose combining the opening angle and the establishment rate of sectors to measure fitness and mutation rates in colonies.

Geometry can also reshape sectors without selection. A rock of radius 1414 cells bends the sectors around it, but over 4040 paired seeds the final front carried 9.0±0.29.0\pm0.2 sectors with the rock and 9.1±0.29.1\pm0.2 without, no measurable difference.

Source: O. Hallatschek and D. R. Nelson · O. Hallatschek

Try it in the Lab

Tell a fitter founder from a lucky one

  1. Choose A faster founder opens a wedge. Predict the half-angle from the geometric rule, then compare it with the wedge you see after the front has grown across the strip.
  2. Set the advantage to zero and replay. Predict how often the gold founder survives now, and explain why it usually disappears.
  3. Choose A big sector by luck alone. Before reading the largest share, predict how far above the fair share of one in twelve it will be, then explain the result without selection.

Exercises with explained answers

  1. A mutant front is 2020 per cent faster than the wild type. What opening angle does the geometric rule predict?

    Show the answer

    cos⁡(Φ/2)=1/1.2=0.833\cos(\Phi/2)=1/1.2=0.833, so Φ/2=33.6∘\Phi/2=33.6^{\circ}, Φ=67∘\Phi=67^{\circ} and tan⁡(Φ/2)=0.663\tan(\Phi/2)=0.663.

  2. For small ss, show that the geometric rule gives tan⁡(Φ/2)≈2s\tan(\Phi/2)\approx\sqrt{2s} when v⋆/v=1+sv^{\star}/v=1+s. How good is this at s=0.1s=0.1?

    Show the answer

    tan⁡(Φ/2)=(1+s)2−1=2s+s2≈2s\tan(\Phi/2)=\sqrt{(1+s)^2-1}=\sqrt{2s+s^2}\approx\sqrt{2s}. At s=0.1s=0.1, 0.2=0.447\sqrt{0.2}=0.447 against the exact 0.4580.458, two per cent low. Small advantages give relatively large angles, which is why angles are a sensitive measure.

  3. In 400400 neutral runs a single founder survives 1313 times. At s=0.1s=0.1 it survives 5454 times. Is the difference significant?

    Show the answer

    The fractions are 0.033±0.0090.033\pm0.009 and 0.135±0.0170.135\pm0.017. The difference is 0.1020.102 with standard error 0.0092+0.0172=0.019\sqrt{0.009^2+0.017^2}=0.019, more than five standard errors. The advantage shows up only when many founders are compared.

Sources and what we recomputed

  1. O. Hallatschek and D. R. Nelson, Life at the front of an expanding population, arXiv:0810.0053v2 (2008)

    Version 2, full text. Eq. (41): a selective advantage gives sector boundaries a mean drift on top of their diffusion. Eq. (42): the drift is tan⁡(Φ/2)\tan(\Phi/2) for opening angle Φ\Phi. Eq. (46): cos⁡(Φ/2)\cos(\Phi/2) is the ratio of wild-type to mutant front speed. Eq. (49): survival probability of a beneficial sector at a linear front; its footnote 4 notes the formal similarity to the Kimura fixation probability in a well-mixed population.

  2. O. Hallatschek, P. Hersen, S. Ramanathan and D. R. Nelson, Genetic drift at expanding frontiers promotes gene segregation, arXiv:0812.2345v1; PNAS 104, 19926 (2007)

    Preprint, version 1, full text. Fig. 2B and the text around it: the mean-square wandering of 92 domain boundaries grows as L2ζL^{2\zeta} with ζ=0.65±0.05\zeta=0.65\pm0.05. The scaling argument after it gives ζ=1+(χ−1)/z\zeta=1+(\chi-1)/z, which is 2/32/3 for the Eden model.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy