Friday, October 21, 2016

A very good introduction for Hardy-Weinberg_by Reed A. Cartwright

Website source: http://www.pandasthumb.org/archives/2004/04/evomath-2-testi-1.html

EvoMath 2: Testing for Hardy-Weinberg Equilibrium

In the first installment of EvoMath, I derived the Hardy-Weinberg Principle and discussed its significance to biology. In the second installment I will demonstrate how to test if a population deviates from Hardy-Weinberg equilibrium.

Recap:

A population is considered to be in Hardy-Weinberg equilibrium if the allele and genotype frequencies are as follows.
GenotypeFrequency
AA
Aa
aa
AlleleFrequency
A
a
evomath1_chart2.png

Test Procedure:

goodness-of-fit test can be used to determine if a population is significantly different from the expections of Hardy-Weinberg equilibrium. If we have a series of genotype counts from a population, then we can compare these counts to the ones predicted by the Hardy-Weinberg model. We conclude that the population is not in Hardy-Weinberg equilibrium if the probability that the counts were drawn under the Hardy-Weinberg model is too small for the deviations to be considered due to random chance. The significance level that is typically used is , i.e. the genotype counts have less than a one in twenty chance of being caused by a population in Hardy-Weinberg equilibrium.
In order to calculate this probability, we will use a test statistic, , which was devised in 1900 by Karl Pearson and has a well characterized distribution. If  are the set of observed counts, and  are the set of expected counts, then
This test statistic has a “chi-square” distribution with  degrees of freedom. Since we are testing Hardy-Weinberg equilibrium with two alleles,  (rationale not shown). Furthermore, it can be shown that if  then . Therefore, if  we will reject the null model and conclude that there is significant statistical support that the population is not in Hardy-Weinberg equilibrium.

Example 1:

Consider the following samples from a population.
GenotypeCount
AA30
Aa55
aa15
AlleleFrequency
A0.575
a0.425
Calculate the  value.
GenotypeObservedExpected(O-E)2/E
AA30330.27
Aa55490.73
aa25180.50
Total1001001.50
Since , we conclude that the genotype frequencies in this population are not significantly different than what would be expected if the population is in Hardy-Weinberg equilibrium.

Example 2:

Race and Sanger (1975) determined the blood groups of 1000 Britons as follows (from Hartl and Clarke 1997).
GenotypeObservedExpected
MM298294.3
MN489496.4
NN213209.3
This results in . As in the previous example, the measured genotype frequencies are not significantly different from the expectations of Hardy-Weinberg equilibrium.

Example 3:

Matthijis et al. (1998) surveyed a group of 54 people suffering from Jaeken syndrome (from Freeman and Herron 2004).
GenotypeObservedExpected
OO1119.44
OR4325.92
RR08.64
This results in . Unlike the previous two examples, the measured genotype frequencies are significantly different from the expectations of Hardy-Weinberg equilibrium. This indicates that one or more of the Hardy-Weinberg conditions are being violated; although, it does not tell us which ones.

Conclusion:

Although to derive the Hardy-Weinberg principle, we assumed that the size of the population was infinite, these statistical tests demonstrate that finite populations can approximately exist in Hardy-Weinberg equilibrium.

References

  • Freeman S and Herron JC (2004) Evolutionary Analysis 3rd ed. Pearson Education, Inc (Upper Saddle River, NJ)
  • Hartl DL and Clarke AG (1997) Principles of Population Genetics 3rd ed. Sinauer Associates, Inc (Sutherland, MA)
  • Matthijis GE et al. (1998) Lack of homozygotes for the most frequent disease allele in carbohydrate-deficient-glycoprotein syndrome type 1A. American Journal of Human Genetics 62: 542-550
  • Race RR and Sanger R (1975) Blood Groups in Man 6th ed. JB Lippincott, Philadelphia

Previous Installments:


Comments: Post a Comment



<< Home

This page is powered by Blogger. Isn't yours?