Simuluesi i Driftit Gjenetik të Popullatës

Dinamika e Driftit Gjenetik përmes Modelit Wright–Fisher

Drifti gjenetik përfaqëson një nga forcave themelore të evolucionit, duke shkaktuar ndryshime të rastësishme në frekuencat e aleleve brenda një popullate përgjatë gjeneratave. Ndryshe nga seleksionimi natyror, ku përhapja e një aleli varet nga vlera e tij përshtatëse, drifti gjenetik është një proces tërësisht neutral. Ai ndodh për shkak të gabimit të mostrimit të rastësishëm gjatë kalimit të gametave nga njëra gjeneratë në tjetrën.

Për të vizualizuar dhe analizuar këtë fenomen, ky simulues përdor modelin neutral diploid Wright–Fisher. Ky model matematik na lejon të ndjekim rrugët e frekuencës së aleleve në popullata të pavarura dhe të vëzhgojmë se si variacioni gjenetik humbet ose fiksohet me kalimin e kohës. Përmes këtij mjeti, përdoruesit mund të kuptojnë se si madhësia e popullatës ndikon drejtpërdrejt në shpejtësinë e humbjes së heterozigozitetit dhe në normat e fiksimit të aleleve.

Parametrat e Konfigurimit të Popullatës

Për të ekzekutuar një simulim, përdoruesi duhet të plotësojë parametrat e kërkuar në ndërfaqen e konfigurimit. Këto të dhëna përcaktojnë strukturën gjenetike dhe kufijtë e popullatës së simuluar:

• Skenarët: Përdoruesit mund të zgjedhin midis skenarëve të paracaktuar si "Popullatë e vogël", "Popullatë e madhe" ose "Alel i rrallë" për të ngarkuar shpejt profile të ndryshme eksperimentale.
• Madhësia efektive e popullatës, Ne: Ky parametër përfaqëson numrin e individëve riprodhues dhe duhet të jetë një numër i plotë nga 1 deri në 100,000.
• Frekuenca fillestare e alelit, p₀: Përpjesëtimi fillestar i alelit të synuar në popullatë, i futur si një proporcion midis 0 dhe 1, duke përfshirë të dy kufijtë.
• Gjeneratat: Numri i gjeneratave përgjatë të cilave do të ndiqet drifti, i cili duhet të jetë një numër i plotë nga 1 deri në 2,000.
• Popullatat e pavarura: Numri i përsëritjeve të pavarura (replikateve) që do të simulohen paralelisht, i cili duhet të jetë një numër i plotë nga 1 deri në 200.
• Vlerë fillestare e rastësishme: Një fushë teksti (deri në 100 karaktere) që shërben si "seed" për të garantuar riprodhueshmërinë e saktë të rrugëve të simuluara. Përdoruesit mund të klikojnë butonin "Vlerë e re fillestare" për të gjeneruar një kod të ri.

Gjatë futjes së të dhënave, sistemi kryen verifikime të sakta. Nëse një fushë nuk plotëson kriteret, shfaqen mesazhe gabimi specifike si:

• ‹field›: “‹token›” nuk është numër.
• ‹field›: vendosni një numër të plotë.
• Madhësia efektive e popullatës duhet të jetë një numër i plotë nga 1 deri në 100,000.
• Frekuenca fillestare e alelit duhet të jetë midis 0 dhe 1, duke përfshirë të dy kufijtë.
• Gjeneratat duhet të jenë një numër i plotë nga 1 deri në 2,000.
• Popullatat e pavarura duhet të jenë një numër i plotë nga 1 deri në 200.
• Vendosni një vlerë fillestare (seed) të rastësishme në mënyrë që simulimi të mund të riprodhohet.
• Mbajeni vlerën fillestare (seed) të rastësishme në 100 karaktere ose më pak.

