Computer simulations on the sympatric speciation modes for the Midas cichlid species complex

Cichlid fishes are one of the best model system for the study of evolution of the species. Inspired by them, in this paper we simulated the splitting of a single species into two separate ones via random mutations, with both populations living together in sympatry, sharing the same habitat. We study the ecological, mating and genetic conditions needed to reproduce the polychromatism and polymorphism of three species of the Midas Cichlid species complex. Our results show two scenarios for the A. Citrinellus speciation process, one with and the other without disruptive natural selection. In the first scenario, the ecological and genetic conditions are sufficient to create two new species, while in the second the mating and genetic conditions must be synchronized in order to control the velocity of genetic drift.

💡 Research Summary

The paper presents a computational investigation of sympatric speciation in the Midas cichlid species complex, focusing on how a single ancestral population can split into two distinct species while sharing the same habitat. The authors first synthesize biological knowledge about the complex’s polychromatism (multiple colour morphs) and polymorphism (morphological variation such as jaw shape), then translate these traits into a set of quantitative parameters for an individual‑based simulation model. The model consists of three interacting modules: (1) an ecological selection module that maps resource distribution and habitat heterogeneity onto a fitness landscape, explicitly incorporating disruptive natural selection that creates two fitness peaks associated with different ecological niches; (2) a genetic variation module that introduces random mutations at a configurable rate (μ) and effect size (σ), modelling the underlying loci as multi‑allelic to capture continuous colour and morphological variation; and (3) a mating/assortative‑pairing module that determines the probability of a mating event based on phenotypic similarity and a sexual‑selection strength parameter (s). The mating module can be toggled between “synchronised” (strict assortative mating within colour morphs) and “random” regimes, allowing the authors to explore how behavioural isolation interacts with ecological forces.

Two principal scenarios are explored. In the first, disruptive natural selection is strong. Even when the initial population is genetically homogeneous, the fitness landscape quickly bifurcates, driving sub‑populations toward opposite peaks. Under these conditions, genetic drift and mutation alone are sufficient to maintain divergence; assortative mating emerges as a by‑product of ecological segregation, and reproductive isolation solidifies without any explicit behavioural constraint. In the second scenario, disruptive selection is weak or absent. Here, ecological forces alone cannot generate stable divergence. The simulation shows that only when sexual selection is simultaneously strong and assortative mating is synchronised do the two colour‑morph clusters develop a reproductive barrier. In this regime, genetic drift plays a crucial role: stochastic fluctuations in allele frequencies allow each morph to accumulate private mutations, and the synchronised mating system prevents gene flow, allowing drift‑driven divergence to be amplified.

The study thus demonstrates that both ecological (disruptive) selection and sexual (assortative) selection can independently, or more realistically, synergistically, drive sympatric speciation. For the Midas cichlids, which exhibit rapid colour and morphological shifts, the model predicts that even a modest ecological gradient can be sufficient if coupled with strong mate choice, whereas strong mate choice can compensate for weak ecological differentiation. The authors also quantify the impact of drift velocity on the timing of reproductive isolation, highlighting that in small populations the stochastic component can be decisive.

Limitations are acknowledged. Parameter values are derived from literature and may not capture the full variability observed in natural lakes; the model collapses a multi‑dimensional environmental space into a single disruptive selection axis, ignoring factors such as temperature, pH, and vertical stratification. Moreover, the simulation is confined to a two‑dimensional spatial grid, omitting depth‑related habitat structuring that is known to influence cichlid behaviour. Future work is suggested to integrate field‑collected genomic and ecological data for parameter calibration, to expand the model to three dimensions, and to incorporate additional environmental variables.

In conclusion, the paper provides a robust computational framework that reproduces the observed polychromatism and polymorphism of three Midas cichlid species and clarifies the conditions under which sympatric speciation can occur. By contrasting scenarios with and without disruptive natural selection, it underscores the importance of synchronised mating systems and genetic drift in shaping the speciation trajectory. The findings have broader implications for understanding rapid adaptive radiations in other freshwater fishes and for designing empirical studies that test the relative contributions of ecological versus sexual selection in natural populations.