Substantial compositional turnover of fungal communities in an alpine ridge-to-snowbed gradient

Experimental evidence that root-associated fungi improve plant growth at high altitude

Unravelling how species communities change along environmental gradients requires a dual understanding: the direct responses of the species to their abiotic surroundings and the indirect variation of these responses through biotic interactions. Here, we focus on the interactive relationships between plants and their symbiotic root-associated fungi (RAF) along stressful abiotic gradients. We investigate whether variations in RAF community composition along altitudinal gradients influence plant growth at high altitudes, where both plants and fungi face harsher abiotic conditions. We established a translocation experiment between pairs of Bistorta vivipara populations across altitudinal gradients. To separate the impact of shifting fungal communities from the overall influence of changing abiotic conditions, we used a root barrier to prevent new colonization by RAF following translocation. To characterize the RAF communities, we applied DNA barcoding to the root samples. Through the utilization of joint species distribution modelling, we assessed the relationship between changes in plant functional traits resulting from experimental treatments and the corresponding changes in the RAF communities. Our findings indicate that RAF communities influence plant responses to stressful abiotic conditions. Plants translocated from low to high altitudes grew more when they were able to associate with the resident high-altitude RAF compared to those plants that were not allowed to associate with the resident RAF. We conclude that interactions with RAF impact how plants respond to stressful abiotic conditions. Our results provide experimental support that interactions with RAF improve plant stress tolerance to altitudinal stressors such as colder temperatures and less nutrient availability.

Mycena species can be opportunist-generalist plant root invaders

Traditional strict separation of fungi into ecological niches as mutualist, parasite or saprotroph is increasingly called into question. Sequences of assumed saprotrophs have been amplified from plant root interiors, and several saprotrophic genera can invade and interact with host plants in laboratory growth experiments. However, it is uncertain if root invasion by saprotrophic fungi is a widespread phenomenon and if laboratory interactions mirror field conditions. Here, we focused on the widespread and speciose saprotrophic genus Mycena and performed (1) a systematic survey of their occurrences (in ITS1/ITS2 datasets) in mycorrhizal roots of 10 plant species, and (2) an analysis of natural abundances of 13 C/15 N stable isotope signatures of Mycena basidiocarps from five field locations to examine their trophic status. We found that Mycena was the only saprotrophic genus consistently found in 9 out of 10 plant host roots, with no indication that the host roots were senescent or otherwise vulnerable. Furthermore, Mycena basidiocarps displayed isotopic signatures consistent with published 13 C/15 N profiles of both saprotrophic and mutualistic lifestyles, supporting earlier laboratory-based studies. We argue that Mycena are widespread latent invaders of healthy plant roots and that Mycena species may form a spectrum of interactions besides saprotrophy also in the field.

First report of the ectomycorrhizal fungal community associated with two herbaceous plants in Inner Mongolia, China

Ectomycorrhizal (EM) fungi play a vital role in ensuring plant health, plant diversity, and ecosystem function. However, the study on fungal diversity and community assembly of EM fungi associated with herbaceous plants remains poorly understood. Thus, in our study, Carex pediformis and Polygonum divaricatum in the subalpine meadow of central Inner Mongolia, China were selected for exploring EM fungal diversity and community assembly mechanisms by using llumina MiSeq sequencing of the fungal internal transcribed spacer 2 region (ITS2). We evaluated the impact of soil, climatic, and spatial variables on EM fungal diversity and community turnover. Deterministic vs. stochastic processes for EM fungal community assembly were quantified using β-Nearest taxon index scores. The results showed that a total of 70 EM fungal OTUs belonging to 21 lineages were identified, of which Tomentella-Thelephora, Helotiales1, Tricholoma, Inocybe, Wilcoxina were the most dominant EM fungal lineages. EM fungal communities were significantly different between the two herbaceous plants and among the two sampling sites, and this difference was mainly influenced by soil organic matter (OM) content and mean annual precipitation (MAP). The neutral community model (NCM) explained 45.7% of the variations in EM fungi community assembly. A total of 99.27% of the β-Nearest Taxa Index (βNTI) value was between -2 and 2. These results suggest that the dominant role of stochastic processes in shaping EM fungal community assembly. In addition, RC_bray values showed that ecological drift in stochastic processes dominantly determined community assembly of EM fungi. Overall, our study shed light on the EM fungal diversity and community assembly associated with herbaceous plants in the subalpine region of central Inner Mongolia for the first time, which provided a better understanding of the role of herbaceous EM fungi.

