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:

Screenshot 2024-12-20 at 11 30 20 AM

Screenshot 2024-12-20 at 11 30 33 AM

Screenshot 2024-12-20 at 11 30 42 AM

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.

Screenshot 2024-12-20 at 11 30 54 AM

Figure 4. Genotype pie chart for the WA population

Going forward

  1. Finish genotyping WA samples
  2. Troubleshoot QC gel discrepancies by rerunning MA CC samples
  3. Redo all Chelex samples?
  4. Identify samples for confirmation sequencing
  5. Determine methods for comparing population responses
  6. Troubleshoot lipid assay protocol
  7. Conduct lipid assay for crabs of interest
Written on December 20, 2024