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.) a 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í.