Wright–Fisherov model i genetički drift
Simulator genetičkog drifta populacije je besplatni online alat koji omogućava korisnicima da simuliraju i vizualiziraju neutralni genetički drift koristeći diploidni Wright–Fisherov model. Postavljanjem simulirane populacije za razmnožavanje, korisnici mogu pratiti i upoređivati reproducibilne putanje frekvencije alela kroz više nezavisnih populacija kako bi uočili kako se genetička varijacija gubi ili fiksira kroz generacije. Alat pomaže u razumijevanju uticaja nasumičnog uzorkovanja na frekvencije alela tokom vremena, prikazujući simuliranu heterozigotnost, fiksaciju i stope gubitka, bez pretenzija da predviđa stvarne evolucijske ishode.
U svakoj generaciji, Wright–Fisherov prijelaz uzorkuje sljedeći broj alela i dijeli ga s dostupnim kopijama gena prema formuli:
Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne)
Prije uzorkovanja, ovi momenti opisuju raspon mogućih sljedećih generacija:
E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne)
Ovaj matematički model pokazuje kako se u svakoj generaciji vrši nasumični odabir iz konačnog skupa genskih kopija, što dovodi do fluktuacija u frekvenciji alela čak i u potpunom odsustvu selektivnog pritiska [i, j].
Efektivna naspram popisne veličine populacije
Prilikom postavljanja parametara u interfejsu pod nazivom "Postavke populacije", ključni parametar je "Efektivna veličina populacije, Ne", koja se unosi u polje "jedinke koje se razmnožavaju". Efektivna veličina populacije (Ne) predstavlja veličinu idealizirane populacije koja gubi genetičku raznolikost istom stopom kao i stvarna populacija koja se proučava [i, j].
U stvarnim populacijama, popisna veličina (ukupan broj jedinki) je gotovo uvijek veća od efektivne veličine [i, j]. Razlozi za ovo odstupanje uključuju:
- Neravnomjeran omjer spolova među jedinkama koje se razmnožavaju.
- Veliku varijansu u reproduktivnom uspjehu pojedinih jedinki.
- Fluktuacije u veličini populacije kroz generacije (uska grla).
- Prisustvo prostorne strukture i ograničen protok gena.
Zbog toga se u simulator unosi isključivo efektivna veličina populacije kako bi se osigurala tačnost teorijskih proračuna drifta.
Dinamika fiksacije i gubitka alela
Kada se pokrene simulacija, alat generiše grafikon pod nazivom "Putanje frekvencije alela" sa ARIA oznakom "Frekvencija alela po generacijama za ponovljene simulirane populacije". Na ovom grafikonu, X-osa nosi naziv "Generacija", dok Y-osa nosi naziv "Frekvencija alela p". Referentna linija označava vrijednost "Očekivano p = p₀".
Tokom generacija, frekvencije alela u nezavisnim replikatima lutaju dok ne dosegnu jednu od dvije apsorbirajuće granice:
- "Fiksirano na p = 1" – alel postaje jedini prisutan u populaciji.
- "Izgubljeno na p = 0" – alel u potpunosti nestaje iz populacije.
U neutralnom modelu, vjerovatnoća da će se alel na kraju fiksirati jednaka je njegovoj početnoj frekvenciji p₀ [i, j]. Ako simulacija završi prije nego što sve populacije dosegnu ove granice, simulator prikazuje udio populacija koje su označene kao "Nije apsorbovano do T". Za populacije koje su doživjele fiksaciju ili gubitak, prikazuje se "Srednja generacija apsorpcije", a ukoliko nijedna populacija ne dosegne granice, prikazuje se vrijednost "nije dostignuto".
Opadanje heterozigotnosti i genetičke raznolikosti
Heterozigotnost (H) predstavlja udio jedinki koje nose različite alele na posmatranom lokusu [i, j]. Genetički drift neumitno smanjuje heterozigotnost u konačnim populacijama jer nasumično uzorkovanje vodi ka fiksaciji jednih i gubitku drugih alela [i, j].
Očekivani udio različitih parova alela opada kroz t generacija na sljedeći način:
H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ [i, j]
Simulator izračunava i upoređuje teorijsku vrijednost "Očekivano H u T" sa stvarnim simuliranim rezultatom pod nazivom "Srednje simulirano H u T". Kroz ovaj prikaz jasno se vidi kako brzina gubljenja genetičke raznolikosti direktno zavisi od efektivne veličine populacije – manje populacije gube heterozigotnost znatno brže od većih [i, j].
Paradoks srednje vrijednosti kroz populacije
Jedan od najzanimljivijih aspekata genetičkog drifta jeste ponašanje prosječne frekvencije alela kroz mnogo nezavisnih populacija [i, j]. Dok pojedinačne populacije divlje driftaju prema fiksaciji ili gubitku, "Srednje konačno p" u svim replikatima ostaje blizu početne vrijednosti p₀.
Ovaj fenomen se naziva paradoksom srednje vrijednosti: na nivou pojedinačne populacije, genetički drift je potpuno nepredvidiv i vodi ka gubitku varijacije [i, j]. Međutim, na nivou metapopulacije (skupa svih nezavisnih populacija), prosječna frekvencija alela ostaje stabilna i predvidiva [i, j].
Metoda, pretpostavke i ograničenja
Ova stranica implementira neutralni diploidni Wright–Fisherov model opisan u referencama Nature Education [j] i Stanford Human Genetics [k]. Svaka nova generacija je nasumični uzorak od 2Ne kopija alela iz prethodne generacije. Izvori provjereni 6. augusta 2026.
Svako binomno izvlačenje koristi tačno Bernoullijevo uzorkovanje; ne koristi se normalna aproksimacija. Proračuni zadržavaju JavaScript dvostruku preciznost. Frekvencije se nalaze na koracima od 1 ÷ (2Ne); prikaz se zaokružuje na najviše šest decimalnih mjesta.
Model se zasniva na sljedećim rigoroznim pretpostavkama:
- Konstantna efektivna veličina populacije.
- Diskretne generacije koje se ne preklapaju.
- Nasumično nezavisno uzorkovanje.
- Selektivna neutralnost (nema prirodne selekcije).
- Jedan nevezani lokus sa dva alela.
- Bez mutacije ili migracije (protoka gena).
Ne koristite ovo kao prognozu kada su selekcija, promjena populacije, neravnoteža spolova ili reproduktivnog uspjeha, preklapanje generacija, ukrštanje u srodstvu, vezanost gena, protok gena ili struktura populacije od značaja. Inženjerski faktor sigurnosti se ne primjenjuje na ovaj stohastički biološki model.
Svi parametri i simulirane populacije obrađuju se lokalno i ostaju u ovom pregledniku. Podaci se ne šalju na vanjske servere.
Primjer proračuna prve generacije
Za vizualizaciju matematičkog procesa, simulator prikazuje supstituciju konkretnih vrijednosti za prvu generaciju:
- "Prva generacija koristi 2 ×
‹population›=‹copies›izvlačenja kopija alela." - "E[X₁] =
‹copies›×‹frequency›=‹mean›kopija A." - "Za prvu generaciju, Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›." - "Prvo iscrtano ponavljanje uzorkovalo je X₁ =
‹sampled›kopija A, pa je p₁ =‹sampled›÷‹copies›=‹frequency›."
Često postavljana pitanja (FAQ)
Trebam li unijeti popisnu veličinu ili efektivnu veličinu populacije?
Unesite efektivnu veličinu populacije, Ne: veličinu idealizirane populacije za razmnožavanje koja bi imala drift istom stopom kao i populacija koju želite proučavati. Popisna veličina može biti mnogo veća kada su omjeri spolova, reproduktivni uspjeh, fluktuacije populacije ili struktura populacije neravnomjerni.
Da li je neutralna vjerovatnoća konačne fiksacije jednaka početnoj frekvenciji?
Da, u ovom neutralnom modelu bez mutacije ili migracije, alel koji počinje s frekvencijom p₀ ima vjerovatnoću konačne fiksacije p₀ i vjerovatnoću gubitka 1 − p₀. Stope prikazane iznad pokrivaju samo uneseni broj generacija i ponavljanja, pa će varirati u zavisnosti od sjemena.
Zašto srednja frekvencija alela može ostati blizu p₀ dok heterozigotnost opada?
Neutralni drift nema preferirani smjer kroz mnogo nezavisnih populacija, pa prosječna frekvencija u očekivanju ostaje p₀. Međutim, pojedinačne populacije se šire prema 0 ili 1, a obje granice imaju nultu heterozigotnost. Varijacija se stoga gubi čak i dok srednja vrijednost među populacijama ostaje gotovo nepromijenjena.
Može li ovo predvidjeti šta će stvarna populacija uraditi?
Ne. Simulator prikazuje ishode pod namjerno idealiziranim neutralnim Wright–Fisherovim modelom. Stvarna predviđanja zahtijevaju opravdane procjene efektivne veličine i bilo koje selekcije, mutacije, migracije, promjene veličine populacije, preklapanja generacija, vezanosti gena i strukture populacije koje su od značaja za dati slučaj.