Zbornik gozdarstva in lesarstva 77, s. 43 - 60 GDK: 165.3:231:174.7 Picea abies (046) Prispelo / Recived: 26. 07. 2005 Izvirni znanstveni lanek Sprejeto / Accepted: 19. 09. 2005 Original scientific paper GENETSKI VIDIK NARAVNE OBNOVE SMREKOVEGA SESTOJA NA STALNI RAZISKOVALNI PLOSKVI ŠIJEC NA POKLJUKI Gregor BOŽI 1 Izvleek Namen raziskave je bil preveriti sposobnost naravnega obnavljanja genetskega vira Picea abies (L.) Karst. na objektu stalne raziskovalne ploskve Šijec na Pokljuki, Slovenija. Genetsko strukturo zrelega sestoja smreke in genetsko strukturo naravnega mladja v razlino velikih pomladitvenih jedrih smo preuili z elektroforetsko analizo 15 polimorfnih izoencimskih genskih lokusov desetih encimskih sistemov (GDH, F-EST, GOT, IDH, LAP, MDH, PGI, PGM, SKDH, 6-PGDH). Pod doloenimi pogoji je bila genetska struktura 101 dreves naravnega pomladka in 64 dreves sestoja, iz katerega je nastal, popolnoma primerljiva. Obnova smrekovega sestoja z naravnim pomlajevanjem zagotavlja trajno ohranjanje genskega vira smreke in je v primeru vejih pomladitvenih jeder tudi zagotovilo za neprekinjeno ohranjanje njegove prilagoditvene sposobnosti. Kljune besede: Picea abies (L.). Karst., naravno pomlajevanje, genetska struktura, izoencimi, ohranjanje gozdnih genskih virov THE GENETIC ASPECT OF THE SPRUCE STAND NATURAL REGENERATION IN THE PERMANENT FOREST RESEARCH PLOT ŠIJEC ON THE POKLJUKA PLATEAU Abstract The aim of the research was to check the ability of natural regeneration for sustainable conservation of Picea abies (L.) Karst. genetic resource in the permanent forest research plot Šijec on the Pokljuka plateau, Slovenia. The genetic structures of spruce stand and its natural regeneration centres of different size have been studied by means of electrophoretic analysis of ten enzyme systems (GDH, F-EST, GOT, IDH, LAP, MDH, PGI, PGM, SKDH, 6-PGDH) for 15 polymorphic isozyme gene loci. Under certain conditions, the genetic structures of 64 adult mother trees of the basic stand compared to 101 young spruce trees from the regeneration centres included in the study were closely similar. The regeneration of spruce stand by natural regeneration permanently conserves the spruce genetic resource and is, in the case of bigger regeneration centres, also an assurance of its continuous adaptability. Key words: Picea abies (L). Karst, natural regeneration, genetic structure, isozymes, forest gene conservation 1 dr. G. B., univ.dipl.inž.gozd., Gozdarski inštitut Slovenije, Vena pot 2, 1000 Ljubljana, SLO gregor.bozic@gozdis.si 44 Zbornik gozdarstva in lesarstva, 77 VSEBINA CONTENTS 1 UVOD ................................................................................................ 45 INTRODUCTION 2 MATERIAL IN METODE.............................................................. 46 MATERIAL AND METHODS 3 REZULTATI .................................................................................... 49 RESULTS 4 RAZPRAVA...................................................................................... 54 DISCUSSION 5 SUMMARY....................................................................................... 56 POVZETEK 6 VIRI ................................................................................................... 57 REFERENCES 7 ZAHVALA ........................................................................................ 60 ACKNOWLEDGEMENTS 45 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … 1 UVOD INTRODUCTION Zaradi dolge življenjske dobe so drevesa nenehno izpostavljena monim vestranskim in spreminjajoim se razmeram življenjskega okolja. Velik genetski potencial je pogoj za njihovo prilagodljivost spremembam življenjskega okolja, razmnoževanju in preživetju sploh. Glede na nepredvidljive spremembe življenjskega okolja v prihodnosti je bistvenega pomena, da to genetsko variabilnost zavarujemo in jo ohranimo. Poznavanje razlinih modelov in konceptov nam daje možnost za izbiro primernih metod izbora in formalnega varstva gozdnih genskih virov. eprav genske vire lahko varujemo s statinimi in dinaminimi postopki, ex situ ali in situ (KOSKI et al. 1997), prevladuje soglasje, da so za ohranitev prilagoditvene sposobnosti populacij dinamini postopki in situ primernejši od vseh drugih postopkov (FINKELDEY 1993). Uinkovito varstvo genetske prilagoditvene sposobnosti populacij pa poleg izbire ustreznih genskih virov, t.j. populacij, ki kažejo doloene ali posebej za prilagajanje pomembne genetske znake, terja tudi preverjanje uresnievanja in uspeha varstvenih ukrepov (GREGORIUS 1991). To med drugim pomeni tudi preverjanje sposobnosti posameznih varovanih populacij za uspešen prenos genetskih informacij z nosilcev na bodoe generacije. Z obnovo gozda se namre oblikujejo dednostne zasnove bodoega gozda, od katere je odvisna stabilnost bodoih sestojev. Preverjanje prilagoditvene sposobnosti populacij v primeru varstva gozdnih genskih virov in situ opravimo z raziskavo genetske strukture zrelega sestoja in njegovega potomstva v posameznih varovanih populacijah. Tako lahko še pravoasno odkrijemo nevarnosti in sprejmemo nadaljnje ukrepe za prepreevanje morebitne izgube genskih virov. To je še zlasti pomembno pri avtohtonih populacijah gozdnih drevesnih vrst, ki so rezultat dolgotrajne naravne selekcije in naravnega razvoja v danem okolju. Gozdni genski viri so namre lahko ogroženi kljub naravni obnovi sestojev, saj lahko antropogeni vplivi neposredno ali posredno vplivajo na biotsko pestrost na ravni ekosistema, vrstne raznolikosti in genetske variabilnosti znotraj vrste (MUHS 1997). Z našo raziskavo želimo predvsem osvetliti genetski vidik postopne obnove zrelega sestoja z naravnim pomlajevanjem in pri tem spoznati doloene lastnosti populacije, ki so pomembne za izdelavo strokovnih izhodiš za ohranjanje obstojee naravne genetske variabilnosti smreke in situ. Raziskava je specifino usmerjena v analizo uspešnosti prenosa genetskih informacij z nosilcev v zrelem sestoju na potomce z naravnim 46 Zbornik gozdarstva in lesarstva, 77 pomlajevanjem. Hipoteza, ki jo preverjamo, je, da naravni pomladek v primernem obsegu lahko zagotavlja ohranitev genetske pestrosti sestoja v pomlajevanju. Hkrati pa množinost mlade generacije v pomladitvenih jedrih še ne pomeni pojavljanja novih kombinacij genetskih zasnov in s tem tudi ne prispeva k veji variabilnosti v populaciji mladja. Cilj raziskave je ugotoviti genetsko strukturo naravnega pomladka, ki je nastajal postopno z veanjem pomladitvenih jeder, in jo primerjati z genetsko strukturo zrelega sestoja na tem rastišu. Ugotoviti želimo, ali so razlike v genetskih strukturah med zrelim sestojem in naravnim mladjem ter ali so razlike tudi med pomladitvenimi jedri. Zanima nas tudi, ali se frekvence genotipov na posameznih genskih lokusih vzorene populacije oziroma njihovih delnih populacij ujemajo s frekvencami genotipov, priakovanih s Hardy-Weinbergovim ravnotežjem, ki nastanejo s pomojo panmiktine reprodukcije. 2 MATERIAL IN METODE MATERIAL AND METHODS 2.1 RAZISKOVALNI OBJEKT RESEARCH OBJECT Genetski vidik postopne obnove sestoja smreke (Picea abies (L.) Karst.) smo preuevali na objektu stalne raziskovalne ploskve Gozdarskega inštituta Slovenije pri visokem barju Šijec na Pokljuki (N 46° 20′ 06″ / E 13° 59′ 42″, 1170 m n.m.). Objekt gradijo zrela drevesa, stara od 130 do 200 let (BOŽI / LEVANI 1998), in pomladek, v katerem smo izbrali 2 razlino veliki pomladitveni jedri (D in E), ki sta bili tudi predmet dosedanjih detajlnih inštitutskih raziskav v okviru projekta Rizosfera, L4-7402 (primerjaj SMOLEJ et al. 2000, DIACI / SMOLEJ / RUPEL 2000, KRAIGHER et al. 2000). Starost pomladka v manjšem jedru D je bila ocenjena na najmanj 15 let, v vejem jedru E na najmanj 20 let. 2.2 VZORENJE SAMPLING METHODOLOGY V vzorcu zrelega sestoja smo na površini 1 ha zajeli 64 nakljuno izbranih vitalnih, nadraslih oziroma sorastlih dreves razline starosti. Z vsakega od poskusnih dreves smo v zimskem obdobju odvzeli vejo s speimi popki, ki smo jih uporabili za ekstrakcijo 47 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … encimov. V vzorcu naravnega mladja pa smo zajeli dve pomladitveni jedri (manjše pomladitveno jedro »D« in veje pomladitveno jedro »E«). Glede na posamezno pomladitveno jedro smo enakomerno zajeli drevesa iz treh višinskih razredov, tj. nizka, srednje visoka in visoka. Z vsakega drevesa smo iz spodnjih vej in stranskih poganjkov odvzeli 5 speih popkov. Velikost vzorca je bila okoli 50 dreves na pomladitveno jedro. 2.3 GENETSKE ANALIZE GENETIC ANALYSES Vzorene subpopulacije smreke smo analizirali z metodo izoencimske horizontalne elektroforeze na škrobnem gelu. Ekstrakcijski pufer smo pripravili po RHODES (1977) (modif. KONNERT / MAURER 1995). Elektroforezo, barvanje gelov in oditavanje elektroforegramov smo opravljali po standardnih metodoloških postopkih za analizo smrekovih vzorcev (KONNERT / MAURER 1995). Analizirali smo 10 encimskih sistemov, ki so kodirani s 15 v analizi polimorfnimi genskimi lokusi, in sicer: glutamat dehidrogenaza (GDH, EC 1.4.1.2, GDH-A), fluorescenne esteraze (FEST, EC 3.1.1.1, FEST-B), glutamat oksalacetat transaminaza (GOT oziroma AAT, EC 2.6.1.1, GOT-A, GOT-B, GOT-C), izocitrat dehidrogenaza (IDH, EC 1.1.1.42, IDH-A, IDH-B), levcin aminopeptidaza (LAP, EC 3.4.11.1, LAP-B), malat dehidrogenaza (MDH, EC 1.1.1.37, MDH-B, MDH-C), fosfoglukoza izomeraza (PGI, EC 5.3.1.9, PGI-B), fosfoglukomutaza (PGM, EC 2.7.5.1, PGM-A), šikimat dehidrogenaza (SKDH, EC 1.1.1.25, SKDH-A), 6- fosfoglukonat dehidrogenaza (6-PGDH, EC 1.1.1.44, 6-PGDH-B, 6-PGDH-C). V oklepajih so podane okrajšave za encimske sisteme, kode E.C. ter analizirane genske lokuse. Kot polimorfen lokus smo upoštevali vsak genski lokus, na katerem smo ugotovili vsaj še en alel, ne glede na njegovo relativno pogostnost v preuevani populaciji. Izoencimske analize smo opravili v genetskem laboratoriju na Bayerische Amt fuer forstliche Saat- und Pflanzenzucht v Teisendorfu, Nemija, pod vodstvom dr. Monike Konnert. Frekvenne porazdelitve alelov na polimorfnih genskih lokusih smo opisali s 3 alelnimi profili, ki smo jih doloili po FINKELDEY (1993) kot: i) nizko stopnjo polimorfizma, e smo na posameznem genskem lokusu zasledili pogosti alel (alelna frekvenca > 80 %) in enega ali ve redkih alelov, ii) visoko stopnjo polimorfizma, e sta na posameznem genskem lokusu lahko opažena najmanj dva alela kot prevladujoa (alelna 48 Zbornik gozdarstva in lesarstva, 77 frekvenca > 20 %), iii) netipini profil, e alelno strukturo na posameznem genskem lokusu ne moremo uvrstiti v 1. ali 2. profil. Alelno pestrost subpopulacij smo ocenili s povprenim številom alelov na lokus (A/L). Za opis genetske variabilnosti smo uporabili merila alelne raznolikosti (ν; GREGORIUS 1978, 1987) in heterozigotnosti (Ha, Hc; GREGORIUS / KRAUHAUSEN / MUELLER- STARCK 1986). Stopnjo dejansko opažene heterozigotnosti (Ha) smo izraunali z relativnim deležem genskih lokusov, na katerih je osebek v populaciji heterozigoten. Pogojeno heterozigotnost (Hc) pa smo izraunali z razmerjem med ugotovljeno hetrozigotnostjo (Ha) in glede na ugotovljene frekvence najvejo možno hetrozigotnostjo (Hmax). Stopnjo genetske diferenciranosti med njimi smo podali z merilom genetske razdalje (d0; GREGORIUS 1974). Statistino primerjavo genetskih struktur med subpopulacijami pomladitvenih jeder oziroma zrelega sestoja z njegovim pomladkom smo napravili z G-testom homogenosti alelnih frekvenc po lokusih. Nielno hipotezo homogenosti alelnih porazdelitev med subpopulacijami zavraamo na nivoju tveganja α = 0,05. Za izraun posameznih parametrov genetske variabilnosti in statistinih primerjav smo uporabili program GSED (GILLET 1998). Za preverjanje, kako se v sestoju ali mladju ugotovljene frekvence genotipov pri posameznih lokusih ujemajo s frekvencami genotipov, priakovanimi po Hardy- Weinbergovem ravnotežju, smo uporabili natanni Hardy-Weinbergov verjetnostni test (HALDANE 1954, GUO / THOMPSON 1992) v programu GENEPOP 3.1d (RAYMOND / ROUSSET 1995). Nielna hipoteza, ki jo testiramo, je, da je zveza gamet nakljuna. Ocena verjetnosti (p) je izraunana z metodo polnega preštevanja po LOUIS / DEMPSTER (1987), ob prisotnosti ve kot 5 alelov na lokus pa z metodo markovske verige po GUO / THOMPSON (1992). Za tvorbo markovskih verig smo uporabili naslednje vhodne vrednosti: vhodno število (1000), število serij (100), število ponovitev na serijo (1000). Priakovano število homozigotov oziroma heterozigotov smo izraunali s korekcijo za male vzorce po LEVENE (1949). Hardy-Weibergovo strukturo preverjamo po lokusih loeno za vsako testno subpopulacijo posebej. Pri tem se ocena verjetnosti (p) nanaša izkljuno na nielno hipotezo (H0), ki je definirana kot nakljuno združevanje gamet; Fis = 0 (testiramo, ali vrednost Fis ponazarja znailni odklon od panmiksije). Ob zavrnitvi nielne hipoteze smo za sprejem alternativne 49 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … hipoteze (H1), ki je definirana kot primankljaj ali presežek števila heterozigotov (Fis > 0 ali Fis < 0) v subpopulaciji, uporabili strožji U-test. Koeficiente inbridinga (Fis) smo ocenjevali ob pomoi Wrightove F- statistike s cenilko po WEIR / COCKERHAM (1984) in s cenilko po ROBERTSON / HILL (1984), ki ima pri nielni hipotezi najmanjšo varianco in je tudi v zelo tesni povezavi z U-testom (ROUSSET / RAYMOND 1995). Ocena koeficienta inbridinga Fis po Wrightovi F-statistiki nakazuje, v kakšnem obsegu se dejanski heterozigotni deleži populacije razlikujejo od priakovanih heterozigotnih deležev, ki nastanejo pri panmiksini oplodnji (Hardy – Weinbergovo ravnotežje). 3 REZULTATI RESULTS Relativne alelne frekvence po posameznih lokusih za sestoj in mladje ter za posamini pomladitveni jedri so navedene v preglednici 1. Od 15 v analizi polimorfnih genskih lokusov 10 lokusov kaže nizko stopnjo alelnega polimorfizma. Ti so GDH-A, FEST-B, GOT-A, GOT-B, IDH-A, IDH-B, MDH-B, MDH-C, PGM-A in SKDH-A. Frekvenca glavnega alela vedno presega 80 %. Visoko stopnjo alelnega polimorfizma odsevajo 4 lokusi, in sicer: GOT-C, PGI-B, 6-PGDH-B in 6-PGDH-C. Lokus LAP-B ima netipien profil. Frekvenca glavnega alela je v razponu od 64,4 % do 72,4 %. Frekvence nad 10 % lahko dosegata še 2 alela, medtem ko se drugi aleli pojavljajo v nizkih frekvencah. Preglednica 1: Alelne frekvence v subpopulacijah zrelega sestoja smreke, njegovega mladja ter pomladitvenih jedrih, s prikazom za 15 polimorfnih genskih lokusov Table 1: Allelic frequencies in subpopulations of adult spruce stand, its young component and regeneration centres, by 15 polymorphic gene loci Lokus Locus Alel Allele Sestoj Stand n = 64 Jedro D Regeneration centre D n = 49 Jedro E Regeneration centre E n = 52 Mladje (D+E) Young component n = 101 IDH-A 2 3 0,016 0,984 0,000 1,000 0,000 1,000 0,000 1,000 IDH-B 2 3 1,000 0,000 0,990 0,010 1,000 0,000 0,995 0,005 MDH-B 1 2 3 0,008 0,984 0,008 0,000 1,000 0,000 0,000 1,000 0,000 0,000 1,000 0,000 (se nadaljuje …) 50 Zbornik gozdarstva in lesarstva, 77 Lokus Locus Alel Allele Sestoj Stand n = 64 Jedro D Regeneration centre D n = 49 Jedro E Regeneration centre E n = 52 Mladje (D+E) Young component n = 101 MDH-C 2 4 0,000 1,000 0,000 1,000 0,010 0,990 0,005 0,995 LAP-B 1 2 3 4 5 6 7 0,055 0,023 0,141 0,680 0,000 0,102 0,000 0,010 0,010 0,122 0,724 0,000 0,122 0,010 0,048 0,010 0,096 0,644 0,038 0,163 0,000 0,030 0,010 0,109 0,683 0,020 0,144 0,005 GOT-A 1 2 0,008 0,992 0,000 1,000 0,000 1,000 0,000 1,000 GOT-B 1 2 0,008 0,992 0,010 0,990 0,000 1,000 0,005 0,995 GOT-C 2 3 4 5 0,344 0,008 0,633 0,016 0,408 0,000 0,561 0,031 0,404 0,000 0,587 0,010 0,406 0,000 0,574 0,020 PGM-A 2 3 0,953 0,047 0,949 0,051 0,942 0,058 0,946 0,054 PGI-B 2 3 0,312 0,688 0,224 0,776 0,394 0,606 0,312 0,688 GDH-A 1 2 0,008 0,992 0,000 1,000 0,000 1,000 0,000 1,000 SKDH-A 1 2 3 5 6 7 0,023 0,047 0,883 0,016 0,008 0,023 0,051 0,010 0,929 0,010 0,000 0,000 0,029 0,010 0,942 0,010 0,010 0,000 0,040 0,010 0,936 0,010 0,005 0,000 6PGDH-B 1 2 3 5 0,000 0,633 0,008 0,359 0,031 0,714 0,000 0,255 0,000 0,654 0,000 0,346 0,015 0,683 0,000 0,302 6PGDH-C 1 2 4 5 0,008 0,602 0,000 0,391 0,000 0,429 0,020 0,551 0,000 0,577 0,000 0,423 0,000 0,505 0,010 0,485 FEST-B 1 2 0,047 0,953 0,031 0,969 0,087 0,913 0,059 0,941 n - število dreves v vzorcu / number of trees sampled V sestoju smo na 15 polimorfnih genskih lokusih odkrili 40 alelov, v mladju pa skupaj 37 razlinih alelov. Razlika v vrednosti parametra A/L za genski sklad izhaja iz števila 51 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … odkritih oziroma neodkritih redkih alelov (preglednica 2). V vzorcu “mladje” nismo odkrili 9 redkih alelov, ki smo jih sicer dobili v vzorcu “sestoj” (IDH-A2, MDH-B1, MDH-B3, GOT-A1, GOT-C3, GDH-A1, SKDH-A7, 6-PGDH-B3, 6-PGDH-C1). Po drugi strani pa v vzorcu “sestoj” nismo odkrili 6 redkih alelov, ki smo jih dobili v vzorcu “mladje” (IDH-B3, MDH-C2, LAP-B5, LAP-B7, 6-PGDH-B1, 6-PGDH-C4). Ker vse navedene alele uvršamo med redke (s frekvenco < 0,050) in veinoma splošne alele v preuevanih populacijah smreke v Sloveniji (BOŽI 2002), menimo, da jih nismo odkrili zaradi velikosti vzorca. Preglednica 2: Primerjava kazalnikov genetske variabilnosti zrelega sestoja smreke (ZR) z njegovim pomladkom (ML): število alelov na lokus (A/L), alelna raznolikost (ν), povprena heterozigotnost (Ha, Hc), genetska razdalja (d0) in vrednost G testa homogenosti genetskih struktur po lokusih, z nivojem znailnosti α = 0,05 (*) Table 2: Comparison of genetic variability parameters between the adult (ZR) and young component of the stand (ML): number of alleles per locus (A/L), allelic diversity (ν), average heterozygosity (Ha, Hc), genetic distance (d0) and the value of G–test of the homogeneity of genetic structures by loci; significance level α = 0,05 (*) Lokus Locus A/L ν Ha (%) Hc (%) d0 (%) G vrednost G value ZR ML ZR ML ZR ML ZR ML IDH-A 2 1 1,03 1,00 3,1 0,0 100,0 - 1,6 3,807 nz IDH-B 1 2 1,00 1,01 0,0 1,0 - 100,0 0,5 0,984 nz MDH-B 3 1 1,03 1,00 3,1 0,0 100,0 - 1,6 3,807 nz MDH-C 1 2 1,00 1,01 0,0 1,0 - 100,0 0,5 0,984 nz LAP-B 5 7 2,02 2,00 50,0 49,5 78,0 78,1 7,0 8,788 nz GOT-A 2 1 1,02 1,00 1,6 0,0 100,0 - 0,8 1,899 nz GOT-B 2 2 1,02 1,01 1,6 1,0 100,0 100,0 0,3 0,104 nz GOT-C 4 3 1,93 2,02 51,6 60,4 70,2 70,9 6,6 3,222 nz PGM-A 2 2 1,10 1,12 6,2 10,9 66,7 100,0 0,8 0,093 nz PGI-B 2 2 1,75 1,75 46,9 40,6 75,0 65,1 0,1 0,000 nz GDH-A 2 1 1,02 1,00 1,6 0,0 100,0 - 0,8 1,899 nz SKDH-A 6 5 1,28 1,14 23,4 8,9 100,0 69,2 6,9 11,207 * 6PGDH-B 3 3 1,89 1,79 46,9 43,6 63,8 68,7 6,5 5,927 nz 6PGDH-C 3 3 1,94 2,04 37,5 49,5 47,1 50,0 10,4 6,779 nz FEST-B 2 2 1,10 1,13 9,4 9,9 100,0 83,3 1,3 0,243 nz Genski sklad Gene pool 2,67 2,47 1,24 1,23 18,9 18,4 77,9 76,1 3,0 nz – ni znailno / non significant 52 Zbornik gozdarstva in lesarstva, 77 Alelna raznolikost osebkov v mladju se kljub doloenim razlikam na posameznih lokusih v povpreju bistveno ne razlikuje od alelne raznolikosti osebkov v zrelem sestoju. Enako velja tudi za stopnjo povprene dejanske in pogojene heterozigotnosti. Najveja možna stopnja pogojene heterozigotnosti je v subpopulaciji zrelega sestoja realizirana 77,9 %, v mladju pa 76,1 %. Da imata sestoj in njegovo mladje zelo podobni, skoraj identini alelni strukturi, nakazuje tudi podatek o povpreni alelni razdalji med njima (d0 = 3 %). Vzoreni populaciji se loita le v 3 % alelnem deležu preuevanega genskega sklada. Najvejo genetsko diferenciacijo smo dobili na lokusih 6-PGDH-C (10,4 %), LAP-B (7,0 %), SKDH-A (6,9 %), GOT-C (6,6 %) in 6-PGDH-B (6,5 %). Na enem od preuevanih lokusov je razlika v alelnih strukturah celo znailna, s stopnjo tveganja α = 0,05. Preglednica 3: Primerjava kazalnikov genetske variabilnosti v mladju na osnovi analize manjšega (ML-D) in vejega (ML-E) pomladitvenega jedra: število alelov na lokus (A/L), alelna raznolikost (ν), povprena heterozigotnost (Ha, Hc), genetska razdalja (d0) in vrednost G testa homogenosti genetskih struktur po lokusih, z nivojem znailnosti α = 0,05 ( *); α = 0,01 (**); Table 3: Comparison of genetic variability parameters in young components of the stand according to smaller (ML-D) and bigger (ML-E) regeneration centres: number of alleles per locus (A/L), allelic diversity (ν), average heterozygosity (Ha, Hc), genetic distance (d0) and the value of G–test of the homogeneity of genetic structures by loci; significance level α = 0,05 (*); α = 0,01 (**); Lokus Locus A/L ν Ha (%) Hc (%) d0 (%) G vrednost G value ML-D ML-E ML-D ML-E ML-D ML-E ML-D ML-E IDH-A 1 1 1,00 1,00 0,0 0,0 - - - - IDH-B 2 1 1,02 1,00 2,0 0,0 100,0 - 1,0 1,452 nz MDH-B 1 1 1,00 1,00 0,0 0,0 - - - - MDH-C 1 2 1,00 1,02 0,0 1,9 - 100,0 1,0 1,332 nz LAP-B 6 6 1,80 2,20 42,9 55,8 77,8 78,4 11,7 10,829 nz GOT-A 1 1 1,00 1,00 0,0 0,0 - - - - GOT-B 2 1 1,02 1,00 2,0 0,0 100,0 - 1,0 1,452 nz GOT-C 3 3 2,07 1,97 63,3 57,7 72,1 69,8 2,5 1,228 nz PGM-A 2 2 1,11 1,12 10,2 11,5 100,0 100,0 0,7 0,044 nz PGI-B 2 2 1,53 1,91 36,7 44,2 81,8 56,1 17,0 6,860 ** GDH-A 1 1 1,00 1,00 0,0 0,0 - - - - SKDH-A 4 5 1,16 1,13 6,1 11,5 42,9 100,0 2,3 1,973 nz 6PGDH-B 3 2 1,74 1,83 36,7 50,0 64,3 72,2 9,1 6,004 * 6PGDH-C 3 2 2,05 1,95 36,7 61,5 40,9 72,7 14,8 6,810 * FEST-B 2 2 1,06 1,19 6,1 13,5 100,0 77,8 5,6 2,961 nz Genski sklad Gene pool 2,27 2,13 1,21 1,24 16,2 20,5 72,5 78,4 4,5 nz – ni statistino znailno / not statistically significant 53 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … Vrednosti v preglednici 2 dovoljujejo ugotovitev, da se sestoj in naravni pomladek v genetski strukturi komaj razlikujeta, to velja za srednje vrednosti genetskih parametrov in tudi za veino genskih lokusov. Hipoteze, da množinost mlade generacije v pomladitvenih jedrih pomeni tudi pojavljanje novih kombinacij genetskih zasnov in prispeva k veji genetski variabilnosti v populaciji, v tem primeru ne moremo potrditi. Preglednica 4: Koeficient inbridinga Fis zrelega sestoja, pomladitvenih jeder in mladja po lokusih; krepki tisk oznauje znailni odmik (p < 0,05) frekvenc genotipov od Hardy-Weinbergovega ravnotežja, podrtane vrednosti oznaujejo znailni deficit heterozigotov. Table 4: Inbreeding coefficient Fis of adult stand, its young component of the stand and regeneration centres by loci. Bold values present the significant deviations (p < 0,05) of genotype frequencies from the Hardy-Weinberg equilibrium. Underlined values present the significance heterozygote deficiency. Lokus Locus Sestoj Stand Jedro D Regeneration centre D Jedro E Regeneration centre E Mladje (D+E) Young components of the stand (D+E) Fis (W&C) Fis (R&H) Fis (W&C) Fis (R&H) Fis (W&C Fis (R&H) Fis (W&C) Fis (R&H) IDH-A -0,008 -0,008 - - - - - - IDH-B - - - - - - - - MDH-B -0,004 -0,000 - - - - - - MDH-C - - - - - - - - LAP-B +0,017 +0,063 +0,047 +0,052 -0,013 -0,018 +0,014 +0,007 GOT-A - - - - - - - - GOT-B - - - - - - - - GOT-C -0,064 -0,029 -0,213 -0,130 -0,161 -0,086 -0,191 -0,111 PGM-A +0,308 +0,311 -0,043 -0,040 -0,052 -0,052 -0,053 -0,053 PGI-B -0,083 -0,084 -0,045 -0,045 +0,084 +0,084 +0,059 +0,060 GDH-A - - - - - - - - SKDH-A -0,078 -0,019 +0,553 +0,269 -0,030 -0,006 +0,279 +0,119 6PGDH-B +0,011 +0,003 +0,143 +0,083 -0,095 -0,096 +0,019 +0,010 6PGDH-C +0,235 +0,119 +0,292 +0,636 -0,252 -0,253 +0,033 +0,507 FEST-B -0,041 -0,042 -0,021 -0,021 +0,158 +0,160 +0,119 +0,120 rtica (-) pomeni, da bodisi ni podatka ali da je bil zabeležen samo en alel ali pa da sta bila zabeležena dva alela, vendar je bil eden od obeh zabeležen le v enem primerku. Fis (W&C): ocena koeficienta inbridinga po WEIR / COCKERHAM (1984). Fis (R&H): ocena koeficienta inbridinga po ROBERTSON / HILL (1984). »-» Means no data or one allele or two alleles presented with single occurrence only Fis (W&C): inbreeding coefficient estimation after WEIR / COCKERHAM (1984). Fis (R&H): inbreeding coefficient estimation after ROBERTSON / HILL (1984). 54 Zbornik gozdarstva in lesarstva, 77 Razlike med pomladitvenima jedroma so opazne na posameznih genskih lokusih in pri srednjih vrednostih. Te, z izjemo parametra genetske pestrosti A/L, vedno dosegajo veje vrednosti v pomladitvenem jedru iz veje sestojne vrzeli, ki se od sestoja tudi bistveno ne razlikuje. Pomladitveni jedri imata vejo stopnjo genetske diferenciranosti, kot e primerjamo zreli sestoj z vzorenim mladjem v celoti. V zrelem sestoju in pomladku se opažene frekvence genotipov pri veini lokusov ujemajo s frekvencami genotipov, priakovanimi po Hardy-Weinbergovem ravnotežju. Znailne odmike frekvenc genotipov od Hardy-Weinbergove strukture smo v sestoju ugotovili le na lokusu 6-PGDH-C, v mladju pa poleg lokusa 6-PGDH-C še na lokusu SKDH-A. Na istih lokusih se pojavlja tudi znailni primanjkljaj heterozigotov. Znailni primanjkljaj heterozigotov v mladju je nastal izkljuno zaradi vkljuitve manjšega pomladitvenega jedra (D) v vzorec, saj pri analizi vejega pomladitvenega jedra (E) teh odmikov nismo ve dobili. Pojava znailnih presežkov heterozigotov nismo ugotovili. 4 RAZPRAVA DISCUSSION Genetski vidik naravne obnove zrelega sestoja smreke (Picea abies (L.) Karst.) in situ smo preuevali na stalnem raziskovalnem objektu Gozdarskega inštituta Slovenije v Triglavskem narodnem parku na Pokljuki. Raziskave genetskih struktur naravnega pomladka je težko pojasniti, ker zajemajo selekcijo fertilnosti (reproduktivni sistem) in tudi razline oblike preživetvenih selekcij sejank (GREGORIUS et al. 1985). e obstaja selekcija, potem z neposredno primerjavo genetskih struktur zrelih sestojev z njihovim pomladkom ni mogoe loiti fertilnostne od preživetvene selekcije (SCHOLZ 1989). To dejstvo ne zmanjšuje pomena takih raziskav. Le z njimi lahko namre ugotovimo, ali in v kakšnem obsegu naravni pomladek zagotavlja ohranitev genetske pestrosti sestoja v pomlajevanju. Genetska primerjava zrelega sestoja z njegovim naravnim potomstvom v dveh razlino velikih pomladitvenih jedrih in starostjo osebkov od najmanj 15 in 20 let je pokazala, da je pod doloenimi pogoji genetska struktura naravnega pomladka in sestoja, iz katerega je 55 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … ta nastal, popolnoma primerljiva. Za obe subpopulaciji je v povpreju nakazano tudi podobno pomanjkanje heterozigotov. Redke razlike se nanašajo predvsem na alele, ki se v populacijah smreke na Pokljuki (BOŽI 2002) pojavljajo v zelo nizkih frekvencah. Znailne razlike v alelnih frekvennih porazdelitvah med subpopulacijama smo dobili samo na 1 od 15 lokusov. Hipoteze, da množinost mlade generacije v pomladitvenih jedrih pomeni tudi pojavljanje novih kombinacij genetskih zasnov in prispeva k veji variabilnosti v populaciji mladja, v našem primeru ne moremo potrditi. Rezultati tudi ne dovoljujejo nobenih jasnih zakljukov glede genetske selekcije. Ni nobenih alelov, katerih nosilci bi bili v prvih življenjskih letih v vladajoem okolju oškodovani ali obogateni. Znailne odmike od Hardy-Weinbergove strukture smo opazili v zrelem sestoju in v njegovem mladju. Pri zrelem sestoju smo dobili statistino znailni odmik od Hardy- Weinbergovega ravnotežja samo na lokusu 6-PGDH-C, v mladju pa tudi na lokusu SKDH-A. Na obeh lokusih je prišlo do odmika od Hardy-Weinbergove strukture izkljuno zaradi primanjkljaja heterozigotov. Pojava znailnih presežkov heterozigotov nismo ugotovili. Rezultati primerjave dveh razlino velikih pomladitvenih jeder so pokazali, da obstajajo med njima veje razlike, kot so med mladjem in zrelim sestojem. Med pomladitvenima jedroma so bile razlike v alelnih porazdelitvah alelov statistino znailne na 3 od 11 v mladju polimorfnih lokusih. Za veje pomladitveno jedro ugotavljamo vejo alelno raznolikost in vejo heterozigotno stopnjo pomladka ter nekoliko manjši presežek homozigotnih osebkov kot v manjšem pomladitvenem jedru. Znailne odmike od Hardy- Weinbergovega ravnotežja smo v mladju dobili le v manjšem, ne pa tudi v vejem pomladitvenem jedru. Vidni homozigotni presežek (oziroma primanjkljaj heterozigotov) v majhnem jedru pojasnjujemo z njegovim nastankom. Homozigoti bi lahko nastali s samoopraševanjem. Kaže se tendenca nenakljunega združevanja gamet in s tem odmik od panmiktine reprodukcije. Glede na preuevane izoencimske genske lokuse menimo, da genetska struktura zrelega sestoja in njegovega pomladka lahko zagotavlja trajno ohranitev genskih virov na rastišu smrekovega sestoja Šijec na Pokljuki. Naravno pomlajevanje zagotavlja celostno, naravno dano genetsko pestrost populacije in je le v primeru v vejih pomladitvenih 56 Zbornik gozdarstva in lesarstva, 77 jedrih tudi zagotovilo za njihovo neprekinjeno prilagoditveno sposobnost na spreminjajoe se razmere njihovega življenjskega okolja. Na osnovi rezultatov naše raziskave bi lahko pri varstvu gozdnih genskih virov smreke in situ upoštevali izhodiše, da mora biti varstvena enota izbrana tako, da lahko zagotavlja možnosti trajnega prenosa genskih virov na naslednjo generacijo z naravnim pomlajevanjem. 5 SUMMARY POVZETEK The genetic aspects of the natural regeneration of a spruce stand by the means of isoenzyme gene markers were studied. The aim of the research was to check the ability of stand natural regeneration for sustainable conservation of forest genetic resources. The studies were carried out at a permanent forest research plot Šijec on the Pokljuka plateau (N 46° 20′ 06″ / E 13° 59′ 42″, 1170 m a.s.l.). The genetic structures of the adult stand and its natural regeneration in two regeneration centres of different size, which have already been included in a detailed study of natural regeneration processes in this plot (DIACI / SMOLEJ / RUPEL 2000, KRAIGHER et al. 2000), were analysed. Isoenzyme analysis was preformed by means of horizontal starch gel electrophoresis on dormant buds extracts of sampled trees. Ten enzyme systems coded by 15 in our analysis of polymorphic gene loci, were utilised: GDH-A, FEST-B, GOT-A, GOT-B, GOT-C, IDH-A, IDH-B, LAP-B, MDH-B, MDH-C, PGI-B, PGM-A, SKDH-A, 6-PGDH-B, 6- PGDH-C. Spruce sub-populations were compared with respect to the average number of alleles per locus (A/L), the heterozygosity (Ha, Hc) and allele diversity (ν), and the genetic distance (d0) between them estimated. The statistical comparison of allele structures was made by a G-test of the homogeneity of frequencies by loci. In order to check the hypothesis that in the stand and its natural regeneration the frequencies of genotypes for different loci correspond to the frequencies of the genotypes in the Hardy- Weinberg equilibrium, the exact probability test was used. For the study of a potential over- or under-quantification of heterozygotes in the sub-population, the U-test was applied. The inbreeding coefficients (Fis) were estimated with Wright's F- statistics with 57 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … the evaluator after WEIR / COCKERHAM (1984) and with the evaluator after ROBERTSON / HILL (1984). All calculations were made with the programmes GSED (GILLET 1998) and GENEPOP 3.1d (RAYMOND / ROUSSET 1995). Results of genetic studies of natural regeneration of native spruce stand at the site of the permanent research plot Šijec on the Pokljuka plateau has shown that natural regeneration enables a sustainable conservation of gene resources. Under certain conditions, the genetic structures of the spruce regeneration and of the basic stand are closely similar. The regeneration of the stand by natural regeneration is, in the case of bigger regeneration centres, also an assurance of their continuous adaptability in the changing environmental conditions. 6 VIRI REFERENCES BOŽI, G., 2002. Genetske raziskave naravnih populacij smreke (Picea abies (L.) Karst.) v Sloveniji: doktorska disertacija. - Ljubljana, Univerza v Ljubljani, Biotehniška fakulteta, oddelek za gozdarstvo in obnovljive gozdne vire, 136 s. BOŽI, G. / LEVANI, T., 1998. Starost in morfološke znailnosti domnevno avtohtone smreke (Picea abies (L.) Karst.) na obmoju visokega barja Šijec na Pokljuki. V: Diaci, J. (ur.). Gorski gozd. Zbornik referatov. XIX. Gozdarski študijski dnevi, Logarska dolina, 26.-27. marec 1998. - Ljubljana, Univerza v Ljubljani, Biotehniška fakulteta, Oddelek za gozdarstvo in obnovljive gozdne vire, s. 243-254. DIACI, J. / SMOLEJ, I. / RUPEL, M., 2000. Raziskave svetlobnih razmer in zakonitosti pomlajevanja smreke na trajni raziskovalni ploskvi Šijec. V: Kraigher, H. (ur.), Smolej, I. (ur.). Rizosfera: raziskave gozdnih tal in rizosfere ter njihov vpliv na nekatere fiziološke parametre gozdnega drevja v izbranih gozdnih ekosistemih, sestojnih tipih in razvojnih fazah gozda. - Strokovna in znanstvena dela, 118, Gozdarski inštitut Slovenije, s. 222-243. FINKELDEY, R., 1993. Die Bedeutung allelischer Profile fuer die Konservierung genetischer Ressourcen bei Waldbaeumen. - Goett. Forstgenet. Bericht, 14, 176 s. 58 Zbornik gozdarstva in lesarstva, 77 GILLET E. M., 1998. GSED – Genetic Structures from Electrophoresis Data. Version 1.1, - Institut fuer Forstgenetik und Forstpflanzenzuechtung, Universitaet Goettingen, 48 s. GREGORIUS H. R., 1974. Genetischer Abstand zwischen Populationen. I. Zur Konzeption der genetischen Abstandsmessung. - Silvae Genetica, 23, s. 22-27. GREGORIUS, H. R., 1978. The concept of genetic diversity and its formal relationship to heterozygosity and genetic distance. - Math. Biosc., 41, s. 253-271. GREGORIUS, H. R. 1987. The relationship between the concepts of genetic diversity and differentiation. - Teor. Appl. Genet., 74, s. 397-401. GREGORIUS, H. R. 1991. Gene conservation and the preservation of adaptability. V: Seitz, A. / Loeschke, V. (ur.). Species Conservation. - A Population Biological Approach, s. 31- 47. GREGORIUS, H.R. / HATTEMER, H.H. / BERGMANN, F. / MUELLER-STARCK, G., 1985. Umwelbelastung und Anpassungsfähigkeit von Baumpopulationen. - Silvae Genetica, 34, 6, s. 230-241. GREGORIUS, H. R. / KRAUHAUSEN, J. / MUELLER-STARCK, G., 1986. Spatial and temporal genetic differentiation among the seed in stand of Fagus sylvatica L. - Heredity, 57, s. 255-262. GUO, W. / THOMPSON, E.A., 1992. Performing the exact test of Hardy-Weinberg proportions for multiple alleles. - Biometrics, 43, s. 805-811. HALDAN, J.B.S., 1954. An exact test for randomness of mating. - Journal of Genetics, 52, s. 631-635. KONNERT, M. / MAURER, W., 1995. Isozymic Investigations on Norway Spruce (Picea abies (L.) Karst.) and European Silver Fir (Abies alba Mill.). A Practical Guide to Separation Methods and Zymogram Evaluation, 79 s. KOSKI, V. / SKROPPA, T. / PAULE, L. / WOLF, H. / TUROK, J. 1997. Tecnical guidelines for genetic conservation of Norway spruce (Picea abies (L.) Karst.). – Rome, EUFORGEN, IPGRI, 42 s. KRAIGHER, H. / SIMONI, P. / SMOLEJ, I. / DIACI, J., 2000. Zakljuki, aplikacije in perspektive raziskav. V: Kraigher, H. / Smolej, I. (ur.). Rizosfera: raziskave gozdnih tal in rizosfere ter njihov vpliv na nekatere fiziološke parametre gozdnega drevja v izbranih gozdnih ekosistemih, sestojnih tipih in razvojnih fazah gozda. - Strokovna in znanstvena dela, 118, Gozdarski inštitut Slovenije, s. 246-262. 59 Boži, G.: Genetski vidiki naravne obnove smrekovega sestoja na … LEVENE, H., 1949. On a matching problem arising in genetics. - Ann. Math. Stat., 20, s. 91-94. LOUIS, E. J. / DEMPSTER, E. R., 1987. An exact test for Hardy-Weinberg and multiple alleles. - Biometrics, 48, s. 361-372. MUHS, H. J., 1997. Characterisation, testing and use of forest genetic resources. V: Jurc, M. / Hoevar, M. (ur.). Znanje za gozd. Zbornik ob 50. obletnici obstoja in delovanja Gozdarskega inštituta Slovenije. – Ljubljana, Gozdarski inštitut Slovenije, s. 331-340. RAYMOND, M. / ROUSSET, F., 1995. GENEPOP (Version1.2): Population genetics software for exact test and ecumenicism. - Journal of Heredity, 86, s. 248-249. ROBERTSON, A. / HILL, W. G., 1984. Deviations from Hardy-Weinberg proportions: sampling variances and use in estimation of inbreeding coefficients. - Genetics, 107, s. 713-718. RHODES, M. J. C., 1977. The extraction and purification of enzymes from plant tissues. - Proceedings of the Biochemical Society, 14, s. 254-248. ROUSSET, F. / RAYMOND, M., 1995. Testing hetrozygote excess and deficiency. - Genetics, 140, s. 1413-1419. SCHOLZ, F., 1989. Impotrance of the genetic structure in tree species for forest ecosystems under the influence of air pollutants. V: Urlich B. (ur.). Internationaler Kongress Waldschadensforschung: Wiessenstand und Perspektiven. Friedrichshafen, s. 478-497. SMOLEJ, I. / ATER, M. / URBANI, M. / SIMONI, P. / KUTNAR, L., 2000. Naravne razmere, preteklo gospodarjenje in stanje gozda na raziskovalnih ploskvah. V: Kraigher, H. / Smolej, I. (ur.). Rizosfera: raziskave gozdnih tal in rizosfere ter njihov vpliv na nekatere fiziološke parametre gozdnega drevja v izbranih gozdnih ekosistemih, sestojnih tipih in razvojnih fazah gozda. - Strokovna in znanstvena dela, 118, Gozdarski inštitut Slovenije, s. 11-32. WEIR, B. S. / COCKERHAM, C. C., 1984. Estimating F-statistic for the analysis of population structure. - Evolution, 38, s. 1358-1370. 60 Zbornik gozdarstva in lesarstva, 77 7 ZAHVALA ACKNOWLEDGEMENTS Prispevek je nastal v okviru Raziskovalne in Programske skupine za gozdno biologijo, ekologijo in tehnologijo (RS/PS0404-003) in raziskovalnega projekta »Razvoj molekularnih in biokemijskih baz podatkov v gozdarstvu« (L4-4450-0404), ki sta ga financirala Javna agencija za raziskovalno dejavnost RS in Ministrstvo za kmetijstvo, gozdarstvo in prehrano RS.