Genetic Insights Into Wild Hybrid House Mice and Speciation

genomic ancestry analysis in wild hybrid house n.w
1 / 61
Embed
Share

Explore genomic ancestry analysis in wild hybrid house mice conducted by Ph.D. student Megan Frayer, delving into the genetics of speciation and the house mouse hybrid zone's role in understanding subspecies differentiation. Discover how this research sheds light on the ongoing process of speciation between Mus musculus domesticus and Mus musculus musculus.

  • Genetics
  • Speciation
  • Hybrid Zone
  • Genomic Analysis
  • House Mice

Uploaded on | 0 Views


Download Presentation

Please find below an Image/Link to download the presentation.

The content on the website is provided AS IS for your information and personal use only. It may not be sold, licensed, or shared on other websites without obtaining consent from the author. If you encounter any issues during the download, it is possible that the publisher has removed the file from their server.

You are allowed to download the files provided on this website for personal or commercial use, subject to the condition that they are used lawfully. All files are the property of their respective owners.

The content on the website is provided AS IS for your information and personal use only. It may not be sold, licensed, or shared on other websites without obtaining consent from the author.

E N D

Presentation Transcript


  1. Genomic Ancestry Analysis in Wild Hybrid House Mice Megan Frayer Ph.D. Student, Laboratory of Genetics UW-Madison HTCondor Week 2022

  2. Genetics of Speciation

  3. Genetics of Speciation

  4. Genetics of Speciation

  5. Genetics of Speciation

  6. Genetics of Speciation

  7. Genetics of Speciation

  8. The house mouse hybrid zone can tell us about how speciation is proceeding between these subspecies M. m. domesticus M. m. musculus

  9. The house mouse hybrid zone can tell us about how speciation is proceeding between these subspecies M. m. domesticus M. m. musculus

  10. The house mouse hybrid zone can tell us about how speciation is proceeding between these subspecies M. m. domesticus M. m. musculus

  11. The house mouse hybrid zone can tell us about how speciation is proceeding between these subspecies M. m. domesticus M. m. musculus

  12. The house mouse hybrid zone can tell us about how speciation is proceeding between these subspecies M. m. domesticus M. m. musculus

  13. The house mouse hybrid zone can tell us about how speciation is proceeding between these subspecies M. m. domesticus M. m. musculus

  14. ATCGTCAGTCAGTCGATCGATACGTAGCATGCAGTACGATGCAGTACGATGATACG TAGCAGTCAGACACGTAGCTATGCATCGTACGTCATGCTACGTCATGCTACTATGC

  15. ATCGTCAGTCAGTCGATCGATACGTAGCATGCAGTACGATGCAGTACGATGATACG TAGCAGTCAGACACGTAGCTATGCATCGTACGTCATGCTACGTCATGCTACTATGC

  16. ATCGTCAGTCAGTCGATCGATACGTAGCATGCAGTACGATGCAGTACGATGATACG TAGCAGTCAGACACGTAGCTATGCATCGTACGTCATGCTACGTCATGCTACTATGC

  17. D O M H E T M U S Centromere Telomere

  18. Parameter grid search

  19. Parameter grid search What is the combination of input parameters with the highest likelihood?

  20. Parameter grid search Parameter defaultRate timeSince Admixture ancestryProp1 ancestralRate1 ancestralRate2 mutation1 mutation2 miscopyRate Miscopy Mutation Values to be tested 0.8 0.86 0.99 1.15 1000 0.4 41000 14000 1E-04 3.4E-05 0.01 0.01 3750 0.5 69250 23650 1E-05 3.4E-06 0.001 0.001 6500 0.6 97500 33290 1E-06 3.4E-07 1E-04 1E-04 9250 12000 14750 20815 1E-07 3.4E-08 1E-05 1E-05 35158 1E-08 3.4E-09 1E-06 1E-06 49500 5.1E-05 5.1E-06 5.1E-07 5.1E-08 5.1E-09

  21. Parameter grid search Parameter defaultRate timeSince Admixture ancestryProp1 ancestralRate1 ancestralRate2 mutation1 mutation2 miscopyRate Miscopy Mutation Values to be tested 0.8 0.86 0.99 1.15 1000 0.4 41000 14000 1E-04 3.4E-05 0.01 0.01 3750 0.5 69250 23650 1E-05 3.4E-06 0.001 0.001 6500 0.6 97500 33290 1E-06 3.4E-07 1E-04 1E-04 9250 12000 14750 20815 1E-07 3.4E-08 1E-05 1E-05 35158 1E-08 3.4E-09 1E-06 1E-06 49500 5.1E-05 5.1E-06 5.1E-07 5.1E-08 5.1E-09 108,000 combinations of parameters to be tested

  22. Parameter grid search Parameter defaultRate timeSince Admixture ancestryProp1 ancestralRate1 ancestralRate2 mutation1 mutation2 miscopyRate Miscopy Mutation Values to be tested 0.8 0.86 0.99 1.15 1000 0.4 41000 14000 1E-04 3.4E-05 0.01 0.01 3750 0.5 69250 23650 1E-05 3.4E-06 0.001 0.001 6500 0.6 97500 33290 1E-06 3.4E-07 1E-04 1E-04 9250 12000 14750 20815 1E-07 3.4E-08 1E-05 1E-05 35158 1E-08 3.4E-09 1E-06 1E-06 49500 5.1E-05 5.1E-06 5.1E-07 5.1E-08 5.1E-09 108,000 combinations of parameters to be tested

  23. Parameter grid search Create input files Run parameter tests Compile and analyze results

  24. Parameter grid search Create input files Run parameter tests Compile and analyze results

  25. parameter_test.dag Create Input Files

  26. parameter_test.dag Examples of files to print: Create Input Files Submit files Executables Input for programs being run Scripts that will need to be run

  27. parameter_test.dag Create Input Files Parameter Test 1

  28. parameter_test.dag Create Input Files Parameter Test 1 Parameter Test 2

  29. parameter_test.dag Create Input Files Parameter Test 1 Parameter Test 2 Parameter Test 3

  30. parameter_test.dag Create Input Files Parameter Test 1 Parameter Test 2 Parameter Test 3 Parameter Test n

  31. parameter_test.dag Create Input Files Parameter Test 1 Parameter Test 2 Parameter Test 3 Parameter Test n Compile results/create summaries

  32. parameter_test.dag Create Input Files SUBDAG_EXTERNAL Parameter Test 1 Parameter Test 2 Parameter Test 3 Parameter Test n Compile results/create summaries

  33. parameter_test.dag Create Input Files SUBDAG_EXTERNAL Parameter Test 1 Parameter Test 2 Parameter Test 3 Parameter Test n Before HTC: 2 hours/test 24.6 years/108,000 tests With HTC: 2 hours/test 10 days/108,000 tests 24.6 years Compile results/create summaries 10 days

  34. Testing with Simulated Chromosomes

  35. Testing with Simulated Chromosomes How well is the program performing?

  36. Testing with Simulated Chromosomes Simulate Chromosomes Determine the true ancestry map Infer ancestry using the method to be tested Compare the true and inferred maps

  37. Testing with Simulated Chromosomes Simulate Chromosomes Determine the true ancestry map Infer ancestry using the method to be tested Compare the true and inferred maps

  38. inference_testing.dag Create Input Files Parameter Set 3 Set 1 Set 1 Set 1 Set 2 Set 1 Set 2 Set 3 Set 1 Set 3 Set n Set 1 Set n Inference Test Inference Test Inference Test Inference Test Compile results/create summaries

  39. inference_testing.dag Create Input Files

  40. inference_testing.dag Create Input Files Inference Test Set 1

  41. inference_testing.dag Create Input Files Inference Test Set 1 Set 1 Inference Test

  42. inference_testing.dag Create Input Files Inference Test Set 1 Set 1 Set 1 Inference Test

  43. inference_testing.dag Create Input Files Set 1 Set 1 Set 1 Set 2 Set 1 Set 2 Inference Test Inference Test

  44. inference_testing.dag Create Input Files Parameter Set 3 Set 1 Set 1 Set 1 Set 2 Set 1 Set 2 Set 3 Set 1 Set 3 Inference Test Inference Test Inference Test

  45. inference_testing.dag Create Input Files Parameter Set 3 Set 1 Set 1 Set 1 Set 2 Set 1 Set 2 Set 3 Set 1 Set 3 Set n Set 1 Set n Inference Test Inference Test Inference Test Inference Test

  46. inference_testing.dag Create Input Files Parameter Set 3 Set 1 Set 1 Set 1 Set 2 Set 1 Set 2 Set 3 Set 1 Set 3 Set n Set 1 Set n Inference Test Inference Test Inference Test Inference Test Compile results/create summaries

  47. inference_testing.dag Create Input Files Parameter Set 3 Set 1 Set 1 Set 1 Set 2 Set 1 Set 2 Set 3 Set 1 Set 3 Set n Set 1 Set n Inference Test Inference Test Inference Test Inference Test Before HTC: 3 hours/test 6.25 days/50 tests With HTC: 3 hours/test 10 hours/50 tests 6.25 days 10 hours Compile results/create summaries

  48. Simulations

  49. Simulations Simulate data and run a script to make a summary

  50. Variables Template Submit Files simulation.dag Replicate 1 Replicate 2 Replicate 3 Replicate n Simulation.config DAGMAN_MAX_JOBS_IDLE = 1000

More Related Content