THE INFLUENCE OF WARMING, SITE CHARACTERESTICS, AND HOST PLANT ON ROOT-ASSOCIATED FUNGAL COMMUNITIES FROM ALEXANDRA FIORD IN THE CANADIAN HIGH ARCTIC by Kei E. Fujimura B.S., University o f Wisconsin, 1991 M.S., Oregon State University, 1999 DISSERTATION SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in NATURAL RESOURCES AND ENVIRONMENTAL STUDIES THE UNIVERSITY OF NORTHERN BRITISH COLUMBIA August 2005 © Kei E. Fujimura, 2005 1^1 Library and Archives Canada Bibliothèque et Archives Canada Published Heritage Branch Direction du Patrimoine de l'édition 395 W ellington Street Ottawa ON K 1A 0N 4 Canada 395, rue W ellington Ottawa ON K 1A 0N 4 Canada Your file Votre référence ISBN: 978-0-494-28458-2 Our file Notre référence ISBN: 978-0-494-28458-2 NOTICE: The author has granted a non­ exclusive license allowing Library and Archives Canada to reproduce, publish, archive, preserve, conserve, communicate to the public by telecommunication or on the Internet, loan, distribute and sell theses worldwide, for commercial or non­ commercial purposes, in microform, paper, electronic and/or any other formats. AVIS: L'auteur a accordé une licence non exclusive permettant à la Bibliothèque et Archives Canada de reproduire, publier, archiver, sauvegarder, conserver, transmettre au public par télécommunication ou par l'Internet, prêter, distribuer et vendre des thèses partout dans le monde, à des fins commerciales ou autres, sur support microforme, papier, électronique et/ou autres formats. The author retains copyright ownership and moral rights in this thesis. Neither the thesis nor substantial extracts from it may be printed or otherwise reproduced without the author's permission. L'auteur conserve la propriété du droit d'auteur et des droits moraux qui protège cette thèse. Ni la thèse ni des extraits substantiels de celle-ci ne doivent être imprimés ou autrement reproduits sans son autorisation. In compliance with the Canadian Privacy Act some supporting forms may have been removed from this thesis. Conformément à la loi canadienne sur la protection de la vie privée, quelques formulaires secondaires ont été enlevés de cette thèse. While these forms may be included in the document page count, their removal does not represent any loss of content from the thesis. Bien que ces formulaires aient inclus dans la pagination, il n'y aura aucun contenu manquant. Canada Abstract Arctic systems are expected to be impacted earlier and more severely by global warming than temperate ecosystems. However, much o f the research on the impact o f warming on arctic ecosystems has centered on plant communities. One objective o f this thesis was to examine how passive warming would impact the root-associated fungal community at Alexandra Fiord, Nunavut. The root-associated fungal community consists mostly o f mycorrhizal, dark-septate and hyaline-septate fungi, which are considered important mutualists in arctic ecosystems. The objective was to compare the fungal community from plots warmed by open-top chambers to ambient plots, using two methodologies: 1) fungal DNA extraeted directly from root tips with terminal restriction fragment length polymorphisms (T-RFLPs) used to estimate variation, and 2) fungal cultures isolated from root tips to which PCR-RFLP techniques were applied to assess variation. T-RFLPs were used to examine the root-associated fungal community on Salix arctica. Differences between the communities were analyzed using canonical correspondence analysis (CCA). Genotype diversity was tested using a 2-way, 2-stage, nested ANOVA. Warming did not significantly change genotype cumulative frequency or diversity o f the root-associated fungal community, but cumulative frequency tended to increase on the warmed plots. Genotype richness was significantly different according to site, which was correlated with differences in soil chemistry. Again site, not warming, was the main factor that distinguished the root-associated fungal community of Salix arctica, Saxifraga oppositifolia, Cassiope tetragona, and Dryas integrifolia based on fungal cultures. Warming did not have a detectable impact on cumulative frequency and diversity, based on CCA and a nested, 3-way ANOVA. Fungal cultures were identified based on sequence analysis and morphology. Phialocephala fortinii ii was the most frequently identified taxon, but almost half o f the fungal isolates remained unknown. The root-associated fungal community was examined along a glacier forefront characterized by a directional, non-replacement primary plant succession pattern. CCA was used to examine genotype frequency; linear regressions were used to test for changes of cumulative frequency and diversity as succession advanced. The fungal community on only one o f the host plants increased in frequency and richness as succession advanced. The darkand hyaline-septate endophyte communities were distinct on different host plants, providing evidence for host specificity and higher diversity than previously reported. Ill Table o f Contents Abstract ...............................................................................................................................................ii Table o f Contents......................................................................................................................... iv List of Tables............................................................................................................................. viii List of Figures.............................................................................................................................. xi Acknowledgements................................................................................................................... xiii 1. Introduction................................................................................................................................. 14 A. Rationale........................................................................................................................... 14 B. Literature review .............................................................................................................. 15 1. Concepts in community ecology....................................................................................15 a) 2. Definition o f community................................................................................................. 15 Biodiversity...................................................................................................................... 18 a) Definitions......................................................................................................................... 18 b) Assessing biodiversity....................................................................................................21 c) Rank/abundance curves..................................................................................................24 3. Methods to assess root-fungal com m unities................................................................28 4. Concepts in mycorrhizal fungal community ecology..................................................32 a) Mycorrhizal fungal community structure o f arctic and alpine System s.................. 36 b) Mycorrhizal fungal succession in relation to plant community succession 39 5. Global climate change: effects on aboveground plant community structure 41 6. Impact o f global climate change on root fungal communities.................................. 48 C. Research objectives and hypotheses............................................................................. 49 D. Thesis organization......................................................................................................... 52 iv IL Impact o f warming on the frequency and diversity o f the root-associated fungal community on Salix arctica from the Canadian High A rctic............................................... 53 A. Introduction.........................................................................................................................53 B. Materials and m ethods...................................................................................................... 54 1. Study S ite ............................................................................................................................54 2. V egetation.......................................................................................................................... 55 3. Field collection...................................................................................................................56 4. Sampling from roots.......................................................................................................... 57 5. DNA extraction and ITS-T-RFLP analysis.................................................................... 58 6. Matching root tips with sporocarps................................................................................. 59 7. Data analysis...................................................................................................................... 60 C. Results................................................................................................................................. 64 1. Soil properties.................................................................................................................... 64 2. Genotype frequency.......................................................................................................... 65 3. D iversity..............................................................................................................................66 4. Community composition.................................................................................................. 70 D. Discussion...........................................................................................................................71 III. How host plant, warming, and site affect the culturable root-associated fungal community from the Canadian High Arctic..................................................................................................82 A. Introduction.........................................................................................................................82 B. Materials and m ethods...................................................................................................... 83 1. Field collection...................................................................................................................83 2. Sampling from roots.......................................................................................................... 84 3. DNA extraction and ITS-RFLP analysis........................................................................ 84 4. Data analysis...................................................................................................................... 85 5. DNA sequencing............................................................................................................... 87 C. Results................................................................................................................................. 91 1. IT S-R FL Ps......................................................................................................................... 91 2. DNA sequencing............................................................................................................... 96 3. Culture m orphology..........................................................................................................97 4. Identifieations based on moleeular techniques and m orphology.............................. 102 D. Discussion.........................................................................................................................105 IV. The root-associated fungal community along a directional, non-replacement succession chronosequence........................................................................................................................ 112 A. Introduetion.......................................................................................................................112 B. Materials and m ethods.................................................................................................... 113 1. Study site ...........................................................................................................................113 2. Field eollection.................................................................................................................114 3. Sampling from roots........................................................................................................ 114 4. Root m icroscopy............................................................................................................. 115 5. DNA Extraction and ITS-T-RFLP analysis.................................................................115 6. Matching root tips with sporocarps...............................................................................117 7. Statistical analysis........................................................................................................... 117 C. Results............................................................................................................................... 119 D. Discussion.........................................................................................................................126 V. Synthesis o f results...................................................................................................................133 VI VI. Conclusion............................................................................................................................... 138 Vil. Literature cited................................................................................................................ 144 V lll. Appendix for chapter 2 ......................................................................................................162 VIII. Appendix for chapter 3 .................................................................................................. 182 IX. Appendix for chapter 4 ............................................................................................................ 199 X. Appendix for chapter 5............................................................................................................. 214 Vll List o f Tables Table 1.1 Summary o f experiments o f warming on arctic plants used for present study........ 47 Table 2.1 2-way, 2-stage nested ANOVA for the effects o f site, treatment, treatment replicate, and plant specimen replicate on genotype frequency...........................................67 Table 2.2 Summary o f Tokeshi’s model used to describe the rank abundance curves at each plot. 1.............................................................................................................................................70 Table 3.1 Number o f successful samples isolated from root tips based on RFLPs.................. 92 Table 3.2 Species richness o f cultures according to host plant, site, and treatment based on R FLPs.......................................................................................................................................... 95 Table 3.3 Species richness o f cultures according to host plant, site, and treatment based on identification made morphologically........................................................................................95 Table 4.1. Number o f successful extractions according to the number o f attempts at extraction and total number of genotypes for each host plant.............................................121 Table A2.1 2-way, nested ANOVA table for effects o f site, treatment, and treatment replicate on soil properties...................................................................................................................... 159 Table A2.2 2-way ANOVA table for the effects o f site, treatment, and treatment replicate on genotype richness.................................................................................................................... 162 Table A2.3 2-way ANOVA table for the effects o f site, treatment, and treatment replicate on genotype evenness.................................................................................................................... 163 Table A2.4 Relative abundance of genotypes from control to OTC for each primer-enzyme combination................................................................................................................................164 Table A3.1 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Phialocephala.............................................182 V lll Table A3.2 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Helotiales.................................................. 184 Table A3.3 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - unknown.................................................... 186 Table A3.4 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Hymenoscyphus......................................... 193 Table A3.5 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - 194 Table A3.6 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Dothideales................................................ 195 Table A3.7 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Phoma..........................................................195 Table A3.8 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Cryptosporiopsis........................................ 196 Table A3.9 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Cadophora.................................................. 197 Table A3.10 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Colispora.....................................................197 Table A 3.11 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Agaricales................................................... 197 Table A3.12 Comparison o f identification (based on molecular techniques) with corresponding morphological identification -Ceratobasidium........................................... 198 Table A4.1 Genotype richness for host plants at chronosequence zone...................................199 IX Table A4.2 Sorensen’s Quantitative Index between host plants and chronosequenee zones200 Table A5.1 Comparison o f LS means o f soil chemistry between succession study (control plot only) and warming studies (T-RFLP-OTC and culture studies)................................217 X List o f Figures Fig. 1.1 Example o f how T-RFLP is obtained in comparison to RFLP..................................... 33 Fig. 2.1. CCA on warming treatments and site with soil properties as biplots....................... 67 Fig. 3.1 Species cumulative frequency for Salix arctica and Saxifraga oppositifolia on different sites............................................................................................................................... 92 Fig. 3.2 Species richness for Salix arctica and Saxifraga oppositifolia on different sites. .. 93 Fig. 3.3 NMS for host plant and warming treatments per site based on phi distances............ 94 Fig. 3.4 Neighbor-joining best tree based on ITS and the Dothideales.................................. 96 Fig. 3.5 Neighbor-joining best tree based on ITS and Helotiales........................................... 98 Fig. 3.6 Neighbor-joining best tree based on LSU and the Agaricales.................................. 99 Fig. 3.7 Neighbor-joining best tree based on LSU and the Ceratobasidiomycetes..............99 Fig. 3.8 Neighbor-joining tree based on LSU and ascomycetes.............................................. 100 Fig 3.9 Abundance o f Phialocephala fortinii based on number of morphologically identified cultures differentiated by host plant and treatment...............................................................101 Fig 4.1 Canonical Correspondence Analysis for effects o f host plant and chronosequence zone............................................................................................................................................. 122 Fig 4.2 Cluster analysis using Ward’s method based on Euclidean distances.......................... 123 Fig. A2.1 NMS for site and warming treatment effect..............................................................165 Fig. A2.2 Cluster analysis o f plots, treatments and plant specimen based on frequency of genotypes....................................................................................................................................169 Fig. A2.3 Least square means for soil properties.........................................................................171 Fig. A2.4 Least square means for genotype frequency for each treatment per plot................173 Fig. A2.5 Rank abundance curves for lowland sites................................................................... 174 XI Figure A2.6 Rank abundance curve for highland dolomitic site ................ 176 Fig. A2.7 Rank abundance curves for highland granitic site..................................................... 178 Fig. A2.8. Least square means o f genotype richness for each treatment per site...................180 Fig. A2.9. Least square means o f genotype evenness measured by Hurlbert’s PIE index.. 181 Fig. A4.1 NMS of host plant and chronosequence zones..........................................................206 Fig A4.2 LS means o f soil property for each chronosequence zone.........................................207 Fig. A4.3 Simple regression for genotype richness for Luzula confusa...................................209 Fig. A4.4 Simple regression o f genotype cumulative frequency of Luzula confusa.............. 210 Fig. A4.5 Rank abundance curve for each chronosequence zone for Luzula confusa 211 Fig. A4.6 Rank abundance curve for each chronosequence zone for Papaver lapponicum. 212 Fig. A4.7 Rank abundance curve for each chronosequence zone for Salix arctica............... 213 Fig. A4.8 Rank abundance curve for each chronosequence zone for Saxifraga oppositifolia ............................................................................................................................................................214 Fig. A4.9 Rank abundance curve for each chronosequence zone for Cassiope tetragona... 215 Fig. A4.10 Rank abundance curve for each chronosequence zone for Dryas integrifolia. ..216 Xll Acknowledgements I would like to thank my advisor, Keith Egger, for all his support, guidance, and patience. 1 would like to thank my supervisory committee members, Randy Currah, Art Fredeen, Kathy Lewis, and Hugues Massicotte for all their support and help. Greg Henry was generous with his field sites and field consultation, and Cecilia Alstrom-Rapaport was invaluable for statistical guidance. 1 am indebted to Randy Currah who graciously accepted my cultures and to Melissa Day for the laborious work o f identifying all the cultures morphologically. 1 would also like to thank Clive Owens o f the Ministry o f Forests Analytical Research Laboratory (Victoria, BC) for processing my soil samples, and Lito Arocena and Yuriko Yano for their insightful soil data interpretation. 1 am grateful to Brent Murray for all his advice in the lab and the following people for their laboratory assistance (without their help I would still be writing my thesis!): Tamara Bereck - exhaustive isolations from root tips Jennifer Catherall - started extractions and RFLPs for cultures Susan Gibson - started sequencing o f cultures, RFLPs Sara Jamieson - preparing roots for extractions Edward Kraay - database creation and technical support Linda Rehaume - help finishing sequences from cultures and helping in lab Rachel Botting, Ld Kraay, Susan & Garth Frizzell, Susan Gibson, Nabla Kennedy, Brian Milakovic, Tim Phaneuf, Jennifer Psyllakis, Linda Rankin, Susan Robertson, Laura Ryser, Linda Tackaberry, Andrew Walker, and members o f the Lgger/Massicotte/Murray lab were wonderful for their insights and camaraderie. I would finally like to thank my parents, Bob and Shigeko Fujimura, for all their support, encouragement, and guidance. X lll I. Introduction A. Rationale The root-associated fungal community includes mycorrhizal fungi, dark and hyaline septate fungi, and possibly parasitic or pathogenic fungi. O f these, the mycorrhizal and darkand hyaline-septate fungi are the most abundant members o f this community. Mycorrhizal fungi are known to have mutualistic relationships with vascular plants, and are important components o f most ecosystems. Dark and hyaline septate fungi have been reported to be both pathogenic and mutualistic. Their role as mutualists is hypothesized to be greater in arctic ecosystems with the absence o f arbuscular mycorrhizae (Bledsoe et al. 1990). Global warming is an important source o f disturbance of arctic ecosystems. The primary effect o f global warming is the increase in mean temperature. Global warming is also significant because o f its secondary effects; it has been linked to the increase o f sea levels, hurricane occurrence, fire and insect outbreak in coniferous forests, and species extinction. These problems are compounded by the release o f greenhouse gases. Arctic environments are opportune ecosystems to examine ecological questions about the root-associated fungal community. These environments have low plant species diversity, which simplifies the complexity o f examining the fungal community found on roots. In addition, arctic ecosystems are expected to be impacted more severely by global warming, and in advance o f other ecosystems. By simulating warming in arctic environments using open-top chambers, insights can be gained into how warming may impact the root-associated fungal community in other environments. Warming indirectly affects arctic ecosystems by causing glaciers to recede. This provides an opportunity to examine changes in the belowground community in response to a 14 unique form of primary plant succession that occurs in the high arctic, directional nonrcplaccmcnt succession, where host plants arc not replaced as succession proceeds so diversity increases along the chronosequence. Alexandra Fiord, Nunavut provides an excellent opportunity to study the effects of climate change on the root-associated fungal community. The biology and autecology of plants in this area have been studied extensively, and passive warming experiments are part o f the International Tundra Experiment (ITEX), which was created to monitor the effects of warming in arctic regions. Alexandra Fiord also hosts a unique type o f primary plant succession, which provides a natural experimental design for studying how the rootassociated fungal community on a common suite o f host plants responds to an increase in plant diversity along a chronosequence. B. Literature review 1. Concepts in community ecology a) Definition of community There has been considerable debate among ecologists regarding the conceptual definition o f “community” . Wilson (1991), for example, questions whether plant communities are really integrated, discrete entities. He argues that plant communities do not exist unless the definition is delimited by a list o f criteria including assembly rules, niche limitation, discreteness, discontinuity, and integratedness (Wilson 1991, Palmer and White 1994). Looijen and van Andel (1999) asserted that the problem with the definition of community is that the term is too ambiguous: 1) it can be applied to different levels o f taxa. 15 2) no objective boundaries can be made, and 3) communities are heterogeneous with respect to species composition. They suggest the following definition: “community may be defined as a set o f individuals o f two or more species that occur in the intersection o f areas occupied by populations o f these species” (Looijen and van Andel 1999). This limits the definition of communities to be used only for coexisting species belonging to a single taxonomic group, such as phyla or class, that has a static boundary (Looijen and van Andel 1999). Parker (2001) refutes Looijen and van Andel’s definition because o f scale limitations and unidentified assumptions. He argues that the scale limitation leads to ambiguity or conflict with respect to what organisms belong to a given community (Parker 2001). He identifies three assumptions from Looijen and van Andel’s model which are often violated in community ecological studies: 1) there must be a unique underlying process, 2) there must be consistency o f processes among replicates and, 3) there must be independence from other communities. Instead, Parker’s definition combines concepts from Brand and Parker (1995) and Pickett et al. (1992): ‘communities are continuous in time and space, and processes underlie composition and dynamics’. His definition includes a conceptual model where the community focuses on a single individual and its interactions with other members o f the community (i.e. consumers, symbionts, pathogens, mutualists, and competitors); each individual o f the community has its own set o f interactions (Parker 2001). Some argue to forgo the conceptual definition (McCune and Grace 2002, Palmer and White 1994) and use an operational definition (Palmer and White 1994), such as a “collection o f organisms found at a specific place and time” (McCune and Grace 2002). This operational definition is similar to Parker’s definition (2001) but does not include his conceptual model, which implies that interactions are necessary in a community. The 16 operational model assumes that variation in species composition is random spatially and temporally (McCune and Grace 2002). With the operational definition, conceptual or theoretical implications are circumvented (McCune and Grace 2002, Palmer and White 1994). To avoid the ambiguous use o f the term community, Fauth et al. (1996) proposed a more restricted set o f definitions. ‘Taxa’ is reserved for phylogenetically-related species regardless o f where they occur; ‘communities’ include all species co-oecurring in one place (i.e. corresponds to the operational definition Palmer and White (1994)); and ‘guilds’ are sets o f species exploiting the same resources (Fauth et al. 1996). Other concepts are used to describe overlapping combinations o f these three terms. ‘Assemblages’ are phylogenetically- restricted groups that occur in a community (i.e. overlap between taxa and communities). When guilds and taxa overlap, then a group o f related species exploit the same resource. ‘Local guilds’ are formed from the intersection of communities and guilds and are groups o f species that share a common resource and occur in the same community (Fauth et al. 1996). When all three overlap, then the term ‘ensemble’ is used; this is a taxonomically restricted group o f species that exploit the same resources, and are located in one place (Fauth et al. 1996). The term that perhaps best applies to this study is the ‘ensemble’ according to Fauth (1996), because this study incorporates phylogenetic, community and guild perspectives. However, operationally, I will use community in the sense o f Palmer and White (1994), which corresponds to Fauth’s definition. 17 2. Biodiversity a) Definitions Biodiversity is a loosely applied term, and has been used to describe diversity from the genetic level to the biome level (Hooper et al. 2005). Biodiversity can be described as the number o f different genotypes, species, ecosystem types, etc. and includes the evenness o f tbeir distribution (Hooper et al. 2005). Species richness is only part o f this definition and refers to the number o f taxonomic units (usually species or genotypes). Ecosystem function is the effect o f the activities o f organisms on the physical and chemical processes o f an environment. A functioning ecosystem is characterized by these processes (Naeem et al. 1999). According to Naeem et al. (1999), ecosystem functioning can be measured by quantifying rates o f movement, such as nutrient transportation, or by measuring growth or production, such as plant stem growth or seed production. Disturbance used to be limited to events that were massively destructive and rare (Rykiel 1985). This definition is no longer acceptable because disturbances are not always catastrophic and can be recurring events in ecosystem. Rykiel (1985) attempted to formulate a general definition by defining disturbance as a physical force, agent, or process that causes a perturbation in an ecological component or system. Disturbance can be abiotic or biotic and: 1) can cause destruction, where biomass is quantitatively reduced; 2) can cause discomposition, where certain populations are eliminated, reduced, added, or expanded; 3) can cause interference, where matter, energy and/or processes are hindered; or 4) can be caused by suppression, where natural disturbances are prevented (Rykiel 1985). The outcome of a disturbance is perturbation, which is a deviation of values that are used to describe the properties o f the ecological component or system (Rykiel 1985). Rykiel’s 18 definition assumes that reference conditions must be known in order to understand disturbance (Pickett et al. 1989). White and Jentsch (2001) argued that one general definition is unachievable for disturbance. They argued that disturbance should incorporate four different areas: 1) variation in disturbance events, which would include the timing and intensity o f the disturbance; 2) variation in the disturbance effects within an ecosystem, which would cover spatial and temporal variation; 3) variation in ecosystem response by including differences in biota and physical environments; this would include the rates of response and species adaptation would be included; and 4) inferences o f scale o f observations and measurements, which are affected by observations, sampling, and analysis by the researcher. These four topics that they argue should be included are covered in the definition proposed by Pickett et al. (1989). The definition by Pickett et al (1989) consists o f identifying the object disturbed, distinguishing what is and is not disturbed, and recognizing the minimal level o f hierarchical organization. The following levels are their recommended hierarchical organization: individual, population, community, ecosystem, and landscape. Disturbance, which is defined as an external force o f a given level, would affect the structure, function, and attributes of that level. For example, at the community level, structural disturbance would include effects on vertical and horizontal patterns, species composition, or functional groups; functional disturbance would include effects on resource levels, competition, or mutualistic interactions; and attribute disturbance would include effects on coexistence, evenness, or dominance. At the ecosystem level, structural disturbance would include effects on functional groups, functional disturbance would include effects on fluxes in the ecosystem, and attribute 19 disturbance would include effects on resistance and resilience of the ecosystem. In this hierarchical organization, temporal and spatial scales are also incorporated and are defined within the context o f each level. Temporal and spatial effects that occur on a broader scale would be called ‘disturbance regime’ such as the fire regime in grasslands, where the reoecurrenee o f the disturbance is every few years. Stability refers to how an ecosystem or community responds to disturbance. For ecosystems, stability would apply to populations or communities and their abiotic environment, such as analysis o f nutrient dynamics (Barbour et al. 1999). For communities, stability is measured by determining how the community composition and diversity responds to disturbance (Barbour et al. 1999). Stability has two components, resilience and resistance, and the overall response to disturbance is determined by the interaction between the two (Barbour et al. 1999). On the ecosystem level, resilience is the ‘ability o f an ecosystem to return to predisturbance conditions’ (Barbour et al. 1999), which may take a long time. Resistance is the ‘ability o f an ecosystem to resist changes in response to disturbance’ (Barbour et al. 1999). These terms can also have community definitions. Community resilience is when the community returns to the same species composition after a disturbance, and community resistance is where species composition does not change due to disturbance (Tokeshi 1999). However, there are problems with the community resistance/resilience concepts. One problem is determining how much change in the community composition is needed before the community is considered ‘disturbed’. Another problem is determining whether it is necessary for the community to return to its ‘original’ composition for stability to return (Tokeshi 1999). Also, community composition changes over time without disturbances (Tokeshi 1999), so these terms may be difficult to apply. 20 Although ecosystems may become unstable, they can still continue to function. How well an ecosystem responds to disturbance depends on its resilience. Processes may be retarded, but if the ecosystem is resilient, then these processes can return to pre-disturbance conditions. b) Assessing biodiversity Assessing biodiversity has become an important issue because o f the increased rate of loss of diversity due to anthropogenic activities. One argument used to support conservation is that preservation o f biodiversity will maintain ecosystem functioning. Maintenance o f biodiversity has become a surrogate for ecosystem function (Naeem 2002). For example, Naeem et al. (1995) found that higher diversity correlated with an increase in community respiration, productivity and nutrient retention in a mesocosm study, and therefore alteration o f biodiversity can affect ecological processes. Tilman et al. (1997a) examined how plant species diversity, functional diversity, and functional composition affect plant productivity, %N in plants, total N in plants, soil NH4, soil NO3, and light penetration. They found that functional diversity, but not plant diversity, significantly impacts these functional variables by positively affecting plant productivity and total N and negatively affecting soil NO3, soil NH4, plant %N, and light penetration. They also found that many species in monocultures have comparatively less biomass than when they are found in multifunctional group plots, which supports the hypothesis that higher diversity increases ecosystem productivity. They concluded that: 1) functional composition and diversity are significant determinants in grassland ecosystem processes and 2) not all plant species are equal, so the loss o f one species may be more deleterious than another. 21 Higher species richness alone is not sufficient to explain the impact o f biodiversity on ecosystem functioning. When examining competition in a resource model for plants, the variance o f the model is explained more by species identification than species richness (Tilman et al. 1997a). Likewise, Hooper and Vitousek (1997) found that the identification o f the functional groups explains more variance than species richness, and that species and combinations o f species, rather than species richness, control yields and nutrients. Another issue eoneeming how biodiversity affects ecosystem functioning is complementary effects versus selection effects (Hooper and Vitousek 1997, Cardinale et al. 2002). Complementary effect theory attempts to explain how resource use by organisms affects ecosystem processes (Cardinale et al. 2002, Loreau and Hector 2001). According to this theory, species diversity can increase while avoiding competition, and species are able to co-exist, especially in environments with limiting resources, by: 1) partitioning resources, where each species can use nutrients, water, or other resources differently instead o f all species competing for or using the same resources; or 2) niche differentiation, where different species avoid using the same resources as other organisms in time and/or space (McKane et al. 2002). Loreau et al. (2001) found that feeding performance by caddisfly larvae improves in the presence o f other taxa. They concluded that the increase in species diversity o f other aquatic arthropods leads to interspecific facilitation. In contrast, selection effect theory is applied when species diversity is correlated with the probability that a dominant species uses most o f the available resources (Cardinale et al. 2002), so the formation o f the community is heavily dependent on these dominant species. In a study where a mathematical model was used to test for complementary effects and selection effects on monocultures and mixed species o f grasses, results explained by the 22 selection effect theory were not as reliable as by the complementary effect theory (Loreau et al. 2001). In complementary effect theory, the performance o f communities can go beyond the additive effects o f individual species (Loreau et al. 2001). Schwartz et al. (2000) criticized studies that link greater biodiversity with an increase o f ecosystem productivity. In both observational and experimental studies, conflicting results have been reported, with some studies having negative or no results and others that are variable through time and space (Schwartz et al. 2000). In addition, other problems with experimental studies include; 1) hidden manipulations such as weeding, which would change the diversity and composition o f the experiment; 2) addition o f species to poor environments, which probably would not occur in nature; and 3) extrapolation o f results to the whole ecosystem when only one trophic level has been included (Schwartz et al. 2000). With theoretical models, the role o f rare species may be missed in stabilizing ecosystems. Also, models may assess stability and function on the wrong scale; models may apply only at a local rather than an ecosystem scale. Schwartz et al. (2000) concluded that where the relationship between biodiversity and ecosystem function are positive, the relationship is not linear as studies suggest, but that the function saturates after a few species or functional groups, creating more o f an asymptotic relationship. For the 23 observational and experimental studies they examined where they could graph biodiversity against ecosystem function, they found that 60% o f these studies produced the asymptotic graph. Although the relationship between biodiversity and ecosystem function is an important question, doubts about this relationship are prevalent because of conflicting results in finding this asymptotic relationship, and also because the number o f species or functional groups to fulfill the functions of an ecosystem is not known. 23 However, as Loreau et al. (2001) suggested, biodiversity may not be as important in the maintenance o f an ecosystem as it is in helping to facilitate changes in the environment. Bengtsson (1998) argued that biodiversity is not mechanically linked to ecosystem functions. Simple measures o f species richness assume that all species are equal in their function. This is unlikely to be true; therefore measuring diversity is pointless unless the function of species is known. He contends that knowing the species and their functions will explain the processes and stability o f an ecosystem. However knowing the functions o f all species is currently impossible. Even though all the functions of species are not known, linking biodiversity to ecosystem functioning remains important because more diverse ecosystems may include redundant species that could fulfill ecosystem functions when dominant species are lost. This would increase the probability o f withstanding or rebounding from disturbances. c) Rank/abundance curves The shape o f rank/abundance plots (a.k.a. dominance/diversity curve or Whittaker plots) is used to determine which species abundance model best describes the data (Magurran 2004). These rank/abundance plots are helpful in that: 1) different patterns o f species richness are easily shown, 2) the relative abundance o f species-poor communities is easily seen, and 3) emphasis is placed on the differences in evenness for contrasting communities (Magurran 2004). Species abundance models can be categorized into two main groups - biological and statistical (Tokeshi 1999). The biological models are also called niche apportionment models and include the following; geometric series, broken stick, MacArthur fraction, dominance pre-emption, random fraction, dominance decay, random assortment, composite, and power 24 fraction (Tokeshi 1999). These biological models are based on the assumption that species divide the niche space among species that live in a community in different ways (Magurran 2004). These models have been criticized for possibly being too simplistic and confusing in terms o f bow the niches are apportioned, but they can be valuable tools for understanding niche differentiation (Magurran 2004). Depending on which model fits the shape o f the rank/abundance plot, interpretations can be made about bow the niche space is divided in the community. Rank/abundance plots that fit a geometric series model often describe communities that are species-poor, such as those found in harsh environments or in early stages o f succession (Magurran 2004). Those communities fitting the Mac Arthur's broken stick model (or random niche boundary hypothesis) are interpreted as having their species competing equally for one resource. However, Mac Arthur’s broken stick model assumes that the niche space is partitioned simultaneously, which probably does not happen in nature (Magurran 2004). Tokeshi developed a set o f niche apportionment models that forgo the assumption of simultaneous niche partitioning o f the broken stick model (Magurran 2004). Two o f his models examine extreme cases when the least or most abundant species are invaded by new species (Tokeshi 1999). The dominance pre-emption model is where new species invade niche spaces occupied by the least abundant species in an existing community (niche fragmentation), or alternatively where a new species takes approximately half o f a new niche space (niche filling) (Tokeshi 1999). In these cases, the dominant species remain so. The dominance decay model is a model o f the other extreme where the largest niche space, instead of the smallest, is appropriated by new species (Tokeshi 1999). 25 Tokeshi has three models that examine how niche apportionment occurs when new species invade all potential niches and not just the ones occupied by the least and most abundant species. The Mac Arthur fraction model is similar to the broken stick model, but it assumes that the niche spaces are invaded sequentially rather than simultaneously; however, the same conclusions can be drawn from both models (Tokeshi 1999). In this model, the probability o f a community being invaded depends on species abundance or niche size, so niche space o f more abundant species will likely be invaded before less abundant species. This model implies a uniform distribution and may be applicable only to small communities with related species (Magurran 2004). In the random fraction model, all species have the same probability o f being invaded by a new species, so the abundance o f species does not influence the chances o f being selected (Tokeshi 1999). This model fits situations where new species compete for niche spaces randomly over an existing niche that is already occupied by an assemblage o f species (Magurran 2004). Magurran (2004) finds this model to be innovative with a wide range o f applications. The power fraction model is similar to the random fraction model, but it is used for species rich assemblages because most o f the niche apportionment models are applicable to communities with small species assemblages (Magurran 2004). The random assortment model assumes that there is no relationship, or a weak one, between niehe apportionment and species abundance (Magurran 2004), so the abundances of species are independent o f each other (Tokeshi 1999). This may be used for situations where communities are in a state a flux from major environmental changes and competition is not limited by species abundance (Magurran 2004). 26 The composite model is achieved by taking two or more o f these niche apportionment models into account to describe how a niche is divided. Tokeshi realized that using only one model may be too simplistic for a community; however, knowing where to set the boundary between more and less abundant species may be problematic (Magurran 2004). Statistical models were initially created so researchers could objectively compare species abundance between communities (Magurran 2004). Even though some o f these models have been labeled as statistical, ecological implications have been drawn from these statistical models. For example, the log normal model is a statistical model, but the ecological implication is that it explains situations where new species come to a niche in a random order rather than in fixed intervals such as in the geometric series model (Magurran 2004). The log normal model has been found to fit many datasets and is commonly used (Magurran 1988). However this model has been criticized because it requires a large number o f species so the log normal distribution may be a mathematical artifact o f a large sample size and so may have few biological implications (Magurran 1988, Tokeshi 1999). These models are useful in assessing plant communities based on species abundance distributions. Arctic ecosystems are considered harsh environments and, as expected, plant communities in the arctic fit the geometric series model (McKane et al. 2002). These species abundance models can be used to assess the root-associated fungal community and test if these communities fit models similar to their plant counterparts. These models will also be helpful in assessing how the root-associated fungal community responds to the direct, non­ replacement succession (see below) found on Alexandra Fiord. 27 3. Methods to assess root-fungal communities The techniques used in this study to examine the root-associated fungal communities included morphotyping, PCR-RPLP, T-RFLP, and DNA sequencing. The term ‘rootassociated fungal communities’ will be used because some o f the techniques do not differentiate mycorrhizal, endophytic, and pathogenic fungi. Because the techniques used are commonly applied for examining mycorrhizal fungal communities, much o f the review will be based on this group o f fungi. Sporocarp collections and morphotyping have been commonly used to assess mycorrhizal fungal communities. However, sporocarp production was found to be an inaccurate estimate o f composition o f ectomycorrhizal fungi found on root tips below ground (Gardes and Bruns 1996a, Dahlberg et al. 1997, Jonsson et al. 1999), and corresponds to only approximately 20% (Jonsson et al. 1999) to 30% (Dahlberg et al. 1997) o f the belowground ectomycorrhizal fungi. Often the most abundant sporocarps do not correspond with the most abundant mycorrhizae (Gardes and Bruns 1996a). Ectomycorrhizal fungal species may rarely, or never, fruit, or may produce small or hypogeous fruiting bodies that are missed in surveys, which could potentially lead to inaccurate estimation of mycorrhizal fungal species found on root tips. Ectomycorrhizal communities described by morphotyping is limiting in that at least half o f the species are completely unknown (Gardes and Bruns 1996a), especially if the morphotype comes from the field (Karen et al. 1997). Morphotyping also requires a high level o f skill (Kârén et al. 1997), and so often takes more time to learn than molecular based techniques (Dahlberg 2001). Better results in identification through morphotyping require more phenotypic characteristics, but then fewer samples can be examined (Horton and Bruns 2001). The efficiency of distinguishing mycorrhizal taxa improves with RPLP analysis over 28 morphotyping, circumventing phenotypic plasticity, where multiple species may be grouped as the same morphotype (Horton and Bruns 2001). Jonsson et al. (1999) distinguished 20 morphotypes from 7152 mycorrhizae but found 42 RFLP-types from 212 root tips that successfully amplified. PCR-RFLP has been a useful tool in researching mycorrhizal fungal communities. The internally transcribed spacer (ITS) region o f the nuclear-encoded ribosomal RNA (nrDNA) gene repeat is often used because; 1) it is readily amplified with fungal-specific primers, allowing amplification o f fungal DNA from mixed genomes, such as plant and fungal DNA in ectomycorrhizal root tips (Gardes and Bruns 1996a) and 2) it is divergent enough for identifying species within a genus (White et al. 1990). The nrDNA is often used for fungal studies because DNA sequences are highly conserved among organisms, variability is high between species and minimal within a species (Egger 1995), and the ribosomal repeat is a multi-copy gene (Gardes and Bruns 1993), making it easier to amplify. Studies involving PCR-RFLP o f the ITS region have focused on comparing above and belowground fungal species composition (Gardes and Bruns 1996a, Dahlberg et al. 1997, Jonsson et al. 1999); describing differences in composition due to treatment or changing environments (Horton et al. 1999, Erland et al. 1999, Kemaghan 2001); or describing changes in composition due to succession (Nara et al. 2003). In addition to the confirmed lack o f correlation between above and belowground fungal composition, these studies find that a few widespread species generally account for most o f the abundance o f mycorrhizal fungi (Gardes and Bruns 1996a, Dahlberg et al. 1997, Jonsson et al. 1999, Horton et al. 1999, Erland et al. 1999, Kemaghan 2001) and that spatial variation is a large determinant for species composition, at least for Swedish forests after a low intensity fire (Dahlberg 2001). 29 Terminal Restriction Fragment Length Polymorphism (T-RFLP) analysis is a relatively rapid and accurate, PCR-based tool to identify taxa (Avaniss-Aghajani et al. 1996, Martfnez-Murcia et al. 1995). It has been used for examining microbial communities found in sludge (Marsh et al. 1998, Liu et al. 1997), termite guts (Liu et al. 1997), aquifer sand from groundwater (Liu et al. 1997), and for determining the effects o f temperature on the microbial community in rice fields (Chin et al. 1999). A few more ectomycorrhizal fungal community studies have been based on T-RFLP analysis, such as analyzing how the increase of CO 2 would change the mycorrhizal fungal community (Klamer et al. 2002) and determining the soil vertical distribution o f ectomycorrhizal fungi (Dickie and Koide 2002). T-RFLPs are helpful when organisms have indistinct morphologies (Avaniss-Aghajani et al. 1994), and this technique avoids creating time-consuming cloning o f organisms, which may not work for all species (Bruce 1997). T-RFLPs are similar to restriction length polymorphisms (Clement et al. 1998) in that they are both used to characterize the ITS region o f the nrDNA for differentiating taxa in describing community diversity. However, the protocol for T-RFLPs differs from RFLPs in that each primer is fluorescently labeled with a different dye, and fragments are separated on a 6 % polyacrylamide gel rather than an agarose gel (see Fig. 1). The fluorescence allows the samples to be detected by automated DNA sequencers/fragment analyzers and because o f the polyacrylamide gel has higher resolution (Avaniss-Aghajani et al. 1996), allowing detection o f fragments that are only 1 or 2 base pairs different (Totsch et al. 1995). Once the ITS region is amplified with the fluorescent dye-labeled primers, restriction endonuclease enzymes are used to digest the region, as with the RFLP methodology. Instead o f visually detecting multiple bands as with RFLP analysis, only one fragment, the terminal fragment, is 30 detected because the laser only detects the digested product with the fluorescent dye (Fig. 1). If both primers are labeled, two sets o f results are generated: one from the forward primer and the other from the reverse primer. Multiple bands in PCR products, which indicate the presence o f more than one organism, are problematic in RFLP analysis because assigning fragments to their respective fungus is difficult or impossible. In RFLPs, individual taxa are indicated by a unique pattern o f multiple bands, but when multiple taxa o f fungi are on the root, fragment patterns become too complex for analysis. T-RFLPs circumvent this problem by detecting only terminal fragments, so in theory, each fragment should represent a unique taxon. Several restriction enzymes may be needed to differentiate taxa that share restriction sites, but multiple restriction enzymes are used for RFLP analyses as well. Clement et al. (1998) list potential problems with T-RFLPs: 1) PCR primers may differentially amplify certain species, therefore, measurements of relative abundance in a community may not be accurate; 2 ) unequal relative abundance may occur due to different optimum annealing temperatures for different species; 3) evidence o f fragments may be limited by electrophoresis technology; and 4) accurate community analysis needs multiple digestive enzymes. The problems listed can also be applied to RFLP analysis as well, and cannot be resolved without more advanced technologies that reduce the number o f samples that are analyzed. T-RFLPs is a valuable tool in mycorrhizal fungal community analyses. Like bacterial systems, mycorrhizal fungi are often morphologieally indistinct. In addition, this method circumvents the phenotypic plasticity o f the mycorrhizal fungi on different hosts. Unlike RFLP analyses, this tool can detect and distinguish multiple mycorrhizal fungi on the same root tip. 31 Although RFLP and T-RFLP analyses are powerful tools for assessing communities, they are limited in identifying taxa. One way to identify taxa using RFLPs and T-RFLPs is to compare fragments with a database that already exists. However, comparing RFLP fragments with those in databases created by other researchers is often not feasible because restriction enzymes and primers are not standardized. DNA sequencing is important in filling this gap, especially when the goal is to identify unknown taxa. Some reasons why sequencing is more successful are: 1) there is a central database (GenBank) where scientists deposit their sequences, and 2 ) sequences are not restricted by choice o f endonucleases or primers; as long as the unknown sequence has an overlapping segment in GenBank, then identification to at least order is plausible. The caveat with comparing sequences from GenBank is that submitted sequences are not checked for accuracy, so the identification of the sequences may not be reliable (Bridge et al. 2003). 4. Concepts in mycorrhizal fungal community ecology Mycorrhizal fungal community ecologists often have the triple task o f describing and interpreting the fungal community structure as well as extrapolating their results to the plant community. Some topics that are addressed in mycorrhizal fungal community ecology are complementary effects (e.g. Koide 2000, Perry et al. 1989), community structure (e.g. Horton and Bruns 1998, Dahlberg et al. 1997, Gardes and Bruns 1996a), and the role o f mycorrhizae in plant communities’ resistance to change (e.g. Horton and Bruns 1998). Mycorrhizal fungal community ecologists have approached the issue o f complementary effects by examining facilitation and niche differentiation (e.g. Dickie and Koide 2002, Helm et al. 1999, Titus and del Moral 1998). Facilitation is an important concept in mycorrhizal community ecology because o f the guild concept, where the diversity o f the mycorrhizal 32 Fluorescent-dyelabeled primer Fluorescent-dyelabeled primer V 18s small subunit 1r ITS 1 5.8s region ITS 2 I Fluorescent tag Cy 5.0 28s large subunit PCR Reaction Fluorescent tag Cy 5.5 ITSl - 5.8s - 1TS2 regions ▼ 1r ITS region cut with an endonuclease enzyme Terminal end that would be fluorescently tagged, e.g. with Cy 5.0, for T-RFLP analysis Appear as on an RFLP gel Terminal fragment shown on automated sequencer Fig. 1.1 Example o f how T-RFLP is obtained in comparison to RFLP 33 fungal community and the plant community can stabilize their plant-soil ecosystem after disturbance or stress (Perry et al. 1989). The benefits o f linkages between plants by mycorrhizal fungi, are generally expressed by bow plant communities may profit rather than by bow the fungal communities may profit. Possible benefits for plants include: 1) seedlings may benefit by linking into a larger fungal network via fungal bypbae (Newman 1988); 2) interplant exchange o f nutrients (Newman 1988); 3) interspecific competition between plants may be altered if nutrients are received from one central network (Newman 1988); 4) competition between plants may be reduced (Newman 1988); 5) nutrients from dying plants may pass directly to living plants (Newman 1988); 6 ) stabilization o f succession patterns because some fungi eould associate with both pre- and post-disturbance plant hosts (Horton and Bruns 1998); and 7) improved plant survival (Trappe and Luoma 1992). Fungal linkages between plants have been demonstrated both in the laboratory and in the field. Phosphorus transfer between Pinus sylvestris and Finns conforta via Suillus bovinus was found in vitro using ^^P, in which the labeled phosphorus did not move only to the plants but throughout the whole fungal network (Finlay and Read 1986). Simard et al. (1997) showed a net transfer o f carbon between Be tula papyrifera and Pseudotsuga menziesii in the field. By shading P. menziesii seedlings and not B. papyrifera, they found that P. menziesii seedlings are sinks for carbon, and that carbon could transfer from a sink to a souree or along a nutrient gradient. They also found that earbon exebange occurs between the two ectomycorrhizal plants {B. papyrifera and P. menziesii), but not with the arbuscular mycorrhizal plant Thujaplicata. Horton et al. (1999) suggested that linkages between ectomycorrhizal and arbutoid plants allow outplanted P. menziesii seedlings to survive in an 34 arbutoid stand o f Arctostaphylos. They conclude that linkages probably do not exist between eeto- and arbuscular mycorrhizal because outplanted P. menziesii seedlings died in arbuscular stands of Adenostoma. Koide (2000) offers two strategies for the role o f complementarity in root colonization by arbuscular mycorrhizae. One strategy is where the fungi are complementary to each other. Although this would allow fungi to coexist on the same root, it does not explain why some antagonism happens between fungi (Koide 2000). The other strategy he proposes is that the function o f the fungus is complementary to those o f the plant, which would lead to redundant species on roots and may explain why some root-fungus relationships have a high level o f speeifieity. The latter is similar to the application o f coexistence theory to ectomycorrhizae, where plants select more beneficial fungi by lengthening or shortening time o f root tip mortality (Hoeksema and Kummel 2003). These concepts may apply to ectomycorrhizal fungi as well, but competition (Wu et al. 1999), different life strategies, such as colonization rate and life spans (Hoeksema and Kummel 2003), must also be considered. Another form o f complementarity is niche differentiation. Although niche differentiation is one o f the oldest explanations for biodiversity, it has not been tested in many ectomycorrhizal fungal communities (Bruns 1995). Recently, niche differentiation was used to explain the vertical distribution o f ectomycorrhizal fungi. Dickie and Koide (2002) used cluster analysis and species diversity measures to differentiate six spatial patterns for fungi. The outcome o f the distribution suggested that niche differentiation explained the vertical distribution o f ectomycorrhizal fungi. 35 Only a few papers link mycorrhizal fungal diversity to ecosystem functioning. Van der Heijden et al. (1998) showed that higher species richness o f arbuscular mycorrhizal fungi increased plant biodiversity, improved plant productivity, lengthened hyphal growth in the soil, increased phosphorus absorption o f plants, and decreased the amount o f phosphorus in the soil. They concluded that there is probably a feedback loop where the plant benefits from the increased amount of phosphorus, and the fungi prosper due to increased carbon, indicated by more hyphal growth. Baxter and Dighton (2001) examined if higher diversity o f ectomycorrhizal fungi would affect plant growth and nutrient acquisition. They concluded that ectomycorrhizal species diversity is more influential in plant biomass and nutrient uptake than the species composition or rate o f colonization. Although a pioneering paper in linking mycorrhizal diversity to ecosystem function, this study examined a community with low diversity o f mycorrhizal fungi (Leake 2001). This may limit applications due to its simplicity o f being an in vitro study. Another confounding factor included using peat and vermiculite for the growth media, which adds excess nutrients. As a result o f these short-comings, it may be premature to draw conclusions about the link between ectomycorrhizal fungi diversity and ecosystem functioning (Leake 2001). a) Mycorrhizal fungal community structure o f arctic and alpine systems Examining mycorrhizal fungal communities in alpine and arctic systems is preferable because confounding factors due to direct human contact, such as logging or fire suppression, that hamper studies in temperate forests have minimal impacts upon these arctic and alpine systems (Trappe 1988). Another advantage is that these systems are thought to be relatively simple in comparison to mycorrhizal fungal communities o f temperate systems (Read 1993). 36 Alpine and arctic systems share similar environmental stresses such as short growing seasons (Trappe 1988, Haselwandter 1987), low air and soil temperatures, and large seasonal and diurnal temperature fluctuations (Haselwandter 1987). Because o f harsher environmental factors found in these systems, traits such as longevity and mycelial spread o f individual fungal genets may be important (Gardes and Dahlberg 1996). Dark-septate endophytes (DSE) are ubiquitous in both alpine and arctic systems (Bledsoe et al. 1990, Cazares 1992). Although Kohn and Stasovski (1990) reported no DSE were found on root tips o f plants from Alexandra Fiord, samples o f DSE from this area were found later (see Chapter 3). Hyaline septate hyphae that are reported and found on several plants may have the same ambiguous function o f being either pathogenic or mutualistic, as has been found with DSE (Jumpponen and Trappe 1998). Arbuscular mycorrhizae (AM) are scarce in the arctic (Bledsoe et al. 1990, Kohn and Stasovski 1990) and higher elevations in alpine systems (Haselwandter 1987). Although, Dalpé and Aiken (1998) found approximately 10% o f Festuca species are associated with arbuscular mycorrhiza in the high Arctic, this is contrary to previous studies where no or very little AM was found (Bledsoe et al. 1990, Kohn and Stasovski 1990). The discrepancy between these may be due to the small sample size used by Bledsoe et al. (1990) and Kohn and Stasovski (1990) (Dalpé and Aiken 1998). Regardless, AM appear to be scarcer in arctic regions and in higher altitudes o f alpine systems. Both AM and ectomycorrhizae have been found on Salix spp. in alpine studies (Trowbridge and Jumpponen 2004) while only ectomycorrhizae have been found in arctic Salix systems (Vare et al. 1992). Because AM are scarce in the arctic in contrast to alpine systems, and DSE are common in both arctic and alpine, Bledsoe et al. (1990) suggested that DSE may replace the 37 functional role o f AM in arctic systems. Jumpponen (1999) suggested that the DSE, particularly Phialocephala fortinii, may allow for transport o f carbohydrates between plants through fungal linkages, as this function has been found for ectomycorrhizal fungi between Betula and Pseudotsuga (Simard et al. 1997). This was suggested because the same genet o f P. fortinii is found on nine different plant species that are classified as ecto-, ericoid, and non-mycorrhizal (Jumpponen 1999). Ericoid mycorrhizae are found on dw arf shrubs in the high arctic and nutritionally stressed alpine plant communities (Haselwandter 1987). Haselwandter (1987) suggested that ericaceous plants are capable o f using more complex nitrogen sources such as proteins or amino acids, which is supported by findings that ericaceous plants take up amino acids in alpine regions (Michelsen et al. 1996). Ericoid mycorrhizal fungi have also been found to access N and P by producing enzymes that break down structural components o f litter such as pectin and hemicellulose in plant cell walls, monophenols, tannins, polyphenols, and lignin (Smith and Read 1997). Although ectomyeorrhizal fungi break down these structures as well, the production o f these enzymes appears to be less common than by ericoid mycorrhizal fungi (Smith and Read 1997). The mycoflora o f alpine and arctic systems are similar in many aspects. Sporocarps of ectomycorrhizal fungi in both systems are sparse in comparison to lower elevation, temperate environments (Trappe 1988), which may be due to climatic factors which strongly influence fruiting (Gardes and Dahlberg 1996). Preliminary data from Alexandra Fiord indicate that several species, such as Russula sp., Cortinarius spp. and Inocybe sp., and Cenococcum geophilum are dominant. Cortinarius spp. are dominant species (comprising 20% o f the abundance) on Salix arctica and Dryas integrifolia roots from the Canadian arctic 38 archipelago (Gardes et al. 2000). Contrary to the findings o f Kohn and Stasovski (1990), Cenococcum geophilum was found in the present study as well as fruiting bodies o f Lycoperdon and Helvella. Lycoperdon spp. are reputed to be myeorrhizal with Picea abies, Pinus nigra, Pinus strobus, Pinus sylvestris, Pseudotsuga menziesii. Eucalyptus spp., and Quercus spp. (Trappe 1962). Helvella aestivalis formed myeorrhizae with Dryas octopetala, and species o f Helvella formed mycorrhizae with Salix reticulata under axenic eonditions (Weidemann et al. 1998). Helvella crispa is reported to form myeorrhizae with Fagus sylvatica and Quercus spp., and H infula with Picea abies (Trappe 1962). b) Mycorrhizal fungal succession in relation to plant community succession Glacier forefronts are commonly used for research on mycorrhizal fungi during primary succession. Primary succession is when pioneer species colonize virgin surfaees (Frankland 1998), and seeondary succession is when the soil is nutrient poor after a disturbance (Smith and Read 1997). Mycorrhizal fungi may improve nutrient-poor conditions for latter species as detected by increasing diversity after volcanic disturbances (Titus and del Moral 1998) and glacial tills (Helm et al. 1999). This facilitative nature o f mycorrhizal fungi is inferred by the successional pattern described for primary succession, which starts with non-mycorrhizal plants, AM plants, then ECM plants (Read 1993) and/or ericoid plants (Cazares 1992). Ectomycorrhizae are thought to colonize in older soils because they have access to nutrients contained in organic residues that are more abundant in later stages after accumulation o f organic matter (Read 1993), such as the increase o f nitrogen and organic matter (Jumpponen et al. 1998) along a chronosequence, which is a sequential change o f related variables in certain properties, from an alpine glacial forefront. 39 Changes in carbohydrates supplied by the host (Dighton and Mason 1985), nitrogen availability (Baar 1996), and soil conditions (Termoshuizen 1991, Kranabetter and Wylie 1998) are some o f the mechanisms suggested for primary succession to progress. For plant communities, primary succession can depend on life history traits, such as seed size and growth rate, maximum height o f the plant, seed rain, and competition. Facilitation and initial site conditions are important for the rate o f change and for species composition and productivity (Chapin et al. 1994). As found with plant communities, primary succession for mycorrhizal fungi probably is not dependent on a single variable. Several researchers found that changes in one variable are not enough to describe fungal succession (Termoshuizen 1991, Baar 1996, Helm et al. 1999, Kranabetter and Wylie 1998). Svoboda and Henry (1987) described three types o f succession; 1) directional, replacement succession with low resistance; 2 ) directional, non-replacement succession in high resistance environments; and 3) non-directional, non-replacement succession in extreme environments. In directional replacement succession, succession goes through serai stages with species replacement until a relatively stable ecosystem is reached. In directional, non­ replacement succession, species are not replaced but live in co-existence with the invading species, which expand slowly. In these systems, in which polar semi-deserts are an example, there is enough space for expansion. Non-directional, non-replacement succession is found in extreme environments, such as polar deserts where few species survive. Several species may invade repeatedly but fail to establish permanently. The lowlands o f Ellesmere Island fit the directional, non-replacement succession description. The mycorrhizal guild system may play an important role in plant competition (Newman 1988, Horton and Bruns 1998) if co-existence between plants is typical and 40 expansion is slow. Although Kropp and Trappe (1982) suggested that pioneer plants may be more host-specific, the case on Ellesmere may be different because succession does not follow the replacement o f plant species but the co-existence o f additional species. Moving away from the glacial forefront, the plant community starts with Papaver lapponicum and Luzula confusa. Salix arctica, Saxifraga oppositifolia, Cassiope tetragona, and Dryas integrifolia eventually appear, and all six species are found not only on the glacier forefront but the rest o f the lowlands o f Alexandra Fiord. Van der Heijden and Vosatka (1999) showed that with AM, the increase o f AM fungal composition and number leads to an increase in plant diversity as well because more variety o f AM fungi allow different plants to establish themselves. Perhaps the increase o f ectomycorrhizal and ericaceous mycorrhizal fungi will have a similar capacity o f increasing plant diversity and stability. Understanding succession on Ellesmere Island will be different from other studies, including those that are conducted in the Arctic. Previous successional studies in the Arctic occurred in the low arctic where trees still grow (Helm et al. 1999, Helm et al. 1996, Brubaker et al. 1995, Chapin et al. 1994) while Ellesmere Island is located in the high arctic where only low shrubs are found. 5. Global climate change: effects on aboveground plant community structure Global warming is a complex type o f disturbance because not only can it have direct effects on an organism or ecosystem, it can also lead to other disturbances. For example, warming has been linked to increase fire frequency (He et al. 2002) and more intense hurricanes (Shen et al. 2000). For this study, the effects o f warming will be examined on the community level even though this disturbance is classified at the ecosystem or landscape 41 level, according to the hierarchical organization of disturbances described by Pickett et al. (1989). Global circulation models predict that arctic systems will not only experience warming before other ecosystems, but also undergo the greatest increase in surface temperatures due to the doubling o f CO 2 (Shaver et al. 1992, Oecbel et al. 1993, Henry and Molau 1997). Climate change will have a more dramatic effect on the arctic than other forms of disturbance, mostly due to its spatial isolation, so findings from the arctic can be used to predict bow other systems may respond (Shaver et al. 1992). By the year 2100, approximately 63% o f biodiversity will be altered due to climate change in the arctic, compared to other human-induced disturbances such as changes in land use (15%), introduction o f exotic species (4%), and changes in atmospheric CO 2 and/or nitrogen deposition (18%) (Chapin et al. 2000). Warming is expected to increase more in winter months (up to 17° C) than during the summer months (~4° C) (Oechel et al. 1993, Edlund 1992), thus lengthening the growing season o f plants (Edlund 1992, Henry and Molau 1997) and altering plant communities through changes in the distribution o f snow in the winter, persistence o f snowbeds, and pattern o f snowmelt (Edlund 1992). Global climate change will likely amplify in arctic regions due to positive feedback loops that include; 1) ice and snow melt that would decrease surface albedo; 2 ) stabilization of the atmosphere that may trap temperature anomalies near the ground surface; 3) cloud dynamics that may amplify change (Overpeck et al. 1997); and 4) the permafrost layer melting sooner (Oechel et al. 1993). Warming in the arctic affects lower latitudes by possibly changing river run-off and the circulation o f the atmosphere, and increasing atmospheric concentrations o f CO 2 and CH4 (Overpeck et al. 1997, Henry and Molau 1997). 42 By understanding how plants established historically, predictions o f how climate change will affect plant species evolutionarily and geographically may be more accurate (Murray 1995). For example, in the early Holocene period, warming increased the number o f shrubs, which parallels the present spread o f dwarf shrubs {Salix spp., Betula nana, and Alnus crispa) in Alaska (Sturm et al. 2001). This is indirect evidence that these regions in Alaska may adapt relatively quickly to climate change (Sturm et al. 2001). Also, fossil records of some species such as Dryas integrifolia and Saxifraga oppositifolia indicate that these plants have existed since the Tertiary period (Murray 1995) and, therefore, have survived temperature fluctuations for at least 1.8 million years. Most o f the present arctic flora established approximately 6000-3000 b.p. (Brubaker et al. 1995) and originated from: 1) survivors from Tertiary forests, northern réfugias from Quaternary glaciation, and Pleistocene migration from Asia; 2) plants that returned during interglacial and post-glacial time from unglaciated areas; and 3) newly evolved species from the Pleistocene and Holocene (Murray 1995). According to Late Quaternary pollen records, species found in the arctic tundra are thought to have expanded southward into much of Canada (Brubaker et al. 1995). Arctic systems are carbon sinks. Current carbon sinks are the wet and moist tussock tundra of arctic systems (Oechel et al. 1993, Shaver et al. 1992). Arctic systems have three times more soil carbon than alpine systems but only 13% o f the plant species richness, which indicates active accumulation o f soil organic matter and little disturbance (Chapin and Komer 1995). Release o f carbon to the atmosphere is predicted to be caused indirectly and not directly from the increase o f temperature (Oechel et al. 1993). Researchers have suggested that the cause o f the loss o f carbon from arctic systems to the atmosphere is due to 43 enhanced drainage and soil aeration, deerease in the water table (Oechel et al. 1993, Shaver et al. 1992, Billings et al. 1983), and increase in respiration, especially from the soil microbial community (Schimel 1995, Billings et al. 1983) o f whieh myeorrhizal fungal hyphae are thought to be a large eontributor (Rygiewicz and Andersen 1994). Change in carbon storage is somewhat constrained by the nitrogen cycle because nitrogen is the primary limiting faetor in aretie systems (Shaver et al. 1992, McKane et al. 1997). With enhaneed drainage and soil aeration, decomposition and release o f carbon will likely occur in systems that have large amounts o f earbon storage such as high latitudinal bogs, and boreal and arctic systems (Oeehel et al. 1993). Billings et al. (1983) found that a 48 ° C warming decreased the net carbon storage in the wet sedge tundra rather than increased net primary production, which they attributed to greater inerease in soil respiration. However, loss from carbon storage may be for the short-term, and eventually increase in above ground plant biomass may compensate for the earbon loss (Oechel et al. 1993). Many factors influence the impacts o f global climate change on above ground plant growth such as water availability, nutrient availability, summer warmth, snowfall, (Edlund 1992, Field et al. 1992), light, and CO 2 levels (Field et al. 1992). For arctie systems, warming may first impact individual plants, indicated by an increase o f vegetation growth (Edlund 1992). Response by plant communities would depend on the eombination of summer warming, snowfall in the winter, possible drought in the summer (Edlund 1992), and resource availability (Field et al. 1992). Global warming may result in major reorganization o f plant communities (Brubaker et al. 1995); however these changes for the plant eommunity may take centuries (Edlund 1992). 44 Several studies have examined the impacts o f global warming on arctic plants in situ by manipulating temperatures with greenhouses (Hobbie and Chapin 1998, Havstrom et al. 1993) or open-top chambers (OTCs) (Henry and Molau 1997, Stenstrom et al. 1997, Jones et al. 1997) placed over plots. Open-top chambers have some advantages over closed greenhouse systems by allowing in more direct solar radiation; lessening the chance of overheating; allowing herbivores and pollinators to the plants; and avoiding decreased relative humidity (Marion et al. 1997). Problems o f both systems include increasing temperature extremes rather than lowering the range o f diurnal temperatures, altering o f wind patterns around the plant, and disturbing the sites (Marion et al. 1997, Hobbie and Chapin 1998). Problems that are unique to OTCs consist o f snow accumulation, disturbance by animals (Marion et al. 1997), and only a small area can be uniformly warmed (Shaver et al. 2000 ). Table 1.1 summarizes experiments o f warming on dominant plants o f arctic systems that will be used in this present study. Experiments using greenhouses to increase air temperature find no significance o f warming on Cassiope tetragona (Hobbie and Chapin 1998, Havstrom et al. 1993), which lead researchers to conclude that perhaps nutrients rather than temperature affect C. tetragona growth (Hobbie and Chapin 1998). Their findings are contrary to what is found when OTCs are used, where warming did increase different factors measuring plant growth (Henry and Molau 1997). This discrepancy may be because major changes in the arctic tundra from warming o f the last glaciation have little similarity in different circumarctic sites (Brubaker et al. 1995), as the two studies are in Alaska and eastern Canada. Another explanation may be that the greenhouse experiments do not allow for enough time for temperature increase to show significant differences as is found with an 45 OTC study on Salix arctica, where Henry and Molau (1997) found no significance after two years but did after four years. 46 Table 1.1 Summary of experiments o f warming on arctic plants used m present study. Plant Reference Location Temperature Plant Growth Manipulation Cassiope Hobbie and Alaska, Greenhouse Biomass decrease. tetragona Chapin (1998) USA Havstrom et al. Spitsbergen, Greenhouse Significant increase from (1993) Norway temperature for leaf indices and shoot growth index. Henry and Ellesmere OTC Significant increase in plant Molau (1997) Island, growth. Canada Stenstrom et Saxifraga OTC Ellesmere Increase in flowering oppositifolia al. (1997) Island, frequency and reproductive Canada success. Salix arctica Jones et al. OTC Increase o f plant growth after Ellesmere (1997) Island, 3-4 years o f warming. IncreaseW in seed production Canada first years o f warming. Stronger impact on male willows than females. Dryas Henry and OTC Ellesmere Significant but not strong integrifolia Molau (1997) Island, increase in vegetative and Canada reproductive phenologies, seed set and weight, germination. NA: not available Predictions of Adaptation Due to Global Warming NA Will probably not migrate south as air temperature increases. NA Will probably be outcompeted by graminoids and forbs. May adjust easily to global warming because it adapts to broad range o f ecosystems. NA 6 . Impact o f global climate change on root fungal communities Much o f the research on global climate change on mycorrhizae has been indirect focusing on the effects o f elevated CO 2. Results have been conflicting showing decreased (Fitter et al. 2000), increased, and no difference for ectomycorrhizal and AM colonization due to elevated CO 2 levels (Fitter et al. 2000, Treseder and Allen 2000). The impact o f elevated CO 2 on mycorrhizal growth and colonization seems to depend on; 1) mycorrhizal fungal species because some species are more sensitive than others to elevated carbon (Fitter et al. 2000); 2) availability o f N, where additional N can negate the effects o f CO 2 on mycorrhizal biomass for some systems (Treseder and Allen 2000); 3) growth rate o f plants because larger plants need more roots (Fitter et al. 2000, Treseder and Allen 2000); and 4) roots lengthening which would increase mycorrhizal colonization (Eissenstat et al. 2000). In a review by Fitter et al. (2000), only two studies that examine the increase o f temperature on AM colonization are published and no studies have currently been published on the effects on ectomycorrhizae. Perry et al. (1990) speculated on the role o f mycorrhizal fungi in climate change in that plant species would migrate during climate change and that sharing mycorrhizal fungi would help with the transition. Although there is a lack o f experiments that examine the impact o f global warming on mycorrhizae, the rhizosphere will probably be an important factor in how ecosystems adjust. The major effects o f global warming may be from its impact on soil processes rather than the increase of biomass o f plants in the tundra (Hobbie and Chapin 1998). Boone et al. (1998) suggested that the rhizosphere would be more sensitive to warmer temperatures than above­ ground plant parts, and that variation in soil respiration is determined by responses o f root respiration and heterotrophs to temperature change. Warmer soil temperature may influence 48 root growth, cell elongation, initiation o f new lateral roots, increase in root respiration and ion uptake, interaction with water and nutrient availability, more N mineralization, less water availability, and earlier initiation o f root growth in the spring (Pregitzer et al. 2000). Although global warming is suspected to have significant impacts on the rhizosphere, much o f implications have been speculative. In addition, information on the impacts on the mycorrhizal community, in particular ectomycorrhiza, is scarce. C. Research objectives and hypotheses Alexandra Fiord provides an opportune site to examine root-associated fungal communities. The role o f these fungi in plant establishment in primary succession increases the understanding of how plants and fungi adapt to nutrient-poor conditions. Future ecological conditions are examined by use o f OTCs to simulate potential global warming scenarios. This study will be one o f the first to examine the impact o f global warming on the root-associated fungal community. Because global warming impacts arctic systems more intensely than temperate environments, this may give insight to the role o f root-associated fungi in facilitating changes to the plant community. Although global warming is suspected to have significant impacts on the rhizosphere, much o f the implications have been speculative. This study examines how warming may impact the root-associated fungal community by using both PCR-based techniques and isolating fungi from root tips. Although there should be some overlap, the PCR-based techniques are more likely to favor mycorrhizal fungi; whereas, fungal isolations would favor faster-growing root endophytes that do not fit the morphological definition o f mycorrhizae. 49 By examining the changes o f the root-associated fungal community in a direct, non­ replacement succession, insight will be gained as to how this community would adjust to a changing plant community. This type o f succession also has the unique characteristic where the increase in biodiversity o f plants happens in vivo while retaining all the same plants. This study can determine if the root-associated fungal community follows a similar trend. Also, this will be the first study to examine the root endophytic community for this type o f succession. Given what the literature indicates about fungal community structure, I expect that the root-associated fungal community, based on direct DNA extraction, will differ according to site and treatment (passive warming versus ambient). The following null hypotheses will therefore be tested: Ho 1.1 The root-associatedfungal community, based on DNA directly extracted from root tips from Salix arctica, will not differ between warmed plots and ambient plots. Ho 1.2 The root-associatedfungal community, based on DNA directly extracted from root tips from Salix arctica, will not differ due to site. Because culture studies may reveal a different perspective than direct DNA amplification studies, I expecte that the root-associated fungal community will differ according to site, treatment (passive warming versus ambient), and host plant. Therefore, the following null hypotheses will be tested: Ho 1.3 The root-associated fungal community, based on cultures isolatedfrom root tips, will not differ between warmed and ambient plots. 50 Ho 1.4 The root-associated fungal community, based on cultures isolated from root tips, will not differ according to the host plants Cassiope tetragona, Dry as integrifolia, Salix arctica, and Saxifraga oppositifolia Ho 1.5 The root-associated fungal community, based on cultures isolatedfrom root tips, will not differ due to site. Because eulture studies only assess culturable fungi, 1 expeet that the fungal communities described by the two methods (direet DNA amplifieation from roots versus culturing) will be different (despite revealing the same patterns aeeording to site, treatment, and host plant). Therefore, the following null hypothesis will be tested: Ho 1.6 The root-associated fungal communities described by the two methods (direct extraction versus culturing) will not differ. The unusual directional, non-replacement sueeession pattern found in high arctic systems permits me to examine how diversity on different host plants varies along a ehronosequenee, without the confounding factor o f host plant replacement. My objective was to examine how the root-associated fungal community changes during a directional, non­ replacement primary plant succession, using Cassiope tetragona, Dryas integrifolia, Luzula confusa, Papaver lapponicum, Salix arctica, and Saxifraga oppositifolia as host plants Therefore, the following hypotheses for this objeetive will be tested: Ho 2.1 The root-associated fungal communities will not differ along a ehronosequenee. Ho 2.2 The root-associated fungal communities will not differ according to host plant. 51 D. Thesis organization This thesis starts with a literature review to provide background on the concepts, theories, and techniques associated with this thesis. This provides more comprehensive information that may not be covered in subsequent chapters. The next two chapters address the first research objective, using two different methods that collectively provide a more complete assessment o f the root-associated fungal community. One method, directly amplifies fungal DNA from root tips, and the second involves isolations o f fungi from roots. These two methods have described different fungal communities found on the same plant host in previous studies. Chapter 2 addresses how warming will impact the root-associated fungal community detected by direct extraction o f fungal DNA from root tips o f one host species and covers hypotheses 1.1-1.2. Chapter 3 addresses the question based upon fungal cultures from several host species, aseptically isolated from plant roots, and will cover hypotheses 1.3-1.5. Chapter 4 covers research objective 2 and chapter 5 is a synthesis o f the research findings. This synthesis will cover hypothesis 1.6 , and will attempt to tie the three studies together. The final chapter is a summary o f the thesis. 52 IL Impact of warming on the frequency and diversity of the root-associated fungal community on Salix arctica from the Canadian High Arctic A. Introduction Global circulation models predict that arctic systems experience will experience a greater effeet o f global warming before other ecosystems due to increased surface temperature and CO 2 levels (Shaver et al. 1992, Oechel et al. 1993, Henry and Molau 1997), decreased surface albedo, alterations in cloud dynamics that may amplify change (Overpeek et al. 1997), and melting o f the permafrost layer (Oechel et al. 1993). Plant communities may ehange due to lengthening o f the growing season (Edlund 1992, Henry and Molau 1997) by inereasing air and soil temperature (Oechel et al. 1993, Edlund 1992) and altering water distribution by changing the dispersal o f snow, increasing the persistence o f snowbeds, and modifying the pattern o f snowmelt (Edlund 1992). The objective o f this study was to examine the impacts of experimental warming on the root-associated fungal community o f Salix arctica in the Canadian high arctic. Although there have been numerous studies on the impact o f experimental warming on arctic plants (e.g. Chapin et al. 1995, Henry and Molau 1997, Hobble and Chapin 1998, Jones et al. 1997, Sturm 2001), no studies have examined how warming may impact the root-associated fungal community even though the rhizosphere may play an important role in plant response (Hobble and Chapin 1998, Boone et al. 1998). To date, only one study has examined the impact of warming on an ectomycorrhizal fungal community from Douglas-fir {Pseudotsuga menziesii Mirb. Franco) seedlings (Rygiewicz et al. 2000), where they found warming to increase species richness. Perry et al. (1990) speculated that the mycorrhizal fungal community may help facilitate migration o f plants. Hobble and Chapin (1998) suggested that 53 the major effects o f global warming will result from the influences o f increased soil temperatures and the subsequent soil processes. Increased soil temperature may influence root growth by increasing cell elongation, initiating new lateral roots, increasing root respiration and ion uptake, interacting with water and nutrient availability, increasing N mineralization, decreasing water availability, and initiating root growth in the spring earlier (Pregitzer et al. 2000). Other studies have indirectly examined the impact o f warming on mycorrhizal fungi by studying how CO 2 fluxes may impact the community (Oechel et al. 1993). Some researchers have suggested that the main cause of the loss o f carbon from arctic systems to the atmosphere is enhanced drainage and soil aeration, decrease in the water table (Oechel et al. 1993, Shaver et al. 1992, Billings et al. 1983) and increase in respiration, especially from the soil microbial community (Schimel 1995, Billings et al. 1983) o f which mycorrhizal fungal hyphae are thought to be a large contributor (Rygiewicz and Andersen 1994). We chose three distinct sites, a lowland site, highland granitic site, and highland dolomitic site, to see how soil type interacts with warming to influence the abundance, composition, and biodiversity o f these communities. In order to address these questions, we used terminal restriction fragment length polymorphism (T-RFLP) analysis, a technique commonly used in prokaryotic systems and but only used recently in mycorrhizal community analysis (Dickie et al. 2002, Klamer et al. 2002). B. Materials and methods 1. Study Site Samples were collected from three sites at Alexandra Fiord on Ellesmere Island, Nunavut, Canada, 78° 53’N, 75° 55’W. One site was located on the lowland (valley bottom) 54 and two on a mountain plateau. The lowland mesic site was at or near sea level and enelosed by 450-700 m-high plateaus to the east and west, a glaeial forefront to the south, and the oeean to the north. The elimate at the lowland site was relatively warm due to frequent sunny skies, relatively warm air masses from the west and south (Labine 1994), and reflection o f sunlight from the surrounding cliffs and oeean (Freedman et al. 1994). The soil was granitic, which is mostly composed o f quartz, feldspar, and mica. The highland sites were located on top o f the western plateau at 450 m. These sites were xeric with sparse vegetation and have desert and semi-desert arctic conditions (Batten and Svoboda 1994), which accounted for the decrease o f diversity o f vascular plants by approximately 40% when compared to lowland areas (Batten and Svoboda 1994). The two highland sites were distinguished by soil type, dolomitic, which is distinguished by high amounts o f calcium magnesium carbonite (CaMg ( € 0 3 )2), and granitic for the other site. 2. Vegetation Alexandra Fiord has been described as a ‘polar oasis’ because it comprises an 8 -km^ pocket of arctic shrubs, mosses, lichens, and sedges nested within vast ice fields. Dryas integrifolia Vahl, Cassiope tetragona (L.) D.Don, and Salix arctica Pall, are the most prominent (Freedman et al. 1994) o f the 92 species o f vascular plants found on the lowlands (Ball and Hill 1994) as well as the less vegetated and less diverse uplands (Batten and Svoboda 1994). The lowland site has been described as a dwarf shrub-cushion plant community (Muc et al. 1994) and is dominated by Salix arctica and Dryas integrifolia, along with Saxifraga oppositifolia L., Cassiope tetragona, Papaver, Pedicularis, sedges, and mosses. The granitic upland site, described as a S. arctica-C. tetragona dominant community (Batten and Svoboda 55 1994) includes S. oppositifolia and D. integrifolia. The upland dolomitic site is a S. oppositifolia-àommstQdi community (Batten and Svoboda 1994) with no D. integrifolia present. To simulate the increased temperatures predicted by global climate models, open-top chambers (OTCs), which covered 0.8 m^ and were 0.3 m high, were used to increase air temperatures by 1-4° C, the predicted temperature range o f global maean increases for the middle of the 21®' century for the Canadian high arctic due to climate change (Henry and Molau 1997). Three 1-m diameter OTCs were placed on each of the three sites; in 1995 (GHR Henry 2000, pers. com.) for the lowland site, and 1993 for the two highland sites (Stenstrom et al. 1997). Plants that were found between 0.5 m and 1.5 m from the OTCs were harvested as controls for a total o f 18 plots. 3. Field collection Two specimens o f Salix arctica were harvested from each o f the plots in August 2000. Plants and surrounding soil were kept in Ziploc® bags (18 x 20 cm) in a permafrost cooler while in the field and in a 4° C refrigerator once back at UNBC until processed. Two 300 gram samples of soil were collected in August 2001 from each plot. Soil was collected no more than 1 m from the harvested plant. Soils were dried and separated in a 2-mm sieve to remove rocks from the samples. One hundred grams from each o f the two replicates were mixed and sent to the Ministry o f Forests, Research Branch Laboratory, Analytical Chemistry Section in Victoria, BC for the following analyses: pH in water, total C and N using combustion elemental analysis, cation exchange capacity and exchangeable cations using 0.1 N barium chloride extraction, available NH 4-N and NO 3-N extracted with 2N KCl, and available P using the Mehlich III protocol. 56 Fungal sporocarps from Alexandra Fiord were collected and tentatively identified to genus; however, no sporocarps were found in any o f the plots per se. A 2- x 2- mm^ piece was extracted from each sporocarp and stored in 50% EtOH for DNA extraction. The remaining sporocarp tissue was dried for storage. 4. Sampling from roots Plant roots were immersed in water for at least 24 hours at 4° C. The root systems were gently cleaned with water and collected in a 0.5 mm sieve (No. 35 USA standard testing sieve, W.S. Tyler, Inc.). Root systems were placed on a numbered grid for random selection o f root tips. Numbers from a random-numbers table were used to select grids for sampling. The root tip that traversed or was closest to the grid was selected, and a 2-cm root section was sampled. Fifteen root tips were randomly selected, described morphologically, placed in a 1.5 ml microcentrifuge tube, and frozen at -40° C until DNA was extracted. Root tips were assessed morphologically (morphotyping) based upon characteristics such as color, tip ramification, and presence and absence o f rhizomorphs, loosely following the techniques o f Agerer (1987-1998) and Goodman et al. (1995). An additional forty root tips were randomly selected, frozen, lyophilized, and stored at -40° C for further DNA extractions if needed. O f these additional tips, approximately 15 were later extracted to increase number o f samples. The remaining root system was frozen and lyophilized. Preparation for DNA extraction was done first for the highland sites because the plants did not appear as robust as at the lowland site. Plant tissue from stems or leaves was also sampled to determine if plant DNA would amplify with the chosen primers. 57 5. DNA Extraction and ITS-T-RELP analysis DNA o f root tips was extracted using the CTAB protoeol o f Gardes and Bruns (1996b), which was modified by excluding the freeze-thaw procedure and by including another purification step o f adding phenolxhloroform-isoamyl aleohol (1:1) (Lee and Taylor 1990). Individual frozen tips were ground in 300 pL o f 2X CTAB buffer (100 mM Tris at pH 8 , 1.4 M NaCl, 20 mM EDTA, 2% CTAB, and 0.2% p-mercaptoethanol) and then incubated at 65° C for one hour. Equal volumes o f phenol: chloroform-isoamyl alcohol (24:1) were added, vortexed, and centrifuged at 13,000 ref for 15 minutes. The supernatant was then transferred to fresh tubes. A second wash o f chloroform-isoamyl alcohol (24:1) was added (v/v), vortexed and eentrifuged for 5 minutes. Again the supernatant was removed and placed in fresh tubes. Nucleic acid was precipitated by adding 500 pL of isopropanol and ineubated for at least 3 hrs at -20° C. This was centrifuged for 15 min, and the isopropanol was removed. The remaining pellet was washed with 70% ethanol and then centrifuged for 10 minutes. The ethanol was removed and allowed to evaporate before the pellet was resuspended in 50 pL o f TE -8 buffer. Plant tissue and sporoearps were extraeted with the same protocol except the extra phenol :chloroform-isoamyl alcohol purification step was not included, and pellets were resuspended in 100 pL o f TE -8 buffer instead o f 50. All extractions were stored in -20° C until use. Amplification o f the nrDNA ITS region was done using lOX Buffer (200 mM TrisHCl [pH 8.4], 500 mM KCl) (Invitrogen), 2X dNTPs, 25 mM MgCb, 0.5 pM o f ITS 1 (White et al. 1990) dye labeled with Cy 5.0, 0.5 pM o f ITS 4 (Gardes and Bruns 1993) dye labeled with Cy 5.5, 5 U o f Platinum Taq DNA Polymerase (Invitrogen), and 0.1 OX DNA template . The following program was used for PCR on a MJ Research thermocycler PTC100: 94 °C (4 min); 48° C (1 min); 72° C (2min); [94 °C (30 sec); 48° C (30 sec); 72° C (1 58 min 30 sec) x 34 cycles]; 72°C (6 min 30 sec). Samples were run on a 0.7% agarose gel to confirm amplification. Unsuccessful amplifications were redone with no dilution and 1:10 dilutions o f the DNA template. Digestions using either Alul or Hin^{ (Invitrogen) and 6 pL o f the PCR produet were completed in 10 pL reactions following the manufacturer’s recommendations and incubated at 37° C for at least 3 hours. One microliter o f the digested samples was mixed with 1.65 pL o f loading dye mixture that contained formamide and two sets o f internal markers at 101 , 200, and 351 bp; one set labeled with Cy 5.0, the other with Cy 5.5. The samples containing the loading dye mixture were denatured at 80° C for two minutes and then quenched on ice. Two microliters were then loaded onto a 6 % polyacrylamide gel with the laser power set at 50%, temperature at 53° C, and current at 1250 V and ran for 60 min. on an OpenGene System Long Tower Sequencer (Bayer International). Fragments for each primer-enzyme combination (i.e. ITS \-Alul, ITS \-HinQ., ITS A-Alul, and ITS 4-7/mfl) were determined using GeneObjects 3.1 fragment analysis software. 6 . Matching root tips with sporocarps Terminal restriction fragments (T-RFs) from root tips were compared to those from sporocarps using TRAMP (T-RFLP analysis matching program) (Dickie et al. 2002) to check for matches. However, because most o f the root tips exhibited multiple fragments, determining T-RFs for individual genotypes was difficult. Identification o f T-RFs for individual genotypes was attempted by matching all possible combinations o f at least three fragments against the known fragments from sporocarps. 59 7. Data analysis Frequencies o f genotypes generated from the T-RPLP analysis were tabulated for each cnzymc-primcr combination. First, genotypes were binned and the average o f the binned numbers was used to identify that genotype, then the number o f genotypes in each bin was entered into frequency tables. Each different primer-enzyme combination was treated as an independent database for determining if site and treatment bad changed the rootassociated fungal communities; therefore all analyses except for the ordination analyses were done in quadruplicate. Non-metric multidimensional sealing (NMS) was used to analyze the frequency tables generated from the four primer-enzyme combinations using PC-ORD, version 4.25 (McCune and Mefford 1999) as an initial exploratory tool to visualize the overall effects of site and treatment on genotype frequency. The four individual primer-enzyme frequency tables were merged into a single table for the ordination, and Beal’s transformation was used to alleviate problems associated with databases containing numerous zeroes (McCune and Grace 2002). The following parameters for NMS were used: Sorensen was used for measuring distance; the configuration for the first run was randomly selected; for all subsequent runs the best configuration from the previous run was used, as recommended by the program; 40 runs were conducted; dimensionality was assessed by examining the scree plot; and 50 randomized runs were used in the Monte Carlo test. Canonical correspondence analysis (CCA) in PC-ORD (McCune and Mefford 1999) was used to examine the influence o f soil properties (i.e. pH, CEC, available NH4, available NO 3, available P, and C:N ratio, and exchangeable cations) in explaining the differences in genotype occurrence among sites. Genotype frequencies were transformed by adding one (to remove zero frequencies), and soil properties were transformed by log 10 except for pH and 60 C:N ratio. Log transformation was used on soil properties because results from the different properties greatly varied so the log transformation compressed high values and spread the low values (McCune and Grace 2002). A Monte Carlo test was used to test for linear relationships between the genotype abundance and soil properties. This test was used because a large number of reiterations were needed to gain a more precise p-value (McCune and Grace 2002). The ordination diagram was based on LC (linear combination) scores, which are linear combinations o f soil properties (McCune and Grace 2002). To further investigate the impacts o f site and treatment on genotype frequency, a nested, two-way analysis o f variance (ANOVA) in Statistica, vers. 6.1 (StatSoft) was used. Because the number of fungi extracted from root tips varied from each plant specimen, genotype frequency was standardized by; total number genotypes per plant specimen^ . . . . . This standardization [1]------------------------------ —— - — ------total number o f root tips used from plant specimen^ r,, assumes that genotypes were not present in unsuccessful amplifications. The statistical model used was: [2] Yijki = p + Si + Tj + Rk(j) + Pi(jk) + (ST)ij + S(ijki)m, where yijki = genotype frequency Si = effect o f site, i = 1,2,3; Tj = effect o f treatment (OTC or control), j = 1,2; Rk(j) = replication effect, k = 1,2,3 nested within theyth treatment. Pi(jk) = plant specimen effect nested within the Ath replicate and jth treatment, 1 = 1,2 . ANOVA was also used to test if warming affected soil properties, which could explain some o f the variation due to site or treatment. The model used was: [3] yijk = p + Si + Tj + Rk(j) + Si*Tj + £(ijk)j., where 61 yyk = the amount (ppm for P and available N ) , concentration (exchangeable cations and CEC), or ratio (C;N ratio, pH) o f soil properties; Si = effect o f site, i = 1,2,3; Tj = effect o f treatment (OTC or control), j = 1,2; Rk(j) = replication effect, k = 1,2,3 nested within the treatment. Genotype diversity was investigated by examining genotype abundance curves, genotype richness, and evenness according to site and treatment. Genotype abundance curves were made by taking the natural log o f genotype frequency and plotting it against the arithmetic ranking o f frequency. Once a model was chosen from the genotype abundance curves, niche apportionment analyses were conducted using PowerNiche, an Excel-based macro that uses niche division algorithms (Drozd and Novotny 1999), which helps determine if abundance is associated with random fraction, power fraction, broken stick, or other niche apportionment models based on the works o f Tokeshi, Sugihara, and Mae Arthur (see Magurran 2004). The selection and division exponents were varied for each test o f the power model with 250 replications. Different models were tested by changing the selection (= k) and division exponent (= m), which is indicative o f which niche apportionment model is chosen, e.g. to test for the broken stick model, one is chosen for both exponents (i.e., k = 1 , m = 1). Other than the broken stick model, the random fraction model (k=0, m =l), the power fraction model (0) to identify fungi that were most closely related to the sequence and to ensure that sequences were fungal. Edited sequences were imported into Mac Vector (Oxford Molecular) and were aligned with the most closely related sequences downloaded from GenBank. Sequences downloaded included Epicoccum nigrum AF455455, Colispora elongata AY148102 Leaf litter Ascomycete AF502763, Epacris microphylla root associate AY268216, Dothideales sp. AY465446, Fungal endophyte MUT 2723 AF373055, Fungal endophyte MUT 585 AF373051, Ascomycete sp. AJ279473, Phoma glomerata AY618248 (Fig. 3.4); Zalerion varium 88 A F169303, Glarea lozoyensis A F169304, Hymenoscyphus epipiphyllus AY348581, Hymenoscyphus monotropae AF 169309, Lachnum bicolor AJ430394, Phialocephala virens AF486132, Xenochalara juniperi AF 184889, Phialophora sp. AY465463, Cadophora sp. AY371512, Phialophora finlandia AF486119, axenic ericoid root isolate AJ430119, ectomycorrhiza cf. Hymenoscyphus ericae (= Rhizoscyphus ericae (Zhang and Zhuang 2004)) AJ430150, Hymenoscyphus sp. AY354244, ectomycorrhizal isolate (Helotiales) AJ430410, Cadophora luteo-olivacea AY249069, Leptodontidium orchidicola A F214578, Cadophora gregata AY249071, Cadophora malorum AY249063, Cadophora sp. AY371506, ectomycorrhiza (cf. Phialocephala fortinii) AJ430214, Phialocephala fortinii crypt sp. AY347405, Phialocephala fortinii AY078138 (Fig 3.5); Inocybe angustispora AY380360, Favolaschia cf. sprucei AF261420, Favolaschia cf. calocera AF261419, Panellus serotinus A F518633, Ripartites metrodii AF042012, Mycena haematopoda AJ406590, Mycena leaiana AF261411, Omphalina rivulicola U66451, Collybia tuberosa AY639884, Clitocybe subvelosa AY647208 (Fig. 3.6); Phlebia lindtneri AF141623, Ceratobasidium goodyerae-repentis AY243523, Uthatobasidium fusisporum A F518664, Thanatephorus cucumeris AF354062, Ceratobasidium sp. AGO AF354094, Ceratobasidium sp. AGL AF354093, Inonotus weirii AY059040, Phlebiopsis gigantea AF141634 (Fig 3.7); Alternaria helianthi A Y \5 4 7 \3 , Ampelomyces sp. AY293794, Ascomycota sp. AJ301960, Byssoascus striatosporus AB040688, Cadophora luteo-olivacea AY249087, Cadophora malorum AY249086, Cenococcum geophilum A Y l 12935, Chalara fungorum AF222462, Chalara kendrickii AF222464, Chalara longispes AF222466, Chalara microchona AF222468, Chalaraparvispora AF222473, Chalara strobilina AF222477, Cryptosporiopsis sp. AY442321, Cryptosporiopsis ericacea AY442323, Cudonia lutea AF433139, ericoid 89 mycorrhizal sp PPO-7 AY599245, ericoid mycorrhizal sp PPO-8 AY599246, ericoid mycorrhizal sp. AY599244, ericoid mycorrhizal sp. PPO-5 AY599243, ericoid mycorrhizal sp. PPO-2 AY599240, euascomycete RFLP type A A F127116, euascomycete RFLP type B AF 127117, euaseomycete RFLP type C AF 127118, Fabrella tsugae AF356694, Graphium rubrum AY266313, Hymenoscyphus ericae (= Rhizoscyphus ericae (Zhang and Zhuang 2004)) AY284122, Hymenoscyphus sp. UBC tra 1436 AY219881, lodosphaeria sp. AF452045, Lachnum cf. bicolor AY544674, Lachnum virgineum AY544646, Lecythophora sp AF353607, Leptosphaeria doliolum U43475, mycorrhizal sp. AF081443, Mycosphaerella mycopappi U43480, Oidiodendron tenuissimum AB040706, Phialocephala fortinii AF326082, Phialocephala dimorphospora AF326081, Phialophora gregata AF222502, Phialophora sp. AF 156922, Phoma herbarum AY293788, Phoma sp. 199 AY293785, Pleomassaria siparia AY004341, Pleospora herbarum U43476, Rhytisma acerinum AF356696, Setomelanomma holmii AF525678, Shiraia bambusicola AB105798, Trematosphaeria heterospora AY016369 (Fig 3.8). Neighbor-joining analysis was used for all trees, with 1000 repetitions for bootstrap values. The Kimura 2-parameter model was used for the distance measure, the transitiomtransversion ratio was estimated, and gaps were distributed proportionally. Bootstrap values greater than 50 were reported. Identifications from samples that were sequenced were applied to those samples with matching RFLP patterns. Those that clustered according to RFLP data were verified by viewing composite gels o f likely samples in Gene Profiler database, vers. 4.02 (Scanalytic, Inc.). These samples were then compared to those identified based on morphological traits. 90 c. Results 1. ITS-RFLPs Out o f 720 isolated cultures, 432 were amplified successfully and digested with all three restriction enzymes for RFLP analysis. For statistical analyses, 30 o f the isolates were deleted for a total o f 402. These isolates were not included for statistical analyses because they were from plants that were not sampled within the experimental design (i.e. they were sequenced to increase the chances o f identification by providing more isolates from the harsher highland sites). The species cumulative frequency for each plot is illustrated in Table 3.1. Isolating fungi from root tips was most successful for Salix arctica and Saxifraga oppositifolia. The average number o f isolates was highest from S. arctica growing in control plots at the granitic site (p=0.02) (Fig 3.1). Although not statistically significant, more isolates from S. oppositifolia were also recovered from the granitic site, averaging 15 isolates from both plant specimens for the granitic site compared to six for the other two (Fig. 3.2). Fungal isolation was particularly unsuccessful for Cassiope tetragona from the warmed plots for all three sites, as well as for the control from the dolomitic site; fungi were only successfully isolated from one plot on the granitic site. In contrast to Salix arctica and Saxifraga oppositifolia, no fungi were successfully isolated from roots o f Dryas integrifolia from the granitic site. Results from NMS showed that the site was the most influential factor in differentiating the root endophytic fungal community, more so than warming or host plant (Fig 3.3). The lowland site had a wider range o f variation when compared to the other two sites. The dolomitic site also had large variation, but the isolates were more similar at this 91 Table 3.1 Number of successful samp es isolated from root tips based on R "LPs. OTC Control Treatment Replicate Replicate 2 Sum 1 Sum 2 3 1 3 (cumulative (cumulative Site frequency) frequency) Lowland 2 C. tetragona 0 0 8 8 8 6 16 21 4 6 7 17 10 5 6 D. integrifolia 10 4 23 4 Salix arctica 9 5 6 15 Saxifraga 5 8 3 16 10 1 6 17 oppositifolia 24 22 64 27 24 Sum (species 18 18 69 abundance) Dolomitic C. tetragona 0 0 0 0 0 0 0 0 D. integrifolia 5 3 5 13 2 8 5 15 2 4 Salix arctica 7 8 17 4 7 15 1 12 Saxifraga 3 8 4 2 17 23 oppositifolia 21 42 Sum (species 15 6 10 14 29 53 abundance) Granitic 0 0 C. tetragona 0 0 10 0 0 10 D. integrifolia 0 0 0 0 0 0 0 0 Salix arctica 17 6 5 24 12 28 16 52 11 Saxifraga 16 16 43 8 13 20 40 oppositifolia Sum (species 17 21 71 42 33 36 25 103 abundance) (a) Salix arctica (b) Saxifraga oppositifolia p = 0 .0 1 7 20.00 - A 15-00 □ Control 10.00 a O TC Lowland DoionntK Granitic □ C ontrol 10.00 OTC Lowland Dolomitic Granitic Site Fig. 3.1 Species cumulative frequency for Salix arctica and Saxifraga oppositifolia on different sites. LSM ± SE 92 (a) Salix arctica p=0.02 12.00 (b) Saxifraga oppositifolia 14.00 10.00 12.00 1.00 10.00 6.00 4.00 2.00 □ C ontrol 8.00 - □ C ontrol HOTC 6.00 HOTC 2.00 ’ 0.00 0.00 Lowland D obm itic Site Granitic Lowland Dolomitic Granitic Site Fig. 3.2 Species richness for Salix arctica and Saxifraga oppositifolia on different sites. LSM ± SE. site than in the community from the lowland site. The genotypes within the granitie site clustered tightly together in the ordination and appeared to have a similar composition. Species richness ranged from 0-13 for Cassiope tetragona, 0-18 for Dryas integrifolia, 11-27 for Salix arctica, and 10-31 for Saxifraga oppositifolia (Table 3.2) for eaeh o f the two replicate plants for each treatment. Beeause fungi were not recovered from all o f the sites for Cassiope tetragona and Dryas integrifolia, ANOVAs were only done for Salix arctica and Saxifraga oppositifolia. There were significantly more different types o f cultures isolated from Salix arctica in the granitic site (p=0.02) (Table 3.2) than the other two sites. Although not statistically significant, more diversity in fungal species was found (Table 3.1, Table 3.3) in the granitic control than the other plots for Salix arctica. For Saxifraga oppositifolia, diversity was greatest (p=0.017) in the granitic site, and a trend showing higher cumulative frequency in the granitic site was found as well. None o f the plant tissues amplified with the primer pair ITS 1/4. 93 Final Configuration, dim ension 1 v s. dim ension 2 2.0 3cb 0.5 Icb 30 # 4cd 0.0 3cd iOb lOb 3cc Icc lcd •0.5 lOd 10c - 1.0 -1.5 -1.5 30c - 1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 Dim ension 1 Fig. 3.3 NMS for host plant and warming treatments per site based on phi distances. The first number represents the site where 1 = lowland, 3 = highland dolomitic, and 4 = highland granitic; the 2nd letter is e for control or O for OTC, and the last letter represents the plant, a = Cassiope tetragona, b = Dryas integrifolia, c = Salix arctica, and d = Saxifraga oppositifolia. 94 Table 3.2 Species richness o f cultures according to host plant, site, and treatment based on RFLPs OTC Control 1 2 Treatment 3 Total for 1 2 3 Total for species species Site Lowland C. tetragona 0 0 6 6 1 6 6 13 2 D. integrifolia 5 7 16 8 5 5 18 5 7 3 15 4 4 S. arctica 5 13 3 S. oppositifolia 4 3 10 10 1 5 16 Dolomitic C. tetragona 0 0 0 0 0 0 0 0 4 3 5 13 2 D. integrifolia 6 5 13 5 2 7 16 4 2 S. arctica 5 11 1 S. oppositifolia 3 6 10 4 1 10 15 Granitie 0 0 0 0 C. tetragona 9 0 0 9 0 0 0 0 D. integrifolia 0 0 0 0 S. arctica 10 6 5 21 11 10 6 27 9 11 31 S. oppositifolia 11 8 7 7 22 Table 3.3 Species richness of cultures according to host plant, site, and treatment based on Treatment Site Lowland C tetragona D. integrifolia S. arctica S. oppositifolia Dolomitic C tetragona D. integrifolia S. arctica S. oppositifolia Granitic C. tetragona D. integrifolia S. arctica S. oppositifolia 1 2 OTC 3 Total for species 0 3 3 8 0 6 4 10 1 1 3 5 1 10 10 23 0 3 5 9 4 4 2 2 3 5 5 10 7 13 12 21 0 4 3 3 0 2 1 3 0 4 3 4 0 10 7 10 0 0 2 5 0 3 4 2 0 4 1 6 0 6 7 12 0 0 8 12 0 0 7 8 0 0 6 8 0 0 21 28 6 0 8 8 0 0 8 14 0 0 7 8 6 0 23 30 1 2 Control 3 Total for species 95 2. DNA sequencing O f the 50 samples that were sequenced, 43 produced high quality results for analysis. Approximately 330 bp were used for phylogenetic analysis o f the ITS region (ITS) and 350 bp o f the nuclear large subunit rRNA gene (LSU). Five trees were formed based on primer pairs and alignment consensus, which were the ITS and the Dothideales; ITS and the Helotiales; LSU and the Ceratobasidiomycetes; LSU and the Agaricales; and LSU and aseomycetes (Figs 3.4-3.8 respectively). Sequencing was able to classify 37.9% o f the samples to at least Order. The ITS - Dothideales tree had two unknown samples (see Fig. 3.4). One o f the unknown samples was affiliated with Colispora elongata although there was <50% bootstrap support, and the other unknown sample was nested in the Dothideales. E picoccum nigrum 1-S3-C -D rin 5.2 C olispora elongata l-S l-C 3 -S a o p 4 .1 81 Epacris m icrophylla root assoc L e a f litter A sco m y cete 87 D oth id eales sp. 97 99 I Fungal endophyte M U T 2723 Fungal endophyte M U T 585 A sc o m y ce te sp. Phom a glom erata l_om 1 Fig. 3.4 Neighbor-joining best tree based on ITS and the Dothideales. Bootstrap values (1000 replications) >75 are reported. 96 All unknown cultures from the ITS - Helotiales (Fig. 3.5) were found to be affiliated with Cadophora, Phialocephala fortinii, or Hymenoscyphus. One sample was assoeiated with species of Cadophora with a bootstrap value o f 100%; three cultures were clustered with P. fortinii supported by a bootstrap value o f 81%; and two samples were affiliated with Hymenoscyphus even though the bootstrap value was <50%. Cultures affiliated with basidiomyeetes were associated with the Agaricales, whieh had <50% bootstrap support, and with Ceratobasidium (77% bootstrap value) and (Fig. 3.6 and 3.7 respeetively). From the LSU -aseomyeete tree (Fig. 3.8), four eultures were affiliated with Cryptosporiopsis (99% bootstrap value); five cultures with Mycosphaerella (85% bootstrap); three with Phoma (84% bootstrap); two with Cadophora (88% bootstrap); three with Hymenoscyphus (90% bootstrap); and four with Phialocephala (95% bootstrap). Identities o f ten eultures were not resolved from this analysis. Not all samples were resolved to genus and are identified by their family names. One sample could not be resolved even to order and was left as an unknown. 3. Culture morphology Culture morphology was based on 1347 eultures (usually two tubes for eaeh isolate). Forty pereent o f the cultures were identified to genus, whieh included Acremonium, Cryptosporiopsis, Geomyces, Leptodontidium, Monodietys, Phialocephala fortinii, Scytalidium, Sebacina, Staphylotrichum, Trichocladium, Trichoderma, Trichosporiella, or Xenosporium. O f these, Phialocephala fortinii was the most abundant taxon isolated, accounting for 85.7% o f the identified eultures. O f the unknowns, cultures with hyaline hyphae were the most abundant (23%) followed by other isolates with dark hyphae other than 97 Xenochalara juniperi Cadophora malorum " Cadophora luteo-olivacea 100 Cadophora 1-S4-0TC l-Saop4.5 Cadophora gregata Cadophora sp. Leptondontidium orchidicola _&L 2-S4-Cl-Saar4.2 100 2-S4-C t-Saar 5.9 {P h ia lo c e p h a la fo r tin ii 2-Sl-Cl-Sær 2.3 {Phialocephala fortinii ) ectomycorrhizal root (cf. P. fortinii) Phialocephala fortinii r l-Sl-C3Saar 1.3 || Phialocephala fortinii crypt sp Phialocephala fortinii Zalerion varium Laehnum bicolor Glarea lozoyensis Hymenoscyphus epiphyllus 74 H yy menoscyphus monotropae 1-Sl-C3-Saar 1.2 Hymenoscyphus 2-S4-OTCl-Saar 1.9 Phialophora finlandia 100 1 Axenic ericoid root isolate Ecto cf. Hymenoscyphus ericae (= Rhizoscyphus ericae) Eeto isolate (Helotiales) Phialocephala virens Phialophora sp. 0.02 Fig. 3.5 Neighbor-joining best tree based on ITS and the Helotiales. Bootstrap values (1000 replications) >75 are reported. Names in parenthesis after samples are identifications based on morphology. 98 In ocyb e angustispora I F avolasch ia cf. sprucei F avolasch ia cf. calocera I Panellus serotinus Ripartites m etrodii M ycen a haem atopoda M ycen a leaiana 1 -S l-C l-D r in 1.6 Orphalina rivu licola ~ C ollyb ia tuberosa C litocyb e su b velosa 0.02 d Fig. 3.6 Neighbor-joining best tree based on LSU and the Agaricales. Phlebia lindtneri 99 C eratobasidium goodyerae-repentis U thatobasidium fusisporum 99 8 9 C eratobasidium sp. A G O S9 94 77 Ceratobasidium sp. A G L Thanatephorus cucum eris l-S l-O T C 2 -S a a r 2 .1 Inonotus w eirii P hleb iop sis gigantea Fig. 3.7 Neighbor-joining best tree based on LSU and the Ceratobasidiomycetes. Bootstrap values (1000 replications) >75 are reported. 99 --------------------- 2-S4-C 3-Saar2.1 Chalara parvispora Cudonialutea Cryptosporiopsis ericacea A scom ycota sp. Cryptosporiopsis sp Helotiales ^ 99 l - S 4 - C 2 - S a o p 3.3 2 -S l-C l-S a o p 4.1 I l-S 4 -C l-S a o p 4.5 l-S4-C 2-Saop 4.5 Ericoid mycorrhizal sp PPO-8 Ericoid mycorrhizal sp PPO-7 Cenococcum geophilum T rematosphaeria heterospora 90 lodosphaeria sp. 98 _| Pleomassaria siparia a s I h ^ y c o sn h a e re lla m v c o n a n n i Euascomycete RFLP type B Phialophora sp. / 1-S4-C l-C ate 5.8 1-Sl-C 2-D rin 2.2 (P h ialoceph alafortin ii) 96 1| l-S 4-O T C l-S aop 4.1 {P h ialoceph alafortim i) Phialocephala U Phialocephala fortinii 1-S I-C l-S aar 4.4 {P h ialoceph alafortim i) 1-S I-C l-C a te 2.3 Lachnum cf. bicolor Mycorrhizal sp. aD2i Fig. 3.8 Neighbor-joining tree based on LSU and ascomycetes. Species in parentheses are based on morphology. Bootstrap values (1000 replications) >75 are reported. 100 p. fortinii (15%). Success o f isolation o f Phialocephala fortinii differed according to host plant species (p < 0.001), site (p < 0.001), and from which site the host plant was taken (p < 0.001). P. fortinii was found on Salix arctica more than other host plants (60.5%), but the number o f successful isolations plummeted for S. arctica on the dolomitic site (Fig 3.9). The cumulative frequency of isolates for the lowland OTC plots was as high as 34, but only 1-2 isolates were found on the dolomitic plots. With Saxifraga oppositifolia as the host, isolation o f P. fortinii was most successful from the granitic site (~ 6 isolates), whereas the number of 1 20 □ Cate control M Cate OTC I ? 15 V % « 10 G 3 Z S Drin Control -1 1 S D rinO T C U ■ Saar OTC B Saop control m TTinUi kû lowland [D Saar control dolomitic M B Saop OTC granitic Site Fig 3.9 Abundance o f Phialocephala fortinii based on number of morphologically identified cultures differentiated by host plant and treatment. Cate = Cassiope tetragona, Drin = Dryas integrifolia, Saar = Salix arctica, and Saop = Saxifraga oppositifolia. OTC = open top chamber, and control refers to ambient conditions. isolates was low (less than 2) on both the lowland and dolomitic sites. Isolation o f P. fortinii from Cassiope tetragona was low from all sites. Isolation from Dryas integrifolia was most 101 successful from the lowland site with 14 isolates from the control and seven from the OTC plots. Isolating P. fortinii from the two highland sites was not very successful, with only two from the control and one from the OTC from the dolomitic site and no isolates from the granitic site. 4. Identifications based on molecular techniques and morphology When combining sequence and morphological data, 193 (44.7%) o f the cultures remained unknown or inconclusive, 63 (14.6%) were identified as Phialocephala fortinii, 44 (10.2%) as members o f the Helotiales, 29 (6.7%) closest to Mycosphaerella, 23 (5.3%) closest to members o f Cryptosporiopsis, 21 (4.9%) closest to Hymenoscyphus, 12 (2.8%) closest to Phoma, 9 (2.1%) closest to Ceratobasidium, 10 (2.3%) closest to Cadophora, 9 (2.1%) as members o f the Dothideales, 7 (1.6%) closest to Geomyces, 4 (0.9%) as members o f the Agaricales, 2 (0.5%) closest to Colispora, 2 (0.5%) closest to Trichoderma, and 1 eaeh (0.2%) closest to Monodictys, Pénicillium, Sebacina, or Trichocladium. Samples whose morphological identifications were incongruent for the two replicates were considered inconclusive. Samples identified based on molecular data and corresponding morphological identification are listed in A3.1-A3.12. Morphological identification as Phialocephala fortinii matched 34 o f 52 samples (65.4%) that were shown to be affiliated with P. fortinii from sequence/RFLP analyses (Table A3.1). An additional five isolates that were affiliated with the Helotiales based on sequence analysis were morphologically identified as P. fortinii (Table A3.2). When isolates with RFLP patterns that matched these five Helotiales were included in the P. fortinii complex, an additional 15 samples were included in the P. fortinii complex. Morphological identification also suggested that 11 samples that could not be identified from 102 sequence/RFLP analyses were P. fortinii (Table A3.3). Each of these 11 samples had a unique RFLP pattern. There were seven different RFLP types based on identification by sequence/RFLP and morphological analyses, which would total 18 different RFLP patterns for all putatively identified P. fortimi. There were some discrepancies between molecular and morphological identifications. P. fortinii was identified morphologically for samples that were affiliated with Mycosphaerella and members in the Dothideales according to the molecular analysis. The RFLP patterns for these samples matched others in the same cluster but were identified differently morphologically. There were also discrepancies between the two identification methods for Hymenoscyphus, Mycosphaerella, and the Dothideales. One sample, identified morphologically as Leptodontidium, was affiliated with Hymenoscyphus based on sequence data, and belonged to an RFLP pattern that had a wide range o f morphological descriptions (Table A3.4). Another discrepancy included a sample that was identified morphologically as Phialocephala fortinii, whose RFLP pattern matched two samples that were affiliated with Hymenoscyphus according to sequence data. In total, four RFLP types were found for this group. Twenty-nine samples were affiliated with Mycosphaerella, which formed five different RFLP patterns (Table A3.5). Two discrepancies between morphological and molecular techniques were found; one sample was identified as Phialocephala fortinii and the other Staphylotrichum, both o f which have RFLP patterns matching several samples associated with Mycosphaerella. Members that were shown to be affiliated with the Dothideales by the molecular methods had some incongruity with the morphological identification (Table A3.6). Two 103 samples were morphologically identified as Phialocephala fortinii (Helotiales) and one as Acremonium, a member o f the Sordariomycetes. For all the isolates that grouped in the Dothideales, there was only one RFLP pattern, and all the samples were isolated only from Saxifraga oppositifolia. Ten samples were affiliated with Phoma according to sequence and RFLP data, with two different RFLP patterns. Morphologically, Monodictys and Trichocladium were identified for four o f these samples (Table A3.7). Phoma belongs to the mitosporic Ascomycetes, while both Monodictys and Trichocladium are teleomorphs that are members o f the Sordariomycetes. Two o f the Monodictys had the same RFLP pattern as Trichocladium and other members o f the Phoma clade, while the third Monodictys had a RFLP pattern that matched the second RFLP type for Phoma. Four samples were associated with Cryptosporiopsis from sequence analysis, but when samples with matching RFLP patterns were included, the number o f samples increased to 23 (Table A3.8). Three o f these samples were supported by morphology. There were two RFLP types for this clade, which had only a difference o f 50 bp for one o f the fragments in MboW', all other fragments in the three enzymes matched. Ten samples were affiliated with Cadophora according to the molecular analysis, with three RFLP types. Three o f these samples were morphologically identified as P. fortinii (Table A3.9). There were only two groups o f cultures associated with Basidiomycetes, the Agaricales and Ceratobasidium (Tables A3.11-A3.12). All the morphological identifications were labeled as ‘hyaline sterile’, with the exception o f one o f the replicates, which was inconclusive. 104 Most o f the samples that were not identified by sequencing/RFLP data remained unknown based on morphology, but a few samples were identified. Seven unknown samples that were morphologically identified as Geomyces were probably the same species because they all had the same RFLP pattern. Although three samples were identified as Trichocladium, each had a distinct RFLP pattern, which may be an artifact o f the small number o f samples. D. Discussion In this study, site was the most influential factor in distinguishing the root endophytic fungal community composition, more so than warming or host plant, which was evident from the NMS analysis. Most of the members consisted o f DSE, o f which Phialocephala fortinii was the most common. This was similar to the findings o f Stoyke et al. (1992) when they isolated fungi from 10 different subalpine host plants, but contrary to Kohn and Stasovski (1990) who reported no DSE from the same study area. In concurrence with previous research (Kohn and Stasovski 1990, Stoyke et al. 1992), hyaline septate fungi were found in the present study. Thirty-one different species were found on the granitic - OTC plot based on RFLP analysis and 30 species with morphological identification, which was comparable to the species richness o f isolated root endophytes found in temperate grasslands (Wilberforce et al. 2003). Differences in species cumulative frequency depended on site and host plant but not due to passive warming. The dolomitic site had significantly fewer isolates than the other two sites. This may be in part due to the soil conditions o f this site, where CEC, pH, Ca, and C:N ratio were higher than from the other two sites (see Chapter 2), the pH substantially higher, pH 8, versus pH 6 for the other two sites. Interestingly the host plants Salix arctica 105 and Saxifraga oppositifolia yielded the most isolates from the granitie site, but at the lowland site, Dryas integrifolia and Cassiope tetragona produced the most isolates.. Although the granitic site had more isolates than the other two sites, it was the most homogenous, as indicated by the tight clustering found in CCA. This may be because isolates from the granitic site were from S. arctica and S. oppositifolia, with the one exception o f isolates from C. tetragona that were found on one o f the plots; whereas, the other sites had isolates from the other host plants. The high variability from the dolomitic site may be due to the low number o f isolates that were found from this site. Even though root tips from Cassiope tetragona were isolated using a protocol that supposedly favors erieaceous plants, the success rate was low when compared to Salix arctica and Saxifraga oppositifolia. Fungal isolates from Dryas integrifolia, a reputed ectomyeorrhizal plant, were also low. The incubation temperature may also have been too high. Although unlikely, the possibility that isolates from D. integrifolia and C. tetragona may favor colder conditions than those found on the other host plants is possible. Species richness tended to be higher for most plots when assessed by RFLP analysis compared to morphological identification, which was in part because cultures morphologically identified as Phialocephala fortinii had 18 different RFLP types. Many o f the cultures that were identified as P. fortinii had unique RFLP types; therefore, the ITS region for these P. fortinii isolates was highly polymorphic, as were strains from Europe and North America (Harney et al. 1997, Grünig et al. 2002a). Our results were contrary to studies that found morphological differences but no or little RFLP differences for the rDNA ITS region (Addy et al. 2000, Jumpponen 1999). Although Addy et al. (2000) covered a larger geographic range, we had 70 isolates o f P. fortinii compared to their 33 isolates and 34 106 for Jumpponen (1999), which may partially explain the wide range o f RFLP types in our study. Our strains may have more characters that were polymorphic as Grünig et al. (2002a) found with their isolates, which had over five times more polymorphic characters than Addy et al. (2000). Phialocephala fortinii is part o f a species complex composed o f cryptic species, which are species that are morphologically similar but have unique genotypes (Grünig et al. 2004, Piercey et al. 2004). Cryptic species are often found in morphologically asexual fungi that are genetically isolated (Taylor et al. 1999), which would likely be the case for these P. fortinii from Ellesmere Island. Grünig et al. (2004) even found different cryptic species on the same root and clusters of cryptic species that were morphologically indistinct but had high diversity at the population level. In addition, genetic variation o f P. fortinii increases at higher latitudes (Ahlich and Sieber 1996, Piercey et al. 2004); which is consistent with the high number o f RFLP types found for P. fortinii at our site. However, some samples may have been misidentified, which is particularly true for those with RFLP patterns that matched other species, and samples with unique RFLP patterns that were not supported by molecular analyses. This inflated the number o f ribotypes and may be an indication o f phenotypic plasticity for culture identification. Most o f the Phialocephala fortinii were isolated from Salix arctica. Although P. fortinii has been reported to be non-host specific, it was isolated more frequently from S. arctica and Saxifraga oppositifolia than from Cassiope tetragona and Dryas integrifolia. This difference may be a combination o f two factors: 1) P. fortinii is found commonly on roots o f ectomyeorrhizal plants (Stoyke et al. 1992, Jumpponen 1999, Addy et al. 2000, Grünig et al. 2001), which would explain its absence on the ericoid C. tetragona, and 2) 107 although D. integrifolia is an ectomyeorrhizal plant, the root systems from Salix arctica and Saxifraga oppositifolia were larger. Also, P. fortinii grew more frequently on the granitie and lowland sites. The number o f isolates decreased significantly from plants grown on the dolomitic site, perhaps due to the different soil chemistry properties found on this site. Surprisingly, isolates from Saxifraga oppositifolia did not follow the same trend as those from Salix arctica. Even though the abundance o f P. fortinii was comparable from the two hosts from the granitie site, there were significantly more isolates found on Salix arctica on the lowland site. Therefore, both host plant and site appear to affect the abundance o f P. fortinii. Twenty fungal isolates were potentially affiliated with Cryptosporiopsis. Only three o f these samples were identified morphologically; the other samples were identified from RELP patterns. These isolates varied in color, suggesting that morphological differences could be due to different stages o f development. Verkley et al. (2003) found that Cryptosporiopsis isolates from Erieaceous plants that grow on MEA and oatmeal agar change colors at different stages o f development; even though our isolates were grown on MMN, this change could occur on this medium as well. Also identification could have been hampered because maeroeonidia are rarely produced in culture (Verkley et al. 2003). Although Verkley et al. (2003) found that isolation o f Cryptosporiopsis is rare after surface sterilization, we were successful even though we surface sterilized the roots, as did Ahlich and Sieber (1996), who isolated C. radicola from Fagus sylvatica, Abies alba, Picea abies, and Pinus sylvestris. Cryptosporiopsis has teleomorphs in Pezicula and Neofabracea. However, our sequences were based on the primer pair for the large subunit o f the rRNA gene and sequences for this region were not available in GenBank for Pezicula and 108 Neofabracea. Therefore, we could not resolve the teleomorph with which our samples were affiliated. Five isolates were affiliated with Mycosphaerella with an 85% bootstrap value and an additional 23 samples were included with matching RFLP patterns. Mycosphaerella has not been reported from roots in Arctic and Antarctic systems, but Cladosporium, which has a Mycosphaerella teleomorph (Wirsel et al. 2002) and is described as an oligotrophic, melanized fungus (Wirsel et al. 2002, Gunde-Cimerman et al. 2003), is commonly found in the Arctic and Antarctica. In particular C. sphaerospermum, C. herbarum (Gunde-Cimerman et al. 2003, Bergero et al. 1999), and C. cladosporoides (Widden and Parkinson 1979, Robinson 2001, Tosi et al. 2002) are found. Cladosporium has been reported to have a wide range o f functions, from saprotrophic, to epiphytic and endophytic in aboveground plant tissues, and pathogenic and mycoparasitic (Held et al. 2005, Wirsel et al. 2002, David 1997). Our samples were probably pathogenic or were soil fungi that were not killed by surface sterilization, because there have been no reported cases o f root endophytic Cladosporium. However, this needs to be investigated further. Other common fungi found in the Arctic include Geomyces and Phoma (Bergero et al. 1999, Robinson 2001, Tosi et al. 2002). We found seven samples that were morphologically identified as Geomyces with matching RFLP patterns. Ten samples were affiliated with Phoma after sequencing and RFLP analyses. Three o f these samples were morphologically identified as Monodictys. This discrepancy may be due to Phoma having Monodictys - like conidia. When Grondona et al. (1997) described Pyrenochaeta dolchi, which belongs to Phoma section Paraphoma, they found this species as having Monodictys like conidia on lateral branches and also stated that the conidia were similar to Phoma cava 109 and Phoma tracheiphla. These fungi were found in soil by previous researchers (Bergero et al. 1999, Robinson 2001) or from moss (Tosi et al. 2002). Because we surface sterilized our roots, these fungi may be endophytes colonizing the root or saprotrophs not adversely affected by the hydrogen peroxide used for sterilization. Our fungal community may be overly representative o f psychrotrophic fungi. Psychrotrophic and psychrophilic fungi can grow at 0 °C, but psychrotrophic fungi can grow above 20°C, while psychrophilic fungi have an optimum growth temperature at 15 °C or below and cannot grow above 20°C (Robinson 2001). We grew our cultures at room temperature, which may have excluded some o f the psychrophilic fungi. Another factor in determining what fungi were successfully isolated is pH, because it influences the growth of cultures (Taber and Taber 1984). Our isolates were grown on MMN, which has a pH o f 5 to 6 . However, we found later that the pH o f our sites was higher at 6 to 8 . The high pH o f 8 was from the dolomitic site, which may help explain the low number o f fungi that successfully isolated from roots. Discrepancies between morphological and PCR-based identification can be the result o f two factors. The morphological identification may be incorrect, or the identity from the sequence database may be incorrect. Bridge et al. (2003) found that approximately 20% of sequences available in GenBank are unreliable for the ITS region o f rDNA for Amanita and Phoma and the rRNA small subunit for members o f the Helotiales. Our findings, however, were probably accurate because many o f the similar taxa downloaded for analysis had more than one sequence for a species. Problems can arise if the unknown sequence is closely matched to an incorrectly identified sequence. Knowledge o f phylogenetic relationships and taxonomy is important to find these possible misidentified sequences. 110 Overall, warming and host plant did not have a noticeable impact on the composition o f the root endophytic fungal community, which was dictated primarily by site. These plots have only been warmed for approximately five years, which may be too short a time to find differences, especially in the arctic where processes are much slower than those found in temperate forests. Examining the fungal communities after another 5 to 10 years may show more effects o f warming on the root endophytic fungal community. Although host specificity was not found for the overall community, Phialocephala fortinii was preferentially isolated from Salix arctica and Saxifraga oppositifolia over other host plants. Our study identified many o f the isolates; however, most o f the isolates remain unidentified. Even samples that have been sequenced remained unknown, which implies that there are many fungi that have yet to be identified in these ecosystems. The role o f these fungi was not determined although many have speculated that root endophytes play an important role for nutrient uptake by plants in harsh environments. Bioassay studies using some o f these isolates would help elucidate some o f the roles o f these fungi. Ill IV. The root-associated fungal community along a directional, non-replacement succession chronosequence from the Canadian High Arctic A. Introduction Facilitation is a type o f complementary effect, a theory to explain how resource use by organisms can direct ecosystem processes (Cardinale et al. 2002, Loreau and Hector 2001). In this theory, species diversity can increase while avoiding competition, and species are able to co-exist especially in environments with limiting resources by: 1) partitioning resources, where each species can use nutrients, water, or other resources differently instead o f all species competing for or using the same resources; or 2 ) niche differentiation, where different species avoid using the same resources as other organisms by means o f time and/or space (McKane et al. 2002). In contrast, selection effect theory is where species diversity is correlated with the chance that a dominant species uses most o f the available resources (Cardinale et al. 2002), so the formation o f the community is heavily dependent on this one central species. Facilitation is thought to play an important role for the primary succession of mycorrhizal fungi in that nutrient poor conditions are improved by pioneer species, as indicated by increasing diversity o f arbuscular mycorrhizal fungi after volcanic disturbances (Titus and del Moral 1998), and o f ectomycorrhizae on glacial tills (Helm et al. 1999), and by the increase o f nitrogen and organic matter along a chronosequence from a sub-alpine glacier (Jumpponen et al. 1998). This facilitative process is also inferred by the successional pattern in the mycorrhizal status of plants, which begins with non-mycorrhizal plants, to arbuscular mycorrhizal, ectomyeorrhizal (Read 1993), and/or ericoid mycorrhizae (Cazares 1992). Given that complementary effect is the main theory to explain resource use by the root112 associated fungal community, selection effect theory cannot be totally discarded. The possiblility that one or two dominant species may use most o f the resources and influence how the community structure will develop subsists. The objectives o f this study were to examine how the mycorrhizal fungal community composition, species diversity and abundance change according to host plant and chronosequence in a directional, non-replacement succession. Directional, non-replacement succession has been used to describe a type o f plant succession, where species found in the youngest plots are also found in older plots; they are not replaced by different plant species (Svoboda and Henry 1987). This type o f succession serves as a natural laboratory where the increase of biodiversity can be examined. This study is the first to our knowledge that examines the root-associated fungal community in this type o f succession. Also, previous research on succession in the arctic occurred in the low arctic where trees are present (Helm et al. 1999, Helm et al. 1996, Brubaker et al. 1995, Chapin et al. 1994). Ellesmere Island is located in the high arctic where only low shrubs grow. B. Materials and methods 1. Study site Samples were collected from the receding western lobe o f the Twin glacier located at the south end o f an 8 -km^ lowland high arctic plant oasis on Ellesmere Island, Nunavut, Canada (78°53’N, 75°55’W) at the end o f July, 2001. The western lobe o f the Twin glacier, which started to advance since the Little Ice Age and began to recede approximately 1960, is diminishing approximately 10 m/yr (GHR Henry, pers. comm.). Plots were placed within zones delineated by the time since exposure: 1990 to present [labeled 1990 plots], 1980-90 [1980 plots], 1970-80 [1970 plots], and 1960-70 [1960 plots]. The area before 1959 was not 113 covered by the western lobe o f the glacier (GHR Henry, pers. eomm.), so plants collected from this area were used as controls. 2. Field collection Plants were seleeted that fulfilled the directional, non-replacement succession model (i.e. once the plant appeared within the chronosequence, it was present for all remaining chronosequence zones including the eontrol). The plants chosen were Luzula confusa Lindeb., Papaver lapponicum (Tolm.) Nordh., Salix arctica Pall., Saxifraga oppositifolia L., Cassiope tetragona (L.) D.Don, and Dryas integrifolia Vahl. Luzula confusa and P. lapponicum were first present in the 1990 plots (Y90), and they were the dominant plants for 1990, 1980, and 1970 plots. S. arctica first appeared in the 1980 plots (Y80), became more abundant and larger in the 1970 plots (Y70), and was dominant in the 1960 plots (Y60). S. oppositifolia was first noted in the 1970 plots, and C. tetragona and D. integrifolia in the 1960 plots. The control (60 C a te Fig 4.1 Canonical Correspondence Analysis for effects o f host plant and chronosequence zone. Rotated 345°. Number preceded by ‘Y ’ represent the earliest time o f exposure from the glacier. 0 .5 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 Rank genotype abundance (b)ITS A-Hinü 3 .5 3 § 'O s t 2 .5 2 OTC 1.5 — Control CD o 1-1 0 .5 0 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 Rank genotype abundance Figure A2.6 Rank abundance curve for highland dolomitic site. Relative abundance is log scaled. 176 (c)ITS \-Alu\ OTC Control 0 t - t - - r : I [ I T T 1 ' 1 ' 1 1 ' 1 1 1 1 ' r i ' f ■ I ' 1 1 1 ■ 1 i 1 i h \ Rank genotype abundance (d)ITS 4-AM 3 .5 -------- — OTC C ontrol 0 .5 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 R a n k g en o ty p e ab u ndance Fig. A2.6 (cont’d) Rank abundance curve for highland dolomitic site. Relative abundance is log scaled. 177 (a) ITS \-HinÛ 3.5 2.5 OTC Control 0.5 1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 Rank genotype abundance (b)ITS 4-//m fl 2 .5 § 13 OTC § Control 0.5 1 3 5 7 9 11 13 15 17 19 21 23 25 2 7 2 9 31 33 35 Rank genotype abundance Fig. A2.7 Rank abundance curves for highland granitic site. Relative abundance is log scaled. 178 (c)ITS \-Alu\ 3.5 aw 3 w 2.5 "9 2 I s OTC ■Control 1.5 0.5 0 1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 Rank genotype abundance (d)ITS A-AM 4 3.5 U 3 fl ft T3 2.5 e i OTC 2 ■control w 1.5 o 1 0.5 0 1 3 5 7 9 11 13 15 17 1921 2 3 2 5 2 7 2931 33 35 3 73 94 1 Rank genotype abundance Fig. A2.7 (cont’d) Rank abundance curves for highland granitic site. Relative abundance is log scaled. 179 (a) ITS \-H inû 35.00 « 30.00 I 25.00 •C 20.00 □ Control *‘ 15.00 10.00 J 5.00 4 BOTC 0.00 Lowland Dolomitic highland Granitic highland Site (b)ITS A-HinÛ 30.00 « 25.00 e X 20.00 h « 15.00 È" o 10.00 □ Control BOTC 5.00 0 00 Lowland Dolomitic h^hland Granitic hi^iland Site 40.00 35.00 I 30.00 4 I 25.00 □ Control ! g, 20.00 - B OTC I " 15.00 J 10.00 5.00 0.00 Lowland Dolomitic h i ^ n d Granitic hi^iland Site (d)ITS 4 -^ M 35.00 « 30.00 I 25.00 •C 2 0 .0 0 I □ C o n tro l' I BOTC 15.00 I 10.00 ^ 5.00 0.00 Lowland Site Fig. A2.8. Least square means o f genotype richness for each treatment per site. Error bar represents one standard error. 180 (a) ITS 0.94 S 0.92 - □ Control HOTC Lowland Dolomitic hi^iland Granitic hi^iland Site (b)ITS A-HinÛ I.GO 0.95 0.90 □ Control B O TC 0.85 0.80 0.75 Lowland Dolomitic hi^iland Granitic hi^iland Site (c)ITS \-Alu\ 1.20 1.00 g 0.80 ; □ Control 0.60 - 0.20 0.00 Lowland Dolomitic h^Jiland Granitic hÿ^tland Site (d)ITS A-AM g 1.00 0.98 ' 0.96 0.94 ; □ Control 0.90 0.86 0.84 0.82 Lowland Dolomitic highland Site Fig. A2.9. Least square means o f genotype evenness measured by Hurlbert’s PIE index. Error bar represents one standard error. 181 VIII. Appendix for chapter 3 Table A3.1 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Phialocephala. Identifications by sequencing and morphology ID based on R FL P and sequence analyses: Phialocephala fortinil RFLP type Ipf 2-S l-C l-D rin3.2 2-C l-S l-S aar2.3 2-S1-OTC 1-Saar 2.1 2-S1-OTC 1-Saar 2.2 RFLP type 2pf ID based on m orphology dark sterile Phialocephala fortinii dark sterile NA 1-Sl-O TCl-Drin 1.1 Phialocephala fortinii l-Sl-O TC2-Saop 5.5 1-Sl-0TC3-Drin 1.3 Phialocephala fortinii NA l-S4-O TC l-Saop4.1 2-Sl-C2-Drin 1.1 Phialocephala fortinii Phialocephala fortinii 2-Sl-C2-Saar 1.1 2-S 1-OTC 1-Saar 1.1 2-S 1-OTC 1-Saar 1.2 2-S 1-OTC 1-Saar 3.1 2-S 1-OTC 1-Saar 4.1 2-Sl-OTC2-Drin2.1 Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii 2-S l-0T C 2-S aar2.2 2-S4-C3-Saar 4.6 RFLP type 3pf 1-Sl-C l-D rin2.1 1-SI-C l-Saar 4.4 l-Sl-O TC 3-C atc4.1A 1-Sl-0TC 3-C atc4.3 1-S4-Cl-Catc3.8 1-S4-C1-Cate 5.7 1-S4-Cl-Catc5.8 l-S4-C2-Saop 4.2 l-S4-O TC l-Saop4.2 2-S4-C1-Saar 2.5 pigmented warty/ dark sterile Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii NA Phialocephala fortinii Phialocephala fortinii NA Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii 182 Table A3.1 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Phialocephala. Identifications by sequencing and morphology ID based on R FL P and sequence analyses : Phialocephala fortinii 2-S4-Saar 5.18 2-S4-C1-Saar 5.4 2-S4-C1-Saar 5.9 2-S4-C2-Saop 4.1a 1-S4-C3-Saar3.4 2-S4-CTL-Cate 2.2 2-S4-OTC3-Saop 1.3 2-S4-OTC3-Saop 5.3 RFLP type 4pf 1-Sl-C2-D rin2.2 1-Sl-C3-Cate 5.3 1-Sl-0TC 2-Saar3.2 1-Sl-0TC2-Saar4.11 1-Sl-0TC2-Saar4.6 l-Sl-O TC2-Saop 5.5a 1-Sl-0TC 3-C ate5.1 RFLP type 5pf 1-Sl-C3-Saar 1.3 1-Sl-C3-Saar 4.2 1-Sl-O TCl-Saar 2.1 1-Sl-0TC 3-C ate3.1 2-Sl-OTC3-Saar5.1 3-Sl-C3-Saar 4.2 RFLP type 6pf 1-Sl-O TC l-Saar 5.1 l-S4-Cl-Saop2.3 2-S4-Cl-Saar 1.2 2-S4-Cl-Saar4.2 ID based on m orphology Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii NA Phialocephala fortinii NA dark sterile Phialocephala fortinii Phialocephala fortinii NA Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii NA Phialocephala fortinii NA Phialocephala fortinii NA Phialocephala fortinii Phialocephala fortinii NA dark septate sterile NA dark sterile Phialocephala fortinii 183 Table A3.2 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Helotiales. Sequence analysis indicated that ID based on morphology ID based on RFLP and sequence analyses: Helotiales RFLP Type 1h 1-Sl-C2-Cate 1.4 1-S4-Ctl-Cate 4.2 1-S4-Ctl-Cate 4.5 NA dark monolioid NA 2-S4-Cl-Saar 1.1a NA 2-S4-C3-Saop 3.2 NA RFLP Type 2 h l-S4-O TCl-Saop 5.9 1-S4-C2-Saar2.7 RFLP Type 3h 1-Sl-C2-Drin5.2 l-S4-OTC2-Saop2.1 2-Sl-OTC2-Saar2.1 2-S4-C3-Saop3.5 2-S4-OTC2-Saar2 2-S4-OTC3-Saar5.2 3-S l-C 3-Saarl.l RFLP Type 4 h l-S4-OTC2-Saop4.2 2-S3-C2-Saop5.1 2-S3-OTCl-Drin 1.2 2-S3-OTCl-Drin 1.4 RFLP Type 5h 1-S4-C2-Saar2.1 1-S4-C2-Saar2.2 RFLP Type 6h 1-S4-Cl-Cate 3.7 l-S4-C l-Saop 1.11 inconclusive inconclusive dark monolioid/ mixed NA NA hyaline sterile NA mixed/ hyaline sterile mixed inconclusive NA dark sterile NA hyaline conidia/mixed NA NA hyaline sterile 184 Table A3.2 Comparison of identification (based on molecular techniques) with corresponding morphological identification - Helotiales. Sequence analysis indicated that ID based on RFLP and sequence analyses: Helotiales RFLP Type ? h 1-Sl-C3-D rin2.4 l-Sl-OTC3-Cate-4.2 l-S4-C l-Saop 5.4 1-S4-C2-Saar5.1 l-S4-C3-Saop5.1 1-S4-0TC1-Saar 4.5 l-S4-O TC l-Saop3.7 l-S4-O TCl-Saop 4.4 1-S4-0TC2-Saar2.1 l-S4-OTC2-Saop4.1 l-S4-OTC2-Saop4.3 2-Sl-OTC2-Drin 1.1 2-S4-C1-Saar 5.20 2-S4-C2-Saar 4.9 2-S4-C2-Saop 3.2 2-S4-C2-Saop 3.3 2-S4-C3-Saop 1.1 2-S4-C3-Saop 3.1 2-S4-OTCl-Saar 1.10 2-S4-OTC1-Saar 2.1 ID based on morphology dark septate monolioid NA NA NA dark septate monolioid dark sterile/ mixed NA NA Phialocephala fortinii dark monolioid Phialocephala fortinii Phialocephala fortinii dark septate monolioid Phialocephala fortinii hyaline sterile inconclusive dark septate monolioid Phialocephala fortinii dark sterile dark septate monolioid 185 Table A3.3 Comparison of identification (based on molecular techniques) with corresponding morphological identification - unknown. Identifications by sequencing were inconclusive. ID based on RFLP and sequence ID based on morphology analyses: Unknowns l-Sl-C 3-Saop 5.4 l-S3-C3-Saop3.3 2-S3-C3-Saar4.1 Ascomycete brown brown 2-S4-C3-Saopl.5 brown 1-S4-0TC2-Saar4.1 Cadophora dark monolioid dark monolioid dark monolioid 1-Sl-C3-Drin 1.3 1-S3-C-Drin 1.1 2-Sl-C2-Saar3.1 2-S3-C2-Saop 5.7 2S4-Cl-Saar 1.16 2-S4-OTC1-Saar 4.4 2-Sl-OTC2-Saop 2.1 l-S4-O TC l-Saop3.3 1-Sl-C3-Drin 4.1 1-Sl-0TC1-Saar5.1 l-Sl-O TC3-Saop 2.2 l-S3-C3-Saop5.2 l-S4-C l-Saop2.3 1-S4-Ctl-Cate 2.7 2-Sl-OTC2-Saar 1.1 2-S3-C3-Drin 4.2 2-S3-C3-Saop4.2 2-S3-OTCl-Saar 1.8 2-S3-OTC3-Drin4.1 2-S4-C2-Saarl.2 1-S4-C2-Saar3.2 1-S4-CTL Cate 2.1 dark monolioid dark monolioid dark monolioid dark septate, few monolioids dark some monolioid dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile dark sterile Geomyces Geomyces 186 Table A3.3 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - unknown. Identifications by sequencing were ID based on RFLP and sequence analyses: Unknowns 1-S4-0TC1-Saar 4.11 1-S4-0TC1-Saar 4.4 l-S4-OTC2-Saop3.1 2-S4-C2-Saop 4.4B 2-S4-C2-Saop 4.8 1-S3-0TC3-Drin 3.3 2-S4-C3-Saar2.1 l-S4-OTC3-Saop 4.3 2-S4-C3-Saar 1.5 3-Sl-C2-Cate-4.1 1-Sl-C3-Cate3.1 1-Sl-C l-D rin2.3 1-S l-C l-D rin2.5 1-Sl-0TC 3-D rin -3.2 1-Sl-0TC 3-Saar3.1 1-Sl-0TC 3-Saar3.2 l-Sl-O TC 3-Saop4.1 1-S3-C2-Drin2.1 1-S3-C-Drin2.1 l-S3-OTC2-Drin2.2B l-S3-OTC3-Saar 1.2 1-S4-C1-Cate 4.1 1-S4-C1-Saar 3.2 l-S4-C l-Saop 1.11 1-S4-C2-Saar4.2 1-S4-C2-Saar4.3 1-S4-Ctl-Cate 4.7 1-S4-0TC1-Saar 3.10 l-S4-OTC2-Saop 4.4 l-S4-OTC3-Saop4.4 ID based on morphology Geomyces Geomyces Geomyces Geomyces Geomyces Hyaline arthroconidia hyaline brown; swollen cells hyaline conidia, monolioid hyaline hyphae, dark spored hyaline monolioid hyaline spherical, ovoid conidia hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile 187 Table A3.3 Comparison o f identification (based on molecular techniques) with corresponding morphological identification. Identifications by sequencing were inconclusive ID based on RFLP and sequence analyses: Unknowns l-S4-OTC3-Saop4.6 l-S4-OTC3-Saop4.9 2-S3-C3-Drin3.3 2-S3-C3-Saop 2.1 2-S3-OTC3-Saar 1.3 2-S4-C1-Saar 3.20 2-S4-C3-Saar 1.3 2-S4-C3-Saop 1.7 2-S4-C3-Saopl.2 2-S4-C3-Saop3.5 2-S4-OTC2-Saop 1.3 2-S4-OTC2-Saop 3.2 2-S4-OTC3-Saar2.1 3-Sl-C2-Cate2.1 3-Sl-OTC3-Drin 1.2 2-S3-OTC3-Saop 4.2 2-S4-OIC3-Saar 2.4 1-Sl-C2-Cate4.1 1-S1-C3-Cate 3.2 1-Sl-C3-D rin2.2 1-Sl-C3-D rin4.2 l-Sl-C 3-Saop2.1 1-S3-Cl-Saar3.4 l-S3-C l-Saop2.2 1-S3-C2-Drin5.5 1-S3-C3-Drin3.1 l-S3-C3-Saop 1.8 l-S3-C3-Saop 4.4 1-S3-0TC1-Saar3.8 1-S4-C1-Cate 4.3 l-S4-C l-Saop 1.7 l-S 4-C l-Saop5.6 1-S4-C2-Saar 1.1 1-S4-0TC1-Saar 4.7 ID based on morphology hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile/ hyaline monolioid hyaline sterile/ hyaline monolioid inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive 188 Table A3.3 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - unknowns. Identifications by sequencing were ID based on RFLP and sequence analyses: Unknowns ID based on morphology l-S4-O TCl-Saop 1.1 l-S4-O TCl-Saop 1.7 l-S4-O TCl-Saop 3.4 l-S4-O TC l-Saop5.5 l-S4-OTC3-Saop4.5 2-S l-C l-D rin4.5 inconclusive inconclusive inconclusive inconclusive inconclusive inconclusive 2-Sl-C l-Saop 3.3 2-S l-C l-S aop4.4 2-S4-Cl-Saar 1.12 2-S4-Cl-Saar 1.5 2-S4-Cl-Saar 1.6 2-Sl-O T C l-Saop2.4 inconclusive inconclusive l-S4-C2-Saop2.3 l-Sl-C 3-Saop 1.2 1-Sl-Cl-Cate2.1 1-Sl-Cl-Catc2.3 1-Sl-C l-Saop 5.1 1-Sl-C2-Cate 4.2 1-Sl-C3-Catc4.1 1-Sl-C3-Catc5.1 1-SI-OTC 1-Saar 1.3 1-Sl-OTCl-Saop 5.2 1-Sl-0TC2-D rin 1.1 1-Sl-0TC 2-D rin4.1 l-Sl-O TC2-Saop 5.4 1-Sl-0TC3-Cate 5.2 1-Sl-0TC3-Catc4.2 1-Sl-0TC3-Cate 5.2 1-Sl-0TC 3-D rin2.1 l-Sl-O TC3-Saar 1.2 1-S3-Cl-D rin4.2 1-S3-C3-Drin 1.3 inconclusive inconclusive inconclusive light brown sterile lightly pigmented sterile Monodictys NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 189 Table A3.3 Comparison of identification (based on molecular techniques) with corresponding morphological identification. Identifications by sequencing were inconclusive ID based on RFLP and sequence analyses: Unknowns l-S3-C3-Saop2.1 l-S3-C3-Saop2.3 l-S3-C3-Saop4.1 l-S3-C3-Saop4.2 l-S3-C3-Saop4.3 l-S3-C3-Saop 5.3 1-S3-C-Drinl.2 1-S3-0TC1-Drin2.1 l-S3-OTC3-Drin2.2 l-S3-OTC3-Saop2 l-S3-OTC3-Saop2.1 l-S3-OTC3-Saop4.5 1-S4-Cl-Cate 3.7 1-S4-Cl-Cate 5.2 1-S4-C2-Saar2.2 1-S4-C2-Saar2.8 l-S4-C2-Saop5.1 1-S4-CTL-Cate 1.6 1-S4-CTL-Cate 2.5 1-S4-CTL-Cate 2.6 1-S4-Ctl-Cate 3.3 1-S4-CTL-Cate 3.5 1-S4-CTL-Cate5.1 1-S4-0TC1-Saar 4.15 l-S4-O TCl-Saop 1.5 l-S4-O TC l-Saop5.2 l-S4-OTC3-Saop 5.2 2-S l-C l-D rin4.1 2-S l-C l-D rin4.2 2-S 1-C l-Saar 2.1 2-S 1-Cl-Saar-3.1 2-Sl-C l-Saop 5.2 ID based on morphology NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 190 Table A3.3 Comparison of identification (based on molecular techniques) with eorresponding morphological identification - unknowns. Identifications by sequencing were ID based on RFLP and sequence analyses: Unknowns 2-S l-C l-Saop5.3 ID based on morphology 2-S3-C3-Saar 3.5 2-S3-C3-Saop 1.1 2-S3-C3-Saop 3.3 2-S3-C3-Saop 5.2 2-S3-OTCl-Saop3.5 2-S3-OTC2-Saar 3.5 2-S3-OTC2-Saar 3.7 2-S3-OTC3-Saar 3.2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 2-S3-OTC3-Saar4.1 2-S3-OTC3-Saar 4.3 2-S3-OTC3-Saop4.1 2-S3-OTC3-Saop 4.4 2-S3-OTC3-Saop 4.5 2-S4-C1-Saar 1.11 2-S4-Cl-Saar 1.14 2-S4-C1-Saar 1.3 NA NA NA NA NA NA NA NA 2-Sl-C l-Saop 5.5 2-Sl-C 2-D rin2.2 2-Sl-C 2-D rin3.2 2-Sl-C2-Saar2.1 2-S3-Cl-D rin3.1 2-S3-C1-Saar 2.1 2-S3-C1-Saar 3.1 2-S3-Cl-Saar 5.1 2-S3-C2-Drin 5.2 2-S3-C2-Drin 1.1 2-S3-C2-Drin 1.2 2-S3-C2-Drin 1.4 2-S3-C2-Drin 5.2 2-S3-C2-Saar 5.2 2-S3-C3-Drin3.1 2-S3-C3-Saar 1 2-S3-C3-Saar 2.5 191 Table A3.3 Comparison o f identification (based on molecular techniques) with corresponding morphological identification. Identifications by sequencing were inconclusive ID based on RFLP and sequence analyses: Unknowns ID based on morphology 2-S4-C1-Saar 2.4 NA NA NA NA NA NA NA 2-S4-C2-Saar 2.4 2-S4-C2-Saop 4 .IB 2-S4-C3-Saar 4.3 2-S4-C3-Saar 4.4 2-S4-C3-Saop 2.3 2-S4-C3-Saop 2.4 2-S4-OTC2-Saar 1.2 2-S4-OTC3-Saar 1.2 2-S4-OTC3-Saop 5.6 3-Sl-C2-Catel.3 3-Sl-OTC3-Drin 1.1 3-Sl-O TC3-Drin5.2 1-S4-OTC1-Saar 4.12 1-Sl-C l-D rin4.1 1-Sl-0TC 2-Saar5.5 1-Sl-0TC 3-D rin5.3 l-S3-OTC3-Saop 4.3 l-S4-O TCl-Saop 4.3 l-S4-OTC3-Saop 4.1 2-Sl-O T C l-Saop2.3 2-Sl-O TC2-Saar3.3 2-S3-OTC3-Saar 4.2 2-S4-C1-Saar 5.1 2-S4-OTC3-Saar 3.4 2-S3-O TCl-D rin4.5 2-Sl-C l-Saop 5.6 2-S4-OTCl-Saar 1.5 2-S4-OTC1-Saar 2.5 2-S3-OTC3-Drin3.1 2-S4-C3-Saar 2.3 2-S4-OTC2-Saar 4.3 NA NA NA NA NA NA Pénicillium Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Phialocephala fortinii Sebacina Trichocladium Trichoderma sporulosum Trichoderma sporulosum white white white 192 Table A3.4 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Hymenoscyphus. Samples showed affiliation ID based on RFLP and sequence analyses: Hymenoscyphus ID based on morphology RFLP type Iny l-S4-C l-Saop2.11 1-S4-0TC1-Saar3.5 dark spherical conidia dark sterile RFLP type 2 hy 1-Sl-C3-Saar 1.2 1-Sl-0TC 1-D rin2.3 intercalary chlamydospore fried egg black type 1-Sl-0TC 1-D rin4.1 1-Sl-0TC 2-D rin3.1 1-S4-CTL-Cate 2.4 inconclusive Leptodontidium 1-S4-0TC2-Saar4.6 2-S l-C l-S aop3.2 2-Sl-C 3-D rin3.2 2-S4-C2-Saar4.10 3-Sl-C3-Saar2.1 RFLP type 3 hy l-S4-O TC l-Saop5.7 2-S4-OTCl-Saar 1.9 RFLP type 4 hy 1-S l-C l-D rin2.2 1-Sl-0TC3-Cate 4.2a 1-S4-Cl-Saar 4.6 2-S4-Cl-Saar 1.16 2-Sl-C 3-D rin3.2 inconclusive inconclusive NA inconclusive NA black type inconclusive inconclusive Phialocephala fortinii NA dark sterile dark monolioid NA 193 Table A3.5 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Mycosphaerella. Samples showed affiliation ID based on RFLP and sequence analyses: Mycosphaerella RFLP type 1m 1-Sl-C2-Cate 5.3 1-S3-0TC1-Drin4.4 ID based on morphology l-S3-O TCl-Saop 3.4 NA inconclusive white l-S3-O TC l-Saop3.6 inconclusive 2-S3-C3-Saop 1.2 white 2-S4-C2-Saar 5.7 2-S4-C2-Saop 2 2-S4-C2-Saop 2.2 2-S4-C2-Saop 2.4 RFLP type 2 m l-S3-C l-Saop4.19 l-S4-OTC3-Saop2.1 RFLP type 3 m l-Sl-C 3-Saop 3.3 1-S3-C Drin 4.9 1-S4-C1-Saar 4.3 1-S4-C2-Saar3.3 1-S4-0TC1-Saar 4.6 inconclusive NA Staphylotrichum inconclusive 2-S4-C2-Saop 44A 2-S4-C2-Saop 4.9 2-S4-OTC1-Saar 4.1 RFLP type 4 m 2-S3-OTCl-Saar 1.9 1 D/G-Saar 3.3 RFLP type 5m 2-S3-C3-Saar 1.7 2-S3-C3-Saar 1.8 2-S3-C3-Saar3.1 2-S3-OTC1-Saar 4.8 2-S3-OTC1-Saar 5.1 2-S3-0TCl-Saar 5.5 1-D/G-Saar 3.4 dark sterile NA inconclusive inconclusive inconclusive Phialocephala fortinii inconclusive dark septate inconclusive NA NA dark sterile NA NA NA NA NA NA dark sterile 194 Table A3 . 6 Comparison of identification (based on molecular techniques) with corresponding morphological identification - Dothideales. Samples showed affiliation to the ID based on morphology ID based on RFLP and sequence analyses: Dothideales l-Sl-C 3-Saop 4.1 1-Sl-O TCl-Saop 1.2 NA 1-Sl-O TCl-Saop 2.1 inconclusive l-Sl-O TC2-Saop 2.2 Phialocephala fortinii l-Sl-O TC2-Saop 5.3 Phialocephala fortinii l-S4-O TC l-Saop3.9 2-S l-C l-Saop 1.6 2-Sl-OTC2-Saop 3.3 hyaline Phialocephala fortinii -like dark sterile Acremonium 2-Sl-OTC2-Saop5.1 black type NA Table A3.7 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Phoma. Samples showed affiliation to Phoma ID based on RFLP and sequence analyses: Phoma RFLP type 1? l-S3-OTC2-Drin2.2A l-S 4-C l-Saop4.7 l-S4-C2-Saop 2.2 l-S4-OTC2-Saop L IB l-S4-OTC2-Saop 2.2 l-S4-OTC3-Saop 3.9 l-S4-OTC3-Saop 4.2 2-S3-OTC3-Saar 2.3 2-S3-OTC3-Saar 3.3 RFLP type 2? l-Sl-O TC 2-Saop3.2 l-Sl-O TC 3-Saop3.2 l-S 3-C l-Saop4.2 ID based on morphology NA NA inconclusive Monodictys Monodictys NA Trichocladium hyaline sterile dark sterile inconclusive Monodictys inconclusive 195 Table A3 . 8 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Cryptosporiopsis. Samples showed affiliation to Cryptosporiopsis according to seq uence and RFLP analyses ID based on RFLP and sequence ID based on morphology analyses: Cryptosporiopsis RFLP type 1l l-Sl-C 3-Saop 3.2 dark sterile NA 1-S4-Cl-Cate 2.2 NA 1-S4-Cl-Cate 2.4 1-S4-C2-Saar3.4 Cryptosporiopsis radicola l-S4-C2-Saop 3.2 brown inconclusive l-S4-C2-Saop3.3 inconclusive l-S4-C2-Saop 3.4 inconclusive l-S4-C2-Saop4.5 inconclusive 1-S4-C3-Saar3.1 1-S4-C3-Saar5.1 l-S4-C3-Saop3.6 1-S4-CTL-Cate 1.2 dark sterile monolioid hyaline sterile hyaline sterile l-S4-OTC2-Saop 3.2 l-S4-OTC3-Saop3.1 l-S4-OTC3-Saop4.11 1-S4-CTL Cate 3.9 2-Sl-C l-Saop 4.1 2-S3-OTC3-Drin 2.2 2-S4-C3-Saar 3.3 2-S4-CTL-Cate 3.6 RFLP type 2 l 1-S3-C2-Drin2.2 1-S3-C2-Drin5.6 inconclusive hyaline sterile inconclusive Cryptosporiopsis radicola inconclusive hyaline sterile white white Cryptosporiopsis radicola NA 196 Table A3.9 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Cadophora. Samples showed affiliation to ID based on morphology ID based on RFLP and sequence analyses; Cadophora RFLP Type Ica 2-S3-C3-Saop 2.2 l-S4-O TC l-Saop4.5 dark sterile NA RFLP Type 2ca 1-SI-C l-Saar 3.6 1-SI-OTC 1-Drin 4.2 1-Sl-0TC 3-C ate4.1 Phialocephala fortinii Phialocephala fortinii inconclusive RFLP Type 3ca 1-S3-0TC2-Drin2.4 1-S4-C2-Saar2.10 1-S4-Ctl Cate 2.1 A 2-S4-CÜ Cate 3.1 2-S4-CÜ Cate 3.3 inconclusive hyaline sterile Phialocephala fortinii inconclusive NA Table A3.10 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Colispora. Samples showed affiliation to Colispora according to sequence anc RFLP analyses ID based on RFLP and sequence ID based on morphology analyses: Colispora 1-S3-C-Drin5.2 l-S4-OTC3-Saop 5.4 NA NA Table A3.11 Comparison o f identification (based on molecular techniques) with corresponding morphological identification - Agarieales. Samples showed affiliation to the ID based on RFLP and sequence analyses: unknown Agarieales 1-SI-C l-D rin 1.6 1-S4-C1-Saar 4.9 2-S4-C1-Saar 3.11 2-S4-C2-Saop 4.3 ID based on morphology hyaline sterile hyaline sterile hyaline sterile hyaline sterile 197 Table A3.12 Comparison o f identification (based on molecular techniques) with corresponding morphological identification -Ceratobasidium. Samples showed affiliation to ID based on RFLP and sequence analyses: Ceratobasidium 1-Sl-0TC 2-Saar2.1 1-Sl-OTC2-Saar2.10 1-S3-C2-Saar 1.3 1-S3-C2-Saar2.3 1-S3-C2-Saar5.3 l-S3-OTC2-Saop2.1 1-S4-Cl-Cate 5.4 2-Sl-C 2-Saop2.2 2-Sl-OTC2-Drin3.1 ID based on morphology hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile hyaline sterile inconclusive hyaline sterile 198 IX. Appendix for chapter 4 Table A4.1 Genotype richness for host plants at chronosequence zone. Number preceded by ‘Y ’ represent the earliest time o f exposure from the glacier. posure Y ear o f exposure (c)Exchangeable Na (d) Exchangeable Mg * 1.6 Y60 g5 5 4 4 3 2 2 Y90 Y 80 Y 70 Y 60 < Y 59 ■ M ean 0 Y90 Y e a r o f exposure Y80 Y70 Y60 4) osu e 0.5 A- 2 .5 ln < Y 5 9 -a - o -Q- a -a- -o 17 19 21 1 23 3 5 Rank genotype abundance 7 9 11 13 15 17 19 21 23 25 27 R a n k g e n o ty p e a b u n d a n c e (c) ITS 4-//m fI (d) ITS A-Alu\ 2 .5 2 6- 1.5 In Y 90 ÿ - • Û - - In Y 9 0 -■ + - - l n Y 8 0 + --ln Y 8 0 0.5 ■a-- - In Y 70 a - * - - l n Y 60 — InY ôO ln < Y 5 9 2 0 .5 — 0 3 5 7 9 Rank genotype abundance 3 5 Fig. A4.6 Rank abundance curve for each chronosequence zone for Papaver lapponicum.. w to 7 9 Rank genotype abundance 13 In Y 7 0 ln < Y 5 9 (a) ITS \-H in ü (b) ITS \-A lu\ 3 2.5 2.5 ■-lnY80 2 ■a B 3 • • -H- - • In Y 7 0 ■■■lnY70 1.5 --ln Y 6 0 - In You ec — h