Het Wright-Fisher-model en genetische drift
De "Simulator voor genetische drift in populaties" is een interactief hulpmiddel voor het simuleren en visualiseren van neutrale genetische drift met behulp van het diploïde Wright–Fisher-model. Binnen de populatiegenetica beschrijft dit model hoe de frequentie van een neutraal allel in een eindige, voortplantende populatie door louter toeval fluctueert van generatie op generatie. Omdat de simulator uitgaat van selectieve neutraliteit, hebben de gemodelleerde allelen geen invloed op de overlevings- of voortplantingskansen van de individuen.
Bij elke generatie trekt de Wright–Fisher-transitie het volgende aantal allelen en deelt dit door de beschikbare genkopieën:
Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne).
Hierin staat Ne voor de effectieve populatiegrootte (ingevoerd in "voortplantende individuen") en pₜ voor de allelfrequentie in generatie t. Omdat elk diploïd individu twee genkopieën draagt, is het totale aantal getrokken allelkopieën per generatie gelijk aan 2Ne. Vóór de steekproef beschrijven deze momenten het bereik van mogelijke volgende generaties:
E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne).
Door deze stochastische overgangen herhaaldelijk te berekenen over opeenvolgende generaties, toont de simulator hoe willekeurige steekproeven leiden tot de fixatie of het verlies van genetische variatie.
Effectieve versus feitelijke populatiegrootte
Bij het configureren van de simulatie voert de gebruiker de parameter "effectieve populatiegrootte" (Ne) in, in plaats van de feitelijke, demografische populatiegrootte (de censusgrootte). De effectieve populatiegrootte is een theoretische waarde die de omvang vertegenwoordigt van een geïdealiseerde populatie die onderhevig is aan dezelfde mate van genetische drift als de werkelijk bestudeerde populatie.
In natuurlijke populaties wijkt de effectieve populatiegrootte vrijwel altijd af van de feitelijke populatiegrootte. Factoren zoals een scheve geslachtsverhouding, sterke variatie in het individuele voortplantingssucces, historische populatieknelpunten (bottlenecks) en fluctuaties in populatiegrootte over de tijd zorgen ervoor dat de genetisch effectieve populatie aanzienlijk kleiner is dan het totale aantal aanwezige individuen. Door Ne als invoerparameter te gebruiken, sluit de simulator direct aan op de wiskundige wetmatigheden van de populatiegenetica.
Dynamiek van fixatie en verlies
Een belangrijk kenmerk van genetische drift in een gesloten populatie zonder mutatie of migratie is dat alleltrajecten uiteindelijk eindigen in een van de twee absorberende grenzen: fixatie (p = 1) of verlies (p = 0). Zodra een allel een van deze grenzen bereikt, kan de frequentie ervan niet meer veranderen, tenzij er nieuwe genetische variatie wordt geïntroduceerd.
In een strikt neutraal model is de kans dat een specifiek allel uiteindelijk wordt gefixeerd exact gelijk aan de startfrequentie (p₀). Een allel met een startfrequentie van 0,2 heeft dus een kans van 20% om uiteindelijk te fixeren en een kans van 80% om verloren te gaan. Omdat simulaties in de praktijk echter over een eindig aantal generaties (T) lopen, zullen niet alle trajecten binnen de looptijd de absorberende grenzen bereiken. De simulator rapporteert daarom na elke run de exacte statistieken voor de gesimuleerde herhalingen.
Heterozygositeit en het verlies van diversiteit
De genetische diversiteit binnen een populatie kan worden uitgedrukt in de verwachte heterozygositeit (H), het aandeel individuen dat naar verwachting heterozygoot is voor de locus. Het verwachte aandeel ongelijke allelenparen neemt over t generaties als volgt af:
H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ.
Deze formule laat zien dat de snelheid waarmee de heterozygositeit afneemt (de mate van coalescentie of genetische verarming) omgekeerd evenredig is met de effectieve populatiegrootte. In kleine populaties verloopt dit verval snel, terwijl grotere populaties hun genetische variatie veel langer behouden.
Er bestaat hierbij een schijnbare paradox: wanneer men het gemiddelde neemt van de allelfrequenties over een groot aantal onafhankelijke populaties, blijft de gemiddelde uiteindelijke frequentie stabiel rond de startfrequentie p₀. Dit komt doordat drift geen inherente richting of voorkeur heeft. Echter, binnen elke individuele populatie daalt de heterozygositeit onherroepelijk naarmate afzonderlijke trajecten naar de grenzen 0 of 1 driften. De gemiddelde heterozygositeit over alle populaties neemt dus af, zelfs wanneer de gemiddelde allelfrequentie constant blijft.
Methode, aannames en limieten
Deze pagina implementeert het neutrale diploïde Wright–Fisher-model zoals beschreven door Nature Education en de Stanford Human Genetics-tekst. Elke nieuwe generatie is een willekeurige steekproef van 2Ne allelkopieën uit de vorige. Bronnen gecontroleerd op 6 augustus 2026.
Modelveronderstellingen
Het model is gebaseerd op de volgende strikte aannames:
- Constante effectieve populatiegrootte over alle generaties.
- Discrete, niet-overlappende generaties.
- Willekeurige en onafhankelijke steekproeftrekking van gameten.
- Volledige selectieve neutraliteit van de bestudeerde allelen.
- Eén niet-gekoppelde locus met maximaal twee allelen.
- Volledige afwezigheid van mutatie en migratie (genenstroming).
Gebruik dit model niet als voorspelling wanneer selectie, populatieverandering, onevenwichtigheid in geslacht of voortplantingssucces, overlappende generaties, inteelt, koppeling, genenstroming of populatiestructuur een rol spelen. Een technische veiligheidsfactor is niet van toepassing op dit stochastische biologische model.
Numerieke precisie en workload
Elke binomiale trekking maakt gebruik van exacte Bernoulli-steekproeven; er wordt geen normale benadering gebruikt. Berekeningen behouden de double-precision van JavaScript. Frequenties liggen op stappen van 1 ÷ (2Ne); de weergave rondt af op maximaal zes decimalen.
Om de prestaties van de browser te waarborgen, geldt er een strikte workloadlimiet van 5.000.000 geplande allel-kopietrekkingen, berekend als 2Ne × generaties × herhalingen. Als een ingestelde configuratie deze limiet overschrijdt, blokkeert de simulator de uitvoering met de foutmelding: "Deze opzet plant ‹draws› allel-kopietrekkingen, wat boven de browserlimiet van ‹limit› ligt. Verklein de populatiegrootte, generaties of herhalingen."
Privacy en gegevensverwerking
Je parameters en gesimuleerde populaties blijven in deze browser. Alle berekeningen, gegevensgeneratie en grafische weergaven vinden lokaal plaats op het apparaat van de gebruiker. Er worden geen invoergegevens, gegenereerde trajecten of unieke seeds geüpload naar externe servers.
Veelgestelde vragen (FAQ)
Moet ik de feitelijke populatiegrootte of de effectieve populatiegrootte invoeren?
Voer de effectieve populatiegrootte, Ne, in: de grootte van een geïdealiseerde voortplantingspopulatie die in hetzelfde tempo zou driften als de populatie die u wilt bestuderen. De feitelijke populatiegrootte kan veel groter zijn wanneer de geslachtsverhoudingen, het voortplantingssucces, populatieschommelingen of de populatiestructuur ongelijkmatig zijn.
Is de neutrale kans op uiteindelijke fixatie gelijk aan de startfrequentie?
Ja, in dit neutrale model zonder mutatie of migratie heeft een allel dat begint met frequentie p₀ een uiteindelijke fixatiekans van p₀ en een verlieskans van 1 − p₀. De hierboven getoonde percentages hebben alleen betrekking op het ingevoerde aantal generaties en herhalingen, dus ze variëren met de seed.
Waarom kan de gemiddelde allelfrequentie rond p₀ blijven terwijl de heterozygositeit daalt?
Neutrale drift heeft geen voorkeursrichting over veel onafhankelijke populaties, waardoor de gemiddelde frequentie naar verwachting p₀ blijft. Individuele populaties verspreiden zich echter richting 0 of 1, en beide grenzen hebben een heterozygositeit van nul. Variatie gaat dus verloren, zelfs terwijl het gemiddelde over de populaties vrijwel onveranderd blijft.
Kan dit voorspellen wat een echte populatie zal doen?
Nee. Het toont uitkomsten onder een bewust geïdealiseerd neutraal Wright–Fisher-model. Echte voorspellingen vereisen verdedigbare schattingen van de effectieve grootte en van eventuele selectie, mutatie, migratie, veranderende populatiegrootte, overlappende generaties, koppeling en populatiestructuur die voor de casus van belang zijn.