Misplaced Pages

Hardy–Weinberg principle: Difference between revisions

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Browse history interactively← Previous editNext edit →Content deleted Content addedVisualWikitext
Revision as of 12:50, 7 April 2014 editClueBot NG (talk | contribs)Bots, Pending changes reviewers, Rollbackers6,439,221 editsm Reverting possible vandalism by 66.114.15.49 to version by Duxwing. False positive? Report it. Thanks, ClueBot NG. (1782866) (Bot)← Previous edit Revision as of 06:21, 11 April 2014 edit undoDicklyon (talk | contribs)Autopatrolled, Extended confirmed users, Rollbackers477,317 edits Reverted to revision 598232216 by XLinkBot (talk): Back to before Duxwing's odd edits. (TW)Next edit →
Line 6: Line 6:
]s: the horizontal axis shows the two ] ''p'' and ''q'' and the vertical axis shows the expected ]. Each line shows one of the three possible genotypes.]] ]s: the horizontal axis shows the two ] ''p'' and ''q'' and the vertical axis shows the expected ]. Each line shows one of the three possible genotypes.]]


The '''Hardy–Weinberg principle''' (also known as the '''Hardy–Weinberg equilibrium, model, theorem''', or '''law''') states that ] and ] frequencies remain constant absent evolutionary influences. These influences include ''non-random ]'', '']'', '']'', '']'', '']'' and '']''. Because one or more of these influences typically are in real populations, the Hardy–Weinberg principle describes an ideal condition against which the effects of these influences can be analyzed. The '''Hardy–Weinberg principle''' (also known as the '''Hardy–Weinberg equilibrium, model, theorem''', or '''law''') states that ] and ] frequencies in a population will remain constant from generation to generation in the absence of other evolutionary influences. These influences include ''non-random ]'', '']'', '']'', '']'', '']'' and '']''. Because one or more of these influences are typically present in real populations, the Hardy–Weinberg principle describes an ideal condition against which the effects of these influences can be analyzed.


In the simplest case of a single locus with two ]s denoted A and a with frequencies {{math|<var>f</var>(A) {{=}} ''p''}} and {{math|<var>f</var>(a) {{=}} ''q''}}, respectively, the expected genotype frequencies are {{math|<var>f</var>(AA) {{=}} ''p''<sup>2</sup>}} for the AA ]s, {{math|<var>f</var>(aa) {{=}} ''q''<sup>2</sup>}} for the aa homozygotes, and {{math|<var>f</var>(Aa) {{=}} 2''pq''}} for the ]s. The genotype proportions ''p''<sup>2</sup>, 2''pq'', and ''q''<sup>2</sup> are called the Hardy-Weinberg proportions. ] of the square of the sum of ''p'' and ''q'', and such a sum, as it represents the total of all possibilities, must be equal to 1. Therefore {{math| (<var>p</var> + <var>q</var>)<sup>2</sup> {{=}} ''p''<sup>2</sup> + 2''pq'' + ''q''<sup>2</sup> {{=}} 1}}. The solution of this equation is {{math|''q'' {{=}} 1 - ''p''}}.] In the simplest case of a single locus with two ]s denoted A and a with frequencies {{math|<var>f</var>(A) {{=}} ''p''}} and {{math|<var>f</var>(a) {{=}} ''q''}}, respectively, the expected genotype frequencies are {{math|<var>f</var>(AA) {{=}} ''p''<sup>2</sup>}} for the AA ]s, {{math|<var>f</var>(aa) {{=}} ''q''<sup>2</sup>}} for the aa homozygotes, and {{math|<var>f</var>(Aa) {{=}} 2''pq''}} for the ]s. The genotype proportions ''p''<sup>2</sup>, 2''pq'', and ''q''<sup>2</sup> are called the Hardy-Weinberg proportions. ] of the square of the sum of ''p'' and ''q'', and such a sum, as it represents the total of all possibilities, must be equal to 1. Therefore {{math| (<var>p</var> + <var>q</var>)<sup>2</sup> {{=}} ''p''<sup>2</sup> + 2''pq'' + ''q''<sup>2</sup> {{=}} 1}}. The solution of this equation is {{math|''q'' {{=}} 1 - ''p''}}.]

Revision as of 06:21, 11 April 2014

