General Instructions: Use the free allele simulation software available from:
http://faculty.washington.edu/herronjc/SoftwareFolder/AlleleA1.html
It does not run on iPad, but you can get versions for either MacOS X or Windows.
Install it and then open it on your system; you should see this:
Macintosh HD:Users:davegray:Desktop:allele main screen.tiff
Familiarize yourself with the basic controls.
Basically the software graphs allele frequency versus time for the “A1” allele (i.e. “p”) at a two allele locus. By clicking the arrow pointing down on the left side of the graph just below “Frequency of allele A1” you can opt to graph the frequency of the A2 allele (i.e. “q”), or genotypic frequencies (which may or may not equal p2 2pq and q2 depending on whether or not your population is meeting the assumptions of Hardy-Weinberg equilibrium, hereafter HW).
Click the run button on the bottom left. You should get a straight line across the screen, and to the right of the graph should then appear “Final Frequencies” of alleles and genotypes. This is after the 500 generations of the X-axis of the graph. You can make the number of generations more or less by clicking the right facing arrow at the end of that axis.
Let’s interpret what happened: p started as 0.5, and after 500 generations was still 0.5, so we would say that there was no evolution at this locus; also note that the genotype frequencies are as we would expect from p2 + 2pq + q2 = 1. HW says this happens when there is (1) no selection, (2) no mutation, (3) no gene flow, (4) no genetic drift, AND (5) random mating. The AND is important, all conditions must be satisfied. What the allele lab allows one to do is systematically violate one or more of those five conditions.
Below is another screenshot. Red oval is where you can change genotypic fitness (really the average fitness of the phenotype produced by that genotype), i.e. selection; blue oval is where you can change mutation rate A1 to A2 and reverse; yellow box is about gene flow, you can change the fraction of migrants each generation and the allele frequency in the source population; black arrow points to where you can control effective population size, which is about genetic drift; grey arrow points to how you control if mating is random or not.
Macintosh HD:Users:davegray:Desktop:allele main screen.tiff
You can run multiple trials for any given set of conditions; you can graph a single line which disappears if you run it again (single), have multiple lines all one color (multiple), or have multiple lines of different colors (auto).
You can clear lines without changing your experimental conditions by clicking clear and you can put all conditions back to the default starting conditions by clocking reset.
Now you are ready to do the lab. Click reset.
Part 1: Selection alone
Selection against a deleterious recessive allele.
Change the Fitness of the A2A2 genotype to 0.5, leaving the others as 1.0. This is equivalent to modeling a disease caused by being homozygous for deleterious recessive alleles. How do you know that the A2 allele is recessive? Because the A1A2 heterozygous genotype is equivalent to the A1A1 homozygous genotype, both have fitness of 1.0. Click run.
What happened?
What is the final frequency of the A1 allele after 500 generations?
Why hasn’t the “bad” A2 allele gone to frequency zero? (hint, after 500 generations what is the frequency of the homozygous recessive A2A2 genotype?).
Maybe it will go to zero if you run it longer; try 1000 generations. What happened?
Selection against a deleterious dominant allele.
Click reset. Then change the genotypic fitness to model selection against a deleterious dominant allele. How would you do this? If it is dominant, then A1A1 fitness and A1A2 fitnesses are equal to each other and lower than the fitness of the A2A2 genotype. So try A1A1 = 0.5, A1A2 = 0.5, and A2A2 = 1.0. Click run.
What happened?
What is the final frequency of the A1 allele after 500 generations?
How fast was the A2 allele lost? If you can’t see it on your graph then you can click clear, change the number of generations to 25, click run, and get a better view.
What if selection were much much much weaker, for example click reset and then try genotypic fitnesses of A1A1 = 0.95, A1A2 = 0.95, and A2A2 = 1.0. Click run.
What happened?
What is the final frequency of the A1 allele after 500 generations?
Selection for a beneficial recessive allele.
Click reset. How would you model selection for a beneficial recessive?
Record the genotypic fitnesses you use to simulate this.
A1A1 =
A1A2 =
A2A2 =
Click run. What happened?
What is the final frequency of the A1 allele after 500 generations?
How fast was the A1 allele lost (if it was)? If you can’t see it on your graph then you can change generation numbers and do it again.
Is modeling selection for a beneficial recessive different from modeling selection against a deleterious dominant? Why or why not?
How do you think the pattern and the outcome of evolution would differ if it you modeled selection for a beneficial dominant allele instead of a beneficial recessive allele? Seriously, what is your prediction?
Selection for a beneficial dominant allele.
OK, now you can test that prediction. Click reset, and then set up a selection scenario to test your idea. What is your scenario?
Record the genotypic fitnesses you use to simulate this.
A1A1 =
A1A2 =
A2A2 =
Click run. What happened?
What is the final frequency of the beneficial allele after 500 generations [probably this will be the A1 allele, but some of you might have set up a scenario that makes the A2 allele the dominant beneficial one]?
How fast was the deleterious allele lost (if it was)? If you can’t see it on your graph then you can change generation numbers and do it again.
Was your prediction supported?
Part 2: Mutation alone
This one is quick and easy. Click reset. Then make the mutation rate of A1 to A2 something that is high but biologically plausible, e.g. 0.0001, and the reverse mutation rate A2 to A1 just leave at zero. This simulates mutation as strongly as is biologically realistic, so we are giving mutation the best chance to show us what it can do. Click run, what happens?
So what do you conclude about the evolutionary power of mutation (all by itself) to change allele frequencies? Is mutation an important evolutionary force (obviously it is incredibly important in providing genetic variation for selection to act upon, but is mutation itself a significant force causing microevolutionary change)?
Part 3: Gene flow alone
Also quick and easy. Click reset. Note the default settings: A1 allele starts in your population at frequency 0.5, with no immigration from some other source population which also has A1 allele frequency of 0.5. Click run and nothing changes; you did that before – it is just HW.
Now change the fraction of migrants in each generation to 0.1, meaning 10% of your new population each generation will be made up of individuals that came from the other source population. Click run and nothing changes, why?
Recall that what gene flow does is to make populations more similar to each other. If your source population is not different in frequency of A1 allele, then Fst = 0 and the frequency of A1 allele in your population won’t change with migration. So now make migrants come from a population that is genetically different from your population.