Thursday, September 4, 2014

LESSON: Selection Pressure, Part Two

In this lesson, we’re going to pick up where my last lesson on selection pressure left off.  Here’s the link to the first lesson if you’d like to familiarize yourself with that first.

Today’s lesson focuses on how selection pressure changes allele frequencies over many generations.  Our sample problem is a follow-up to the scenario in my first lesson:

Despite your good work in Capitol City last time, the strange epidemic persisted and all individuals who lack Factor G are now sterile.  Individuals who lack Factor G are genotype gg; individuals who produce Factor G are GG or Gg and have normal levels of fertility.  How many generations will pass until the recessive g allele has a frequency of less than 1%?

In my last lesson, we figured out that when gg individuals cannot reproduce, the frequency of g drops from 0.3 in the parent population to 0.231 in the first generation after selection.  The simple formula for figuring out g’ if gg individuals are not reproducing is g’ = g/(1+g).

Here, g’ refers to the frequency of g in the next generation.  (And if you want to see the derivation of that simple formula, see my last lesson.)

If we do the math, we see that after one generation of selection, g decreased by 23%:

Slide 1 cropped

So how do we find the generation at which allele g is less than 1%?

Before we jump to those calculations, let’s consider what these numbers really mean.  After the first round of selection, the frequency of g is 23.1% or 0.231.  This number includes the individuals who are gg, even though they are sterile.  Why is that?  Because when Gg heterozygotes mate, their children will be:

1 GG: 2 Gg: 1 gg

The fact that Gg individuals have normal fertility means that there will be gg children, even though those children are not able to reproduce.

Also, let’s consider the number of Gg heterozygotes in the population after one round of selection pressure.  Again, the frequency of g is 0.231.  The number of Gg hets = 2pq.

2pq = 2(0.769)(0.231) = 0.355 or 35.5%

35% of the population here is Gg, which means that g is being propagated by about a third of the population (a substantial fraction!).

Let’s move on to the calculations to determine when the frequency of g will be less than 1%.

I doubt that you’d be asked to solve a problem like this on an exam without being able to use Excel or a similar program.  The calculations are easy on a spreadsheet but very tedious to do by hand.  I used Excel to solve this problem.  Here’s what my spreadsheet looked like to set up the calculations:

Slide 2

Here are the entries I put into Excel to set up this spreadsheet: 

Cell B1: 0.3 (original value of g before selection)

Cell B6: =B1/(1+B1)

(this is the calculation for g after each round of selection: g’ = g/(1+g))

Cell B7: =B6/(1+B6)

Cells B8-B14: Copy/Paste from cell B7 down to B14.

From the image above, you can see that after 9 rounds of selection (Row 14), the frequency of g is 0.081 or 8.1%.  So we’re not done, and we can continue our calculations by copy/pasting Cells A14 and B14 until we reach…

Slide 3

…Cell 102!

We’re looking for the generation at which q drops below 0.01, which happens to occur at the 97th generation of selective pressure.  What’s so interesting to me as a geneticist is that g (or q, in Hardy-Weinberg terms) can persist for so long when gg individuals are sterile.  This example illustrates why recessive deleterious alleles are not easily eliminated from the population.  As long as Gg heterozygotes are healthy, then g will be in the allele pool for a long time.

So that was a lesson on selective pressure and Microsoft Excel techniques.  Got questions, comments, or something else to say?  Tell me in the comments below! 

Thanks for reading!

A Little Something Inspirational

As your new semester gets underway, remember that…

“Education is not filling a bucket, but lighting a fire.” ~William Yates

Happy learning to you!

Monday, September 1, 2014

LESSON: Selection Pressure and Changing Allele Frequencies

We’ve used Hardy-Weinberg equilibria to determine phenotype frequencies in a population.  The Hardy-Weinberg equations assume that a population is stable.  In other words, allele frequencies are not changing—there’s no selection pressure.

But what happens when selection is acting upon a population?  Theoretically, we understand that selection will favor certain alleles or allele combinations over others.  This will shift the balance of alleles toward a new equilibrium.  If we know the magnitude of the selection pressure, we can calculate the effect of selection on allele frequencies.

Let’s work through a sample problem to unpack this set of questions. 

You and your medical team are summoned to Capitol City, where a strange epidemic has rendered 9% of the population sterile. Working rapidly, you discover that there’s a strong correlation between Factor G, a protein found in blood, and the fertile residents: all the fertile residents test positive for Factor G, but all the infertile residents test negative for Factor G. Later genetic and biochemical tests reveal that the population of Capitol City carries two alleles for a gene that is necessary for the production of Factor G such that G is dominant to g. All the fertile residents are genotypically GG or Gg. All gg individuals are now sterile. The original frequencies of G and g were as follows:

G = 0.7

g = 0.3

After the epidemic, what are the allele frequencies of G and g in the next generation?

Let’s define a new term first.

* Fitness, W.  Fitness is a measure of reproductive success and should be a value between 0 and 1.  If there is no selection pressure reducing the reproductive success of a genotype, then W = 1.  If a genotype cannot reproduce (as in our sample problem above, where gg individuals are sterile), then W = 0.

So for our sample problem:

For GG, W = 1

For Gg, W = 1

For gg, W = 0.

We’ll make use of these values below. 

Now we need to dive deeper into the math to connect s to Hardy-Weinberg equilibria.

Slide1 cropped

Slide2

We could have predicted that if q2 = 0, then p2 + 2pq = 1.  So that set of calculations confirms our prediction, but we still have no idea what p and q are after selection.  We need another set of equations for that task.

Slide3

Slide4 cropped

Slide5 cropped

Slide6

Slide 7 redo cropped

Slide 8 cropped

And there you have it!  Now it’s your turn: what if the epidemic, rather than making gg individuals sterile, reduced their fertility by 50%  What would the frequency of G and g be after one round of selection pressure?  (I’ll provide or confirm the answer in the comments when someone asks for it.) 

References:

Comprehensive Genetics coursepack, 2014 version, published by Dr. John Ellison, Texas A&M University