A Wright–Fisher-modell és a genetikai drift működése
A genetikai drift (véletlen genetikai sodródás) a populációgenetika egyik alapvető evolúciós tényezője, amely az allélgyakoriságok generációról generációra történő véletlenszerű ingadozását írja le véges populációkban. A Populációgenetikai drift szimulátor egy ingyenes online eszköz, amely lehetővé teszi a semleges genetikai drift szimulációját és vizualizációját a diploid Wright–Fisher-modell segítségével. Egy szimulált szaporodó populáció beállításával a felhasználók reprodukálható allélgyakorisági útvonalakat követhetnek nyomon és hasonlíthatnak össze több független populációban, bemutatva, hogyan vész el vagy fixálódik a genetikai variabilitás a generációk során.
A szimulátor matematikai alapját a diploid Wright–Fisher-modell adja, amelyben minden egyes generáció az előző generációból vett 2Nₑ allélmásolat véletlenszerű mintája. A modell feltételezi, hogy a generációk elkülönülnek és nem fedik át egymást, a populációméret állandó, a mintavétel teljesen véletlenszerű és független, a vizsgált lókusz szelekciós szempontból semleges, kétallélos és nem kapcsolt, valamint nincs mutáció vagy migráció.
A szimuláció során a következő matematikai összefüggések írják le a folyamatot:
- Átmeneti valószínűség: A Wright–Fisher-átmenet minden generációban mintavétellel meghatározza a következő allélszámot, és elosztja azt a rendelkezésre álló génmásolatokkal: Xₜ₊₁ ~ Binomial(2Nₑ, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Nₑ).
- Várható érték és variancia: A mintavétel előtt ezek a momentumok írják le a lehetséges következő generációk tartományát: E[Xₜ₊₁] = 2Nₑ · pₜ; Var(Xₜ₊₁) = 2Nₑ · pₜ(1 - pₜ); Var(pₜ₊₁) = pₜ(1 - pₜ) ÷ (2Nₑ).
- Heterozigozitás csökkenése: Az eltérő allélpárok várható aránya t generáción keresztül a következőképpen csökken: H₀ = 2p₀(1 - p₀); Hₜ = H₀(1 - 1 ÷ (2Nₑ))ᵗ.
Populáció beállítása és paraméterek
A szimuláció elindításához a felhasználóknak meg kell adniuk a kiindulási paramétereket a „Populáció beállítása” felületen. Az alábbi beállítások állnak rendelkezésre:
- Forgatókönyvek: Előre beállított konfigurációk a gyors indításhoz, mint például a „Kis populáció”, „Nagy populáció” és „Ritka allél” [Forgatókönyvek, Kis populáció, Nagy populáció, Ritka allél].
- Effektív populációméret, Ne: A szaporodó egyedek száma [Effektív populációméret, Ne, szaporodó egyedek]. Az értéknek 1 és 100 000 közötti egész számnak kell lennie.
- Kiindulási allélgyakoriság, p₀: Az allél kezdeti aránya a populációban [Kiindulási allélgyakoriság, p₀, arány, 0–1]. Értéke 0 és 1 között adható meg, a határértékeket is beleértve.
- Generációk: A szimulálni kívánt generációk száma [Generációk, generáció]. 1 és 2000 közötti egész számnak kell lennie.
- Független populációk: A szimulált ismétlések száma [Független populációk, ismétlés]. 1 és 200 közötti egész számnak kell lennie.
- Véletlenszerű mag: Egy legfeljebb 100 karakterből álló szöveges mező, amely biztosítja a szimulációs útvonalak pontos reprodukálhatóságát [Véletlenszerű mag, A véletlenszerű mag legfeljebb 100 karakterből állhat.]. Az „Új mag” gombbal véletlenszerű új érték generálható.
A böngésző védelme érdekében a rendszer egy maximum 5 000 000 tervezett allélmásolat-húzást engedélyez, amelyet a 2Nₑ × generációk × ismétlések képlet alapján számol ki. Ha a beállítások túllépik ezt a korlátot, a szimuláció nem fut le.
Szimulációs eredmények és grafikonok
A futtatás után a szimulátor kirajzolja az „Allélgyakorisági útvonalak” grafikont, amelynek ARIA-címkéje: „Allélgyakoriság generációnként ismételt szimulált populációkban” [Allélgyakorisági útvonalak, Allélgyakoriság generációnként ismételt szimulált populációkban]. A grafikon tengelyei a következők:
- X-tengely: „Generáció”
- Y-tengely: „p allélgyakoriság”
- Referenciavonal: „Várható p = p₀”
A grafikon alatt a „‹shown› / ‹total› útvonal megjelenítése; az összes ismétlés szerepel az összegzésben.” felirat jelzi a megjelenített adatok arányát.
A szimuláció részletes statisztikai összegzést is nyújt az alábbi mutatókkal:
- Fixálódott p = 1 értéknél: Azon populációk aránya, ahol az allél fixálódott.
- Elveszett p = 0 értéknél: Azon populációk aránya, ahol az allél teljesen eltűnt.
- Nincs abszorbeálva a T időpontig: Azon populációk aránya, amelyek a szimuláció végéig nem érték el a határértékeket.
- Átlagos végső p: Az összes ismétlés átlagos allélgyakorisága a szimuláció végén.
- Várható H a T időpontban: A heterozigozitás elméleti értéke a generációk végén.
- Átlagos szimulált H a T időpontban: A szimulált populációkban mért tényleges átlagos heterozigozitás.
- Átlagos abszorpciós generáció: Az a generáció, amelyben a populációk elérték a fixációt vagy a veszteséget (ha ez nem történt meg, a kijelzőn a „nem érte el” felirat látható) [Átlagos abszorpciós generáció, nem érte el].
Képlet, behelyettesítés és pontosság
A „Képlet, behelyettesítés és pontosság” szakasz bemutatja a konkrét futtatás első generációs számításait:
- Az első generáció 2 ×
‹population›=‹copies›allélmásolat-húzást használ. - E[X₁] =
‹copies›×‹frequency›=‹mean›A másolat. - Az első generációra: Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›. - Az első ábrázolt ismétlés mintája X₁ =
‹sampled›A másolat, így p₁ =‹sampled›÷‹copies›=‹frequency›.
Minden binomiális húzás pontos Bernoulli-mintavételt használ; nincs normális eloszlással való közelítés. A számítások megőrzik a JavaScript dupla pontosságát. A gyakoriságok az 1 ÷ (2Ne) lépésközein alapulnak; a kijelző legfeljebb hat tizedesjegyre kerekít.
Módszer, feltételezések és korlátok
Ez az oldal a Nature Education és a Stanford Human Genetics tankönyv által leírt semleges diploid Wright–Fisher-modellt valósítja meg. A források ellenőrizve: 2026. augusztus 6.
A modell szigorú elméleti feltételezéseken alapul. Valós populációk esetében a cenzusméret (a ténylegesen megszámolható egyedszám) gyakran jelentősen eltér az Nₑ effektív populációmérettől az ivararányok eltolódása, a szaporodási siker egyenetlensége vagy a populációingadozások miatt. Mivel a modell semleges, nem számol a természetes szelekcióval, a mutációval, a migrációval (génáramlással), a beltenyésztéssel, a kapcsoltsággal vagy a populációszerkezettel. Emiatt az eszköz kizárólag oktatási és elméleti demonstrációs célokat szolgál, valós evolúciós folyamatok előrejelzésére nem alkalmazható. Erre a sztochasztikus biológiai modellre nem vonatkozik mérnöki biztonsági tényező.
Adatkezelési tájékoztató
A felhasználói adatok védelme érdekében a szimuláció teljes egészében a helyi eszközön fut [privacy]. A megadott paraméterek és a szimulált populációk ebben a böngészőben maradnak, semmilyen adat nem kerül feltöltésre vagy külső szerverre továbbításra [privacy].
Gyakran Ismételt Kérdések
A cenzusméretet vagy az effektív populációméretet kell megadnom?
Az Nₑ effektív populációméretet adja meg: egy olyan idealizált szaporodó populáció méretét, amely ugyanolyan ütemben sodródna (drift), mint a vizsgálni kívánt populáció. A cenzusméret (tényleges egyedszám) sokkal nagyobb lehet, ha az ivararányok, a szaporodási siker, a populációingadozások vagy a populációszerkezet egyenetlenek.
A végső fixáció semleges valószínűsége megegyezik a kiindulási gyakorisággal?
Igen, ebben a mutáció és migráció nélküli semleges modellben a p₀ kiindulási gyakoriságú allél végső fixációjának valószínűsége p₀, az elvesztésének valószínűsége pedig 1 - p₀. A fent látható arányok csak a megadott generáció- és ismétlésszámra vonatkoznak, így a magtól függően változnak.
Miért maradhat az átlagos allélgyakoriság p₀ közelében, miközben a heterozigozitás csökken?
A semleges driftnek nincs kitüntetett iránya a sok független populációban, így az átlagos gyakoriság várható értéke p₀ marad. Az egyes populációk azonban elmozdulnak a 0 vagy az 1 felé, és mindkét határértéknél nulla a heterozigozitás. A variabilitás tehát elvész, miközben a populációk közötti átlag szinte változatlan marad.
Képes ez megjósolni egy valós populáció viselkedését?
Nem. Az eszköz egy szándékosan idealizált, semleges Wright–Fisher-modell szerinti eredményeket mutat be. A valós előrejelzésekhez az effektív méret, valamint az eset szempontjából fontos szelekció, mutáció, migráció, változó populációméret, átfedő generációk, kapcsoltság és populációszerkezet megalapozott becslése szükséges.