Cold Acclimation Green Crab Experiment Part 30
Hardy-Weinberg equilibrium for MA and WA populations
Alright, time to determine if my populations are in Hardy-Weinberg equilibrium! Generally, you’d calculate the proportion of the two alleles in each population:
p = (# CC * 2) + (# CT) / (# crabs * 2) q = (# TT * 2) + (# CT) / (# crabs * 2)
For sanity check, p + q should equal 1! Then, you’d use those allele frequencies to calculate the expected proportion of each genotype you’d see in the population:
Expected proportion of CC = p^2 Expected proportion of CT = 2pq Expected proportion of TT = q^2
Finally, you’d compare the expected proportions versus observed proportions using a Chi-squared test. If a population is in Hardy-Weinberg equilibrium, the observed genotype proportions would not be significantly different than the expected genotype proportions.
Of course online calculators now exist for this, so instead of doing everything by hand, I plugged in the genotype counts into this website.
Table 1. Hardy-Weinberg equilibrium P-values
Population | P-value |
---|---|
MA | 0.0000071191801649826 |
New Bedford | 0.02291149613108 |
Eel Pond | 6.6098836727653e-7 |
WA | 0.00169236548011 |
So…none of my populations are in Hardy-Weinberg equilibrium. With the MA populations, it’s interesting because there’s such a low instance of TT crabs any way you slice the population:
Figures 1-3. Genotype pie charts for the full MA population, New Bedford, and Eel Pond
With the WA population, I have a low instance of CC crabs, but I haven’t finished genotyping so the jury is still out. If I remember correctly, Julia’s crabs were in Hardy-Weinberg equilibrium, but she also picked her crabs based on genotype to have a balanced experimental design.
Figure 4. Genotype pie chart for the WA population
Going forward
- Finish genotyping WA samples
- Troubleshoot QC gel discrepancies by rerunning MA CC samples
- Redo all Chelex samples?
- Identify samples for confirmation sequencing
- Determine methods for comparing population responses
- Troubleshoot lipid assay protocol
- Conduct lipid assay for crabs of interest