This article has multiple issues. Please help improve it or discuss these issues on the talk page. (Learn how and when to remove these messages)
This article includes a list of general references, but it lacks sufficient corresponding inline citations. Please help to improve this article by introducing more precise citations. (May 2011) (Learn how and when to remove this message)
This article may be too technical for most readers to understand. Please help improve it to make it understandable to non-experts, without removing the technical details. (February 2012) (Learn how and when to remove this message)
(Learn how and when to remove this message)
Hardy–Weinberg proportions for two alleles: the horizontal axis shows the two allele frequencies p and q and the vertical axis shows the expected genotype frequencies. Each line shows one of the three possible genotypes.

The Hardy–Weinberg principle (also known as the Hardy–Weinberg equilibrium, model, theorem, or law) states that allele and genotype frequencies in a population will remain constant from generation to generation in the absence of other evolutionary influences. These influences include non-random mating, mutation, selection, genetic drift, gene flow and meiotic drive. Because one or more of these influences are typically present in real populations, the Hardy–Weinberg principle describes an ideal condition against which the effects of these influences can be analyzed.

In the simplest case of a single locus with two alleles denoted A and a with frequencies f(A) = p and f(a) = q, respectively, the expected genotype frequencies are f(AA) = p for the AA homozygotes, f(aa) = q for the aa homozygotes, and f(Aa) = 2pq for the heterozygotes. The genotype proportions p, 2pq, and q are called the Hardy-Weinberg proportions.

If union of gametes to produce the next generation is random, it can be shown that the new frequency f′ satisfies f ( A ) = f ( A ) {\displaystyle \textstyle f'({\text{A}})=f({\text{A}})} and f ( a ) = f ( a ) {\displaystyle \textstyle f'({\text{a}})=f({\text{a}})} . That is, allele frequencies are constant between generations.

This principle was named after G. H. Hardy and Wilhelm Weinberg, who first demonstrated it mathematically.

Derivation

Consider a population of monoecious diploids, where each organism produces male and female gametes at equal frequency, and has two alleles at each gene locus. Organisms reproduce by random union of gametes (the “gene pool” population model). A locus in this population has two alleles, A and a, that occur with initial frequencies f0(A) = p and f0(a) = q, respectively. The allele frequencies at each generation are obtained by pooling together the alleles from each genotype of the same generation according to the expected contribution from the homozygote and heterozygote genotypes, which are 1 and 1/2, respectively:

f t ( A ) = f t ( AA ) + 1 2 f t ( Aa ) {\displaystyle f_{t}({\text{A}})=f_{t}({\text{AA}})+{\frac {1}{2}}f_{t}({\text{Aa}})} 1
f t ( a ) = f t ( aa ) + 1 2 f t ( Aa ) {\displaystyle f_{t}({\text{a}})=f_{t}({\text{aa}})+{\frac {1}{2}}f_{t}({\text{Aa}})} 2

The different ways to form genotypes for the next generation can be shown in a Punnett square, where the proportion of each genotype is equal to the product of the row and column allele frequencies from the current generation.

Table 1: Punnett square for Hardy–Weinberg
Females
A (p) a (q)
Males A (p) AA (p) Aa (pq)
a (q) Aa (qp) aa (q)

The sum of the entries is p + 2pq + q = 1, as the genotype frequencies must sum to one.

Note again that as p + q = 1, the binomial expansion of (p + q) = p + 2pq + q = 1 gives the same relationships.

Summing the elements of the Punnet square or the binomial expansion, we obtain the expected genotype proportions among the offspring after a single generation:

f 1 ( AA ) = p 2 = f 0 ( A ) 2 {\displaystyle f_{1}({\text{AA}})=p^{2}=f_{0}({\text{A}})^{2}} 3
f 1 ( Aa ) = p q + q p = 2 p q = 2 f 0 ( A ) f 0 ( a ) {\displaystyle f_{1}({\text{Aa}})=pq+qp=2pq=2f_{0}({\text{A}})f_{0}({\text{a}})} 4
f 1 ( aa ) = q 2 = f 0 ( a ) 2 {\displaystyle f_{1}({\text{aa}})=q^{2}=f_{0}({\text{a}})^{2}} 5

These frequencies define the Hardy–Weinberg equilibrium. It should be mentioned that the genotype frequencies after the first generation need not equal the genotype frequencies from the initial generation, e.g. f1(AA) ≠ f0(AA). However, the genotype frequencies for all future times will equal the Hardy-Weinberg frequencies, e.g. ft(AA) = f1(AA) for t > 1. This follows since the genotype frequencies of the next generation depend only on the allele frequencies of the current generation which, as calculated by equations (1) and (2), are preserved from the initial generation:

