Human activity has more than doubled the amount of nitrogen entering the global nitrogen cycle, and the boreal forest biome is a nitrogen-limited ecosystem sensitive to nitrogen load perturbation. Although bryophyte-associated microbes contribute significantly to boreal forest ecosystem function, particularly in carbon and nitrogen cycling, little is known about their responses to anthropogenic global change. Amplicon pyrosequencing of the ITS2 region of rDNA was used to investigate how fungal communities associated with three bryophyte species responded to increased nitrogen loads in a long-term fertilization experiment in a boreal Picea abies forest in southern Norway. Overall, OTU richness, community composition, and the relative abundance of specific ecological guilds were primarily influenced by host species identity and tissue-type. Although not the primary factor affecting fungal communities, nitrogen addition did impact the abundance of specific guilds of fungi and the resulting overall community composition. Increased nitrogen loads decreased ectomycorrhizal abundance, with Amphinema, Cortinarius, Russula, and Tylospora OTUs responding negatively to fertilization. Pathogen abundance increased with fertilization, particularly in the moss pathogen Eocronartium. Saprophytic fungi were both positively and negatively impacted by the nitrogen addition, indicating a complex community level response. The overshadowing of the effects of increased nitrogen loads by variation related to host and tissue-type highlights the complexity of bryophyte-associated microbial communities and the intricate nature of their responses to anthropogenic global change.

Although bryophytes are a dominant vegetation component of boreal and alpine ecosystems, little is known about their associated fungal communities. HPLC assays of ergosterol (fungal biomass) and amplicon pyrosequencing of the ITS2 region of rDNA were used to investigate how the fungal communities associated with four bryophyte species changed across an elevational gradient transitioning from conifer forest to the low alpine. Fungal biomass and OTU richness associated with the four moss hosts did not vary significantly across the gradient (p>0.05), and both were more strongly affected by host and tissue-type. Despite largely constant levels of fungal biomass, distinct shifts in community composition of fungi associated with Hylocomium, Pleurozium, and Polytrichum occured between the elevation zones of the gradient. This likely is a result of influence on fungal communities by major environmental factors such as temperature, directly or indirectly mediated by, or interacting with, the response of other components of the vegetation (i.e. the dominant trees). Fungal communities associated with Dicranum were an exception, exhibiting spatial autocorrelation between plots, and no significant structuring by elevation. Nevertheless, the detection of distinct fungal assemblages associated with a single host growing in different elevation zones along an elevational gradient is of particular relevance in light of the ongoing changes in vegetation patterns in boreal and alpine systems due to global climate warming.