Slovník | Vyhledávání | Mapa webu
 
Analýza a modelování dynamických biologických datDiskrétní deterministické modely Aplikace Populační genetika Analýza rovnice Aplikace (75) v autonomním případě

Logo Matematická biologie

Analýza rovnice Aplikace (75) v autonomním případě

Předpokládejme, že poměr počtu zygot určitého genotypu v časovém intervalu a počtu gamet vyprodukovaných jedincem téhož genotypu v čase je stejný pro každý čas V takovém případě jsou podle Aplikace (73) reprodukční zdatnosti jednotlivých genotypů a konstantní, jsou to nezáporné parametry rovnice Aplikace (75). Budeme předpokládat, že alespoň jeden z parametrů je kladný. V opačném případě by totiž uvažovaná populace vymizela hned v první filiální generaci.

Nejprve uvažujme kladné reprodukční zdatnosti heterozygotů, Zavedeme relativní zdatnosti homozygotů vzhledem ke zdatnosti heterozygotů

Zlomek na pravé straně rovnice Aplikace (75) vykrátíme parametrem a upravíme. Dostaneme rovnici

(81)

se dvěma nezápornými parametry. Poněvadž stavová proměnná vyjadřuje relativní frekvenci (tj. pravděpodobnost), je stavovým prostorem uzavřený interval

Při analýze rovnice Aplikace (81) rozlišíme tři případy.

(i) reprodukční zdatnosti všech genotypů jsou stejné. V tomto případě je rovnice Aplikace (81) tvaru

která má jedině konstantní řešení. Pokud tedy reprodukce nezávisí na genotypu, frekvence alel v populaci se nemění. To je Hardyho-Weinbergův zákon.

(ii.) reprodukční zdatnost homozygota genotypu je stejná, jako reprodukční zdatnost heterozygota a reprodukční zdatnost homozygota genotypu je jiná. Jinak řečeno, jedinci genotypů a mají stejný fenotyp, jedinci genotypu mají fenotyp jiný. To odpovídá situaci, že alela je dominantní a alela  je recesivní. V tomto případě má rovnice Aplikace (81) tvar

(82)

Najdeme její rovnovážné body a vyšetříme jejich stabilitu.

Označme na chvíli Rovnice má dva kořeny a Platí

To znamená (viz obr. Aplikace 4), že pro každé řešení rovnice Aplikace (82) platí

Pokud výběr preferuje fenotyp určený recesivní alelou, pak dominantní alela z populace vymizí; pokud výběr preferuje fenotyp určený dominantní alelou, pak recesivní alela z populace vymizí.

Obr. 4. Grafické řešení rovnice Aplikace (81) a jeho průběh pro různé hodnoty parametrů K, k. Tj. vývoj relativní frekvence alely A v populaci pro různé relativní reprodukční zdatnosti homozygotů vzhledem k heterozygotům.

(iii) každá alela nějak přispívá k reprodukční zdatnosti, žádná z nich není dominantní. Příspěvek alel k fenotypu je aditivní, alely jsou semidominantní.

I v tomto případě najdeme rovnovážné body rovnice Aplikace (81) a vyšetříme jejich stabilitu. Nyní označíme

Pak rovnice tj. kubická rovnice

má kořeny 0, 1 a Kořen je rovnovážným bodem rovnice Aplikace (82) pouze tehdy, když což nastane právě tehdy, když nebo (stále totiž předpokládáme ). Dále platí

To znamená (viz obr. Aplikace 4), že pro je rovnovážný bod 0 asymptoticky stabilní a pro je nestabilní. Podobně pro je rovnovážný bod 1 asymptoticky stabilní a pro je nestabilní. Pokud je rovnovážným bodem rovnice Aplikace (81), pak je asymptoticky stabilní v případě a nestabilní v případě ; snadno totiž ověříme, že při označení

platí pro a

což znamená, že

Podívejme se ještě na jeden speciální případ rovnice Aplikace (82), a to takový, když V tomto případě je

a to znamená, že reprodukční zdatnost heterozygota je geometrickým průměrem zdatností jednotlivých homozygotů. Rovnice Aplikace (82) nabude tvar

