Il-Mudell ta' Wright–Fisher u d-Drift Ġenetiku
Id-drift ġenetiku huwa wieħed mill-mekkaniżmi fundamentali tal-evoluzzjoni li jiddeskrivi l-bidliet każwali fil-frekwenzi tal-alleli minn ġenerazzjoni għal oħra. B'differenza mis-selezzjoni naturali, fejn l-alleli jinfirxu jew jintilfu minħabba l-vantaġġ jew l-iżvantaġġ adattiv tagħhom, id-drift ġenetiku jopera b'mod kompletament newtrali.
Il-Simulatur tad-Drift tal-Ġenetika tal-Popolazzjoni jimplimenta l-mudell diploid newtrali ta' Wright–Fisher deskritt minn Nature Education u t-test Stanford Human Genetics. Dan il-mudell jassumi li kull ġenerazzjoni ġdida hija kampjun każwali ta' 2Ne kopji ta' alleli mill-ġenerazzjoni preċedenti. Permezz ta' dan il-proċess ta' kampjunar każwali, il-frekwenzi tal-alleli jvarjaw b'mod stokastiku tul iż-żmien, u dan iwassal biex xi alleli jiġu ffissati fil-popolazzjoni filwaqt li oħrajn jintilfu għal kollox.
Issettjar tal-Popolazzjoni u Parametri tal-Input
Biex tibda simulazzjoni, l-utent għandu jikkonfigura l-parametri fl-interface ta' "Issettjar tal-popolazzjoni". Din is-sistema tippermetti l-analiżi ta' allelu newtrali wieħed, daqs wieħed ta' popolazzjoni diplojde, u bosta ripetizzjonijiet indipendenti. Il-parametri li ġejjin huma meħtieġa:
- Xenarji (Presets): L-utent jista' jagħżel malajr bejn konfigurazzjonijiet lesti bħal "Popolazzjoni żgħira", "Popolazzjoni kbira", jew "Allelu rari".
- Daqs effettiv tal-popolazzjoni, Ne: Dan jirrappreżenta n-numru ta' individwi li jissenslu. L-input irid ikun numru sħiħ minn 1 sa 100,000.
- Frekwenza tal-alleli tal-bidu, p₀: Din hija l-frekwenza inizjali tal-allelu fil-popolazzjoni, imdaħħla bħala proporzjon, 0–1. Il-valur irid ikun bejn 0 u 1, inklużi ż-żewġ limiti.
- Ġenerazzjonijiet: In-numru ta' ġenerazzjonijiet li tulhom se tiġi ssimulata l-popolazzjoni. L-input irid ikun numru sħiħ minn 1 sa 2,000.
- Popolazzjonijiet indipendenti: In-numru ta' repliki indipendenti li għandhom jiġu ssimulati fl-istess ħin. L-input irid ikun numru sħiħ minn 1 sa 200.
- Seed każwali: Qasam tat-test użat biex tiġi żgurata r-riproduċibbiltà tal-mogħdijiet simulati. Is-seed trid tkun ta' 100 karattru jew inqas. L-utent jista' wkoll jikklikkja fuq il-buttuna "Seed ġdida" biex jiġġenera valur każwali ġdid.
Il-brawżer jinfurza limitu ta' xogħol biex jipprevjeni l-iffriżar tal-applikazzjoni. Jekk l-estrazzjonijiet ippjanati ta' kopji ta' alleli (ikkalkulati bħala 2Ne × ġenerazzjonijiet × repliki) jaqbżu l-5,000,000, is-simulazzjoni tiġi mblukkata u jintwera l-messaġġ ta' żball korrispondenti.
Metodu, Assunzjonijiet, u Limiti
Il-mudell matematiku użat f'dan is-simulatur huwa bbażat fuq assunzjonijiet speċifiċi li jiddefinixxu popolazzjoni idealizzata. Dawn l-assunzjonijiet huma:
- Daqs effettiv tal-popolazzjoni kostanti tul il-ġenerazzjonijiet kollha.
- Ġenerazzjonijiet diskreti li ma jidħlux fuq xulxin.
- Kampjunar każwali indipendenti ta' gameti.
- Newtralità selettiva sħiħa (l-ebda allelu ma jipprovdi vantaġġ ta' sopravivenza jew riproduzzjoni).
- Locus wieħed mhux marbut b'żewġ alleli.
- L-ebda mutazzjoni u lanqas migrazzjoni.
Minħabba dawn l-assunzjonijiet, dan il-mudell m'għandux jintuża bħala tbassir għal popolazzjonijiet reali fejn fatturi bħas-selezzjoni, il-bidla fil-popolazzjoni, l-iżbilanċ fis-sess jew fis-suċċess riproduttiv, il-ġenerazzjonijiet li jidħlu fuq xulxin, l-inbreeding, ir-rabta, il-fluss tal-ġeni, jew l-istruttura tal-popolazzjoni huma importanti. Barra minn hekk, fattur ta' sikurezza tal-inġinerija ma japplikax għal dan il-mudell bijoloġiku stokastiku.
Kull estrazzjoni binomjali fis-simulatur tuża kampjunar eżatt ta' Bernoulli minflok approssimazzjoni normali. Il-kalkoli kollha jżommu l-preċiżjoni doppja ta' JavaScript, u l-frekwenzi kollha possibbli jinsabu fuq inkrementi ta' 1 ÷ (2Ne). Il-wiri tal-interface jqarrab il-valuri kollha għal massimu ta' sitt postijiet deċimali.
Il-parametri kollha u l-popolazzjonijiet simulati tiegħek jibqgħu f'dan il-brawżer. Il-kalkoli kollha jsiru lokalment fuq it-tagħmir tal-utent, u l-ebda data ma tiġi mtellgħa jew mibgħuta lil servers esterni.
Formula, Sostituzzjoni, u Preċiżjoni
Il-matematika wara t-tranżizzjonijiet tal-ġenerazzjonijiet u l-istatistika tas-simulazzjoni hija deskritta permezz tal-formuli li ġejjin:
F'kull ġenerazzjoni, it-tranżizzjoni ta' Wright–Fisher tieħu kampjun tal-għadd tal-alleli li jmiss u tiddividih bil-kopji tal-ġeni disponibbli: Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne)
Qabel il-kampjunar, dawn il-mumenti jiddeskrivu l-firxa ta' ġenerazzjonijiet li jmiss possibbli: E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne)
Is-sehem mistenni ta' pari ta' alleli mhux l-istess (eterozigożità) jonqos tul t ġenerazzjonijiet kif ġej: H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ
Meta titħaddem simulazzjoni, is-sistema turi s-sostituzzjonijiet speċifiċi għall-ewwel ġenerazzjoni bbażati fuq il-parametri mdaħħla:
- L-ewwel ġenerazzjoni tuża 2 ×
‹population›=‹copies›estrazzjoni ta' kopji ta' alleli. - E[X₁] =
‹copies›×‹frequency›=‹mean›kopja ta' A. - Għall-ewwel ġenerazzjoni, Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›. - L-ewwel replika mpinġija ħadet kampjun ta' X₁ =
‹sampled›kopja ta' A, għalhekk p₁ =‹sampled›÷‹copies›=‹frequency›.
Interpretazzjoni tar-Riżultati u l-Mogħdijiet
Wara li tikklikkja fuq "Ibda s-simulazzjoni tad-drift", l-interface turi chart bit-titlu "Mogħdijiet tal-frekwenza tal-alleli" u l-ARIA label "Frekwenza tal-alleli skont il-ġenerazzjoni għal popolazzjonijiet simulati ripetuti". L-assi X huwa mmarkat bħala "Ġenerazzjoni" u l-assi Y bħala "Frekwenza tal-alleli p". Linja ta' referenza orizzontali turi l-valur "p mistenni = p₀".
Taħt iċ-chart, l-indikatur tal-mogħdijiet juri: "Qed jintwerew ‹shown› minn ‹total› mogħdija; kull replika hija inkluża fis-sommarju."
It-tabella tas-sommarju tipprovdi l-metriċi li ġejjin dwar il-popolazzjonijiet simulati kollha:
- Fissat f'p = 1: Il-proporzjon ta' repliki fejn l-allelu laħaq fissazzjoni sħiħa.
- Mitluf f'p = 0: Il-proporzjon ta' repliki fejn l-allelu ntilef għal kollox.
- Mhux assorbit sa T: Il-proporzjon ta' repliki li baqgħu jvarjaw mingħajr ma laħqu l-ebda limitu (0 jew 1) sal-aħħar ġenerazzjoni.
- p finali medja: Il-frekwenza medja tal-allelu madwar ir-repliki kollha fl-aħħar ġenerazzjoni.
- H mistennija f'T: L-eterozigożità teoretika kkalkulata bil-formula ta' tnaqqis.
- H simulata medja f'T: L-eterozigożità medja osservata fir-repliki simulati fl-aħħar ġenerazzjoni.
- Ġenerazzjoni medja ta' assorbiment: Il-ġenerazzjoni medja li fiha r-repliki laħqu l-fissazzjoni jew it-telf (turi "ma nlaħaqx" jekk l-ebda popolazzjoni ma ġiet assorbita).
Messaġġi ta' Żball u Validazzjoni
Is-sistema tivvalida l-inputs kollha qabel ma tibda s-simulazzjoni. Jekk jinstabu valuri invalidi, jintwerew dawn il-messaġġi ta' żball speċifiċi:
- Jekk l-input mhuwiex numru:
‹field›: “‹token›” mhijiex numru. - Jekk jiddaħħal numru deċimali fejn huwa meħtieġ numru sħiħ:
‹field›: daħħal numru sħiħ. - Daqs tal-popolazzjoni barra mil-limiti:
Id-daqs effettiv tal-popolazzjoni jrid ikun numru sħiħ minn 1 sa 100,000. - Frekwenza tal-bidu barra mil-limiti:
Il-frekwenza tal-alleli tal-bidu trid tkun bejn 0 u 1, inklużi ż-żewġ limiti. - Ġenerazzjonijiet barra mil-limiti:
Il-ġenerazzjonijiet iridu jkunu numru sħiħ minn 1 sa 2,000. - Repliki barra mil-limiti:
Il-popolazzjonijiet indipendenti jridu jkunu numru sħiħ minn 1 sa 200. - Seed nieqsa:
Daħħal seed każwali biex is-simulazzjoni tkun tista' tiġi rriprodotta. - Seed twila wisq:
Żomm is-seed każwali għal 100 karattru jew inqas. - Limitu ta' xogħol maqbuż:
Dan l-issettjar jippjana ‹draws› estrazzjonijiet ta' kopji ta' alleli, li huwa ogħla mil-limitu tal-brawżer ta' ‹limit›. Naqqas id-daqs tal-popolazzjoni, il-ġenerazzjonijiet, jew ir-ripetizzjonijiet.
Mistoqsijiet Frekwenti (FAQ)
Għandi ndaħħal id-daqs taċ-ċensiment jew id-daqs effettiv tal-popolazzjoni?
Daħħal id-daqs effettiv tal-popolazzjoni, Ne: id-daqs ta' popolazzjoni tat-tnissil idealizzata li jkollha drift bl-istess rata bħall-popolazzjoni li tixtieq tistudja. Id-daqs taċ-ċensiment jista' jkun ferm akbar meta l-proporzjonijiet tas-sessi, is-suċċess riproduttiv, il-fluttwazzjonijiet fil-popolazzjoni, jew l-istruttura tal-popolazzjoni ma jkunux ibbilanċjati.
Il-probabbiltà newtrali ta' fissazzjoni eventwali hija ugwali għall-frekwenza tal-bidu?
Iva, f'dan il-mudell newtrali mingħajr mutazzjoni jew migrazzjoni, allelu li jibda bi frekwenza p₀ għandu probabbiltà ta' fissazzjoni eventwali ta' p₀ u probabbiltà ta' telf ta' 1 − p₀. Ir-rati murija hawn fuq ikopru biss in-numru mdaħħal ta' ġenerazzjonijiet u ripetizzjonijiet, għalhekk ivarjaw skont is-seed.
Għaliex il-frekwenza medja tal-alleli tista' tibqa' qrib p₀ waqt li l-eterozigożità tonqos?
Id-drift newtrali m'għandu l-ebda direzzjoni ppreferuta fuq bosta popolazzjonijiet indipendenti, għalhekk il-frekwenza medja tibqa' p₀ bħala stennija. Madankollu, il-popolazzjonijiet individwali jinfirxu lejn 0 jew 1, u ż-żewġ limiti għandhom eterozigożità ta' żero. Għalhekk, il-varjazzjoni tintilef anke waqt li l-medja bejn il-popolazzjonijiet tibqa' kważi l-istess.