in silico
Phase 1 : In Silico Validation: SLiM Simulation of Bottleneck, Drift, and Introgression
- Objectif
- To determine whether the proposed mechanism (stochastic fixation via drift in Ne<50 followed by parasexual introgression) is plausible under realistic parameter ranges, and to falsify the hypothesis if simulations show that selection (s<-0.01) or migration (m<0.001) dominate drift under all tested conditions.
- Coût estimé
- €800-1,500 (cloud compute time + software licenses, if any)
- Durée estimée
- 4-6 weeks (2 weeks coding + 2 weeks simulation runs + 2 weeks analysis)
- Critères de réussite
- Fixation probability for Ne=30, initial freq=0.01, s=0 · ≥ 0.15 (i.e., at least 15 out of 100 replicates show fixation) · (Count of replicates where allele frequency reaches 1.0 in founder population within 50 generations)
- Introgression frequency in progenitor at generation 100 for m=0.05 · ≥ 0.05 (i.e., allele frequency in progenitor ≥ 5%) · (Mean allele frequency across 100 replicates at generation 100)
- ASTRAL concordance factor for resistance-linked loci · < 0.6 (i.e., significant discordance) · (Mean concordance factor across 100 replicates for the 100-kb window containing the ERG11 locus)
- On continue si
- All three success criteria are met (fixation probability ≥0.15, introgression frequency ≥0.05, concordance factor <0.6).
- On arrête si
- Fixation probability <0.05 for Ne=30, s=0 (i.e., drift alone cannot fix the allele even without selection), OR introgression frequency <0.01 for m=0.05 (i.e., parasexual mating is too rare to cause detectable introgression).
- On réoriente si
- Fixation probability is ≥0.15 but introgression frequency is between 0.01 and 0.05: proceed to Phase 2 but increase the number of replicates or extend the experimental timeline to 150 generations.
- Risques
- SLiM simulations may not accurately model parasexual mating (e.g., chromosome loss rate, mating type switching frequency).Probabilité : mediumParade : Use published parameter ranges from DOI: 10.1371/journal.pgen.1003665 and DOI: 10.1128/EC.00351-07. Perform sensitivity analysis by varying mating rate by 10x.
- Computational time for 100 replicates × 200 generations × 10 parameter combinations may exceed cluster limits.Probabilité : mediumParade : Use a factorial design with 3 levels per variable (Ne: 10, 30, 100; initial freq: 0.001, 0.01, 0.05; s: -0.05, 0, 0.05; m: 0.001, 0.01, 0.05) = 81 combinations × 100 replicates = 8,100 runs. Each run takes ~5 minutes on a single core. Use 16 cores for 42 hours total. Pre-allocate cloud compute budget.
- ASTRAL concordance factor may be uninformative if gene tree resolution is low (short branches).Probabilité : lowParade : Use 100-kb windows to ensure sufficient phylogenetic signal. Also compute D-statistic (ABBA-BABA) as a complementary test.