< Back to Article
Epigenetic switching as a strategy for quick adaptation while attenuating biochemical noise
Fig 4
Genetic potential increases around the bistable region at high nonlinearity.
(A) The f(A) function for some solution genotypes are shown (see row and column titles) exemplifying cases where (1) both 
and 
are monostable (i.e. the associated f(A) and γ ⋅ A intersect only once; all examples with nH = 1, and {nH = 16, KD = 1}), (2) a bistable solution genotype with 
({nH = 16, KD = 50}, f(A) and γ ⋅ A intersect both in A = 20 and A = 80, with an unstable steady state between these steady states), and (3) only 
is bistable, and not a solution for the LOW environment ({nH = 16, KD = 120}). (B) Solution genotypes 
per environment were calculated (Eq 8). (C) The maximum genetic potential for different values of KD and nH (i.e. 
) is shown in the colormap. As reference, the bistable solution genotypes for each environment are delimited by the green lines: 
as dark green, and 
as light green. Additionally, we show as an example the population average KD and nH per cycle for ten simulations over 10000 generations with k = 80, nH = 1, and KD = 10 as the initial genotype (red circle, θ0), and evolutionary parameters N = 10000, u = 0.03, M = 1.1, ν = 0.10, and st = 40 (black lines), and the equivalent CONTROL simulations (i.e. without biochemical noise; gray lines). The colorbar shows the one-mutation distance corresponding to each value of M.
doi: https://doi.org/10.1371/journal.pcbi.1007364.g004