Biophysics and population size constrains speciation in an evolutionary model of developmental system drift

Developmental system drift is a likely mechanism for the origin of hybrid incompatibilities between closely related species. We examine here the detailed mechanistic basis of hybrid incompatibilities for a genotype-phenotype map for developmental system drift under stabilising selection, where the organismal phenotype is conserved, but the underlying molecular phenotypes and genotype can drift. This leads to number of emergent phenomenon not obtainable by modelling genotype or phenotype alone. Our results show that: 1) speciation is more rapid at smaller population sizes with a characteristic, Orr-like, power law, but at large population sizes slow, characterised by a sub-diffusive growth law; 2) the molecular phenotypes under weakest selection contribute to the earliest incompatibilities; and 3) pairwise incompatibilities dominate over higher order, contrary to previous predictions that the latter should dominate. Our results indicate that biophysics and population size provide a much stronger constraint to speciation than suggested by previous models.