in silico
Phase 1: In Silico Validation: Power Analysis, Simulation, and Data Mining
- Objective
- To computationally confirm that the proposed experimental design has sufficient statistical power, that the DGRP panel contains relevant genetic variation for period modulation, and that the introgression scheme is feasible without excessive linkage drag.
- Estimated cost
- €500-1,500
- Estimated duration
- 4-6 weeks
- Success criteria
- Statistical power for Phase 2 period shift detection · ≥80% power to detect Δτ = 0.5 h at α = 0.01/40 · (Monte Carlo simulation output (proportion of significant tests))
- Genetic diversity in DGRP panel for clock modifiers · At least 5 clock modifier genes with MAF > 5% in the DGRP panel · (MAF calculation from DGRP SNP data)
- Introgression feasibility · Simulated residual heterozygosity < 0.1% in background after 10 generations · (Simulation output from AlphaSimR)
- GWAS detection limit · Minimum detectable QTL effect size ≤ 0.75 h (i.e., within the predicted 0.5–1.5 h range) · (genpwr output for linear mixed model)
- Go if
- All four success criteria are met. Specifically: power ≥80%, ≥5 modifier genes with MAF>5%, residual heterozygosity <0.1%, detectable QTL effect ≤0.75 h.
- No-go if
- Any criterion fails: power <70%, <3 modifier genes with MAF>5%, residual heterozygosity >0.5%, detectable QTL effect >1.0 h.
- Pivot if
- If power is borderline (70-80%) but genetic diversity is high, increase sample size to n=24 flies/NIL in Phase 2. If introgression simulation shows >0.5% residual heterozygosity, extend backcrossing to 12 generations.
- Risks
- DGRP panel may lack sufficient genetic variation at clock modifier loci (low MAF).Probability: mediumMitigation: Use complementary wild-type strains from the DSPR (Drosophila Synthetic Population Resource) or inbred lines from the Drosophila Genetic Reference Panel (DGRP) plus additional wild-caught lines. Phase 1 will identify the most diverse strains.
- Simulated power may be overestimated due to unmodeled variance (e.g., batch effects, sex differences).Probability: mediumMitigation: Use conservative variance estimate (σ² = 0.3 h²) in simulations. Include random batch effect in power model.
- Linkage drag around the Per locus may be larger than 500 kb, confounding background effects.Probability: lowMitigation: Simulate with 1 Mb flanking region. If drag persists, use CRISPR/Cas9 to knock in PerS directly into each background (alternative, higher cost).