Does Sex Induce a Phase Transition?

We discovered a dynamic phase transition induced by sexual reproduction. The dynamics is a pure Darwinian rule with both fundamental ingredients to drive evolution: 1) random mutations and crossings which act in the sense of increasing the entropy (or diversity); and 2) selection which acts in the opposite sense by limiting the entropy explosion. Selection wins this competition if mutations performed at birth are few enough. By slowly increasing the average number m of mutations, however, the population suddenly undergoes a mutational degradation precisely at a transition point mc. Above this point, the “bad” alleles spread over the genetic pool of the population, overcoming the selection pressure. Individuals become selectively alike, and evolution stops. Only below this point, m < mc, evolutionary life is possible. The finite-size-scaling behaviour of this transition is exhibited for large enough “chromosome” lengths L. One important and surprising observation is the L-independence of the transition curves, for large L. They are also independent on the population size. Another is that mc is near unity, i.e. life cannot be stable with much more than one mutation per diploid genome, independent of the chromosome length, in agreement with reality. One possible consequence is that an eventual evolutionary jump towards larger L enabling the storage of more genetic information would demand an improved DNA copying machinery in order to keep the same total number of mutations per offspring.

💡 Research Summary

The paper investigates a dynamical phase transition that emerges in sexually reproducing populations when the balance between random mutations (which increase genetic diversity) and natural selection (which reduces diversity by eliminating deleterious alleles) is altered. The authors construct a minimal computational model of diploid organisms, each carrying two chromosomes of length L bits (0/1). In each generation, parental genomes undergo random crossover and recombination, after which an average of m mutations are introduced uniformly across the genome. Mutations are modeled as bit flips, thereby increasing the system’s entropy. Fitness is defined as a decreasing function of the number of “bad” alleles present; individuals with higher fitness have a larger probability of contributing offspring to the next generation, embodying the selective pressure that opposes entropy growth.

The central control parameter is the average number of mutations per offspring, m. By slowly raising m from low to high values, the authors monitor population‑level observables such as mean fitness, genetic entropy, and the distribution of allelic states. They discover a sharp threshold mc at which the system undergoes a qualitative change: for m < mc, selection successfully suppresses the spread of deleterious alleles, the population retains a broad distribution of genotypes, and evolutionary exploration continues. When m exceeds mc, the influx of new mutations overwhelms selection; “bad” alleles proliferate, the population collapses into a nearly homogeneous genetic state, mean fitness drops dramatically, and further adaptive evolution essentially ceases. This abrupt shift is identified as a phase transition in the statistical‑physics sense.

A striking result is that the transition curves become independent of chromosome length L once L is sufficiently large. In other words, the critical mutation load mc does not shift with increasing genomic size; the relevant quantity is the total number of mutations per diploid genome, not the mutation rate per base pair. Moreover, the authors demonstrate that population size N has negligible impact on the location or sharpness of the transition, indicating that the phenomenon is not a finite‑size artifact but a genuine thermodynamic‑like transition.

To characterize the critical behavior, the authors perform finite‑size scaling analyses. Near mc, fluctuations in fitness and entropy grow, and the system exhibits scaling exponents reminiscent of second‑order phase transitions. The transition is thus not merely a smooth crossover but a true singularity in the limit of large L (and, implicitly, large N).

Biologically, the findings provide a quantitative explanation for the empirical observation that most multicellular organisms experience roughly one new mutation per diploid genome per generation. The model predicts that life cannot remain evolutionarily viable when the average mutation load exceeds this value, regardless of how much genetic information is stored. Consequently, any evolutionary expansion of genome size must be accompanied by improvements in DNA replication fidelity or enhanced DNA‑repair mechanisms to keep m near or below mc. This insight links the evolution of genome complexity directly to the evolution of molecular machinery that safeguards genetic integrity.

In summary, the paper presents a clear, mathematically grounded demonstration that sexual reproduction, combined with a simple Darwinian rule set, can generate a mutation‑driven phase transition. Below the critical mutation load, selection maintains diversity and enables ongoing evolution; above it, the system locks into a low‑fitness, genetically uniform state, effectively halting evolution. The work bridges concepts from statistical physics and evolutionary biology, offering a fresh perspective on why mutation rates are tightly constrained in nature and highlighting the co‑evolutionary pressure on replication fidelity as genomes expand.