f 1 ( A ) = f 1 ( AA ) + 1 2 f 1 ( Aa ) = p 2 + p q = p ( p + q ) = p = f 0 ( A ) {\displaystyle f_{1}({\text{A}})=f_{1}({\text{AA}})+{\frac {1}{2}}f_{1}({\text{Aa}})=p^{2}+pq=p\left(p+q\right)=p=f_{0}({\text{A}})}
f 1 ( a ) = f 1 ( aa ) + 1 2 f 1 ( Aa ) = q 2 + p q = q ( p + q ) = q = f 0 ( a ) {\displaystyle f_{1}({\text{a}})=f_{1}({\text{aa}})+{\frac {1}{2}}f_{1}({\text{Aa}})=q^{2}+pq=q\left(p+q\right)=q=f_{0}({\text{a}})}

For the more general case of dioecious diploids that reproduce by random mating of individuals, it is necessary to calculate the genotype frequencies from the nine possible matings between each parental genotype (AA, Aa, and aa) in either sex, weighted by the expected genotype contributions of each such mating. Equivalently, one considers the six unique diploid-diploid combinations:

[ ( AA , AA ) , ( AA , Aa ) , ( AA , aa ) , ( Aa , Aa ) , ( Aa , aa ) , ( aa , aa ) ] {\displaystyle \left}

and constructs a Punnett square for each, so as to calculate its contribution to the next generation's genotypes. These contributions are weighted according to the probability of each diploid-diploid combination, which follows a Multinomial distribution with k = 3. For example, the probability of the mating combination (AA,aa) is 2 ft(AA)ft(aa) and it can only result in the Aa genotype: . Overall, the resulting genotype frequencies are calculated as:

[ f t + 1 ( AA ) , f t + 1 ( Aa ) , f t + 1 ( aa ) ] = f t ( AA ) f t ( AA ) [ 1 , 0 , 0 ] + 2 f t ( AA ) f t ( Aa ) [ 1 / 2 , 1 / 2 , 0 ] + 2 f t ( AA ) f t ( aa ) [ 0 , 1 , 0 ] + f t ( Aa ) f t ( Aa ) [ 1 / 4 , 1 / 2 , 1 / 4 ] + 2 f t ( Aa ) f t ( aa ) [ 0 , 1 / 2 , 1 / 2 ] + f t ( aa ) f t ( aa ) [ 0 , 0 , 1 ] = [ ( f t ( AA ) + f t ( Aa ) 2 ) 2 , 2 ( f t ( AA ) + f t ( Aa ) 2 ) ( f t ( aa ) + f t ( Aa ) 2 ) , ( f t ( aa ) + f t ( Aa ) 2 ) 2 ] = [ f t ( A ) 2 , 2 f t ( A ) f t ( a ) , f t ( a ) 2 ] {\displaystyle {\begin{aligned}&\left\\&\quad =f_{t}({\text{AA}})f_{t}({\text{AA}})\left+2f_{t}({\text{AA}})f_{t}({\text{Aa}})\left+2f_{t}({\text{AA}})f_{t}({\text{aa}})\left\\&\quad \quad +f_{t}({\text{Aa}})f_{t}({\text{Aa}})\left+2f_{t}({\text{Aa}})f_{t}({\text{aa}})\left+f_{t}({\text{aa}})f_{t}({\text{aa}})\left\\&\quad =\left\\&\quad =\left\end{aligned}}}

As before, one can show that the allele frequencies at time t+1 equal those a time t, and so, are constant in time. Similarly, the genotype frequencies depend only on the allele frequencies, and so, after time t=1 are also constant in time.

If in either monoecious or dioecious organisms, either the allele or genotype proportions are initially unequal in either sex, it can be shown that constant proportions are obtained after one generation of random mating. If dioecious organisms are heterogametic and the gene locus is located on the X chromosome, it can be shown that if the allele frequencies are initially unequal in the two sexes , f′(a) in the heterogametic sex ‘chases’ f(a) in the homogametic sex of the previous generation, until an equilibrium is reached at the weighted average of the two initial frequencies.

Deviations from Hardy–Weinberg equilibrium

The seven assumptions underlying Hardy–Weinberg equilibrium are as follows:

  • organisms are diploid
  • only sexual reproduction occurs
  • generations are non overlapping
  • mating is random
  • population size is infinitely large
  • allele frequencies are equal in the sexes
  • there is no migration, mutation or selection

Violations of the Hardy–Weinberg assumptions can cause deviations from expectation. How this affects the population depends on the assumptions that are violated.

  • Random mating. The HWP states the population will have the given genotypic frequencies (called Hardy–Weinberg proportions) after a single generation of random mating within the population. When the random mating assumption is violated, the population will not have Hardy–Weinberg proportions. A common cause of non-random mating is inbreeding, which causes an increase in homozygosity for all genes.

