Open Access

Global climate change is occurring all over the world, but in the Arctic the climate is changing more rapidly and drastically than in many other parts of our planet. Many species that are already at their climatic limit need to adapt to recent climate conditions or migrate in order to not go extinct. The possibilities of adaption include phenotypic plasticity and adaptation to various extents. This is also the case for white spruce P. glauca, which belongs to the conifers and thus in the largest group of gymnosperms still living today. Among the approx. 600 extant conifer species white spruce is one of the most widespread trees in North American boreal forests. Its range extends from 69° N in the Canadian Northwest Territories to the Great Lakes at about 44° N, where it occurs from sea level to an altitude of about 1520 m (Burns and Honkala, 1990). Site related, climate-dependent differences in white spruce reproduction can be seen as a strategy to survive under the harsh climatic conditions at Alaska's treelines: Besides sexual reproduction, the vegetative propagation occurs in the white spruce as an additional reproductive mechanism. This can be realized by "layering" when the lower branches of the tree crown touch the ground and develop roots to later grow as a separate individual with or without a connection to the mother tree. Known as other mechanisms of vegetative propagation are also the rooting of fallen trees which were not completely uprooted, and the "root suckering", in which new shoots sprout from the roots of the tree. However, the latter was not yet observed in the genus Picea. With the help of short, repetitive, non-coding sequences in the genome, which are therefore not subject to selection and are called microsatellites, these clones can be determined by genotyping. For this purpose, using different polymorphic microsatellites, an individual multilocus genotype is created for each tree, by means of which it can be compared with all other trees of the same species. In the first part of this work (article I), the occurrence of clones in three study areas at Alaskan treelines are examined and the reasons for their appearance in variable numbers are discussed. For this purpose, 2571 white spruces (P. glauca) were genotyped and their position was determined via differential GPS in the field. The percentage of clonal trees is higher in areas with harsh climatic conditions and correlates with the height of the lowest branches of the tree crown. This suggests that the vegetative propagation of white spruce is a backup strategy for times when climatic conditions hamper sexual reproduction. The correlation between clone numbers and tree crown height suggests "layering" as the main mechanism for cloning whereas selection for vegetative reproduction seems to be very unlikely shown by the results for genetic differentiation between the clonal and the singleton trees in this study. In the second part of this work (articles II and III), the influence of environmental factors and phenotypic traits on the mycobiome of the needles (including all fungi living on (epiphytic) and in (endophytic) the needles) in our study areas in Alaska was investigated. The mycobiome of the white spruce needles was chosen as a proxy for the parasite infection rate by fungi and thus serves as a fitness parameter. For this purpose, all epiphytic and endophytic fungal species were analyzed by a metabarcoding analysis. In article II, 48 trees of one study area at Alaska’s northern treeline (Brooks Range) were examined for differences in mycobiome due to genetic differentiation, phenotypic characteristics and / or habitat characteristics. The trees used for this study were sampled from two adjacent plots on a south-facing mountain slope with an elevation gradient from 875 to 950 meters above sea level. It could be shown that, in contrast to the trees genotype, the height above sea level, the mountain slope, as well as the height and age of the trees have a significant impact on the mycobiome. The genetic differentiation between the tree individuals, however, showed no significant effect. Based on article II we examined the mycobiome composition of a total of 96 trees in 2 plots (16 trees each) at three sites in Alaska over a distance of 500 kilometers. Additionally, we sampled needles of two different ages for each tree (current year and three years old needles) summing up to 192 samples in total. The incentive of this study (article III) was to investigate the influence of origin and age of spruce needles on their mycobiome and if there is a genetic predisposition that is related to the fungal species community. In addition, the sampling design was improved by collecting needles from all four orientations (North, South, East and West) and sampling trees at a standardized distance to each other to avoid systematic errors. Comparable to article II the influence of the trees genetics on the species community of the epiphytic and endophytic fungi of the white spruce needles seems to be very unlikely. In contrast, a significant influence of the geographic origin and the needle age on the species structure of the needle inhabiting fungal species was found. The phenotypic tree traits height and dbh (diameter at breast height) had only minor influence and did in fact explain less than 2% of the mycobiome variance. Using Illumina sequencing, 10.2 million reads from the nucleotide sequence between the internal transcribed spacer (ITS) genes could be obtained, which yielded in 1575 ribotypes (called operational taxonomic unit, OTU) for the fungi. These were compared with a reference database to compare and assign them to known fungal species. For example, 942 OTUs with >95% similarity could be identified as known species, with 1975 samples identified on genus level and 2683 when determined to family level. The most pronounced difference between the two studies (article II and III) were due to the fungal species of the class of Pucciniomycetes, more specifically the genus Chrysomyxa which belongs to the rust fungi and is plant pathogenic. In the study of article II (sampling in 2012), Pucciniomycetes accounted for only a minor portion of the assigned DNA sequences. In the second study (article III, sampling in 2015) they accounted for more than half of all basidiomycetes found, which in turn contain 20.0% of all DNA sequences, the second largest phylum found beside Ascomycetes (51.4%).

White spruce (Picea glauca (Moench) Voss) is one of the most common conifers in Alaska and various treelines mark the species distribution range. Because treelines positions are driven by climate and because climate change is estimated to be strongest in northern latitudes, treeline shifts appear likely. However, species range shifts depend on various species parameters, probably most importantly on phenotypic plasticity, genetic adaptation and dispersal. Due to their long generation cycles and their immobility, trees evolved to endure a wide variety of climatic conditions. In most locations, interannual climate variability is larger than the expected climate change until 2100. Thus treeline position is typically thought of as the integrated effect of multiple years and to lag behind gradual climate change by several decades. Past dendrochronological studies revealed that growth of white spruce in Alaska can be limited by several climatic variables, in particular water stress and low temperatures. Depending on how the intensity of climate warming, this could result in a leading range edge at treelines limited by low temperatures and trailing treelines where soil moisture is or becomes most limiting. Climate-growth correlations are the dendrochronological version of reaction norms and describe the relationship between an environmental variable and traits like tree-ring parameters (e.g. ring width, wood density, wood anatomy). These correlations can be used to explore potential effects of climate change on a target species. However, it is known that individuals differ with respect to multiple variables like size, age, microsite conditions, competition status or their genome. Such individual differences could be important because they can modulate climate-growth relationships and consequently also range shifts and growth trends. Removing individual differences by averaging tree-ring parameters of many individuals into site chronologies could be an oversimplification that might bias estimates of future white spruce performance. Population dynamics that emerge from the interactions of individuals (e.g. competition) and the range of reactions to the same environmental drivers can only be studied via individual tree analyses. Consequently, this thesis focuses on factors that might alter individual white spruce’ climate sensitivity and methods to assess such effects. In particular, the research articles included explore three topics: 1. First, clones were identified via microsatellites and high-frequency climate signals of clones were compared to that of non-clonal individuals. Clonal and non-clonal individuals showed similar high-frequency climate signals which allows to use clonal and non-clonal individuals to construct mean site chronologies. However, clones were more frequently found under the harsher environmental conditions at the treelines which could be of interest for the species survival strategy at alpine treelines and is further explored in the associated RESPONSE project A5 by David Würth. 2. In the second article, methods for the exploration and visualization of individual-tree differences in climate sensitivity are described. These methods represent a toolbox to explore causes for the variety of different climate sensitivities found in individual trees at the same site. Though, overlaying gradients of multiple factors like temperature, tree density and/or tree height can make it difficult to attribute a single cause to the range of reaction norms (climate growth correlations). 3. Lastly, the third article attempts to disentangle the effect of age and size on climate-growth correlations. Multiple past studies found that trees of different Ages responded differently to climatic drivers. In contrast, other studies found that trees do not age like many other organisms. Age and size of a trees are roughly correlated, though there are large differences in the growth rate of trees, which can lead to smaller trees that are older than taller trees. Consequently, age is an imperfect Proxy for size and in contrast to age, size has been shown to affect wood anatomy and thus tree physiology. The article compares two tree-age methods and one tree-size method based on cumulative ring width. In line with previous research on aging and Wood anatomy, tree size appeared to be the best predictor to explain ontogenetic changes in white spruce’ climate sensitivity. In particular, tallest trees exhibited strongest correlations with water stress in previous year July. In conclusion, this thesis is about factors that can alter climate-growth relationships (reaction norms) of white spruce. The results emphasize that interactions between climate variables and other factors like tree size or competition status are important for estimates of future tree growth and potential treeline shifts. In line with previous studies on white spruce in Alaska, the results of this thesis underline the importance of water stress for white spruce. Individuals that are taller and that have more competitors for water appear to be most susceptible to the potentially drier future climate in Alaska. While tree ring based growth trends estimates of white spruce are difficult to derive due to multiple overlaying low frequency (>10 years) signals, all investigated treeline sites showed highest growth at the treeline edge. This could indicate expanding range edges. However, a potential bottleneck for treeline advances and retreats could be seedling establishment, which should be explored in more detail in the future.

Myxomycetes are protists belonging to the super-group Amoebozoa. The traditional taxonomic system, which is now largely outdated by molecular studies, recognizes five orders: Liceales, Trichiales, Physarales, Stemonitales and Echinosteliales. Molecular phylogenies revealed two basal clades: Physarales and Stemonitales (the so-called dark-spored myxomycetes) are the first; the other above-mentioned orders form the second (the bright-spored myxomycetes). However, except for Echinosteliales none of the traditional orders appears to be monophyletic in the traditionally used delimitation. The dark-spored myxomycetes encompass the majority of the described morphospecies. Due to the high genetic divergence in DNA sequences between the bright- and dark-spored myxomycetes, only the latter are considered in this dissertation. Historically myxomycetes have been described as fungi, due to their macroscopically visible fructifications which, though considerably smaller, resemble those of fungi. These fruit bodies provide enough morphological traits to support a morphological species concept with currently ca. 1000 species described. Therefore diversity studies of myxomycetes have been conducted for over 200 years and a substantial body of data on ecology and distribution of these fructifications exist. From these studies myxomycetes are known to form often distinct communities across terrestrial ecosystems with highly specific habitat requirements, such as snowbanks (nivicolous), herbivore dung (coprophilous) or decaying wood (xylophilous). However knowledge on the myxamoebae – the trophic life stage of the myxomycetes – is very scarce. Only recent advances in molecular techniques such as direct species identification based on DNA sequences from environmental samples (ePCR), have made studies of myxamoebae (and other microbes) possible. From these first molecular based studies myxomycetes are currently estimated to account for between 5 to almost 50% of all soil amoebae, and have been shown to be present in a wide variety of soils. To fully take advantage of these new methods, a molecular DNA marker needs to be established as well as a reference sequence database. The usability of a DNA marker gene depends on its ability to separate species by a distinction between intra- and interspecific divergence between sequences of the same and related species, the so-called ‘barcoding gap’. The first part of this thesis (article I and II) deals with the subject of establishing such a DNA marker and database, and in doing so touches upon the subject of ‘what is a myxomycetes species?’ A total of 1 200 specimens were compiled into a reference database (the largest database to date of dark-spored myxomycetes). The genetic distance from sequence-to-sequence was used to assess genetic clade structures within morphospecies and putative biospecies (sexually isolated linages) were identified. The result was an estimate of hidden diversity, exceeding that of described morphospecies by 99%. The optimum sequence similarity threshold for OTU-picking (genetic species differentiation, denoted Operational Taxonomic Unit) with the used SSU marker was identified as 99.1% similarity. The second part of this thesis (article III and IV) presents ecological studies conducted with NGS (ePCR) in which the established threshold and database are applied and are demonstrated to provide reliable and novel insights into the soil myxamoebae community. It is investigated whether the occurrence of fruit bodies reflects the distribution of soil myxamoebae, and the research questions ‘do myxomycetes show broader realized niches as soil amoebae than as fructifications?’ and ‘are myxamoebae distributions correlated to potential prey organisms (fungi and bacteria)?’ are investigated. In the ecological study presented in article III parallel metabarcoding of bacteria, fungi and dark-spored myxomycete was used for the first time in a joint approach to analyze the communities from an elevational transect in the northern limestone German Alps (48 soil samples). Illumina sequencing of the soil samples revealed 1.68 Mio sequences of a section of the rRNA gene, which were assigned to 578 operational taxonomic units (OTU) from myxomycetes. These show a high similarity (>98%) to 42 different morphospecies (the respective figures for bacteria and fungi were 2.16/5710/215 and 3.68/6133/260, respectively). Multivariate analyses were carried out to disentangle microbial interplay and to identify the main environmental parameters determining the distribution of myxamoebae and thus setting the boundaries for their ecological niches. Potential interactions between the three target organisms were analysed by integrating community composition and phylogenetic diversity with environmental parameters. We identified niche differentiation for all three communities (bacteria, fungi and myxamoebae) which was strongly driven by the vegetation. Bacteria and fungi displayed similar community responses, driven by symbiont species and plant substrate quality. Myxamoebae showed a more patchy distribution, though still clearly stratified among genera, which seemed to be a response to both structural properties of the habitat and specific bacterial taxa. In addition we find an altitudinal species turn-over for all three communities, most likely explained by adaptation to harsh environmental conditions. Finally a high number of myxomycetes OTUs (associated with the genus Lamproderma) not currently represented in our reference database were found, representing potentially novel species. This study is the first to report niche differentiation between the guild of nivicolous (“snowbank”) myxomycetes and thus fine-scale niche differentiation among a predatory soil protist; identifying both potential food preferences and antagonistic interactions with specific bacterial taxa. Finally, the second ecological study (article IV) focuses on comparing the distribution of myxamoebae revealed by ePCR of soil samples with fructifications collected from the same area (714 specimens determined to 30 morphospecies, which form 70 unique ribotypes that can be assigned to 45 ribotype clusters using a 99.1% similarity threshold). The study found a strong coherency between the two inventories, though with species specific relative differences in abundance, which can in part be attributed to the visibility of the fructifications. In addition, a year to year comparison of fructification records gives support to the hypothesis that the abundance of fructifications depends strongly on the onset of snowfall in the previous autumn and the soil temperature regime throughout the winter.

Die mitteleuropäischen Flachbärlappe (Gattung Diphasiastrum) sind in Deutschland alle hochgradig gefährdet und können ohne geeignete Artenhilfsmaßnahmen hier nicht dauerhaft überleben. In der vorliegenden Arbeit werden die Grundlagen für ein Artenhilfsprogramm geschaffen, indem die Reproduktionsbiologie untersucht worden ist und die ökologischen Ansprüche und Gefährdungsursachen ermittelt wurden, um entsprechende Hilfsmaßnahmen für die Arten zu entwickeln. Um Rückschlüsse auf das Reproduktionssystem der Eltern- und Hybridarten zu erhalten, wurde die genetische Diversität ermittelt. Dabei kam das fingerprinting-Verfahren AFLP zum Einsatz, womit Arten und genetisch verschiedene Individuen voneinander abgegrenzt werden können. Verwendet wurden die beiden Primer-Kombinationen EcoRI-AAG / VspI-CT und EcoRI-ACT / VspI-CAG. Die größte genetische Diversität der Elternarten weist D. complanatum auf, die überwiegend durch echte Fremdbefruchtung (outcrossing) entsteht, während die genetische Diversität der beiden anderen Elternarten D. alpinum und D. tristachyum gering ist und nur knapp oberhalb einer definierten Fehlerrate liegt. Die Proben der Hybridarten unterscheiden sich so stark voneinander, dass davon ausgegangen werden muss, dass dies immer wieder neu entstehende F1-Hybriden sind, wenngleich die Unterschiede bei D. oellgaardii vergleichsweise gering aufgrund der geringen genetischen Diversität der Elternarten D. alpinum und D. tristachyum ist. Die Sporenproduktion in den Sporenständen wurde unter anderem direkt durch Zählung von Sporen in den Sporangien unterm Stereomikroskop und Zählung der Sporangien in den Sporenständen ermittelt. Sporangien von L. clavatum enthalten demnach 27735 (± 7492) Sporen pro Sporangium, 105,2 (± 3,7) Sporangien pro Sporenstand und hochgerechnet etwa 2,1 bis 3,8 Millionen Sporen pro Sporenstand. Die terminale Fallgeschwindigkeit liegt für L. clavatum bei 2,16 (± 0,11) cm*s-1 und für D. complanatum bei 2,25 (± 0,10) cm*s-1. Die Sporen haben einen Durchmesser von 29,5 (± 2,0) μm bzw. 32,3 (± 2,5) μm. Die gemessene Geschwindigkeit liegt deutlich unter der theoretischen und lässt sich damit erklären, da Sporen keine perfekten Kugeln sind und aufgrund ihrer stark reliefierten Oberfläche Turbulenzen erzeugt werden, die den Fall verlangsamen. Die Anzahl der in Entfernungen bis 200 m zu einer sporenbildenden Population fliegenden Sporen wurde mithilfe von vertikalen klebenden Sporenfallen bei D. complanatum, D. tristachyum und L. clavatum bestimmt. Nur für die Population von L. clavatum mit 11358 reifen Sporenständen auf kleiner Fläche wurden weitere Berechnungen durchgeführt. Folgende Funktion beschreibt die Anzahl der durch die Luft fliegenden Sporen in einer Höhe von etwa 40 cm über dem Boden in Abhängigkeit zur Entfernung x: f(x) = 45878*x-2,302 (R² = 0,9979). Es konnten selbst in 200 m Entfernung noch einzelne Sporen an den Sporenfallen nachgewiesen werden. Da die maximale theoretische Ausbreitungsdistanz bei nur knapp 130 m liegt, selbst wenn ein konstant horizontal wehender Wind von 100 km/h angenommen wird, müssen aufwärtsgerichtete Luftströmungen eine entscheidende Rolle bei der Fernausbreitung spielen. Die Ansiedlungsversuche wurden im Thüringer Schiefergebirge am Grünen Band bei Brennersgrün (D. alpinum und D. tristachyum) und im Pöllwitzer Wald (D. complanatum) durchgeführt. Als Vergleichsart wurde wieder L. clavatum verwendet. Die Sprossverpflanzungen verliefen insgesamt erfolgreich mit einer Überlebensrate von 8% für D. alpinum, 17% für D. complanatum, 8% für D. tristachyum und 22% für L. clavatum. Der jährliche Rhizomzuwachs liegt bei 0,5 cm, 13,3 (± 1,8) cm, 7,5 cm bzw. 9,1 (± 4,0) cm für die entsprechenden Arten. Die Vegetationsbedeckung vorher abgeplaggter Flächen liegt zwischen 39 und 53% nach zwei Jahren, jedoch mit großen Unterschieden selbst zwischen benachbarten Flächen und wird hauptsächlich durch ein schnelles Mooswachstum bestimmt. Nach etwa fünf Monaten sind keine Sporen auf sterilem Nährstoffmedium gekeimt, obwohl diese unterschiedlich behandelt wurden, zum Beispiel mit Rauchgas, Hitze, konzentrierter Schwefelsäure oder durch Mörsern. Auch die Keimungsversuche an den Wuchsorten waren nach 2,5 Jahren erfolglos. Eine Erklärung kann ein Dormanzstadium unbekannter Dauer vor der Keimung sein. Der Anzahl der jährlich gebildeten vertikalen Sprossbüschel wurde indirekt für je eine Population für D. zeilleri (2,5/Jahr) und D. issleri (2,0/Jahr) bestimmt, indem der Quotient aus der Anzahl der Sprossbüschel an der längsten Rhizomverbindung und dem bekannten Alter des Standorts ermittelt worden ist. Die Zugehörigkeit von Rhizomstücken zu einem Klon wurde mit der AFLP-Methode abgesichert. Die meisten Populationen in Deutschland werden demnach mehrere Jahrzehnte alt, jedoch ohne beobachtete Verjüngung über Prothallien. Eine Erklärung könnten die immer noch sehr geringen pH-Werte in den Unterböden von 3,6 (± 0,24) sein, die giftige Al3+-Ionen pflanzenverfügbar machen.

Leaf-inhabiting fungi are a hyperdiverse group of microbiota found in all terrestrial habitats. Comparative studies targeting the drivers of endophytic fungal biodiversity are rare and identified multiple effectors, such as plant chemistry, climate and seasonal attributes. Our project aimed to study the pattern of the leaf-associated mycobiome of European beech (Fagus sylvatica) at altitudinally distinct sites to reveal diversity, composition and seasonal dynamics of fungal endophytes by a combination of metabarcoding, cultivation and subsequent ecological analyses. This thesis also intended to study the fungal relationship with biotic and abiotic factors: elevation, local site conditions, leaf biochemistry and leaf status. Metabarcoding and cultivation were applied for same leaf samples to trace both environmental drivers and method-dependent signals of the detected fungi. An experimental field site consisting of 100 (2-years old beech) trees was established called ‘beech phytometer’ system at two altitudes (517 and 975 m a.s.l.) in a German mountain forest. Beech trees were randomly selected from both sites as well as from neighboring beech trees. Ten trees from each site were chosen and 10 leaves per tree were sampled. Climatic and leaf biochemistry (Chlorophyll, flavonoid and nitrogen) data were seasonally (Autumn, Spring and Summer) investigated for two continuous years (Oct 2013 to Oct 2015) at these two elevations. In the first year (autumn, 2013) of the project (chapter 3.1), the leaf-inhabiting fungi of natural beech trees were investigated by using high-throughput sequencing (metabarcoding) at three altitudinally distinct sites (with timberline at 1381 m a.s.l.) in the German Alps. This paper focuses on a detailed description and evaluation of metabarcoding amplicon library preparation and a subsequent analytical workflow. Fungal diversity and community composition were compared as a function of different elevated sites and leaf status (i.e., vital or senescent). However, three investigation sites resulted in 969 OTUs (operational taxonomic units) from 820441 sequences. Taxonomic compositions (order) of beech fungi differed strongly among the three sites but were less distinct between the vital and the senescence leaves. Fungal community composition at valley site clearly differed from those of mountain and timberline where differences between mountain and timberline were less prominent. Vital and senescence leaf differed in fungal community structures indicating a strong dynamics of leaf fungi in autumn. Elevation and leaf status were found to be the main explaining factors, which affected the fungal richness and compositions. Another survey (Chapter 3.2) was conducted just after the establishment of the ‘beech phytometer’ trees in the same period (autumn 2013) where leaf mycobiome of the phytometer trees (trees originally came from Northern Germany and grown in nursery) were compared with the fungi of surrounding natural beech habitat at valley (517 m a.s.l.) and mountain site (975 m a.s.l.) in the same location “Untersberg”. Fungal diversity was lowest in the managed habitat in the nursery and was highest in natural habitat. Fungal diversity and compositions significantly associated the origin of the trees. Under natural conditions, the fungi were more diverse at lower altitudes than at higher altitudes. Additionally, leaf chlorophyll and flavonoid contents showed negative correlations with fungal richness in natural stands. In the second year (autumn 2014), another survey (chapter 3.3) was conducted on leaf endophytes of phytometer trees with metabarcoding and cultivation approaches to trace the environmental drivers and method-dependent indications. Metabarcoding resulted in 597 OTUs from 170480 curated ITS1 reads and cultivation revealed 70 OTUs from 438 culture-based Sanger sequences. Both approaches resulted in non-overlapping community compositions and pronounced differences in taxonomic classification and trophic stages. However, both methods revealed similar correlations of the fungal communities with local environmental conditions. Our results indicate undisputable advantages of metabarcoding over cultivation in terms of representation of the major functional guilds, rare taxa and diversity signals of leaf-inhabiting fungi. This stressed out the importance of cultivation for complementing sequence databases with good quality reference data and encouraged the use of both approaches in future microbial biodiversity assessment studies. Phytometer and natural trees were intensively investigated in this study (chapter 3.4) to assess the influence of site characteristics (altitudes, local microclimate), seasonality, leaf biochemistry and leaf age on fungal diversity and composition. In total, our analytical Illumina workflow resulted in 15703599 demultiplexed and ITS1 reads from 165 samples. Clustering at 97% similarity resulted in 1199 OTUs. Climatic parameters were significantly differed between valley and mountain on daily basis but were insignificantly differed on monthly basis. The compositional difference between phytometer and natural mycobiome was significant for combined data as well as for the seasonal data (Oct 2013-Oct 2014). We observed a strong seasonal turnover in phyllosphere fungi in both habitats over the two years of investigation, suggesting that the plant-fungal system not only responds to cyclic climatic conditions but depends as well on various parameters, e.g., geographic position, substrates age and surrounding vegetation. A side (chapter 3.5) study was done to see the connection between the foliar endophytes and foliar phenolic compounds of European aspen (Populus tremula) in the presence and absence of specialist beetles (Chrysomela tremula). A distinct pattern of the leaf endophytes was found to be associated with aspen genotype and chemotype, but this specificity disappeared in the presence of herbivorous beetles. This suggested that leaf endophytes responded to the herbivory in aspen. In general, the altitudinal difference is the most important explaining factor for fungal community differences, which shapes many dependent abiotic and biotic habitat factors. Regarding cost and time per sequence, metabarcoding is superior to cultivation approaches and offers surprisingly profound insights by yielding much more data, allowing to test at once multiple hypotheses in fungal ecology.