Microbial diversity patterns in the root zone of two Meconopsis plants on the Qinghai-Tibet Plateau

In the extreme alpine climate of the Qinghai-Tibet Plateau (QTP), plant growth and reproduction are limited by extremely cold temperatures, low soil moisture, and scarce nutrient availability. The root-associated microbiome indirectly promotes plant growth and plays a role in the fitness of plants on the QTP, particularly in Tibetan medicinal plants. Despite the importance of the root-associated microbiome, little is known about the root zone. This study used high-throughput sequencing to investigate two medicinal Meconopsis plants, M. horridula and M. integrifolia, to determine whether habitat or plant identity had a more significant impact on the microbial composition of the roots. The fungal sequences were obtained using ITS-1 and ITS-2, and bacterial sequences were obtained using 16S rRNA. Different microbial patterns were observed in the microbial compositions of fungi and bacteria in the root zones of two Meconopsis plants. In contrast to bacteria, which were not significantly impacted by plant identity or habitat, the fungi in the root zone were significantly impacted by plant identity, but not habitat. In addition, the synergistic effect was more significant than the antagonistic effect in the correlation between fungi and bacteria in the root zone soil. The fungal structure was influenced by total nitrogen and pH, whereas the structure of bacterial communities was influenced by soil moisture and organic matter. Plant identity had a greater influence on fungal structure than habitat in two Meconopsis plants. The dissimilarity of fungal communities suggests that more attention should be paid to fungi-plant interactions.

Both abundant and rare fungi colonizing Fagus sylvatica ectomycorrhizal root-tips shape associated bacterial communities

Ectomycorrhizal fungi live in close association with their host plants and form complex interactions with bacterial/archaeal communities in soil. We investigated whether abundant or rare ectomycorrhizal fungi on root-tips of young beech trees (Fagus sylvatica) shape bacterial/archaeal communities. We sequenced 16S rRNA genes and fungal internal transcribed spacer regions of individual root-tips and used ecological networks to detect the tendency of certain assemblies of fungal and bacterial/archaeal taxa to inhabit the same root-tip (i.e. modularity). Individual ectomycorrhizal root-tips hosted distinct fungal communities associated with unique bacterial/archaeal communities. The structure of the fungal-bacterial/archaeal association was determined by both, dominant and rare fungi. Integrating our data in a conceptual framework suggests that the effect of rare fungi on the bacterial/archaeal communities of ectomycorrhizal root-tips contributes to assemblages of bacteria/archaea on root-tips. This highlights the potential impact of complex fine-scale interactions between root-tip associated fungi and other soil microorganisms for the ectomycorrhizal symbiosis.

Host generalists dominate fungal communities associated with alpine knotweed roots: a study of Sebacinales

Bistorta vivipara is a widespread herbaceous perennial plant with a discontinuous pattern of distribution in arctic, alpine, subalpine and boreal habitats across the northern Hemisphere. Studies of the fungi associated with the roots of B. vivipara have mainly been conducted in arctic and alpine ecosystems. This study examined the fungal diversity and specificity from root tips of B. vivipara in two local mountain ecosystems as well as on a global scale. Sequences were generated by Sanger sequencing of the internal transcribed spacer (ITS) region followed by an analysis of accurately annotated nuclear segments including ITS1-5.8S-ITS2 sequences available from public databases. In total, 181 different UNITE species hypotheses (SHs) were detected to be fungi associated with B. vivipara, 73 of which occurred in the Bavarian Alps and nine in the Swabian Alps-with one SH shared among both mountains. In both sites as well as in additional public data, individuals of B. vivipara were found to contain phylogenetically diverse fungi, with the Basidiomycota, represented by the Thelephorales and Sebacinales, being the most dominant. A comparative analysis of the diversity of the Sebacinales associated with B. vivipara and other co-occurring plant genera showed that the highest number of sebacinoid SHs were associated with Quercus and Pinus, followed by Bistorta. A comparison of B. vivipara with plant families such as Ericaceae, Fagaceae, Orchidaceae, and Pinaceae showed a clear trend: Only a few species were specific to B. vivipara and a large number of SHs were shared with other co-occurring non-B. vivipara plant species. In Sebacinales, the majority of SHs associated with B. vivipara belonged to the ectomycorrhiza (ECM)-forming Sebacinaceae, with fewer SHs belonging to the Serendipitaceae encompassing diverse ericoid-orchid-ECM-endophytic associations. The large proportion of non-host-specific fungi able to form a symbiosis with other non-B. vivipara plants could suggest that the high fungal diversity in B. vivipara comes from an active recruitment of their associates from the co-occurring vegetation. The non-host-specificity suggests that this strategy may offer ecological advantages; specifically, linkages with generalist rather than specialist fungi. Proximity to co-occurring non-B. vivipara plants can maximise the fitness of B. vivipara, allowing more rapid and easy colonisation of the available habitats.

Glacier retreat in the High Arctic: opportunity or threat for ectomycorrhizal diversity?

Climate change causes Arctic glaciers to retreat faster, exposing new areas for colonization. Several pioneer plants likely to colonize recent deglaciated, nutrient-poor areas depend on fungal partners for successful establishment. Little is known about general patterns or characteristics of facilitating fungal pioneers and how they vary with regional climate in the Arctic. The High Arctic Archipelago Svalbard represents an excellent study system to address these questions, as glaciers cover ∼60% of the land surface and recent estimations suggest at least 7% reduction of glacier area since 1960s. Roots of two ectomycorrhizal (ECM) plants (Salix polaris and Bistorta vivipara) were sampled in eight glacier forelands. Associated ECM fungi were assessed using DNA metabarcoding. About 25% of the diversity was unknown at family level, indicating presence of undescribed species. Seven genera dominated based on richness and abundance, but their relative importance varied with local factors. The genus Geopora showed surprisingly high richness and abundance, particularly in dry, nutrient-poor forelands. Such forelands will diminish along with increasing temperature and precipitation, and faster succession. Our results support a taxonomical shift in pioneer ECM diversity with climate change, and we are likely to lose unknown fungal diversity, without knowing their identity or ecological importance.

Soil P reduces mycorrhizal colonization while favors fungal pathogens: observational and experimental evidence in Bipinnula (Orchidaceae)

Little is known about the soil factors influencing root-associated fungal communities in Orchidaceae. Limited evidence suggests that soil nutrients may modulate the association with orchid mycorrhizal fungi (OMF), but their influence on non-mycorrhizal fungi remains unexplored. To study how nutrient availability affects mycorrhizal and non-mycorrhizal fungi associated with the orchid Bipinnula fimbriata, we conducted a metagenomic investigation within a large population with variable soil conditions. Additionally, we tested the effect of phosphorus (P) addition on fungal communities and mycorrhizal colonization. Soil P negatively correlated with the abundance of OMF, but not with the abundance of non-mycorrhizal fungi. After fertilization, increments in soil P negatively affected mycorrhizal colonization; however, they had no effect on OMF richness or composition. The abundance and richness of pathotrophs were negatively related to mycorrhizal colonization and then, after fertilization, the decrease in mycorrhizal colonization correlated with an increase in pathogen richness. Our results suggest that OMF are affected by soil conditions differently from non-mycorrhizal fungi. Bipinnula fimbriata responds to fertilization by altering mycorrhizal colonization rather than by switching OMF partners in the short term, and the influence of nutrients on OMF is coupled with indirect effects on the whole fungal community and potentially on plant's health.

Habitat specialisation controls ectomycorrhizal fungi above the treeline in the European Alps

Alpine habitats are one of the most vulnerable ecosystems to environmental change, however, little information is known about the drivers of plant-fungal interactions in these ecosystems and their resilience to climate change. We investigated the influence of the main drivers of ectomycorrhizal (EM) fungal communities along elevation and environmental gradients in the alpine zone of the European Alps and measured their degree of specialisation using network analysis. We sampled ectomycorrhizas of Dryas octopetala, Bistorta vivipara and Salix herbacea, and soil fungal communities at 28 locations across five countries, from the treeline to the nival zone. We found that: (1) EM fungal community composition, but not richness, changes along elevation, (2) there is no strong evidence of host specialisation, however, EM fungal networks in the alpine zone and within these, EM fungi associated with snowbed communities, are more specialised than in other alpine habitats, (3) plant host population structure does not influence EM fungal communities, and (4) most variability in EM fungal communities is explained by fine-scale changes in edaphic properties, like soil pH and total nitrogen. The higher specialisation and narrower ecological niches of these plant-fungal interactions in snowbed habitats make these habitats particularly vulnerable to environmental change in alpine ecosystems.

GlobalFungi, a global database of fungal occurrences from high-throughput-sequencing metabarcoding studies

Fungi are key players in vital ecosystem services, spanning carbon cycling, decomposition, symbiotic associations with cultivated and wild plants and pathogenicity. The high importance of fungi in ecosystem processes contrasts with the incompleteness of our understanding of the patterns of fungal biogeography and the environmental factors that drive those patterns. To reduce this gap of knowledge, we collected and validated data published on the composition of soil fungal communities in terrestrial environments including soil and plant-associated habitats and made them publicly accessible through a user interface at https://globalfungi.com . The GlobalFungi database contains over 600 million observations of fungal sequences across > 17 000 samples with geographical locations and additional metadata contained in 178 original studies with millions of unique nucleotide sequences (sequence variants) of the fungal internal transcribed spacers (ITS) 1 and 2 representing fungal species and genera. The study represents the most comprehensive atlas of global fungal distribution, and it is framed in such a way that third-party data addition is possible.

Alpine soil microbial ecology in a changing world

Climate change has a disproportionally large impact on alpine soil ecosystems, leading to pronounced changes in soil microbial diversity and function associated with effects on biogeochemical processes at the local and supraregional scales. However, due to restricted accessibility, high-altitude soils remain largely understudied and a considerable heterogeneity hampers the comparability of different alpine studies. Here, we highlight differences and similarities between alpine and arctic ecosystems, and we discuss the impact of climatic variables and associated vegetation and soil properties on microbial ecology. We consider how microbial alpha-diversity, community structures and function change along altitudinal gradients and with other topographic features such as slope aspect. In addition, we focus on alpine permafrost soils, harboring a surprisingly large unknown microbial diversity and on microbial succession along glacier forefield chronosequences constituting the most thoroughly studied alpine habitat. Finally, highlighting experimental approaches, we present climate change studies showing shifts in microbial community structures and function in response to warming and altered moisture, interestingly with some contradiction. Collectively, despite harsh environmental conditions, many specially adapted microorganisms are able to thrive in alpine environments. Their community structures strongly correlate with climatic, vegetation and soil properties and thus closely mirror the complexity and small-scale heterogeneity of alpine soils.

The Soil Microbiome of GLORIA Mountain Summits in the Swiss Alps

While vegetation has intensively been surveyed on mountain summits, limited knowledge exists about the diversity and community structure of soil biota. Here, we study how climatic variables, vegetation, parent material, soil properties, and slope aspect affect the soil microbiome on 10 GLORIA (Global Observation Research Initiative in Alpine environments) mountain summits ranging from the lower alpine to the nival zone in Switzerland. At these summits we sampled soils from all four aspects and examined how the bacterial and fungal communities vary by using Illumina MiSeq sequencing. We found that mountain summit soils contain highly diverse microbial communities with a total of 10,406 bacterial and 6,291 fungal taxa. Bacterial α-diversity increased with increasing soil pH and decreased with increasing elevation, whereas fungal α-diversity did not change significantly. Soil pH was the strongest predictor for microbial β-diversity. Bacterial and fungal community structures exhibited a significant positive relationship with plant communities, indicating that summits with a more distinct plant composition also revealed more distinct microbial communities. The influence of elevation was stronger than aspect on the soil microbiome. Several microbial taxa responded to elevation and soil pH. Chloroflexi and Mucoromycota were significantly more abundant on summits at higher elevations, whereas the relative abundance of Basidiomycota and Agaricomycetes decreased with elevation. Most bacterial OTUs belonging to the phylum Acidobacteria were indicators for siliceous parent material and several OTUs belonging to the phylum Planctomycetes were associated with calcareous soils. The trends for fungi were less clear. Indicator OTUs belonging to the genera Mortierella and Naganishia showed a mixed response to parent material, demonstrating their ubiquitous and opportunistic behaviour in soils. Overall, fungal communities responded weakly to abiotic and biotic factors. In contrast, bacterial communities were strongly influenced by environmental changes suggesting they will be strongly affected by future climate change and associated temperature increase and an upward migration of vegetation. Our results provide the first insights into the soil microbiome of mountain summits in the European Alps that are shaped as a result of highly variable local environmental conditions and may help to predict responses of the soil biota to global climate change.

A single ectomycorrhizal plant root system includes a diverse and spatially structured fungal community

Although only a relatively small proportion of plant species form ectomycorrhizae with fungi, it is crucial for growth and survival for a number of widespread woody plant species. Few studies have attempted to investigate the fine scale spatial structure of entire root systems of adult ectomycorrhizal (EcM) plants. Here, we use the herbaceous perennial Bistorta vivipara to map the entire root system of an adult EcM plant and investigate the spatial structure of its root-associated fungi. All EcM root tips were sampled, mapped and identified using a direct PCR approach and Sanger sequencing of the internal transcribed spacer region. A total of 32.1% of all sampled root tips (739 of 2302) were successfully sequenced and clustered into 41 operational taxonomic units (OTUs). We observed a clear spatial structuring of the root-associated fungi within the root system. Clusters of individual OTUs were observed in the younger parts of the root system, consistent with observations of priority effects in previous studies, but were absent from the older parts of the root system. This may suggest a succession and fragmentation of the root-associated fungi even at a very fine scale, where competition likely comes into play at different successional stages within the root system.

Sequence clustering threshold has little effect on the recovery of microbial community structure

Analysis of microbial community structure by multivariate ordination methods, using data obtained by high-throughput sequencing of amplified markers (i.e., DNA metabarcoding), often requires clustering of DNA sequences into operational taxonomic units (OTUs). Parameters for the clustering procedure tend not to be justified but are set by tradition rather than being based on explicit knowledge. In this study, we explore the extent to which ordination results are affected by variation in parameter settings for the clustering procedure. Amplicon sequence data from nine microbial community studies, representing different sampling designs, spatial scales and ecosystems, were subjected to clustering into OTUs at seven different similarity thresholds (clustering thresholds) ranging from 87% to 99% sequence similarity. The 63 data sets thus obtained were subjected to parallel DCA and GNMDS ordinations. The resulting community structures were highly similar across all clustering thresholds. We explain this pattern by the existence of strong ecological structuring gradients and phylogenetically diverse sets of abundant OTUs that are highly stable across clustering thresholds. Removing low-abundance, rare OTUs had negligible effects on community patterns. Our results indicate that microbial data sets with a clear gradient structure are highly robust to choice of sequence clustering threshold.

Alpine bistort (Bistorta vivipara) in edge habitat associates with fewer but distinct ectomycorrhizal fungal species: a comparative study of three contrasting soil environments in Svalbard

Bistorta vivipara is a widespread arctic-alpine ectomycorrhizal (ECM) plant species. Recent findings suggest that fungal communities associated with B. vivipara roots appear random over short distances, but at larger scales, environmental filtering structure fungal communities. Habitats in highly stressful environments where specialist species with narrower niches may have an advantage represent unique opportunity to test the effect of environmental filtering. We utilised high-throughput amplicon sequencing to identify ECM communities associated with B. vivipara in Svalbard. We compared ECM communities in a core habitat where B. vivipara is frequent (Dryas-heath) with edge habitats representing extremes in terms of nutrient availability where B. vivipara is less frequent (bird-manured meadow and a nutrient-depleted mine tilling). Our analysis revealed that soil conditions in edge habitats favour less diverse but more distinct ECM fungal communities with functional traits adapted to local conditions. ECM richness was overall lower in both edge habitats, and the taxonomic compositions of ECM fungi were in line with our functional expectations. Stress-tolerant genera such as Laccaria and Hebeloma were abundant in nutrient-poor mine site whereas functional competitors genera such as Lactarius and Russula were dominant in the nutrient-rich bird-cliff site. Our results suggest that ECM communities in rare edge habitats are most likely not subsets of the larger pool of ECM fungi found in natural tundra, and they may represent a significant contribution to the overall diversity of ECM fungi in the Arctic.

Diversity and distribution patterns of root-associated fungi on herbaceous plants in alpine meadows of southwestern China

The diversity of root-associated fungi associated with four ectomycorrhizal herbaceous species, Kobresia capillifolia, Carex parva, Polygonum macrophyllum and Potentilla fallens, collected in three sites of alpine meadows in southwestern China, was estimated based on internal transcribed spacer (ITS) rDNA sequence analysis of root tips. Three hundred seventy-seven fungal sequences sorted to 154 operational taxonomical units (sequence similarity of ≥ 97% across the ITS) were obtained from the four plant species across all three sites. Similar taxa (in GenBank with ≥ 97% similarity) were not found in GenBank and/or UNITE for most of the OTUs. Ectomycorrhiz a made up 64% of the fungi operational taxonomic units (OTUs), endophytes constituted 4% and the other 33% were unidentified root-associated fungi. Fungal OTUs were represented by 57% basidiomycetes and 43% ascomycetes. Inocybe, Tomentella/Thelophora, Sebacina, Hebeloma, Pezizomycotina, Cenococcum geophilum complex, Cortinarius, Lactarius and Helotiales were OTU-rich fungal lineages. Across the sites and host species the root-associated fungal communities generally exhibited low host and site specificity but high host and sampling site preference. Collectively our study revealed noteworthy diversity and endemism of root-associated fungi of alpine plants in this global biodiversity hotspot.

Temporal variation of Bistorta vivipara-associated ectomycorrhizal fungal communities in the High Arctic

Ectomycorrhizal (ECM) fungi are important for efficient nutrient uptake of several widespread arctic plant species. Knowledge of temporal variation of ECM fungi, and the relationship of these patterns to environmental variables, is essential to understand energy and nutrient cycling in Arctic ecosystems. We sampled roots of Bistorta vivipara ten times over two years; three times during the growing-season (June, July and September) and twice during winter (November and April) of both years. We found 668 ECM OTUs belonging to 25 different ECM lineages, whereof 157 OTUs persisted throughout all sampling time-points. Overall, ECM fungal richness peaked in winter and species belonging to Cortinarius, Serendipita and Sebacina were more frequent in winter than during summer. Structure of ECM fungal communities was primarily affected by spatial factors. However, after accounting for spatial effects, significant seasonal variation was evident revealing correspondence with seasonal changes in environmental conditions. We demonstrate that arctic ECM richness and community structure differ between summer (growing-season) and winter, possibly due to reduced activity of the core community, and addition of fungi adapted for winter conditions forming a winter-active fungal community. Significant month × year interactions were observed both for fungal richness and community composition, indicating unpredictable between-year variation. Our study indicates that addressing seasonal changes requires replication over several years.

Arctic fungal communities associated with roots of Bistorta vivipara do not respond to the same fine-scale edaphic gradients as the aboveground vegetation

Soil conditions and microclimate are important determinants of the fine-scale distribution of plant species in the Arctic, creating locally heterogeneous vegetation. We hypothesize that root-associated fungal (RAF) communities respond to the same fine-scale environmental gradients as the aboveground vegetation, creating a coherent pattern between aboveground vegetation and RAF. We explored how RAF communities of the ectomycorrhizal (ECM) plant Bistorta vivipara and aboveground vegetation structure of arctic plants were affected by biotic and abiotic variables at 0.3-3.0-m scales. RAF communities were determined using pyrosequencing. Composition and spatial structure of RAF and aboveground vegetation in relation to collected biotic and abiotic variables were analysed by ordination and semi-variance analyses. The vegetation was spatially structured along soil C and N gradients, whereas RAF lacked significant spatial structure. A weak relationship between RAF community composition and the cover of two ECM plants, B. vivipara and S. polaris, was found, and RAF richness increased with host root length and root weight. Results suggest that the fine-scale spatial structure of RAF communities of B. vivipara and the aboveground vegetation are driven by different factors. At fine spatial scales, neighbouring ECM plants may affect RAF community composition, whereas soil nutrients gradients structure the vegetation.

Different bacterial communities in ectomycorrhizae and surrounding soil

Several eukaryotic symbioses have shown to host a rich diversity of prokaryotes that interact with their hosts. Here, we study bacterial communities associated with ectomycorrhizal root systems of Bistorta vivipara compared to bacterial communities in bulk soil using pyrosequencing of 16S rRNA amplicons. A high richness of Operational Taxonomic Units (OTUs) was found in plant roots (3,571 OTUs) and surrounding soil (3,476 OTUs). The community composition differed markedly between these two environments. Actinobacteria, Armatimonadetes, Chloroflexi and OTUs unclassified at phylum level were significantly more abundant in plant roots than in soil. A large proportion of the OTUs, especially those in plant roots, presented low similarity to Sanger 16S rRNA reference sequences, suggesting novel bacterial diversity in ectomycorrhizae. Furthermore, the bacterial communities of the plant roots were spatially structured up to a distance of 60 cm, which may be explained by bacteria using fungal hyphae as a transport vector. The analyzed ectomycorrhizae presents a distinct microbiome, which likely influence the functioning of the plant-fungus symbiosis.