Për të shmangur mbingarkesën e shfletuesit, ekziston një kufi maksimal prej 5,000,000 tërheqjesh të planifikuara të kopjeve të aleleve, i llogaritur si 2Ne × gjeneratat × përsëritjet. Nëse kalohet ky kufi, shfaqet gabimi: Ky konfigurim planifikon ‹draws› tërheqje të kopjeve të aleleve, mbi kufirin ‹limit› të shfletuesit. Zvogëloni madhësinë e popullatës, gjeneratat ose përsëritjet.

Interpretimi i Rezultateve dhe Grafikëve

Pas ekzekutimit të simulimit, ndërfaqja përditësohet me statusin përkatës. Gjatë procesit shfaqet mesazhi Duke mostruar gjeneratat e ardhshme…, i cili pasohet nga konfirmimi U simuluan ‹replicates› popullata përgjatë ‹generations› gjeneratash duke përdorur vlerën fillestare “‹seed›”. ose udhëzimi Ndryshoni skenarin e popullatës ose ekzekutoni përsëri këtë shembull me vlerë fillestare.

Grafiku kryesor mban etiketën ARIA Frekuenca e aleleve sipas gjeneratës për popullatat e simuluara të përsëritura. Boshti X përfaqëson gjeneratat (Gjenerata), ndërsa boshti Y tregon frekuencën e alelit (Frekuenca e alelit p). Një vijë horizontale referuese tregon vlerën p e pritshme = p₀. Nën grafik, treguesi i rrugëve specifikon: Po shfaqen ‹shown› nga ‹total› trajektore; çdo përsëritje është e përfshirë në përmbledhje.

Tabela e statistikave përmbledhëse ofron të dhënat e mëposhtme për të gjitha popullatat e simuluara:

• Fiksuar në p = 1: Përpjesëtimi i popullatave ku aleli ka arritur fiksimin e plotë.
• Humbur në p = 0: Përpjesëtimi i popullatave ku aleli është zhdukur tërësisht.
• Jo e absorbuar deri në T: Popullatat që mbeten polimorfe (ku frekuenca nuk ka arritur ende në 0 ose 1) në gjeneratën përfundimtare.
• p mesatare përfundimtare: Frekuenca mesatare e alelit në të gjitha replikatet në fund të simulimit.
• H e pritshme në T: Heterozigoziteti teorik i parashikuar për gjeneratën përfundimtare.
• H mesatare e simuluar në T: Heterozigoziteti mesatar i vëzhguar në popullatat e simuluara.
• Gjenerata mesatare e absorbimit: Gjenerata mesatare në të cilën popullatat arritën fiksimin ose humbjen (shfaq nuk u arrit nëse asnjë popullatë nuk është absorbuar).

Metoda Matematikore dhe Saktësia

Simuluesi bazohet në modelin neutral diploid Wright–Fisher. Në këtë model, çdo gjeneratë e re formohet përmes një mostrimi binomial të saktë të kopjeve të aleleve nga gjenerata paraardhëse, pa përdorur përafrime normale. Llogaritjet kryhen me saktësi të dyfishtë (double precision) në JavaScript, ku frekuencat e mundshme qëndrojnë në rritje prej 1 ÷ (2Ne). Shfaqja e vlerave në ndërfaqe rrumbullakoset në një maksimum prej gjashtë shifrash dhjetore.

Seksioni Formula, zëvendësimi dhe saktësia detajon ekuacionet e përdorura për llogaritjet:

1. Tranzicioni Wright–Fisher: Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne).
2. Momentet e gjeneratës së ardhshme: E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne).
3. Rënia e heterozigozitetit: H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ.

Për të ilustruar zëvendësimin matematik në gjeneratën e parë, mjeti shfaq vlerat e mëposhtme të llogaritura:

• Gjenerata e parë përdor 2 × ‹population› = ‹copies› tërheqje të kopjeve të aleleve.
• E[X₁] = ‹copies› × ‹frequency› = ‹mean› kopje A.
• Për gjeneratën e parë, Var(X₁) = ‹copies› × ‹frequency› × (1 − ‹frequency›) = ‹variance›; SD(p₁) = ‹sd›.
• Përsëritja e parë e grafikuar mostroi X₁ = ‹sampled› kopje A, kështu që p₁ = ‹sampled› ÷ ‹copies› = ‹frequency›.