If a population violates one of the following four assumptions, the population may continue to have Hardy–Weinberg proportions each generation, but the allele frequencies will change over time.

  • Selection, in general, causes allele frequencies to change, often quite rapidly. While directional selection eventually leads to the loss of all alleles except the favored one, some forms of selection, such as balancing selection, lead to equilibrium without loss of alleles.
  • Mutation will have a very subtle effect on allele frequencies. Mutation rates are of the order 10 to 10, and the change in allele frequency will be, at most, the same order. Recurrent mutation will maintain alleles in the population, even if there is strong selection against them.
  • Migration genetically links two or more populations together. In general, allele frequencies will become more homogeneous among the populations. Some models for migration inherently include nonrandom mating (Wahlund effect, for example). For those models, the Hardy–Weinberg proportions will normally not be valid.
  • Small population size can cause a random change in allele frequencies. This is due to a sampling effect, and is called genetic drift. Sampling effects are most important when the allele is present in a small number of copies.

Sex linkage

Where the A gene is sex linked, the heterogametic sex (e.g., mammalian males; avian females) have only one copy of the gene (and are termed hemizygous), while the homogametic sex (e.g., human females) have two copies. The genotype frequencies at equilibrium are p and q for the heterogametic sex but p, 2pq and q for the homogametic sex.

For example, in humans red–green colorblindness is an X-linked recessive trait. In western European males, the trait affects about 1 in 12, (q = 0.083) whereas it affects about 1 in 200 females (0.005, compared to q = 0.007), very close to Hardy–Weinberg proportions.

If a population is brought together with males and females with a different allele frequency in each subpopulation (males or females), the allele frequency of the male population in the next generation will follow that of the female population because each son receives its X chromosome from its mother. The population converges on equilibrium very quickly.

Generalizations

The simple derivation above can be generalized for more than two alleles and polyploidy.

Generalization for more than two alleles

Consider an extra allele frequency, r. The two-allele case is the binomial expansion of (p + q), and thus the three-allele case is the trinomial expansion of (p + qr).

( p + q + r ) 2 = p 2 + q 2 + r 2 + 2 p q + 2 p r + 2 q r {\displaystyle (p+q+r)^{2}=p^{2}+q^{2}+r^{2}+2pq+2pr+2qr\,}

More generally, consider the alleles A1, ..., An given by the allele frequencies p1 to pn;

( p 1 + + p n ) 2 {\displaystyle (p_{1}+\cdots +p_{n})^{2}\,}

giving for all homozygotes:

f ( A i A i ) = p i 2 {\displaystyle f(A_{i}A_{i})=p_{i}^{2}\,}

and for all heterozygotes:

f ( A i A j ) = 2 p i p j {\displaystyle f(A_{i}A_{j})=2p_{i}p_{j}\,}

Generalization for polyploidy

The Hardy–Weinberg principle may also be generalized to polyploid systems, that is, for organisms that have more than two copies of each chromosome. Consider again only two alleles. The diploid case is the binomial expansion of:

( p + q ) 2 {\displaystyle (p+q)^{2}\,}

and therefore the polyploid case is the polynomial expansion of:

( p + q ) c {\displaystyle (p+q)^{c}\,}

where c is the ploidy, for example with tetraploid (c = 4):

Table 2: Expected genotype frequencies for tetraploidy
Genotype Frequency
AAAA p 4 {\displaystyle p^{4}}
AAAa 4 p 3 q {\displaystyle 4p^{3}q}
AAaa 6 p 2 q 2 {\displaystyle 6p^{2}q^{2}}
Aaaa 4 p q 3 {\displaystyle 4pq^{3}}
aaaa q 4 {\displaystyle q^{4}}

Depending on whether the organism is a 'true' tetraploid or an amphidiploid will determine how long it will take for the population to reach Hardy–Weinberg equilibrium.

Complete generalization

For n {\displaystyle n} distinct alleles in c {\displaystyle c} -ploids, the genotype frequencies in the Hardy–Weinberg equilibrium are given by individual terms in the multinomial expansion of ( p 1 + + p n ) c {\displaystyle (p_{1}+\cdots +p_{n})^{c}} :

( p 1 + + p n ) c = k 1 , , k n   N : k 1 + + k n = c ( c k 1 , , k n ) p 1 k 1 p n k n {\displaystyle (p_{1}+\cdots +p_{n})^{c}=\sum _{k_{1},\ldots ,k_{n}\ \in \mathbb {N} :k_{1}+\cdots +k_{n}=c}{c \choose k_{1},\ldots ,k_{n}}p_{1}^{k_{1}}\cdots p_{n}^{k_{n}}}

