MultSim: a program to simulate founder effects

To evaluate various colonization models, we wrote a computer program to simulate founder events from the allele frequency distributions observed in the source populations and allele frequency differences between source and sink populations. In short, the program randomly selected alleles from the source population based on their relative frequencies and begins a new population with a specified number of haploid colonists (thus approximating chromosomes rather than individuals). The program then calculates the maximum and minimum number of colonists that could have founded the observed sink population under several scenarios. For example, if the number of colonists was very large, virtually no difference in allele frequency should be observed between the source and sink populations. If the number was very small, there is a high probability of getting a rather large difference in allele frequencies. Our program samples from the source population allele frequencies multiple times and created a frequency distribution of possible genetic divergences, measured as 1-Psa (Bowcock et al 1994). It uses this distribution to derive the 95% confidence intervals of numbers of colonists founding the sink population under each scenario for each locus.

Secondarily, if no alleles are lost in the founding of the new population, one can also perform simulations as above to determine the number of colonists at which there is a 95% chance of losing one or more alleles under each colonization scenario for each locus. This provides a second estimate of the minimum number of colonists that founded the sink population.

Three types of colonization events can be simulated. First, a founding event followed by an instantaneous increase to very large population size (FAST). Thus, all the allele frequency differences between the source and sink populations result from the initial founder event. Second, a founding event followed by an increase to a very large population size (>5000) with increases by a factor of X per generation (SLOW). This model may be more realistic, as genetic drift occurs in the early stages of the colonization after the initial founding event. Finally, the sink population is founded after a series of Y bottlenecks of the specified number of haploid colonists, thus simulating a stepping-stone process in the spread of the species (Y-STEP).

There are some assumptions inherent in these simulations. Our method uses a haploid model rather than a diploid model. Hence, some evolutionary dynamics may not be captured, such as skewed sex ratios or high variance in male mating success. Also, the model assumes that no genetic drift or selection has occurred in the source population subsequent to the migrants leaving to the sink population. We urge potential users of this model to keep these assumptions in mind when performing their calculations.

If you're STILL interested, the files are below. Please cite the following if you use this software:

Noor, M. A. F., M. Pascual, and K. R. Smith. 2000. Genetic variation in the spread of Drosophila subobscura from a nonequilibrium population. Evolution, 54: 696-703.
See Abstract.

Program (Multsim.ZIP), version 1.1, 55K
includes PsaCalc.EXE, a program for calculating 1-proportion shared alleles (1-Psa), 31K

 

Bug fixes

Version 1.1 fixes the request for critical Psa value to the critical (1-Psa) value

 

Reference : Bowcock, A. M., A. Ruiz-Linares, J. Tomfohrde, E. Minch, J. R. Kidd, and L. L. Cavalli-Sforza. 1994. High resolution of human evolutionary trees with polymorphic microsatellites. Nature 368:455-457.

Link back to Noor lab software page.