Phase transition in the genome evolution favours non-random distribution of genes on chromosomes

We have used the Monte Carlo based computer models to show that selection pressure could affect the distribution of recombination hotspots along the chromosome. Close to critical crossover rate, where genomes may switch between the Darwinian purifying selection or complementation of haplotypes, the distribution of recombination events and the force of selection exerted on genes affect the structure of chromosomes. The order of expression of gene s and their location on chromosome may decide about the extinction or survival of competing populations.

💡 Research Summary

The paper presents a computational investigation into how natural selection can shape the distribution of recombination hotspots along chromosomes, thereby influencing the non‑random placement of genes. Using Monte Carlo simulations, the authors construct two competing virtual populations, each composed of diploid individuals carrying a linear chromosome with a fixed number of loci. Every locus is assigned a selection coefficient reflecting its functional importance and a temporal expression order. The central variable is the crossover rate (the probability of a recombination event per megabase), which the authors vary across a low, critical, and high regime.

In the low‑crossover regime, the model exhibits “haplotype complementation”: the two homologous chromosomes retain complementary sets of functional alleles, allowing the population to survive despite a high load of deleterious mutations. As the crossover rate approaches a critical threshold—identified in the simulations as roughly 0.01–0.05 cM/Mb—the system undergoes a sharp phase transition. Below this threshold, complementation dominates; above it, recombination breaks up haplotypes so frequently that purifying selection becomes the primary evolutionary force, efficiently eliminating harmful alleles. This transition mirrors a classic phase‑change phenomenon, where a small change in a control parameter (crossover rate) produces a qualitative shift in population dynamics.

A key insight concerns the spatial arrangement of genes. When recombination events are biased toward chromosome ends, core genes with high selection coefficients that are centrally located experience reduced exposure to crossover‑induced disruption. Consequently, populations with central clustering of essential genes and terminal hotspots display higher survival probabilities. Conversely, placing essential genes near telomeres dramatically raises extinction risk because recombination more readily fragments these loci. The simulations also reveal that genes positioned close to one another benefit from co‑selection: their joint inheritance is reinforced, enhancing the stability of beneficial haplotypes.

The authors further explore how chromosome length, gene number, and the distribution of selection coefficients shift the critical crossover value. Larger chromosomes with many loci tend to have a lower critical threshold, making them more susceptible to the transition toward purifying selection, whereas smaller genomes require higher crossover rates to trigger the same shift. This scaling behavior suggests that the observed patterns of recombination hotspot localization in real organisms—such as the enrichment of hotspots in subtelomeric regions of mammals and plants—may be an adaptive response to balance the protection of vital genes with the generation of genetic diversity.

Overall, the study argues that recombination is not a purely stochastic process; rather, it is subject to evolutionary optimization. The non‑random distribution of recombination hotspots emerges as a consequence of selection pressures that favor chromosome architectures capable of preserving essential functions while still permitting adaptive variation. These findings have broad implications for our understanding of chromosome evolution, the formation of gene clusters, and the maintenance of complex phenotypes. Future work should integrate empirical recombination maps and environmental stressors to validate and extend the model, potentially informing breeding strategies and the management of genetic disease risk in natural and agricultural populations.