I am in Brazil this week to give a talk about speciation in continuous populations, so I figure I will save a little effort by summarizing some of the stuff that is in that talk. By the way, Sao Carlos is a wonderful town, and if you can come up with an excuse to come to Brazil I strongly recommend it. Of course, if you are a vegetarian, while it SHOULD be a good place for you, man do these folks like their meat!
So, on to speciation. This is a project I did with Maggie Eppstein, currently chair of the University of Vermont Computer Science department, and Josh Payne, then a computer science graduate student and now a postdoc at the University of Zurich. In other words, it was big on computer science, and maybe a bit less big on biology, but the implications for biology are important. (Payne, J.L., Eppstein, M.J., and Goodnight, C.J. “Sensitivity of Self-Organized Speciation to Long Distance Dispersal”, Proceedings of the 2007 IEEE Symposium on Artificial Life (Alife’07), pp. 1-7, 2007. Eppstein, M.J., Payne, J. L., and Goodnight, C.J., “Underdominance, Multiscale Interactions, and Self-Organizing Barriers to Gene Flow”, J. Artificial Evolution and Applications (special issue on Biological Applications), Volume 2009, Article ID 725049, 13 pages, 2009
This project stemmed out of a discussion I had with several people at NECSI (http://necsi.edu/). So, in talking to them I came to realize that when you had two types that were incompatible, they were distributed randomly on a plain, interactions were local, and the types spread contagiously that there would be this interesting process of coarsening. That is if you started out with the types randomly distributed the locally common type would have an advantage and increase, whereas the locally rare type would have a disadvantage and decrease. The net result would be your random distribution would devolve into regions that were primarily one type, and other regions that were the opposite type. I also found out that the boundary between these regions would wander about randomly across the plain, but it would get trapped by any sort of boundary
In any case the question becomes whether speciation this process of coarsening in populations with local interaction lead to speciation. We decided to look at the very simplest case, that of simple underdominance. That is, can we take a population with within locus underdominance, and have it self organize into two species.
By the way, as with any good theoretician, I will use Mayr’s biological species concept as my definition of a species. I know that this is really not a great definition, but heck, I got to meet Professor Mayr on several occasions and he said nice things about my research, so I stand by my decision. . .
First off, the trivial case is where the underdominance is so extreme that the heterozygote is lethal:
Genotype A1A1 A1A2 A2A2
Fitness 1 0 1
Obviously, such a population will immediately split into two reproductively isolated “species”. Yes, this works, but it is biologically totally unconvincing.
I won’t go into the details, however, what we can show is that it is quite easy to introduce mild underdominance into a viscous population. Thus, we would be hardly surprised to see underdominance such as:
Genotype A1A1 A1A2 A2A2
Fitness 1 0.9 1
The problem is that this is not an effective barrier to gene flow. Indeed we were able to show that it takes a complete reproductive isolation to prevent the movement of a neutral gene across the boundary.
Thus, we are only going to allow mild underdominance, yet we want the population to divide into two completely reproductively isolated populations. It turns out that it is not that hard if we allow multiple loci with mild underdominance. To illustrate this consider two loci with underdominance of 50%, so that the double heterozygote is lethal.
In other words, simply by having localized mating and underdominance we are half way there. We get the coarsening, and regions of the two species, but in all cases there is a pathway through a viable hybrid zone between any two regions.
This is where we had to introduce some epistasis. In our next iteration we added a bit of epistasis so that the two opposite corner genotypes were favored:
Two underdominant loci. Double homozygotes of the same numerical value have the highest fitness, single heterozygotes (yellow) have lower fitness, double heterozygote (white) is lethal.
This now results in coarsening and eventually speciation:
This is actually very interesting, because you will note that in the early stages all four homozygous forms are formed, but the single locus heterozygous boundaries wander around randomly. When two such come into contact they coalesce into a single stronger boundary, eventually leading to speciation. Importantly, this is not limited to only two loci. It turns out that this coarsening and coalescence of leaky boundaries will continue regardless of how many loci are involved. Thus, even extremely mild underdominance at a large number of loci will eventually lead to speciation in this model:
This is a very simple simulation, yet it makes the important point that when there are genetic incompatibilities speciation can easily occur. Indeed, this implies that the speciation may be the expected outcome for a widespread species with limited gene flow.