Multilocus model of sympatric speciation. III. Computer simulations.

In the previous papers some analytical results were obtained for the limit stages of sympatric speciation. The present paper aims at finding the scope of validity for these results. The Monte Carlo computer model of this process was created and studied. We deal with two aspects of the speciation process: the development of reproductive isolation between the forming species and the extinction of the intermediate individuals. It has been shown that the advantage of the allele of reproductive isolation increases with the growth of its frequency. The extinction of the intermediates goes differently with various numbers of loci involved in speciation. If reproductive isolation is due to differences in two or four loci, then the completion of extinction of the intermediates requires the strongest disruptive selection, so that the necessary conditions for speciation found previously are also proved to be sufficient. But with eight and probably with a larger number of loci, the selection required to promote speciation in a population that is far from it is considerably stronger than selection needed at the last stage of speciation. Consequently, under some intensities of disruptive selection the final state of a population depends initial state. The conditions under which the stationary state of a population is characterized by bimodal distribution of phenotypes are also found.