Applications

The Hardy–Weinberg principle may be applied in two ways, either a population is assumed to be in Hardy–Weinberg proportions, in which the genotype frequencies can be calculated, or if the genotype frequencies of all three genotypes are known, they can be tested for deviations that are statistically significant.

Application to cases of complete dominance

Suppose that the phenotypes of AA and Aa are indistinguishable, i.e., there is complete dominance. Assuming that the Hardy–Weinberg principle applies to the population, then q {\displaystyle q} can still be calculated from f(aa):

q = f ( aa ) {\displaystyle q={\sqrt {f({\text{aa}})}}}

and p {\displaystyle p} can be calculated from q {\displaystyle q} . And thus an estimate of f(AA) and f(Aa) derived from p 2 {\displaystyle p^{2}} and 2 p q {\displaystyle 2pq} respectively. Note however, such a population cannot be tested for equilibrium using the significance tests below because it is assumed a priori.

Significance tests for deviation

Testing deviation from the HWP is generally performed using Pearson's chi-squared test, using the observed genotype frequencies obtained from the data and the expected genotype frequencies obtained using the HWP. For systems where there are large numbers of alleles, this may result in data with many empty possible genotypes and low genotype counts, because there are often not enough individuals present in the sample to adequately represent all genotype classes. If this is the case, then the asymptotic assumption of the chi-squared distribution, will no longer hold, and it may be necessary to use a form of Fisher's exact test, which requires a computer to solve. More recently a number of MCMC methods of testing for deviations from HWP have been proposed (Guo & Thompson, 1992; Wigginton et al. 2005)

Example χ test for deviation

This data is from E.B. Ford (1971) on the Scarlet tiger moth, for which the phenotypes of a sample of the population were recorded. Genotype-phenotype distinction is assumed to be negligibly small. The null hypothesis is that the population is in Hardy–Weinberg proportions, and the alternative hypothesis is that the population is not in Hardy–Weinberg proportions.

Table 3: Example Hardy–Weinberg principle calculation
Phenotype White-spotted (AA) Intermediate (Aa) Little spotting (aa) Total
Number 1469 138 5 1612

From which allele frequencies can be calculated:

p = 2 × o b s ( AA ) + o b s ( Aa ) 2 × ( o b s ( AA ) + o b s ( Aa ) + o b s ( aa ) ) = 1469 × 2 + 138 2 × ( 1469 + 138 + 5 ) = 3076 3224 = 0.954 {\displaystyle {\begin{aligned}p&={2\times \mathrm {obs} ({\text{AA}})+\mathrm {obs} ({\text{Aa}}) \over 2\times (\mathrm {obs} ({\text{AA}})+\mathrm {obs} ({\text{Aa}})+\mathrm {obs} ({\text{aa}}))}\\\\&={1469\times 2+138 \over 2\times (1469+138+5)}\\\\&={3076 \over 3224}\\\\&=0.954\end{aligned}}}

and

q = 1 p = 1 0.954 = 0.046 {\displaystyle {\begin{aligned}q&=1-p\\&=1-0.954\\&=0.046\end{aligned}}}

So the Hardy–Weinberg expectation is:

E x p ( AA ) = p 2 n = 0.954 2 × 1612 = 1467.4 E x p ( Aa ) = 2 p q n = 2 × 0.954 × 0.046 × 1612 = 141.2 E x p ( aa ) = q 2 n = 0.046 2 × 1612 = 3.4 {\displaystyle {\begin{aligned}\mathrm {Exp} ({\text{AA}})&=p^{2}n=0.954^{2}\times 1612=1467.4\\\mathrm {Exp} ({\text{Aa}})&=2pqn=2\times 0.954\times 0.046\times 1612=141.2\\\mathrm {Exp} ({\text{aa}})&=q^{2}n=0.046^{2}\times 1612=3.4\end{aligned}}}

Pearson's chi-squared test states:

χ 2 = ( O E ) 2 E = ( 1469 1467.4 ) 2 1467.4 + ( 138 141.2 ) 2 141.2 + ( 5 3.4 ) 2 3.4 = 0.001 + 0.073 + 0.756 = 0.83 {\displaystyle {\begin{aligned}\chi ^{2}&=\sum {(O-E)^{2} \over E}\\&={(1469-1467.4)^{2} \over 1467.4}+{(138-141.2)^{2} \over 141.2}+{(5-3.4)^{2} \over 3.4}\\&=0.001+0.073+0.756\\&=0.83\end{aligned}}}