Supozimet dhe Kufizimet e Modelit

Ky simulues është një model teorik i idealizuar dhe nuk duhet të përdoret për të parashikuar ecurinë e popullatave reale ku faktorët jo-neutralë luajnë rol. Modeli mbështetet në supozimet e mëposhtme:

• Madhësi efektive konstante e popullatës përgjatë gjeneratave.
• Gjenerata diskrete dhe pa mbivendosje.
• Mostrim i rastësishëm dhe tërësisht i pavarur i gametave.
• Neutralitet selektiv (nuk ka përzgjedhje natyrore).
• Një lokus i vetëm i palidhur me saktësisht dy alele.
• Mungesë e plotë e mutacionit dhe migrimit (rrjedhës së gjeneve).

Mos e përdorni këtë mjet si parashikim kur në popullatën tuaj studimore ka rëndësi seleksionimi, ndryshimi i madhësisë së popullatës, çekuilibri gjinor ose i suksesit riprodhues, gjeneratat e mbivendosura, kryqëzimi i afërt (inbreeding), lidhja gjenetike, rrjedha e gjeneve ose struktura e popullatës. Gjithashtu, vlen të theksohet se një faktor sigurie inxhinierik nuk zbatohet për këtë model biologjik stokastik.

Privatësia dhe Përpunimi i të Dhënave

Të gjitha llogaritjet dhe simulimet kryhen lokalisht në pajisjen tuaj. Parametrat tuaj dhe popullatat e simuluara qëndrojnë në këtë shfletues dhe asnjë e dhënë nuk ngarkohet apo dërgohet në serverë të jashtëm.

---

Pyetje të Shpeshta (FAQ)

A duhet të vendos madhësinë e regjistruar apo madhësinë efektive të popullatës?
Vendosni madhësinë efektive të popullatës, Ne: madhësinë e një popullate riprodhuese të idealizuar që do të pësonte drift gjenetik me të njëjtën shkallë si popullata që synoni të studioni. Madhësia e regjistruar (census size) mund të jetë shumë më e madhe kur raportet gjinore, suksesi riprodhues, luhatjet e popullatës ose struktura e popullatës janë të pabarabarta.

Pse frekuenca mesatare e alelit mund të mbetet afër p₀ ndërsa heterozigoziteti bie?
Drifti gjenetik neutral nuk ka drejtim të preferuar nëpër shumë popullata të pavarura, kështu që frekuenca mesatare mbetet p₀ si pritshmëri. Megjithatë, popullatat individuale shpërndahen drejt 0 ose 1, dhe të dy kufijtë kanë heterozigozitet zero. Prandaj, variacioni humbet edhe pse mesatarja ndërmjet popullatave mbetet pothuajse e pandryshuar.

A është probabiliteti neutral i fiksimit përfundimtar i barabartë me frekuencën fillestare?
Po, në këtë model neutral pa mutacion apo migrim, një alel që fillon me frekuencën p₀ ka probabilitet përfundimtar fiksimi p₀ dhe probabilitet humbjeje 1 − p₀. Shkallët e shfaqura më sipër mbulojnë vetëm numrin e futur të gjeneratave dhe përsëritjeve, kështu që ato do të ndryshojnë sipas vlerës fillestare (seed).

A mund të parashikojë kjo se çfarë do të bëjë një popullatë reale?
Jo. Ai tregon rezultatet sipas një modeli neutral të idealizuar qëllimisht të Wright–Fisher. Parashikimet reale kërkojnë vlerësime të mbrojtshme të madhësisë efektive dhe të çdo seleksionimi, mutacioni, migrimi, ndryshimi të madhësisë së popullatës, gjeneratave të mbivendosura, lidhjes gjenetike dhe strukturës së popullatës që kanë rëndësi për rastin.