1 / 16

Lab 3 : Exact tests and Measuring of Genetic Variation

Lab 3 : Exact tests and Measuring of Genetic Variation. χ 2 - test. Where,. Chi-square Assumptions : Finite # of observations. Observations are independent. Samples collected randomly. Large sample size (>20; >50).

jagger
Download Presentation

Lab 3 : Exact tests and Measuring of Genetic Variation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Lab 3 : Exact tests and Measuring of Genetic Variation

  2. χ2 - test Where,

  3. Chi-square Assumptions : • Finite # of observations. • Observations are independent. • Samples collected randomly. • Large sample size (>20; >50)

  4. Example: Suppose you caught 5 Bluegill fish and detected two alleles (A1 and A2) and observed that all 5 fish were A1A2 heterozygotes. Calculate allele frequencies and do a χ2 – test to determine whether the population is in HWE. Conclusion: Reject Ho at α = 0.05 , because calculated χ2-value (=5) is more than critical χ2- value with 1 d.f. (≈ 3.84) i.e. Bluegill population is not in HWE.

  5. Why is the previous conclusion is not reliable? Because, it violates the assumption of large sample size. As a rule of thumb, the Chi-square test should not be used when the expected number for any genotype class is less than 5.

  6. Exact Test • Calculate the probability of observing N11=0, N12=5, N22=0 under HWE using the multinomial probability equation . 2. Generate all possible permutations of 5 A1 alleles and 5 A2 alleles into 3 genotypes i.e. 10! =3,628,800. 3. Calculate probability of observing each of these samples under HWE using multinomial probability equation. 4. Determine proportion of samples, whose probability is ≤ 0.0313. 5. If proportion (p-value) is less than 0.05, then reject Ho at α = 0.05.

  7. p – value = Sample with probability ≤ 0.0313 / Total # of sample = 3/30 = 0.10 Conclusion: Because p- value is more than 0.05, hence we can't reject Ho i.e. The bluegill population is in HWE at α = 0.05.

  8. Generation of all possible samples and calculation of probability for each sample is computationally intensive. It will require too much time and is practically impossible for large samples. In practice, exact tests are done by sampling a distribution generated from a Markov Chain (beyond the scope of this course).

  9. Measures of Genetic Variation • Heterozygosity (Gene diversity). • Number of alleles (Allele diversity). • Effective number of alleles. • Percentage of polymorphic loci.

  10. 1. Heterozygosity (Gene diversity) • - Most commonly used measure of genetic variation. • Observed heterozygosity (HO) = Proportion of heterozygotes in a sample. • Expected heterozygosity(HE) = Heterozygosity expected under HWE. = Expected homozygosity under HWE = p12 + p22 + P32 + …….+ pn2 For small sample size(< 50), unbiased HEcan be calculated by :

  11. 2. Number of Alleles(n): - Number of alleles present at a locus in a population. - Also called allele diversity (A). - Strongly influenced by sample size. 3. Effective number of Alleles(ne): The number of alleles a population would have if all alleles were at equal frequency 4. Proportion of polymorphic loci (P): - Not so useful for highly variable loci like Microsatellites. - Locus selection bias

  12. GenAlEx

  13. GenAlEx

  14. Problem 1. Use GenAlEx to perform the following analyses based on the human SSR data: • Calculate the genetic variation measures HO, HE, Na, and Ne for all loci in all populations. Include the estimated values of these measures for all loci in a population you will be assigned during the lab. What can you conclude about the allele frequencies of the 10 loci by comparing Na to Ne? • Calculate the average HO and HE across loci for your assigned population. Can you predict anything about the test of HWE based on these values? • c.GRADUATE STUDENTS ONLY: For your assigned population, calculate the average HO and HE weighted by the number of individuals for which data were available. • d.Perform a Chi-square test of HWE for all loci in all populations. Include a summary of the test for your assigned population in the lab report. How do you interpret the results of this test?

  15. Problem 2. Perform an exact test of HWE for all loci and all populations using Arlequin. Include the results for your assigned population in the lab report. How do these results compare to those from the Chi-square test and why? Which test do you trust more? • Go to the data worksheet and select the Export Data submenu of GenAlEx. • Select Arlequin and click OK on the export parameters window. • Save the file with an extension “.arp” in the class data folder. You could, for example, name your file human_ssr.arp

More Related