There is 1 degree of freedom (degrees of freedom for test for Hardy–Weinberg proportions are # genotypes − # alleles). The 5% significance level for 1 degree of freedom is 3.84, and since the χ value is less than this, the null hypothesis that the population is in Hardy–Weinberg frequencies is not rejected.

Fisher's exact test (probability test)

Fisher's exact test can be applied to testing for Hardy–Weinberg proportions. Since the test is conditional on the allele frequencies, p and q, the problem can be viewed as testing for the proper number of heterozygotes. In this way, the hypothesis of Hardy–Weinberg proportions is rejected if the number of heterozygotes is too large or too small. The conditional probabilities for the heterozygote, given the allele frequencies are given in Emigh (1980) as

prob [ n 12 | n 1 ] = ( n n 11 , n 12 , n 22 ) ( 2 n n 1 , n 2 ) 2 n 12 , {\displaystyle \operatorname {prob} ={\frac {n \choose n_{11},n_{12},n_{22}}{2n \choose n_{1},n_{2}}}2^{n_{12}},}

where n11, n12, n22 are the observed numbers of the three genotypes, AA, Aa, and aa, respectively, and n1 is the number of A alleles, where n 1 = 2 n 11 + n 12 {\displaystyle n_{1}=2n_{11}+n_{12}} .

An example Using one of the examples from Emigh (1980), we can consider the case where n = 100, and p = 0.34. The possible observed heterozygotes and their exact significance level is given in Table 4.

Table 4: Example of Fisher's Exact Test for n = 100, p = 0.34.
Number of heterozygotes Significance level
0 0.000
2 0.000
4 0.000
6 0.000
8 0.000
10 0.000
12 0.000
14 0.000
16 0.000
18 0.001
20 0.007
22 0.034
24 0.067
26 0.151
28 0.291
30 0.474
32 0.730
34 1.000

Using this table, one must look up the significance level of the test based on the observed number of heterozygotes. For example, if one observed 20 heterozygotes, the significance level for the test is 0.007. As is typical for Fisher's exact test for small samples, the gradation of significance levels is quite coarse.

However, a table like this has to be created for every experiment, since the tables are dependent on both n and p.

Inbreeding coefficient

The inbreeding coefficient, F (see also F-statistics), is one minus the observed frequency of heterozygotes over that expected from Hardy–Weinberg equilibrium.

F = E ( f ( Aa ) ) O ( f ( Aa ) ) E ( f ( Aa ) ) = 1 O ( f ( Aa ) ) E ( f ( Aa ) ) {\displaystyle F={\frac {\operatorname {E} {(f({\text{Aa}}))}-\operatorname {O} (f({\text{Aa}}))}{\operatorname {E} (f({\text{Aa}}))}}=1-{\frac {\operatorname {O} (f({\text{Aa}}))}{\operatorname {E} (f({\text{Aa}}))}}} ,

where the expected value from Hardy–Weinberg equilibrium is given by

E ( f ( Aa ) ) = 2 p q {\displaystyle \operatorname {E} (f({\text{Aa}}))=2pq}

For example, for Ford's data above;

F = 1 138 141.2 = 0.023 {\displaystyle F=1-{138 \over 141.2}=0.023} .

For two alleles, the chi-squared goodness of fit test for Hardy–Weinberg proportions is equivalent to the test for inbreeding, F = 0.

The inbreeding coefficient is unstable as the expected value approaches zero, and thus not useful for rare and very common alleles. For: E = 0, O > 0, F = −∞ and E = 0, O = 0, F is undefined.

History

Mendelian genetics were rediscovered in 1900. However, it remained somewhat controversial for several years as it was not then known how it could cause continuous characteristics. Udny Yule (1902) argued against Mendelism because he thought that dominant alleles would increase in the population. The American William E. Castle (1903) showed that without selection, the genotype frequencies would remain stable. Karl Pearson (1903) found one equilibrium position with values of p = q = 0.5. Reginald Punnett, unable to counter Yule's point, introduced the problem to G. H. Hardy, a British mathematician, with whom he played cricket. Hardy was a pure mathematician and held applied mathematics in some contempt; his view of biologists' use of mathematics comes across in his 1908 paper where he describes this as "very simple".

To the Editor of Science: I am reluctant to intrude in a discussion concerning matters of which I have no expert knowledge, and I should have expected the very simple point which I wish to make to have been familiar to biologists. However, some remarks of Mr. Udny Yule, to which Mr. R. C. Punnett has called my attention, suggest that it may still be worth making...
Suppose that Aa is a pair of Mendelian characters, A being dominant, and that in any given generation the number of pure dominants (AA), heterozygotes (Aa), and pure recessives (aa) are as p:2q:r. Finally, suppose that the numbers are fairly large, so that mating may be regarded as random, that the sexes are evenly distributed among the three varieties, and that all are equally fertile. A little mathematics of the multiplication-table type is enough to show that in the next generation the numbers will be as (p+q):2(p+q)(q+r):(q+r), or as p1:2q1:r1, say.
The interesting question is — in what circumstances will this distribution be the same as that in the generation before? It is easy to see that the condition for this is q = pr. And since q1 = p1r1, whatever the values of p, q, and r may be, the distribution will in any case continue unchanged after the second generation

The principle was thus known as Hardy's law in the English-speaking world until 1943, when Curt Stern pointed out that it had first been formulated independently in 1908 by the German physician Wilhelm Weinberg. William Castle in 1903 also derived the ratios for the special case of equal allele frequencies, and it is sometimes (but rarely) called the Hardy–Weinberg–Castle Law.

Derivation of Hardy's equations

The derivation of Hardy's equations is illustrative. He begins with a population of genotypes consisting of pure dominants (AA), heterozygotes (Aa), and pure recessives (aa) in the relative proportions p:2q:r with the conditions noted above, that is,

p + 2 q + r = 1 {\displaystyle p+2q+r=1} .

Rewriting this as (p + q) + (q + r) = 1 and squaring both sides yields Hardy’s result:

p 1 + 2 q 1 + r 1 = ( p + q ) 2 + 2 ( p + q ) ( q + r ) + ( q + r ) 2 = 1 {\displaystyle p_{1}+2q_{1}+r_{1}=(p+q)^{2}+2(p+q)(q+r)+(q+r)^{2}=1}

Hardy’s equivalence condition is

E 1 = q 1 2 p 1 r 1 = [ ( p + q ) ( q + r ) ] 2 ( p + q ) 2 ( q + r ) 2 = 0 {\displaystyle {\begin{aligned}E_{1}&={q_{1}^{2}-p_{1}r_{1}}\\&={^{2}-(p+q)^{2}(q+r)^{2}=0}\end{aligned}}}

for generations after the first.

Now to prove Hardy's statement: the successor generation's proportions equal, p : 2 q : r = p 1 : 2 q 1 : r 1 {\displaystyle \textstyle p:2q:r=p_{1}:2q_{1}:r_{1}} , if and only if the equivalence condition q 2 = p r {\displaystyle \textstyle q^{2}=pr} holds. If p : 2 q : r = p 1 : 2 q 1 : r 1 {\displaystyle \textstyle p:2q:r=p_{1}:2q_{1}:r_{1}} , then

q 2 = q 1 2 = [ ( p + q ) ( q + r ) ] 2 = ( p + q ) 2 ( q + r ) 2 = p 1 r 1 = p r {\displaystyle {\begin{aligned}q^{2}&=q_{1}^{2}\\&=\left^{2}\\&=\left(p+q\right)^{2}\left(q+r\right)^{2}\\&=p_{1}r_{1}\\&=pr\end{aligned}}} ,

so the equivalence condition holds. This shows equal proportions imply the equivalence condition. Now suppose the equivalence condition q 2 = p r {\displaystyle \textstyle q^{2}=pr} holds.

p = p ( p + 2 q + r ) = p 2 + 2 p q + p r = p 2 + 2 p q + q 2 = ( p + q ) 2 = p 1 {\displaystyle {\begin{aligned}p&=p\left(p+2q+r\right)\\&=p^{2}+2pq+pr\\&=p^{2}+2pq+q^{2}\\&=\left(p+q\right)^{2}\\&=p_{1}\end{aligned}}}

By symmetry, r = r 1 {\displaystyle r=r_{1}} . Finally,

q = q ( p + 2 q + r ) = q p + 2 q 2 + q r = q 2 + q ( p + r ) + q 2 = q 2 + q ( p + r ) + p r = ( p + q ) ( q + r ) = q 1 {\displaystyle {\begin{aligned}q&=q\left(p+2q+r\right)\\&=qp+2q^{2}+qr\\&=q^{2}+q(p+r)+q^{2}\\&=q^{2}+q(p+r)+pr\\&=\left(p+q\right)\left(q+r\right)\\&=q_{1}\end{aligned}}} .

This shows p = p 1 {\displaystyle \textstyle p=p_{1}} , q = q 1 {\displaystyle \textstyle q=q_{1}} , and r = r 1 {\displaystyle \textstyle r=r_{1}} . Thus, proportions p : 2 q : r = p 1 : 2 q 1 : r 1 {\displaystyle \textstyle p:2q:r=p_{1}:2q_{1}:r_{1}} equal. This shows the equivalence condition implies the proportions equal. Therefore, the proportions equal if and only if the equivalence condition holds, as was claimed.

Since the equivalence condition always holds for the second generation, all succeeding generations have the same proportions.

Numerical example

An example computation of the genotype distribution given by Hardy's original equations is instructive. The phenotype distribution from Table 3 above will be used to compute Hardy's initial genotype distribution. Note that the p and q values used by Hardy are not the same as those used above.

sum = o b s ( AA ) + 2 × o b s ( Aa ) + o b s ( aa ) = 1469 + 2 × 138 + 5 = 1750 {\displaystyle {\begin{aligned}{\text{sum}}&={\mathrm {obs} ({\text{AA}})+2\times \mathrm {obs} ({\text{Aa}})+\mathrm {obs} ({\text{aa}})}={1469+2\times 138+5}\\&=1750\end{aligned}}}
p = 1469 1750 = 0.83943 2 q = 2 × 138 1750 = 0.15771 r = 5 1750 = 0.00286 {\displaystyle {\begin{aligned}p&={1469 \over 1750}=0.83943\\2q&={2\times 138 \over 1750}=0.15771\\r&={5 \over 1750}=0.00286\end{aligned}}}

As checks on the distribution, compute

p + 2 q + r = 0.83943 + 0.15771 + 0.00286 = 1.00000 {\displaystyle p+2q+r=0.83943+0.15771+0.00286=1.00000\,}

and

E 0 = q 2 p r = 0.00382 {\displaystyle E_{0}=q^{2}-pr=0.00382\,}

For the next generation, Hardy's equations give,

q = 0.15771 2 = 0.07886 p 1 = ( p + q ) 2 = 0.84325 2 q 1 = 2 ( p + q ) ( q + r ) = 0.15007 r 1 = ( q + r ) 2 = 0.00668 {\displaystyle {\begin{aligned}q&={0.15771 \over 2}=0.07886\\\\p_{1}&=(p+q)^{2}=0.84325\\2q_{1}&=2(p+q)(q+r)=0.15007\\r_{1}&=(q+r)^{2}=0.00668\end{aligned}}}

Again as checks on the distribution, compute

p 1 + 2 q 1 + r 1 = 0.84325 + 0.15007 + 0.00668 = 1.00000 {\displaystyle p_{1}+2q_{1}+r_{1}=0.84325+0.15007+0.00668=1.00000\,}

and

E 1 = q 1 2 p 1 r 1 = 0.00000 {\displaystyle E_{1}=q_{1}^{2}-p_{1}r_{1}=0.00000\,}

which are the expected values. The reader may demonstrate that subsequent use of the second-generation values for a third generation will yield identical results.

Graphical representation

It is possible to represent the distribution of genotype frequencies for a bi-allelic locus within a population graphically using a de Finetti diagram. This uses a triangular plot (also known as trilinear, triaxial or ternary plot) to represent the distribution of the three genotype frequencies in relation to each other. Although it differs from many other such plots in that the direction of one of the axes has been reversed.

The curved line in the above diagram is the Hardy–Weinberg parabola and represents the state where alleles are in Hardy–Weinberg equilibrium.

It is possible to represent the effects of Natural Selection and its effect on allele frequency on such graphs (e.g. Ineichen & Batschelet 1975)

The de Finetti diagram has been developed and used extensively by A. W. F. Edwards in his book Foundations of Mathematical Genetics.

References

  1. The term frequency usually refers to a number or count, but in this context, it is synonymous with probability.
  2. http://www.mun.ca/biology/scarr/2900_HW_for_dioecious.html
  3. Hartl DL, Clarke AG (2007) Principles of population genetics. Sunderland, MA: Sinauer
  4. ^ Emigh, Ted H. (1980). "A Comparison of Tests for Hardy–Weinberg Equilibrium". Biometrics. 36 (4): 627–642. doi:10.2307/2556115. JSTOR 2556115.
  5. Crow, James F. (1999). "Hardy, Weinberg and language impediments". Genetics. 152 (3): 821–825. PMC 1460671. PMID 10388804.
  6. Stern, Curt (1962). "Wilhelm Weinberg". Genetics. 47: 1–5.

Bibliography

External links

Population genetics
Key concepts
Selection
Effects of selection
on genomic variation
Genetic drift
Founders
Related topics

Categories: