key: cord-0837147-c4dopubu authors: Fujihara, Ichiro; Furusawa, Mitsuru title: A differential equation, deduced from a DNA-type genetic algorithm with the lagging-strand-biased mutagenesis date: 2022-03-26 journal: Heliyon DOI: 10.1016/j.heliyon.2022.e09155 sha: d036e7b856b1ddd0e62cfb98741e13e98820a384 doc_id: 837147 cord_uid: c4dopubu Recent evidence indicated that the significant fidelity difference existed between the leading and lagging-strand. Using a 2-D DNA-type genetic algorithm (GA), we have shown that the lagging-stand-biased mutagenesis (disparity mutagenesis) has a significant advantage in the promotion of evolution, compared with the traditional parity mutagenesis model. Our aim of the present study was to deduce a differential equation which well reflected the result of simulations. The analytical solution of the differential equation obtained was in good agreement with the results of the simulation experiment. Comparing the results of the disparity mutagenesis with those of the parity one, the characteristics of the disparity mutagenesis were discussed in terms of relative mutation rates between the lagging and leading strands. Conditions of the extinction of a species were also discussed. When any two species are arbitrarily chosen, it is believed that they have evolved from a common ancestral species. The ancestral species, however, had been extinct ages ago in most cases. The memory of the mutations might be vanished away during long evolutionary time. Accordingly, it would be more difficult to ascertain the key mutations that caused this divergence. Thirty years ago, one of the present authors, MF, thought of the idea that experimental acceleration of evolution using a living organism might reveal the process of species differentiation. To realize the acceleration of evolution, we proposed a specific mutant, "disparity mutator", in which the lagging-strand-biased mutagenesis was achieved [1, 2] . Using disparity mutators, we could experimentally accelerate evolution in bacteria [3, 4] , yeasts [5] , malaria parasites [6] and mice [7] . By the way, as a plausible way to simulate the evolutionary process and to approach general optimization problems, J. Holland proposed a genetic algorithm (GA) in 1975 [8] . Subsequently, GA was mainly used for the solution of complex optimization problems. Unexpectedly, however, the contribution of GA to evolutionary genetics was very small. For instance, to our knowledge there has been no review or report on the relationship between GA and biological evolution in the last five years. E-mail address: ichiro@las.osaka-sandai.ac.jp (I. Fujihara). Living things evolve as a population (species). The population is formed by the proliferation of each individual member, which is definitely controlled by DNA. However, studies of conventional population genetics have been carried out in the past without paying attention to the mechanism of DNA replication. Instead, average mutation rates per DNA (or cell) replication, and the cell death by deleterious mutations have mostly been considered. As the result, a specific value, "errorthreshold", was worked out. Namely, it was stated that a cell dies when mutation number/cell/replication >1 [9, 10] . On the contrary, we firstly developed a DNA-type genetic algorithm (GA), which mimicked the molecular structure of DNA and its replication machinery. We simulated the process of evolution using this DNA-type GA [11] . The distinctive characteristic of our study is that the fidelity difference between the leading and lagging strand was considered, based on the evidence of the fidelity difference in living things [12, 13, 14] . Our GA consisted of two complemental strand algorithms as if DNA was replicated. From the results of simulations, here we obtained a differential equation for our disparity mutagenesis model. In the present study, the lagging-stand-biased mutagenesis was applied as a matter of practical convenience. The effects of the disparity mutagenesis can be summarized as follows. 1) When the fidelity of leading DNA strand is enough high, the genotype of a parent DNA might be transmitted to either daughter DNAs as it is. Therefore, the ancestral genotype is forever transmitted in terms of a replicon, even when the mutation rates of the lagging strand are high enough. 2) As a result, the error-threshold of the DNA population rises beyond 1. The GA used in the present study was basically the same one as previously published, which was the 2-D DNA-type GA [15] . This digital DNA takes into consideration gene-to-gene interactions and semiconservative replication using a leading and lagging strand according to the DNA of living things. The aim of the present study was to deduce a differential equation from the results of simulations using this GA. Two differential equations were deduced from the conventional parity mutagenesis and from the disparity mutagenesis, respectively. The solutions of these equations agree well with the results of the simulations. We examine the mathematical characteristics of these equations and discuss the biological meanings of the disparity mutagenesis model. Here we summarize the characteristics of the disparity systems derived from the results of the previous paper. 1. Most importantly, the disparity mutagenesis significantly increased the error-threshold. This evidence would clearly explain the increased adaptability shown by the living disparity mutators. 2. Furthermore, the increase of the error-threshold in the disparity mutagenesis would introduce an evolutionary jump when a drastic environmental change occurs. The parity mutagenesis also could produce an evolutionary jump, but only within a very limited range of mutation rates. 3. In a long stable state of environments, evolution eventually reached a plateau in the disparity mutagenesis due to the high fidelity of the leading strand, even when the average mutation rates were sufficiently high. It might provide us the logic why living fossils exist. 4 . In short, the leading strand with a high fidelity guarantees the stability of species. On the other hand, the lagging strand with a low fidelity makes it possible to explore a much broader fitness landscape. After all, although it seems to be contradictory, the disparity mutagenesis realizes the expansion of diversity and the increase of stability of species at the same time. In the disparity system, the lagging strand introduces more mutations to promote evolution, but it is also more likely to cause death of a cell. Therefore, the leading strand introduces as few mutations as possible to preserve the species. The rate of evolution and the stability of the species seem to depend on the balance of the number of mutations entering two strands. On the other hand, in the parity system mutations are equally distributed in both strands. For the simulation of evolution processes, the DNA-type 2-D GA was used, the genome of which consisted of double-stranded DNA, and semi-conservatively replicated using the leading and lagging strand. The genome consisted of fifty functional genes (FC). A single FC gene had three unique fitness landscapes. The two hundred non-functional junk genes (JK) were distributed evenly in number between FC. Two neighbor genes including JK are interacting, that means all genes are interacting. This model GNA is expressed as (FC=50, JK=200). The details of the contents of the GA and conditions for simulations were described in the-previous study [15] . In our previous report [15] we regarded DNA as a sequence of connected genes on a two-dimensional information lattice and assumed that neighboring genes interacted with each other. At the initial state of the simulation, the population consists of 100 cells (or genomes) and the coordinates of the (fitness score) of each gene are set to "0" on the horizontal axis in every individual cell. This Figure legend: The vertical axis represents the order of the genes, and the horizontal axis represents the scores for each gene. If the distance between two vertically adjacent genes is more than "1" in the coordinates, the gene interaction is broken and the cell dies (e.g., the enlarged image on the left, where the two genes in the middle are connected by a dotted line indicating the break). The dark squares show the positions at which genes get points. means that from the start of the simulation, the early homogeneous cell population had to survive in a completely new environment. Fig. 1 shows an example of how genes are arranged on a twodimensional lattice and how the cell's alive/dead was determined. In this model, genetic variation is indicated by a shift in the coordinates of the gene. It is assumed that the survival of a cell requires the interaction of all neighboring genes, whether JK or FC genes. If the distance between adjacent genes is less than or equal one, they are regarded interacting. When the gene interaction was interrupted at any position by mutations, the chromosome immediately died. Under these assumptions, we carried out Monte Carlo simulations to investigate how the DNA sequence changes during cell division by repeatedly introducing mutations in the leading and lagging strands of DNA. The main purpose of the calculations was to clarify the difference of between the parity and disparity mutagenesis. For the disparity systems, even if the fidelity of the lagging chain is low, the robustness is sufficient if the fidelity of the leading chain is high. Acceleration or deceleration of evolution can only be brought about by environmental changes. Fig. 1 illustrates the survival conditions when a gene undergoes mutation on the information lattice and how the gene gets scores. The vertical axis represents the order of the genes, and the horizontal axis represents the scores for each gene. In the previous report [15] , the simulation was performed under the following conditions. 1) The mutation rate of the lagging strand was an integer and was normally distributed around (mutation rate of lagging strand). 2) The mutation rate of the leading strand was real and 0 ≤ < 1 ( ; mutation rate of leading strand). See the previous report [15] for more details on the model. In the previous study, it was shown that the effects of and on cell growth in the early stages of the simulation were rapid and rather ) with different mutation rates in the disparity system. Figure legend: The genome consisted of fifty functional genes (FC) and two hundred junk genes (JK). Fitness scores are shown when =0 and =2-7. Y axis: fitness score. X axis: generation. The initial population was 100. When the population size became more than 2,000, truncation selection was applied to adjust the size to about 2,000. values: black line ( = 2), red (3), orange (4), green (5), blue (6), purple (7). complicated. Therefore, in the present study, we focused on the early stage of simulation and changed the simulation conditions as follows: 1') The mutation rate of the lagging strand is simply . 3') There is no bias in favor of or against mutation. Left and right movements were made with equal probability. Fig. 4(a) , the final and the rate of increase in at =0 were smaller than those at =0.025. In general, the final placement of genes and the speed of rise depends on random numbers employed. Simulations were performed using several different random number sequences under the same conditions of =4 and =0. As a result, it was found that fluctuates about ±5 with different random number sequences. This may be the reason the final and the rate of rise in at =0 is smaller than those at =0.025. ,ext, were observed at =0.191 for =4 and =0.060 for =5. These simulation results suggest the high fidelity of the leading strand is an essential condition for a fast rise and high score of . Fig . 5 shows the dependence of for the parity system, when (= ). As shown in this figure, the rate of increases of and the final height of became higher as increased from 0.2. It was in a wide range of (= ) = 0.3-0.7 that the was able to rise faster and reach a higher point. In a previous paper [15] , we evaluated the stability of cells as a species with Shannon-Wiener ′ . The stability of the cells as a species was evaluated in the same way as before. The cells were stable as a species only in a small region where the mutation rate was about (= ) = 0.7. The speed of rise and the highest value of began to decrease rapidly when (= ) was beyond 0.8. To quantitatively determine the magnitude of and its generational change, the following two factors appear to be particularly important. They are 1) The change in the cell number before and after cell division, Δ , and 2) the number of mutations introduced into the surviving cells after cell division, . To calculate these quantities in parity and disparity systems, it is needed to evaluate the probability of survival of cells from the leading and lagging strands after mutation and the expected values of mutations entering the surviving cells. Let be the probability of a cell's death when a single mutation is introduced. The probability that a cell survives after the insertion of a single mutation is (1-). The leading strand of both parity and disparity systems receives a (<1.0) mutation per cell division. Let and be the probability that a progeny born from the leading strand will be dead or alive, respectively. They are also expected values of dead and alive cells from a single cell. They are expressed as follows: In the disparity system, a single cell division introduces mutations in the lagging strand and mutations in the leading strand. In the case of leading strands, the probability that a progeny born from the leading strand will be alive or dead, is the same with those for the parity system, as expressed by the equations (1) and (2). In a single cell division, mutations are inserted in the lagging strand. Let and be the probability that a progeny born from the lagging strand will be alive or dead, respectively. They are expressed as follows: In the parity system, a cell's DNA can be considered to replicate using leading and lagging strands with the same probability of mutation (= ). The lagging strands can, then, be regarded as leading strands in parity system. The increase in the number of cells due to cell division, Δ , is obtained by subtracting 1 from the expected value of the number of living cells from the two leading strands. After cell division, if two cells are born alive from two leading strands (parity system), then Δ =1. Conversely, if no cells remain, then Δ = -1. For the disparity system, the expected number of surviving cells after division is the sum of and which are represented by equations (2) and (3), respectively. The increase in the number of cells due to cell division, Δ , is obtained by subtracting 1 from the expected value of living cells from the leading and lagging strands. It is represented as follows: The equation (6) indicates that as increases, the difference in the number of cells before and after cell division, Δ , decreases almost exponentially. However, it also indicates that the number of cells is never reduced if = 0. In the disparity system, the fidelity of the leading strand is particularly important to maintain the cell number. In the leading strand, the probability of mutation is and the probability of survival is (1-) . Thus, the average number of mutations introduced into the surviving leading strand, , , is given by: Then, in the parity system, the expected number of mutations inserted into the two leading strands of the next generation, , is written as follows: In the disparity system, mutations are introduced into the lagging strand and the surviving probability is (1 − ) . Thus, the number of mutations remaining in the surviving lagging strand, , , is given by: Then, in the disparity system, the expected number of mutations inserted into the leading and lagging strands of the next generation can be written by using equations (7) and (9) as follows: As with Δ , when selection pressure is applied to the number of cells according to the value, the larger the , the more favorable the increase in . As described in Section 4, the larger Δ , the more daughter cells there are after cell division and the larger , the more mutations are inserted in the surviving daughter cells. The product of Δ and is an excellent indicator of the rate of evolution of a species and the final value that can be reached. Then, we assumed that the evolutional activity coefficients, (Δ , ), which expresses the speed of rise of , can be expressed as follows: Δ and for the parity systems are expressed by equations (5) and (8), respectively, and those for the disparity systems are expressed by equation (6) and (10) . In the present study, , is a constant which is related to the pattern of the fitness landscape shown in Fig. 1 . However, the dependence of on the landscape pattern has not tested at present. (Δ , ) for the parity system can be regarded as a function of (= ). While (Δ , ) for the disparity system can be regarded as a function of and . From here on, instead of (Δ , ), we use the notation pari ( ) for parity systems and disp ( , ) for disparity systems. The evolutionary activity coefficient for parity system, pari ( ), is obtained using equations (5) and (8), and is expressed as follows: And the evolutionary activity coefficient for disparity system, disp ( , ), is obtained using equations (6) and (10), and is expressed as follows: In a disparity system, the smaller both and are, the larger the value of disp ( , ) will be. For example, the sequence of the magnitudes of disp ( ,0) is disp (2,0)> disp (3,0)> disp (4,0)> disp (5,0)> and so on. Fig. 2 shows that the lower the , the faster the rise in a disparity system. Then, we assumed the following relations hold between ∕ and pari ( ) for the parity systems and ∕ and disp ( , ) for the disparity systems. ∕ ∝ ( ( )) (parity system) ∕ ∝ ( ( , )) (disparity system) (15) where, is generation and the functions ( pari ( )) and ( disp ( , )) are a monotonically increasing function of pari ( ) and disp ( , ). Equations (14) and (15) imply that the larger pari ( ) and disp ( , ) are, the faster the rise of rises. is the sum of the fitness scores obtained when 50 functional genes in DNA come to a scoring point on the information lattice as shown in Fig. 1 . As can be seen in Fig. 1 , the adjacent genes must be linked to each other, so cannot be larger than a certain . Then, in addition to equation (14) or equation (15), we must also consider that the final value of will be less than . As shown in Figs. 2-5, when evolutionary stage was early and was low, seemed to move rapidly to higher position. The rate of increase in was assumed to be proportional to the difference between the highest value and the current . This situation can be expressed by a differential equation as follows: where, is the height of the plateau that the cell population can reach and is dependent on for parity systems and dependent on and for disparity systems. However, as shown in Fig. 2 , when = 2 to 5 in the disparity system, the final plateau height was almost the same, if = 0. Even at =6 and =7, the still seems to be rising after 150,000 generations. Then, of the disparity system is assumed to be when =0, and is considered to be zero in the extinction limit ,ext. ( , ) should be closely related to the evolutionary activity coefficient as well as the rate of increase in . We assume here that the dependence of ( , ) on and is as follows: Similarly, we assume that the dependence of ( ) on for parity system is assumed as follows: where disp ( , ) and pari ( ) are the normalized evolutional activity coefficients expressed as follows: and Fig. 6 shows the dependences of pari ( ) on for the parity and the dependence of disp ( , ) on and for the disparity systems ( =2 to 5). In the disparity systems, ( , ) is 1 at =0 Fig. 6 . The dependence of pari ( ) on for the parity and disp ( , ) on and for the disparity systems. Probabilities of a cell dying due to a single mutation, , were 0.45 for both parity and disparity systems. Figure legend: The relation between and pari ( ) or disp ( , ) black line (parity, = e), red (disparity, =2), green (disparity, =3), light blue (disparity, =4), and purple (disparity, =5). and decreases monotonously with increase in . The negative value of disp ( , ) indicates is over the extinction limit, ,ext. On the other hand, in the parity system, pari ( ) draws like a parabola with respect to , and is highest when is between 0.5 and 0.6. These data suggest that the plateau is highest when is zero for the disparity system, and when is between 0.5 to 0.6 for the parity system. If disp ( , ) or pari ( ) is less than zero, the cell population becomes extinct. Considering the equations (15), (16) and (17), the generation dependence of of the disparity system can be expressed by the following differential equation. By transforming this equation, we obtained the following separable differential equation. Solving the equation (24), we obtain the following equation. where ( , ) is given by equation (19) . Similarly, the generation dependence of of the parity system can be expressed using equations (14) , (16) and (18), as follows: Next, we consider the specific expressions of ( ( , )) for the disparity systems and ( ( )) for the parity systems. As mentioned earlier ( ( , ))and ( ( )) are monotonically increasing functions of disp ( , ) and pari ( ), respectively. First, for simplicity, we assumed the relationships ( ( )) = ( ) and ( ( , )) = ( , ). Then, for the parity system, using equations (12) and (20), equation (26) becomes as follows: For the disparity system, using equations (13) and (19), equation (25) becomes as follows: Fig. 7 . The generation dependence of for when =0, calculated for the disparity system calculated from equation (28). Figure legend: values: black line ( =2), red (3), orange (4), green (5), blue (6), purple (7). As mentioned earlier, the initial value of in the simulations was 35 and increased to about 140 in the highest case. In the present calculations using equations (27) and (28), the initial values were set to 0 and the highest plateau values were set to 100. To draw following figures from equation (28) and equation (27), we needed the probability of the cell dying due to a single insertion of mutation, . It was set at 0.45 which was evaluated from the relation between ln( ,ext) and which will be explained in the following Section 7. The integral constant was determined based on the simulation results of the disparity system ( =2, =0) to satisfy that the calculated value of is zero at 0th generation and reaches the plateau value of 100 ( ) after 20,000 generations. In this paper, the same integration constant was used for all the calculations for both the disparity and parity systems. Fig. 7 shows the generation dependence of the for the disparity systems at =0, calculated using equation (28). Comparing with Fig. 2 , the dependence of in the disparity system with =0 is well represented by equation (28). When =0, the smaller the value of , the faster the rises. This is consistent with the results presented in Fig. 2 . As we can estimate from equations (6) and (10), the number of cells that survives and the number of mutations introduced are higher, when =2. Both quantities became smaller as increased. Figs. 8 and 9 show dependences of of disparity systems for =2 to =5 calculated from equation (28). As shown in these figures, the dependence of in the disparity system also seems to be well represented by equation (28). In the early stages, rises rapidly like a logarithmic function: the smaller the , the faster the rise in . In the disparity system, of about 2 would be sufficient for performing rapid evolution under the present conditions. As equation (6) indicates, Δ decreases rapidly as increases. If the number of cells does not increase sufficiently, the truncation selection pressure does not work effectively and consequently the increase in is reduced. When increases at any , the cell growth decreases and eventually reaches zero. The negative Δ in equation (6) indicate the cells become extinct. The graphs in Figs. 8 and 9 illustrate well the characteristics of the graphs in Figs. 3 and 4 , respectively. Fig. 10 shows the generation dependence of for the parity system calculated using equation (27). As stated above, in the early stages of the disparity system, both the rate of increase of and the magnitude Fig. 8(a) . The generation dependence of for =2 of the disparity system at different values calculated from equation (28). of decreased as increased from zero at any . However, in the parity system, as (= ) increases, the rate of rise in and the height of the plateau increases. Though the curves for =0.5 and =0.6 are too close to distinguish, both the rate of increase and the magnitude of are maximum around =0.5-0.6. Both the rate of increase and the magnitude of begin to decrease beyond >0.6. This is because in the parity system, when is small, Δ is large but is small. Conversely, when is large, Δ is small but is large. Therefore, the behavior of depends on the balance between Δ and . The fast increase rate and large value of were observed as well as in the disparity system. However, they are limited in a narrow region around = 0.5 to 0.6. This behavior is well represented by the dependence of pari ( ) shown in Fig. 6 . Comparing Fig. 10 with Fig. 5 , we can see that the equation (27) is able to represent the dependence of of the parity systems as well as the disparity systems. However, for =0.2, the calculated values seem to be slightly lower than the simulated values. Here we also consider the relationship between the extinction of population and the mutation parameters; and . The probability Table 1 are plotted. ,ext: filled black circles (FC=50, JK=50), filled red triangles (FC=50, JK=100), filled green squares (FC=50, JK=200) and filled purple diamonds (FC=50, JK=450). The linear approximation line for each system is represented by the same colors. that a cell will die when a single mutation is introduced is treated as a constant. The effect of differences in the fitness landscape on the FS values has not yet been clarified. The distribution of the fitness points in each FC gene and the number of JK genes between FC genes may affect the speed of rise and the height of the plateau through . In this paper, we confirmed the effect of the number of JK genes on the FS values. The ,ext (the extinction limit) in the disparity systems, which consisted of 50 FC genes and different number of JK genes, is shown in Table 1 . The logarithm of ,ext against is plotted in Fig. 11 . By the way, we have already dealt with this relationship in equation (6) . From equation (6), if Δ <0, the number of cells in the disparity system goes down with increasing generations, so the ,ext becomes Taking the logarithm of both sides of equation (29) gives the following equation: From equation (30), we see that the slope of the graph is ln(1-). The probability of surviving by the insertion of a single mutation into a cell, (1-), was estimated from the slope of the line with ranging from 2 to 5. In the present system (FC=50, JK=200), the value of (1-) was about 0.55. Therefore, was set to 0.45 in this study. As shown in Fig. 11 , the dependence of ln( ,ext) is similar to each other for all ratios of FC and JK. In addition, as shown at the bottom of Table 1 , becomes slightly larger when the number of JK genes increases, but the difference is small. This suggests that the number of JK genes between FC genes does not have much influence on the results. The present treatment shows that , the number of mutations introduced into a cell, and Δ , the changes in cell number due to cell division, play an important role in the evolution of the present genetic algorithm. The probability that a cell will die from the introduction of a single mutation, , plays an important role in this process. If =0, cells from the leading strands never die. In such case of disparity system, cells can just barely maintain their numbers in the next generation, even if many mutations are introduced at one time into the lagging strand. However, in a parity system, the rate of increase and magnitude in can only be expected within a narrow range of (= ). We conclude that the disparity system is more favorable to evolution than the parity system. Concerning the general implications of asymmetric replication of DNA, see the review [13, 14] . We proposed here two differential equations for evolution (the parity model and disparity model), based on the molecular mechanism of semiconservative and asymmetrical DNA replication. The biological implications of equation (27) for the parity model and equation (28) for the disparity model were summarized as follows: 1) Evolution is a function of the number of generations: Judging from the present results, it might appear that evolution could be the function of time. This is because a synchronized growth system is used. However, one thing is for sure -evolution is the function of the number of generations. 2) Increase of error threshold: Equation (27) (conventional parity model) led to catastrophe when the error rate of the leading and the lagging strands ( = ) was beyond a bit 1. However, equation (28) (disparity model) did not lead to catastrophe even when =15, if was significantly small; meaning that the evolution ability of this population was still maintained. 3) Evolutionary jump: At the start of evolution experiments in the present study, a genetically homogeneous population was exposed to an absolutely novel environment. The fast rise of (fitness score) observed in the parity and disparity models may mean that the environmental shock induces an evolutionary jump. We showed in the previous study that drastic environmental changes artificially introduced in a stable state of evolution induced an evolutionary jump without increasing mutation rates [15] . This jump was followed by the explosive evolution that was reminiscent of the Cambrian explosion. The quick response to the transmissible cancer observed in Tasmanian devil might be explained by assuming a disparity mutagenesis [16] . 4) Stable evolution: To maintain an appropriate and to keep low genetic diversity would be indispensable conditions for a stable species. These conditions are satisfied only at very narrow mutation rates ( = =0.5-0.7) in the parity model. In the disparity model, the error threshold is considerably increased. This is because once acquired must be guaranteed forever by the very low values which are close to zero in the broader range of mutation rates. This evidence may mean that it is easy to form an evolutionary plateau. In other words, the existence of living fossil or dead-end of evolution would enable us to consider this at least theoretically. 5) Pandemic of COVID-19 and disparity mutagenesis: Genome of SARS-CoV-2 is not double-stranded DNA but single-stranded RNA. The all-RNA molecules appearing in an infected cell are + or -stranded RNA. These RNA molecules act as a template when replicating. In other words, a single-stranded RNA is "conservatively" replicated. Therefore, the template molecule must be error-less, but the newly-synthesizing strand is error-prone due to the low fidelity of the viral RNA polymerase. Thus, the mode of the RNA replication can be regarded as a typical disparity mutagenesis. The equation (28), =0, may provide a simple explanation of the strong infectivity of single-stranded RNA viruses including SARS-CoV-2. Ichiro Fujihara: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Wrote the paper. Mitsuru Furusawa: Analyzed and interpreted the data; Wrote the paper. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Data will be made available on request. The authors declare no conflict of interest. No additional information is available for this paper. Promotion of evolution: disparity in the frequency of strandspecific misreading between the lagging and leading strands enhances disproportionate accumulation of mutations Asymmetrical DNA replication promotes evolution: disparity theory of evolution Preferential replication-dependent mutagenesis in the lagging DNA strand in Escherichia coli A conspicuous adaptability to antibiotics in the Escherichia coli mutator strain dnaQ49 Isolation of thermos-tolerant mutants by using proofreading-deficient DNA polymerase as an effective mutator in Saccharomyces cerevisiae Generation of rodent malaria parasites with a high mutation rate by destructing proofreading activity of DNA polymerase Germline mutation rates and the long-term phenotypic effects of mutation accumulation in wild-type laboratory mice and mutator mice Adaptation in Natural and Artificial Systems The molecular quasispecies Evolutionary Dynamics: Exploring the Equations of Life, President and Fellows of A neo-Darwinian algorithm: asymmetrical mutations due to semiconservative DNA-type replication promote evolution The inherent asymmetry of DNA replication Implications of fidelity difference between the leading and the lagging strand of DNA for the acceleration of evolution The nature of mutations induced by replicationtranscription collisions Disparity mutagenesis model possesses the ability to realize both stable and rapid evolution in response to changing environments without altering mutation rates Rapid evolutionary response to a transmissible cancer in Tasmanian devils The authors thank Dr. F. Rueker for his critical reading of the manuscript, and Dr. N. Gouda for his helpful advice and comments.