což je Bevertonova-Holtova rovnice, jejíž řešení je podle Riccatiho a Bernoulliova rovnice dáno formulí

Opět tedy platí

Je-li pak je možný genetický polymorfismus. Ten je ale v případě nestabilní; záleží na počáteční genetické struktuře populace, která z alel převládne. Trvalý genetický polymorfismus je možný jen v případě, že tj. reprodukční zdatnost každého z homozygotů je menší než reprodukční zdatnost heterozygota.

 

Z dosavadních úvah jsme vyloučili případ tj. možnost, že heterozygoti nejsou schopni reprodukce. V takovém případě má autonomní rovnice Aplikace (75) tvar

(83)

Pokud pak a alela by z populace vymizela hned v prvním časovém kroku. Takovou situaci bychom mohli interpretovat tak, že alela představovala nějakou škodlivou (smrtící) mutaci. Pokud pak a z populace by bezprostředně vymizela alela Dále budeme předpokládat tj. že ani alela ani alela nepředstavuje smrtící mutaci.

Konstantní posloupnosti a jsou evidentně řešením rovnice Aplikace (83) pro jakékoliv hodnoty parametrů vyjadřují možnost, že v modelované populaci má uvažovaný gen jedinou alelu. Budeme hledat další řešení rovnice Aplikace (83), které není identicky nulové ani jednotkové. Uvažujme proto rovnici spolu s počáteční podmínkou

(84)

Substituce

převede počáteční úlohu Aplikace (83), Aplikace (84) na počáteční úlohu pro lineární rovnici

která má řešení

Zpětnou substitucí dosatneme řešení původní počáteční úlohy Aplikace (83), Aplikace (84) ve tvaru

Posloupnost s obecným členem je vybraná z geometrické posloupnosti s počátečním členem a s kvocientem Hodnota kvocientu určuje chování řešení. Je-li

pak

a řešení úlohy Aplikace (83), Aplikace (84) je konstantní, Je-li resp. pak

přitom je řešení rovnice Aplikace (83) monotonní.

Pokud je tedy existuje rovnovážný stav který je repulsivní, a dva asymptoticky stabilní rovnovážné stavy 0 a 1. Ke kterému ze stabilních rovnovážných stavů bude řešení konvergovat, závisí na počáteční podmínce. Kvalitativně se tedy jedná o stejnou situaci jako na obr. Aplikace 4, případ

Závěr: Autonomní rovnice Aplikace (75) s nezápornými parametry uvažovaná na stavovém prostoru má rovnovážná řešení v krajních bodech tj. řešení a

Pokud jsou všechny parametry stejné, pak má rovnice pouze konstantní řešení; každý bod stavového prostoru je rovnovážný.

Pokud přičemž alespoň jedna z nerovností je ostrá, pak rovnice Aplikace (75) nemá uvnitř stavového prostoru   žádné rovnovážné body.

Pokud nebo pak má autonomní rovnice Aplikace (75) izolovaný rovnovážný bod

uvnitř stavového prostoru


Dostatečné podmínky pro asymptotickou stabilitu nebo repulsivitu izolovaných rovnovážných bodů jsou:

  • Je-li nebo pak je stacionární řešení asymptoticky stabilní.
     
  • Je-li nebo pak je stacionární řešení repulsivní.
     
  • Je-li nebo pak je stacionární řešení asymptoticky stabilní.
     
  • Je-li nebo pak je stacionární řešení repulsivní.
     
  • Je-li pak je stacionární řešení asymptoticky stabilní.
     
  • Je-li pak je stacionární řešení repulsivní.

Biologickou evoluci lze chápat jako změnu frekvencí alel v průběhu času. Představme si, že celá populace má na příslušném lokusu alelu a v počátečním čase se u nějakého jedince objeví její mutace, alela Pokud v takovém případě bude přičemž alespoň jedna z těchto nerovností je ostrá, pak se bude mutovaná alela v populaci šířit a nakonec v ní převládne; pokud pak budou obě alely v populaci dlouhodobě koexistovat, populace se stane geneticky polymorfní.

 
vytvořil Institut biostatistiky a analýz Lékařské fakulty Masarykovy univerzity