Hermir fyrir erfðaflökt í stofnerfðafræði
Erfðaflökt (e. genetic drift) er einn af grunnkröftum þróunarinnar. Í litlum sem stórum þýðum valda tilviljanakenndar sveiflur við úrtakstöku kynslóða því að samsætur annaðhvort festast eða tapast með tímanum. Þessi hermir keyrir hið hlutlausa tvílitna Wright–Fisher líkan sem lýst er af Nature Education og í Stanford Human Genetics kennslubókinni. Hver ný kynslóð er slembiúrtak af 2Ne samsætuafritum úr þeirri fyrri. Heimildir voru skoðaðar 6. ágúst 2026.
Með því að nota þetta tól er hægt að fylgjast með og bera saman endurskapanlega ferla samsætu tíðni yfir margar kynslóðir. Hermirinn sýnir hvernig tilviljun ein og sér stýrir erfðasamsetningu þegar önnur öfl eins og val, stökkbreytingar eða fólksflutningar eiga sér ekki stað. Allir útreikningar og hermingar fara fram staðbundið í vafra notandans. Gildi þín og hermd þýði verða eftir í þessum vafra og engin gögn eru send annað.
Uppsetning þýðis og breytur
Í viðmótinu „Uppsetning þýðis“ getur notandinn skilgreint upphafsskilyrði hermunarinnar. Hægt er að velja á milli þriggja tilbúinna sviðsmynda:
- Stórt þýði
- Lítið þýði
- Sjaldgæf samsæta
Til að sérsníða hermunina eru eftirfarandi breytur slegnar inn:
| Heiti reits í viðmóti | Lýsing og eining | Leyfileg gildi |
|---|---|---|
| Virk þýðisstærð, Ne | Fjöldi undaneldisdýra í hinu hermda þýði | Heil tala frá 1 til 100.000 |
| Upphafstíðni samsætu, p₀ | Upphafshlutfall samsætunnar í þýðinu | Hlutfall á bilinu 0 til 1 (báðir endapunktar meðtaldir) |
| Kynslóðir | Fjöldi kynslóða sem hermt er yfir | Heil tala frá 1 til 2.000 |
| Sjálfstæð þýði | Fjöldi endurtekna þýða í sömu keyrslu | Heil tala frá 1 til 200 |
| Slembiútsæði | Textastrengur sem tryggir endurskapanleika | Allt að 100 stafir |
Hægt er að smella á hnappinn „Nýtt útsæði“ til að fá slembigildi fyrir útsæðið. Ef slegið er inn ógilt gildi mun hermirinn birta viðeigandi villuboð:
- Ef gildi er ekki tala:
‹field›: „‹token›“ er ekki tala. - Ef slegið er inn aukastafatala þar sem krafist er heillar tölu:
‹field›: sláðu inn heila tölu. - Ef virk þýðisstærð er utan marka:
Virk þýðisstærð verður að vera heil tala frá 1 til 100.000. - Ef upphafstíðni er utan marka:
Upphafstíðni samsætu verður að vera á milli 0 og 1, að báðum endamörkum meðtöldum. - Ef kynslóðafjöldi er utan marka:
Kynslóðir verða að vera heil tala frá 1 til 2.000. - Ef fjöldi sjálfstæðra þýða er utan marka:
Sjálfstæð þýði verða að vera heil tala frá 1 til 200. - Ef útsæði vantar:
Sláðu inn slembiútsæði svo hægt sé að endurskapa hermunina. - Ef útsæði er of langt:
Hafðu slembiútsæðið 100 stafir eða færri.
Reiknigeta og takmarkanir vafrans
Til að koma í veg fyrir að vafrinn frjósi við stórar hermanir er hámarksfjöldi áætlaðra útdrátta takmarkaður við 5.000.000 samsætuafrit. Þessi fjöldi er reiknaður sem:
Áætlaðir útdrættir = 2Nₑ × kynslóðir × sjálfstæð þýði
Ef þessi mörk eru brotin birtast villuboðin: Þessi uppsetning gerir ráð fyrir ‹draws› útdráttum á samsætuafritum, sem er yfir ‹limit› takmörkunum vafrans. Minnkaðu þýðisstærð, kynslóðir eða endurtekningar.
Formúla, gildisinnsetning og nákvæmni
Sérhver tvíkostaútdráttur notar nákvæma Bernoulli-úrtakstöku; engin normalnálgun er notuð. Útreikningar halda tvöfaldri nákvæmni í JavaScript. Tíðnir liggja á hlutfallinu 1 ÷ (2Ne); birtingin námundar að hámarki við sex aukastafi.
Stærðfræðilegt líkan hermisins byggir á eftirfarandi jöfnum:
- Wright–Fisher umbreytingin: Við hverja kynslóð dregur Wright–Fisher umbreytingin slembiúrtak fyrir næsta samsætufjölda og deilir honum með tiltækum genaafritum: Xₜ₊₁ ≈ Binomial(2Nₑ, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Nₑ)
- Væntigildi og dreifni: Fyrir úrtakstöku lýsa þessi væntigildi og dreifni sviði mögulegra næstu kynslóða: E[Xₜ₊₁] = 2Nₑ · pₜ; Var(Xₜ₊₁) = 2Nₑ · pₜ(1 - pₜ); Var(pₜ₊₁) = pₜ(1 - pₜ) ÷ (2Nₑ)
- Arfblendni (Heterozygosity): Vænt hlutfall ólíkra samsætupara minnkar yfir t kynslóðir á eftirfarandi hátt: H₀ = 2p₀(1 - p₀); Hₜ = H₀(1 - 1 ÷ (2Nₑ))ᵗ
Dæmi um gildisinnsetningu í fyrstu kynslóð
Þegar hermirinn er keyrður reiknar hann út gildi fyrir fyrstu kynslóðina og birtir eftirfarandi skref:
- Fyrsta kynslóð notar 2 ×
‹population›=‹copies›útdrætti á samsætuafritum. - E[X₁] =
‹copies›×‹frequency›=‹mean›A-afrit. - Fyrir fyrstu kynslóðina er Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›. - Fyrsta teiknaða endurtekningin fékk X₁ =
‹sampled›A-afrit í úrtaki, svo p₁ =‹sampled›÷‹copies›=‹frequency›.
Niðurstöður og myndræn framsetning
Eftir að smellt hefur verið á „Keyra flökthermun“ sýnir hermirinn stöðuna með skilaboðunum Dregur úrtak fyrir næstu kynslóðir… og síðan Hermdi ‹replicates› þýði í gegnum ‹generations› kynslóðir með því að nota útsæðið „‹seed›“. Ef breytingar eru gerðar birtast skilaboðin Breyttu sviðsmynd þýðisins eða keyrðu þetta útsæðisdæmi aftur.
Ferlar samsætu tíðni
Línuritið sýnir þróun samsætutíðninnar yfir tíma. Á x-ásnum er „Kynslóð“ og á y-ásnum er „Samsætu tíðni p“. Á ritinu er einnig sýnd viðmiðunarlína sem merkt er „Vænt p = p₀“. Undir ritinu stendur: Sýnir ‹shown› af ‹total› ferlum; öll endurtekningasöfn eru tekin með í samantektinni.
Samantekt tölfræði
Eftirfarandi tölfræðilegu mælikvarðar eru reiknaðir út frá öllum endurteknum þýðum:
- Fest við p = 1: Hlutfall þýða þar sem samsætan náði festingu.
- Tapast við p = 0: Hlutfall þýða þar sem samsætan tapaðist alveg.
- Ekki runnið saman við T: Hlutfall þýða þar sem tíðnin er enn á milli 0 og 1 við lok kynslóða (T).
- Meðaltal loka-p: Meðaltal samsætutíðni allra þýða í lok hermunarinnar.
- Vænt H við T: Fræðileg væntanleg arfblendni eftir T kynslóðir.
- Meðaltal hermdrar H við T: Raunveruleg meðalarfblendni sem mældist í hermdu þýðunum við kynslóð T.
- Meðalkynslóð við samruna: Meðaltal þeirra kynslóða þar sem samsætan annað hvort festist eða tapaðist (birtir
ekki náðef engin samruni átti sér stað).
Aðferð, forsendur og takmarkanir
Líkanið gerir ráð fyrir ákveðnum einföldunum sem mikilvægt er að hafa í huga við túlkun niðurstaðna:
- Stöðug virk þýðisstærð: Stærðin Nₑ breytist ekki á milli kynslóða.
- Stakar kynslóðir sem skarast ekki: Allir foreldrar deyja áður en afkvæmi æxlast.
- Slembið sjálfstætt úrtak: Engin makavalshlutdrægni eða landfræðileg takmörkun.
- Hlutlaust val: Samsætan hefur engin áhrif á hæfni (survival eða æxlun).
- Eitt ótengt tveggja samsætu sæti: Engin tengsl við önnur genasæti.
- Engar stökkbreytingar eða fólksflutningar: Engin ný gen flæða inn í kerfið.
Ekki nota þetta sem spá þegar val, þýðisbreytingar, ójafnvægi kynja eða æxlunarárangurs, skarast kynslóðir, skyldleikarækt, tengsl, genaflæði eða þýðisgerð skipta máli. Verkfræðilegur öryggisstuðull á ekki við um þetta slembna líffræðilega líkan.
Algengar spurningar
Á ég að slá inn raunverulegan einstaklingafjöldi eða virka þýðisstærð?
Sláðu inn virka þýðisstærð, Ne: stærð á einfölduðu undaneldisþýði sem myndi flökta á sama hraða og þýðið sem þú ætlar að rannsaka. Raunverulegur einstaklingafjöldi getur verið mun meiri þegar kynjahlutföll, æxlunarárangur, sveiflur í þýði eða þýðisgerð eru ójöfn.
Eru hlutlausu líkurnar á endanlegri festingu jafnar upphafstíðninni?
Já, í þessu hlutlausa líkani án stökkbreytinga eða fólksflutninga eru líkurnar á endanlegri festingu samsætu sem byrjar í tíðninni p₀ jafnar p₀ og líkurnar á tapi 1 − p₀. Hlutföllin sem sýnd eru hér að ofan ná aðeins yfir þann fjölda kynslóða og endurtekninga sem sleginn er inn, svo þau munu breytast eftir útsæðinu.
Af hverju getur meðaltíðni samsætu haldist nálægt p₀ á meðan arfblendni minnkar?
Hlutlaust flökt hefur enga ákveðna stefnu yfir mörg sjálfstæð þýði, þannig að meðaltíðnin helst p₀ í væntingu. Einstök þýði dreifast hins vegar í átt að 0 eða 1, og á báðum endamörkum er arfblendni engin. Breytileiki tapast því jafnvel þó að meðaltal þýða haldist nánast óbreytt.
Getur þetta spáð fyrir um hvað raunverulegt þýði mun gera?
Nei. Hún sýnir niðurstöður samkvæmt viljandi einfölduðu og hlutlausu Wright–Fisher líkani. Raunverulegar spár krefjast áreiðanlegs mats á virkri stærð og hvers kyns vali, stökkbreytingum, fólksflutningum, breytilegri þýðisstærð, skarast kynslóðum, tengslum og þýðisgerð sem máli skipta í hverju tilfelli.