Molekulární přístupy v behaviorální ekologii nPárovací systémy (mating systems) nExtra-pair fertilizations nConspecific brood parasitism nMate choice (pre- and postcopulatory) nSocial breeding (relatedness) nManipulace poměrem pohlaví nAdaptive sex ratio nSex ratio conflicts nSex-biased dispersal (pop-gen methods) nForaging nIdentifying prey, individual food preferences n n Parentage analysis infidelity Msat paternita - barvy Nevera002 Určení otce, matky nebo obou Jones AG, Ardren WR, 2003, Mol Ecol (invited review) Proč zjišťovat paternitu? nmimopárové oplození – otec vychovatel vs. genetický otec nstanovení počtu reprodukčně-aktivních jedinců v populaci nstudium faktorů ovlivňujících fitness (reprodukční úspěšnost) nvícenásobná paternita – analýzy párovacího systému (detekce promiskuity aj.) n Paternita x maternita x rodičovský pár •Známe jednoho rodiče Neznáme žádného rodiče Historie •již první studie genetického polymorfismu ® genetické techniky mohou rozlišit nejisté otcovství ® jedno z hlavních témat molekulární ekologie •chromosomal polymorphism •allozyme electrophoresis •DNA fingerprinting – overturn of existing paradigms in behavioural ecology (birds) •statistical techniques for single-locus polymorphism (allozymes) – departure from practice (i.e. DNA fingerprinting) • Historie 1 2 3 4 5 6 7 8 9 10 allozymes DNA fingerprinting, RAPD absence of theoretical models low variation multi-locus genotype Využití mikrosatelitů (dnes nejčastěji používané markery pro stanovení paternity) •Spojení teorie (statistické modely) a praxe (vysoce polymorfní single-locus markery) •Tandemová opakování krátkých motivů •Např. (CTTT)n nebo(CA)n •Vysoce polymorfní •Jednoduchá mendelovská dědičnost - kodominantní Alely se liší délkou – analýza dvou jedinců CTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTT CTTTCTTTCTTTCTTTC CTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTT CTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTT + - 1 2 Příklad analýzy jednoho lokusu – fragmentační analýza PCR produktů (kodominantní markery = jednoduchá mendelovská dědičnost) Msat paternita - barvy Genotyp (bp) Matka: 125/131 Otec: 131/134 Potomek 1: 125/134 Potomek 2: 131/137 Sledovaný otec mohl zplodit potomka 1, ale zcela jistě není otcem potomka 2 + - ? Proteinový fingerprinting •Hohol severní Andersson & Åhlund 2000 • •Proteiny z bílku vajec • •isoelectric focusing in immobilized pH gradients • •Vnitrodruhový hnízdní parazitismus (více něž polovina hnízd) • •Zanáší si zřejmě příbuzné samice • Bucephala clangula Statistical methods of parentage analysis • Exclusion (vyloučení) nekompatibilita mezi rodičem a potomkem ® zamítnutí rodičovství • • Categorical and fractional likelihood nelze vyloučit více jedinců – pravděpodobnost závisí na frekvencích alel v populaci • • Parental reconstruction z genotypů potomků z jedné rodiny rekonstrukce genotypu rodičů •srovnání se skupinou potenciálních rodičů nebo srovnání navzájem mezi potomky k určení „multiple paternity“ Určení paternity – exclusion principle Matka Otec Mládě 1 Mládě 2 Fragmentační analýza 3LOKCELE tools_1a Je možné analyzovat několik lokusů najednou. Princip vyloučení („exclusion principle“) Msat paternita - barvy Genotyp (bp) Matka: 125/131 Otec: 131/134 Potomek 1: 125/134 Potomek 2: 131/137 Sledovaný otec mohl zplodit potomka 1, ale zcela jistě není otcem potomka 2 + - ? Kolik je třeba lokusů? •Nutný počet není absolutní • •Dáno: üpočtem alel üheterozygotností üčetností alel üvelikostí vzorku ü •ideální - hypervariabilní lokusy stačí třeba jen 3 • •Tedy nejdříve je nutná pilotní studie dospělců v populaci • •→ výpočet Exclusion probability • • Nielsen et al. 2001 různý počet alel, různý počet možných rodičů Exclusion probability •Pravděpodobnost, že genotyp náhodně vybraného jedince nebude pasovat na potomka •Závisí na počtu lokusů a alel •Jiné hodnoty pro prvního a druhého rodiče •Dobré jsou hodnoty okolo 0,99 •Vypočítá například Cervus • •od v. 3.0 výše už počítá „the combined non-exclusion probabilities across all loci“ = these represent the average probability of not excluding a single randomly-chosen unrelated individual from parentage at one or more loci. (dobré jsou hodnoty pod 0,01) Exclusion - komplikace •Mendelian rules of inheritance: genotyping errors, null alleles, mutations ® false exclusions in strict conditions • „Null alleles“ a „genotyping error“ • lokus 1 lokus 2 • null alleles genotyping error • •Matka 100 150 300 350 • •Samec 1 100 100 300 367 • •Mládě 150 150 350 365 • • • • • • Samec 1 je vždy opravdovým otcem, ale jednoduchá „exclusion“ metoda ho vždy vyloučí Exclusion - komplikace •Mendelian rules of inheritance: genotyping errors, null alleles, mutations ® false exclusions in strict conditions •Extended family structure (příbuznost potenciálních rodičů) – nejhorší jsou sourozenci •Linked loci – decrease of variation •Loci on sex chromosomes •Problémy narůstají s rostoucím množstvím jedinců a lokusů •„Exclusion“ je velmi užitečná metoda např. v experimentech, kde jsou všichni rodiče předem známi a zgenotypováni • Př.: Faktory ovlivňující celoživotní fitness u hořavky duhové 10 x stejné uspořádání 10 x stejné uspořádání Reichard et al., Mol Ecol 2008, Reichard et al., Evolution 2009 - 3 páry hořavek v každé nádrži - přes 4000 embryí (odebírány v pravidelných intervalech, 5 mikrosatelitů - velikost těla - kondice - zbarvení duhovky - velikost gonád - parazitofauna Př.: Faktory ovlivňující celoživotní fitness u hořavky duhové Reichard et al., Mol Ecol 2008, Reichard et al., Evolution 2009 „oportunity for sexual selection“ „velikost těla u samců rozhoduje“ Nahloučená distribuce = selekce na velikost gonád Rovnoměrná distribuce = selekce na velikost těla Mutace Ibarguchi et al. 2004 •Obvykle se zanedbávají • •Mohou však být i poměrně časté – vysoká mutační rychlost mikrosatelitů (4.5 x 10-2 - 5.1 x 10-6 na lokus) • •Navíc častější u samců (2 až 6x) • •→ Vyloučení na základě jednoho lokusu může být špatně !!! • Řešení •Použít více lokusů • •Probability of resemblance (PR) Podívám se na ostatní lokusy a vypočítám s jakou pravděpodobností mohou alely sdílet dva nepříbuzní jedinci. • Uria lomvia •22% mimopárová paternita bez korekce na mutace • •7% mimopárová paternita s korekcí soulad s jinou použitou technikou Ibarguchi et al. 2004 Categorical x fractional likelihood •„likelihood“ skóre je vypočítáno na základě genotypů potomka a všech nevyloučených rodičovských genotypů (na základě frekvencí alel v populaci) • •Categorical l.: potomek jako jednotka přiřazen otci, biologicky validní • •Fractional l.: potomek přiřazen všem kompatibilním otcům. Statisticky přesnější a proto někdy výhodnější (např. lepší pro srovnání reprodukčního úspěchu různých kategorií samců) • Categorical allocation a)hledáme otce b)hledáme matku •c) hledáme oba rodiče •d) hledáme oba rodiče, ale neznáme pohlaví jednotlivých kandidátů • •... více viz program CERVUS 1) Fractional allocation •přiřadí určitou frakci (mezi 0 a 1) každého potomka všem nevyloučeným kandidátním rodičům • •the portion of an offspring allocated to a particular candidate parent is proportional to its likelihood of parenting the offspring compared to all other nonexcluded candidate parents Megaptera novaeangliae Nielsen et al. 2001 •Až 25 samců se pokouší pářit se samicí •Různé role samců - dominantní společník a vyzyvatel, sekundární společníci •6 mikrosatelitových lokusů •Fractional likelihood paternity method •Dominantní samci mají asi 3x více mláďat (je to ale neprůkazné) Categorical and fractional likelihood - nevýhody •male’s likelihood of paternity is positively related to his homozygosity (homozygotes have higher likelihood to have the same allele with the progeny) ® be carefull when estimating variance in reproductive success • •categorical – overestimate the reproductive success of homozygotes and underestimate that of heterozygotes • •fractional – a priori equal probability of paternity (required by the technique) underestimates variance in reproductive success • Genotypic (parental) reconstruction •rekonstrukce parentálních genotypů z velkého množství genotypů potomků (full- or halfsibs) •exhaustive algorithm •extremely computationally intensive •nejčastěji při analýzách paternity u ryb • Choosing the method •sampling of unknown parents (all, some, none) •sampling large groups of full- or halfsibs (yes, no) •a priori known parents (both, one, neither) • preferred technique no sufficient markers alternative technique choosing computer program Počítačové programy •Cervus, Newpat, Probmax, Kinship, Famoz, Papa, Parente, Patri, Gerud • •available online – see Jones & Ardren (2003) • •the most important differences – handling sources of difficulty (null alleles, mutations, scoring errors) Co je nejdůležitější? •Kvalitní data z terénu či experimentů! • •Nejlépe všichni dospělci z populace, behaviorální a jiná fenotypická data, mláďata přiřazena do rodin (matka + sociální partner) • •Špatná data = špatné výsledky Stane se, že nemohu získat vzorky všech dospělých jedinců z populace (potenciálních rodičů) nebo nelze říci, která mláďata jsou z jedné rodiny •Nemusí se zdařit najít rodiče všem potomkům •I tak ale mohu zjistit ledacos zajímavého (multiple paternity) • Příklad genotypy matky a embryí • lokus 1 lokus 2 • •Matka 100 150 300 350 • •1. embryo 100 115 300 320 •2. embryo 150 120 350 310 •3. embryo 120 100 350 365 • • • • • • Alely od matky • lokus 1 lokus 2 • •Matka 100 150 300 350 • •1. embryo 100 115 300 320 •2. embryo 150 120 350 310 •3. embryo 120 100 350 365 • • •Na druhém lokusu více než dvě alely, které nejsou od matky → více otců • • • ... nebo mutace – data nutno korigovat – viz Uria lomvia Apodemus agrarius mláďata až od tří samců AF „Multiple paternity“ závisí na sezóně (populační hustotě dostupných samců?) - v létě již prakticky neexistují vrhy zplozené jedním samcem - výrazné mezidruhové rozdíly v rámci rodu Apodemus Bryja et al., Behav Ecol Sociobiol 2008 Orangutan na Sumatře (paternity) •Nápadný dimorfismus (vystouplé tváře a hrdelní vaky u samců) • •Předpoklad: Samci bez nápadných znaků hormonálně suprimováni a nemnoží se • •Mikrosatelity – 50 % mláďat je od samců bez nápadných znaků → alternativní strategie Lepomis (maternity) •Samec (bourgeois male) hlídá hnízdo s jikrami • •L. marginatus, L. punctaus, L. auritus samec má v hnízdě jikry od několika samic • •L. macrochirus tři typy samců (alternativní reprodukční strategie): bourgeois (> 7 let, staví hnízda) mladí samci (vjíždějí do hnízd a vypouštějí spermie) nespárovaní staří samci (napodobují samice) 20% potomků není od bourgeois L. macrochirus Skladování spermií •Dny u savců týdny u ptáků nebo hmyzu měsíce u mloků roky u hadů a želv • •Chrysemys picta mikrosatelity → po 3 roky mláďata od stejného otce, opakované páření se stejným otcem je nepravděpodobné Chrysemys picta Varan – Zoo Liberec – mláďata po dvou letech od úmrtí samce ES Ptáci •Řada druhů považována za monogamní (pěvci) • •Ale u 75 % druhů ptáků mimopárová mláďata • •Skutečně monogamní pěvci – jen 14% (Phylloscopus) • • • • • Sialia sialis 8-35 % mimopárových mláďat Emberiza schoeniclus až 55% mimopárových mláďat Malurus cyaneus až 72% mimopárových mláďat Odhad % mimopárových mláďat •Dostatečný vzorek pro odhad • •Více než 200 Halichoerus grypus tuleň kuželozubý Wilmer et al. 1999 •Samci z centra kolonie úspěšnější (až 30x) •Samci reprodukčně aktivní 10 let i déle •Pro více než polovinu mláďat nenalezeni otci → role páření ve vodě • •Kolonie, dimorfismus→ polygynie •Dvě skotské kolonie North Rona a Isle of May •Vzorky ze zhruba deseti let •9 mikrosatelitových lokusů •IDENTITY, NEWPAT, CERVUS • Nejčastěji používané programy • •CERVUS •NEWPAT •FAMOZ •a další Jones AG, Ardren WR (2003) Methods of parentage analysis in natural populations. Molecular Ecology 12:2511-2523. •využívá výpočtu „likelihood“ • •přidává váhu jednotlivým alelám na základě jejich frekvence • •možnost nastavení míry genotypizační chyby, bere v úvahu nulové alely CERVUS CERVUS •PRVNÍ VERZE •Marshall TC, Slate J, Kruuk LEB, Pemberton JM (1998) Statistical confidence for likelihood-based paternity inference in natural populations. Molecular Ecology 7 (5) , 639–655 •AKTUÁLNÍ VERZE 3.0 •Kalinowski ST, Taper ML, Marshall TC (2007) Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Molecular Ecology 16: 1099-1006. • •DOWNLOAD - www.fieldgenetics.com •Možnost použití • • a) hledáme otce • b) hledáme matku • c) hledáme oba rodiče • d) hledáme oba rodiče, ale neznáme pohlaví jednotlivých kandidátů CERVUS Typický příklad: Maternita a paternita u netopýrů M M M M M M M F F F F F F F F F F F F F F J J J J J J J J J Pipistrellus nathusii – jižní Čechy (Jahelková, Hulva, Bryja) Omezení v programu CERVUS •pouze diploidní data • •pouze kodominantní lokusy (např. mikrosatelity) • •lokusy nesmí být ve vazbě (test v Genepop 4.0) • •jen malé odchylky od HW rovnováhy Vstupní soubor •2 možnosti • – vytvoření csv souboru v MS Excell – – konverze z formátu jiného programu (Genepop, Kinship, Genetix) Vytvoření csv souboru Input file – konverze z jiného formátu •Příklad – konverze z Genetix funkce CONCATENATE Uložit jako „Text oddělený tabulátory“ Cesta v Genetix: Fichier ® Importer ® vybrat soubor, označit Text avec séparateur pokračovat volbou Fichier ® Enregister sous ® Načtení dat do Cervus Otevření souboru: .CSV File ® View text file ® vybrat soubor Z Genetix Tools ® Convert genotype file ® Genetix to Cervus ® vybrat soubor (Source file) a vybrat umístění a název konvertovaného souboru (Save as) Analýza dat 1) Frekvence alel lépe počítat pouze pro dospělce, bez mláďat Analýza dat Locus k N HObs HExp PIC NE-1P NE-2P NE-PP NE-I NE-SI HW F(Null) loc1 4 6 0.667 0.773 0.652 0.721 0.555 0.385 0.141 0.431 NS +0.0167 loc2 5 6 0.833 0.818 0.708 0.658 0.481 0.300 0.105 0.401 * -0.0355 loc3 5 6 0.667 0.803 0.692 0.674 0.498 0.315 0.114 0.410 *** +0.0050 počet alel počet zgenotypovaných jedinců pozorovaná heterozygotnost očekávaná heterozygotnost míra informativnosti daného polymorfismu non-exlusion probabilities p pro test na nulové alely p pro test odchylky od HW rovnováhy Analýza dat 2) Simulace lépe počítat pouze pro dospělce, bez mláďat Analýza dat – výstup simulace Analýza dat •3) Vlastní určení paternity (a/nebo maternity) Kroky před a)vytvoření souboru s mláďaty b) b) vytvoření souboru s potencionálními otci (resp. matkami) Analýza dat 3) Vlastní určení paternity Analýza dat Analýza dat výběr zobrazení výsledků analýzy Výstup • před otevřením souboru s výstupem je třeba opět vyměnit čárky za středník • • mismatch – nesedící alely, nemožné zdědit po rodičích (může být ale způsoben nějakou chybou) • • non-exclusion probability – pravděpodobnost, že nepříbuzný kandidát není vyloučen jako rodič Výstup • • LOD skóre – přirozený logaritmus pravděpodobnosti - menší než 0 – spíše není rodičem - rovný 0 – může být a nemusí být rodičem - větší než 0 – může být rodičem - • hodnota Delta – rozdíl mezi dvěma nejbližšími LOD skóre. Její kritická hodnota je vypočítávána při simulacích • • znak * označuje otce/matky určené s 95% pravděpodobností • znak + označuje otce/matky určené s 80% pravděpodobností • znak – označuje nejpravděpodobnějšího otce/matku, který ale není určen jako rodič VÝSTUP Výstup PŘÍKLADY matka, otec, mládě, alela, LOD skóre…. PŘÍKLAD 1 •máme 5 mláďat, známe jejich matky • •existuje 10 potencionálních otců, každý z kandidátů může být otcem každého z mláďat • •soubory: genotypy, mláďata, otcové •máme 5 mláďat, neznáme jejich matky ani otce, zajímají nás jen matky • •existuje 10 potenciálních matek pro všechna mláďata • •soubory: jen genotypy2, ostatní soubory je třeba vytvořit PŘÍKLAD 2 •10 mláďat, neznáme ani otce, ani matku • •k dispozici genotypy 10 samců a 15 samic • •soubory – jen genotypy3, zbytek je třeba vytvořit PŘÍKLAD 3 •opět 5 mláďat, známe jejich matky • •pro každé mládě existují 4 kandidátní otcové (označení kand, číslo mláděte, písmeno a-d) • •soubory – pouze genotypy4, zbytek třeba vytvořit PŘÍKLAD 4