The transition to speciation in the finite genome Derrida-Higgs model: a heuristic solution

TL;DR

This paper develops an analytical and heuristic framework to understand the critical genome size needed for sympatric speciation in a finite genome Derrida-Higgs model, accounting for genetic similarity and population parameters.

Contribution

It introduces a new analytical theory for similarity distribution and a heuristic method to determine the critical genome size for speciation in finite genomes.

Findings

Derived a condition for speciation in the infinite genome limit.

Developed an analytical theory for similarity distribution without mating restrictions.

Proposed a heuristic to compute critical genome size, matching simulation results.

Abstract

The process of speciation, where an ancestral species divides in two or more new species, involves several geographic, environmental and genetic components that interact in a complex way. Understanding all these elements at once is challenging and simple models can help unveiling the role of each factor separately. The Derrida-Higgs model describes the evolution of a sexually reproducing population subjected to mutations in a well mixed population. Individuals are characterized by a string with entries $\pm1$ representing a haploid genome with biallelic genes. If mating is restricted by genetic similarity, so that only individuals that are sufficiently similar can mate, sympatric speciation, i.e. the emergence of species without geographic isolation, can occur. Only four parameters rule the dynamics: population size $N$, mutation rate $\mu$, minimum similarity for mating $q_{min}$ and genome size $B$. In the limit $B\rightarrow\infty$, speciation occurs if the simple condition $q_{min}>(1+4\mu N)^{-1}$ is satisfied. However, this condition fails for finite genomes, and speciation does not occur if the genome size is too small. This indicates the existence of a critical genome size for speciation. In this work, we develop an analytical theory of the distribution of similarities between individuals, a quantity that defines how tight or spread out is the genetic content of the population. This theory is carried out in the absence of mating restrictions, where evolution equations for the mean and variance of the similarity distribution can be derived. We then propose a heuristic description of the speciation transition which allows us to numerically calculate the critical genome size for speciation as a function of the other model parameters. The result is in good agreement with the simulations of the model and may guide further investigations on theoretical conditions for species formation.