Skip to main content
Figure 1 | BMC Evolutionary Biology

Figure 1

From: Error-threshold exists in fitness landscapes with lethal mutants

Figure 1

Relationship between the error-threshold and various parameters. Lines are calculated from Eq. (4) of the model without back mutation. The notation and parameters are as follows (unless stated otherwise): The broken line is for κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ 1 = 0.05; the solid line is for κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ 1 = 0.5; the dotted line is for κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ 1 = 1.0; f0/f1 = 0.99-1. Points (circles and plus signs) are calculated through the comparison of the first and second greatest eigen values of the model with back mutation [see "How high must lethality be to introduce an error-threshold" and Fig. 4(A) for details]. κ = 0.05 (the same as in Fig. 3) and f d = 0.99d, unless stated otherwise. Note that Eq. (4) is the necessary condition of the error-threshold without back mutation, and thus, it overestimates 1 - qmin compared to the calculation through the comparison between the eigen values. (A) Relationship between the error-threshold and sequence length. The graph shows that the error-threshold becomes smaller as ν increases. Note the faster decrease shown by the circles compared to that by the dashed line. Although there is a scaling between ν and κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ for the lines, as seen from κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ (ν - 1) term in Eq. (4), this is merely due to the fact that only G MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFge=raaa@382C@ 0 and G MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFge=raaa@382C@ 1 are considered in Eq. (4). (B) Relationship between the maximum tolerable genomic mutation rate per replication [i.e., (1 - qmin)ν] and ν. The graph shows that (1 - qmin)ν increases as ν increases if f0/f1 is kept constant (the thick lines and the circles), but also shows that this does not hold if f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ is kept constant (the thin lines and the plus signs). f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ represents the relative reproductive advantage of the fittest genotype compared to the perfectly randomized genotype which is not lethal. The thin lines are obtained by setting f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ to a constant value such that it coincides with the corresponding thick line at ν = 20; i.e., f 0 / f ν 2 = 0.99 − 20 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaGaeGimaadabeaakiabc+caViabdAgaMnaaBaaaleaadaWcaaqaaGGaciab=17aUbqaaiabikdaYaaaaeqaaOGaeyypa0JaeGimaaJaeiOla4IaeGyoaKJaeGyoaKZaaWbaaSqabeaacqGHsisldaWcaaqaaiabikdaYiabicdaWaqaaiabikdaYaaaaaaaaa@3D2B@ , and f0/f1 is calculated as 0.99 − 20 ν MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqaIWaamcqGGUaGlcqaI5aqocqaI5aqodaahaaWcbeqaaiabgkHiTmaalaaabaGaeGOmaiJaeGimaadabaacciGae8xVd4gaaaaaaaa@3547@ by assuming a multiplicative fitness landscape. The circles are for a constant f0/f1. The plus signs are obtained by setting f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ to a constant value such that it coincides with the circles at ν = 20 [ f 0 / f ν 2 = 0.99 − 20 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaGaeGimaadabeaakiabc+caViabdAgaMnaaBaaaleaadaWcaaqaaGGaciab=17aUbqaaiabikdaYaaaaeqaaOGaeyypa0JaeGimaaJaeiOla4IaeGyoaKJaeGyoaKZaaWbaaSqabeaacqGHsisldaWcaaqaaiabikdaYiabicdaWaqaaiabikdaYaaaaaaaaa@3D2B@ ; thus, f d = 0.99(20/ν)d]. (The lines for κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ 1 = 0.5 are not shown for visibility.) (C) Relationship between (1 - qmin)ν and the intensity of selection f0/f1. The graph shows that for a given value of f0/f1, a longer sequence can tolerate a greater (1 - qmin)ν. In other words, for a given value of (1 - q)ν, a longer sequence can be maintained (i.e., 1 - q is below the error-threshold) for a smaller value of f0/f1. The thick lines and circles are for ν = 80. The thin lines and plus signs are for ν = 40. For all plots, f d = (f0/f1)-d. (The lines for κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ 1 = 0.5 are not shown for visibility) (D) Relationship between (1 - qmin)ν and f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ . The graph shows that for a given value of f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ , (1 - qmin)ν does not differ between different sequence length (this is comparable to the results of [3, 14], and note a parallelism between f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ and the superiority parameter σ defined there). For κ = 0.05 (dashed lines), there is still a noticeable difference, but this is due to the fact that only G MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFge=raaa@382C@ 0 and G MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFge=raaa@382C@ 1 are considered to obtain them although (1 - qmin)ν is large. [compare this result with the circles and plus signs, which are obtained by considering all G MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFge=raaa@382C@ d (0 ≤ d ≤ ν) for almost the same value of κ]. For all plots, f d = (f0/ f ν 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzdaWgaaWcbaWaaSaaaeaaiiGacqWF9oGBaeaacqaIYaGmaaaabeaaaaa@30EE@ )(-2/ν)d. The notation and the rest of parameters are the same as those in (C). (the lines for ν = 40 and ν = 80 are almost on top of each other for κ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacaaaa@2E74@ 1 = 0.5 and 1.0.)

Back to article page