Das Wright–Fisher-Modell der neutralen Gendrift
Der Simulator für Gendrift in der Populationsgenetik basiert auf dem klassischen, neutralen diploiden Wright–Fisher-Modell. Dieses mathematische Fundament beschreibt, wie Genvarianten (Allele) in einer endlichen, idealisierten Fortpflanzungsgemeinschaft allein durch den Zufall der Gametenstichprobenziehung von Generation zu Generation schwanken. Da keine Selektion auf das untersuchte Allel wirkt, ist der Prozess rein stochastisch.
Bei jeder Generation zieht der Wright–Fisher-Übergang die nächste Allelanzahl als Stichprobe und teilt sie durch die verfügbaren Genkopien:
Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne)
Vor der Stichprobenziehung beschreiben diese Momente den Bereich möglicher nächster Generationen:
E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne)
Jede binomiale Ziehung nutzt exakte Bernoulli-Stichproben; es wird keine Normalverteilungsapproximation verwendet. Berechnungen behalten die JavaScript-Double-Precision bei. Häufigkeiten liegen auf Schritten von 1 ÷ (2Ne); die Anzeige rundet auf maximal sechs Dezimalstellen.
Effektive versus tatsächliche Populationsgröße
Für die Simulation ist die Eingabe der effektiven Populationsgröße, Ne (gemessen in "fortpflanzungsfähige Individuen") entscheidend. Diese unterscheidet sich oft stark von der tatsächlichen Anzahl lebender Individuen (Zensusgröße) einer realen Population.
Die effektive Populationsgröße beschreibt die Anzahl an Individuen einer idealisierten Population, die denselben Grad an genetischer Drift aufweist wie die untersuchte reale Population. In der Natur senken ungleiche Geschlechterverhältnisse, extreme Schwankungen der Populationsgröße über Generationen hinweg sowie eine hohe Varianz im individuellen Fortpflanzungserfolg die effektive Größe weit unter die Zensusgröße. Im Simulator muss eine ganze Zahl von 1 bis 100.000 eingegeben werden.
Dynamik von Fixierung und Verlust
Ein zentrales Merkmal der neutralen Gendrift ist, dass Allele im Laufe der Zeit entweder dauerhaft in der Population verankert werden oder vollständig verschwinden. Sobald eine Flugbahn eine der beiden absorbierenden Grenzen erreicht – die vollständige Fixierung bei p = 1 oder den vollständigen Verlust bei p = 0 –, kann sie diesen Zustand ohne Mutation oder Migration nicht mehr verlassen.
Die mathematische Wahrscheinlichkeit, dass ein neutrales Allel schlussendlich fixiert wird, entspricht exakt seiner Ausgangshäufigkeit p₀. Ein seltenes Allel mit einer geringen Ausgangshäufigkeit geht daher mit hoher Wahrscheinlichkeit verloren, während ein bereits weit verbreitetes Allel eine entsprechend höhere Fixierungswahrscheinlichkeit besitzt. Der Simulator erfasst diese Dynamiken über die eingegebenen Generationen hinweg und gibt die Raten unter den Bezeichnungen "Fixiert bei p = 1", "Verloren bei p = 0" sowie "Nicht absorbiert bei T" aus. Falls alle simulierten Verläufe vor dem Ende der Generationen fixiert oder verloren gehen, wird die "Mittlere Absorptionsgeneration" berechnet; andernfalls zeigt das System "nicht erreicht" an.
Zerfall der Heterozygotie und genetische Vielfalt
Die genetische Drift führt in endlichen Populationen unweigerlich zu einem Verlust an genetischer Variabilität. Dieser Verlust lässt sich über die Heterozygotie (H) quantifizieren, welche den Anteil mischtbiger (heterozygoter) Individuen in einer Population darstellt.
Der erwartete Anteil ungleicher Allelpaare nimmt über t Generationen wie folgt ab:
H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ
Dieser mathematische Zusammenhang verdeutlicht, dass die Geschwindigkeit des Diversitätsverlusts direkt von der effektiven Populationsgröße Ne abhängt. In kleinen Populationen schwindet die Heterozygotie rapide, während große Populationen ihre genetische Vielfalt über wesentlich mehr Generationen hinweg bewahren. Der Simulator stellt dieser theoretischen Erwartung ("Erwartetes H bei T") den tatsächlich gemessenen Wert der Simulationen gegenüber ("Mittleres simuliertes H bei T").
Das Paradoxon des Mittelwerts
Bei der Betrachtung vieler unabhängiger Populationen (Wiederholungen) zeigt sich ein scheinbares Paradoxon: Während jede einzelne Population einen unvorhersehbaren Zufallspfad einschlägt und schließlich die genetische Vielfalt durch Fixierung oder Verlust komplett einbüßt, bleibt die mittlere Allelhäufigkeit über alle Replikate hinweg im Erwartungswert stabil bei p₀.
Der Mittelwert der simulierten Pfade ("Mittleres finales p") weicht durch die begrenzte Anzahl an Wiederholungen stochastisch leicht ab, bestätigt jedoch das theoretische Prinzip, dass neutrale Drift keine gerichtete Evolution verursacht. Sie verändert lediglich die Verteilung der Varianz – weg von der Variabilität innerhalb einzelner Populationen, hin zu einer maximalen genetischen Differenzierung zwischen den Populationen.
Grenzen und Annahmen des Modells
Das im Simulator genutzte Wright–Fisher-Modell basiert auf strikten theoretischen Annahmen, die in realen biologischen Systemen selten vollständig erfüllt sind. Zu diesen Annahmen gehören:
- Eine konstante effektive Populationsgröße über alle Generationen hinweg.
- Diskrete, nicht-überlappende Generationen.
- Zufällige, unabhängige Stichprobenziehung der Gameten (Panmixie).
- Vollständige selektive Neutralität des untersuchten Locus.
- Ein einzelner, ungekoppelter Locus mit genau zwei Allelen.
- Das vollständige Fehlen von Mutation und Migration (Genfluss).
Verwenden Sie dieses Tool daher nicht als Vorhersage für reale evolutionäre Prozesse, wenn Selektion, Populationsänderungen, Ungleichgewichte im Geschlechterverhältnis oder Fortpflanzungserfolg, überlappende Generationen, Inzucht, Kopplung, Genfluss oder Populationsstrukturen eine Rolle spielen. Ein technischer Sicherheitsfaktor ist auf dieses stochastische biologische Modell nicht anwendbar.
Funktionsweise und technische Limits
Um eine Überlastung des Endgeräts zu verhindern, erzwingt die Anwendung ein browserseitiges Limit von maximal 5.000.000 geplanten Allelkopie-Ziehungen. Dieser Wert berechnet sich aus dem Produkt:
2Ne × Generationen × Wiederholungen
Sollte dieses Produkt überschritten werden, blockiert das System die Ausführung mit der Fehlermeldung: „Dieses Setup plant ‹draws› Allelkopie-Ziehungen, was das Browser-Limit von ‹limit› überschreitet. Reduzieren Sie die Populationsgröße, die Generationen oder die Wiederholungen.“.
Alle Parameter und simulierten Populationen verbleiben in diesem Browser. Die Verarbeitung erfolgt vollständig lokal auf Ihrem Gerät; es werden keine Daten hochgeladen oder extern gespeichert.
Häufig gestellte Fragen (FAQ)
Sollte ich die Zensusgröße oder die effektive Populationsgröße eingeben?
Geben Sie die effektive Populationsgröße Ne ein: die Größe einer idealisierten Fortpflanzungsgemeinschaft, die mit derselben Rate driften würde wie die Population, die Sie untersuchen möchten. Die tatsächliche Individuenzahl (Zensusgröße) kann viel größer sein, wenn Geschlechterverhältnisse, Fortpflanzungserfolg, Populationsschwankungen oder die Populationsstruktur ungleichmäßig sind.
Entspricht die neutrale Wahrscheinlichkeit einer schlussendlichen Fixierung der Ausgangshäufigkeit?
Ja, in diesem neutralen Modell ohne Mutation oder Migration hat ein Allel, das mit der Häufigkeit p₀ startet, eine schlussendliche Fixierungswahrscheinlichkeit von p₀ und eine Verlustwahrscheinlichkeit von 1 − p₀. Die oben gezeigten Raten decken jedoch nur die eingegebene Anzahl an Generationen und Wiederholungen ab, sodass sie je nach Startwert variieren.
Warum kann die mittlere Allelhäufigkeit nahe p₀ bleiben, während die Heterozygotie sinkt?
Die neutrale Drift hat über viele unabhängige Populationen hinweg keine bevorzugte Richtung, sodass die mittlere Häufigkeit im Erwartungswert bei p₀ bleibt. Einzelne Populationen driften jedoch in Richtung 0 oder 1, und beide Grenzen weisen eine Heterozygotie von null auf. Daher geht die Variation verloren, selbst wenn der Mittelwert über die Populationen hinweg nahezu unverändert bleibt.
Kann dies vorhersagen, wie sich eine reale Population verhalten wird?
Nein. Sie zeigt Ergebnisse unter einem bewusst idealisierten, neutralen Wright–Fisher-Modell. Reale Vorhersagen erfordern belastbare Schätzungen der effektiven Größe sowie aller für den Fall relevanten Faktoren wie Selektion, Mutation, Migration, sich ändernde Populationsgrößen, überlappende Generationen, Kopplung und Populationsstruktur.