Community dynamics of soil-borne fungal communities along elevation gradients in neotropical and palaeotropical forests.

Because of their steep gradients in abiotic and biotic factors, mountains offer an ideal setting to illuminate the mechanisms that underlie patterns of species distributions and community assembly. We compared the composition of taxonomically and functionally diverse fungal communities in soils along five elevational gradients in mountains of the Neo- and Palaeotropics (northern Argentina, southern Brazil, Panama, Malaysian Borneo and Papua New Guinea). Both the richness and composition of soil fungal communities reflect environmental factors, particularly temperature and soil pH, with some shared patterns among neotropical and palaeotropical regions. Community dynamics are characterized by replacement of species along elevation gradients, implying a relatively narrow elevation range for most fungi, which appears to be driven by contrasting environmental preferences among both functional and taxonomic groups. For functional groups dependent on symbioses with plants (especially ectomycorrhizal fungi), the distribution of host plants drives richness and community composition, resulting in important differences in elevational patterns between neotropical and palaeotropical montane communities. The pronounced compositional and functional turnover along elevation gradients implies that tropical montane forest fungi will be sensitive to climate change, resulting in shifts in composition and functionality over time.

The Compositional Turnover of Grapevine-Associated Plant Pathogenic Fungal Communities Is Greater Among Intraindividual Microhabitats and Terroirs than Among Healthy and Esca-Diseased Plants.

Grapevine trunk diseases (GTDs) are a major threat to the wine industry, causing yield loss and dieback of grapevines. While the increasing damage caused by GTDs in recent decades have spurred several studies on grapevine-associated pathogenic fungi, key questions about the emergence and severity of GTDs remain unanswered, including possible differences in plant pathogenic fungal communities in asymptomatic and symptomatic grapevines. We generated fungal DNA metabarcoding data from soil, bark, and perennial wood samples from asymptomatic and symptomatic grapevines sampled in three terroirs. We observed larger compositional differences in plant pathogenic fungi among different plants parts within grapevine plants than among individual grapevines. This is driven by the dominance of GTD-associated fungi in perennial wood and non-GTD pathogens in soil, as well as by the lack of significant differences among asymptomatic and Esca symptomatic grapevines. These results suggest that fungi generally associated with Esca disease belong to the core grapevine microbiome and likely are commensal endophytes and/or latent saprotrophs, some of which can act as opportunistic pathogens on stressed plants. In addition, we found significant compositional differences among sampling sites, particularly in soil, which suggest a certain influence of local edaphic and mesoclimatic factors on plant pathogenic fungal communities. Furthermore, the observed differences among terroirs in plant pathogenic fungal communities in grapevine woody parts indicate that environmental factors likely are important for the development of Esca disease and further studies are needed to investigate the abiotic conditions on fungal compositional dynamics in Esca-affected plants.

Botrytis cinerea expression profile and metabolism differs between noble and grey rot of grapes.

Botrytis cinerea, the pathogen causing grey rot (GR) with important economic losses in fruit crops, can also cause noble rot (NR) of grape berries under certain environmental conditions, leading to metabolic and physical changes necessary for producing highly regarded botrytized wines. The functional genes involved in biochemical processes in these harmful vs. beneficial berry rot types are still scarcely understood. We generated and analyzed transcriptomic data from healthy (H), NR and GR grape berries collected in the Tokaj wine region in Hungary. Our study shows that B. cinerea is most active in NR, followed by GR and H berries. In addition, expression profiles differed qualitatively between NR and GR, and to a smaller extent between months. Several functional genes expressed during NR were linked to well-known physico-chemical changes in botrytized grape berries, including berry skin degradation and the formation of metabolites favorable for botrytized wine production. In addition, we found that B. cinerea appeared to express genes involved in the biosynthesis of antimicrobials during NR, but not in GR, which likely contributes to the dominance of this fungus during NR.

Evolution of endemism on a young tropical mountain.

Tropical mountains are hot spots of biodiversity and endemism, but the evolutionary origins of their unique biotas are poorly understood. In varying degrees, local and regional extinction, long-distance colonization, and local recruitment may all contribute to the exceptional character of these communities. Also, it is debated whether mountain endemics mostly originate from local lowland taxa, or from lineages that reach the mountain by long-range dispersal from cool localities elsewhere. Here we investigate the evolutionary routes to endemism by sampling an entire tropical mountain biota on the 4,095-metre-high Mount Kinabalu in Sabah, East Malaysia. We discover that most of its unique biodiversity is younger than the mountain itself (6 million years), and comprises a mix of immigrant pre-adapted lineages and descendants from local lowland ancestors, although substantial shifts from lower to higher vegetation zones in this latter group were rare. These insights could improve forecasts of the likelihood of extinction and 'evolutionary rescue' in montane biodiversity hot spots under climate change scenarios.

Can arbuscular mycorrhizal fungi from non-invaded montane ecosystems facilitate the growth of alien trees?

It is generally assumed that recruitment and expansion of alien species along elevation gradients are constrained by climate. But, if plants are not fully constrained by climate, their expansion could be facilitated or hindered by other factors such as biotic interactions. Here, we assessed the composition of arbuscular mycorrhizal fungi (AMF) in soils along an elevation gradient (i.e. 900 m, 1600 m, 2200 m and 2700 m a.s.l.) through a fungal DNA meta-barcoding approach. In addition, we studied in the greenhouse the effects of AMF on growth and phosphorous (P) nutrition of seedlings of the alien trees Gleditsia triacanthos, Ligustrum lucidum and Pyracantha angustifolia cultivated in soils from those elevations, spanning the elevation at which they already form monospecific stands (below 1450 m a.s.l.) and higher elevations, above their current range of distribution in montane ecosystems of Central Argentina. For comparison, we also included in the experiment the dominant native tree Lithraea molleoides that historically occurs below 1300 m a.s.l. Arbuscular mycorrhizal fungal community composition showed strong community turnover with increasing elevation. The effects of these AMF communities on plant growth and nutrition differed among native and alien trees. While P nutrition in alien species' seedlings was generally enhanced by AMF along the whole gradient, the native species benefited only from AMF that occur in soils from the elevation corresponding to its current altitudinal range of distribution. These results suggest that AMF might foster upper range expansion of these invasive trees over non-invaded higher elevations.

Soil lead pollution modifies the structure of arbuscular mycorrhizal fungal communities.

The impact of lead (Pb) pollution on native communities of arbuscular mycorrhizal fungi (AMF) was assessed in soil samples from the surroundings of an abandoned Pb smelting factory. To consider the influence of host identity, bulk soil surrounding plant roots soil samples of predominant plant species (Sorghum halepense, Bidens pilosa, and Tagetes minuta) growing in Pb-polluted soils and in an uncontaminated site were selected. Molecular diversity was assessed by sequencing the 18S rDNA region with primers specific to AMF (AMV4.5NF/AMDGR) using Illumina MiSeq. A total of 115 virtual taxa (VT) of AMF were identified in this survey. Plant species did not affect AMF diversity patterns. However, soil Pb content was negatively correlated with VT richness per sample. Paraglomeraceae and Glomeraceae were the predominant families while Acaulosporaceae, Ambisporaceae, Archaeosporaceae, Claroideoglomeraceae, Diversisporaceae, and Gigasporaceae were less abundant. Acaulosporaceae and Glomeraceae were negatively affected by soil Pb, but Paraglomeraceae relative abundance increased under increasing soil Pb content. Overall, 26 indicator taxa were identified; four of them were previously reported in Pb-polluted soils (VT060; VT222; VT004; VT380); and five corresponded to cultured spores of Scutellospora castaneae (VT041), Diversispora spp. and Tricispora nevadensis (VT060), Diversispora epigaea (VT061), Glomus proliferum (VT099), and Gl. indicum (VT222). Even though AMF were present in Pb-polluted soils, community structure was strongly altered via the differential responses of taxonomic groups of AMF to Pb pollution. These taxon-specific differences in tolerance to soil Pb content should be considered for future phytoremediation strategies based on the selection and utilization of native Glomeromycota.

Vegetation on mesic loamy and sandy soils along a 1700-km maritime Eurasia Arctic Transect.

QUESTIONS: How do plant communities on zonal loamy vs. sandy soils vary across the full maritime Arctic bioclimate gradient? How are plant communities of these areas related to existing vegetation units of the European Vegetation Classification? What are the main environmental factors controlling transitions of vegetation along the bioclimate gradient? LOCATION: 1700-km Eurasia Arctic Transect (EAT), Yamal Peninsula and Franz Josef Land (FJL), Russia. METHODS: The Braun-Blanquet approach was used to sample mesic loamy and sandy plots on 14 total study sites at six locations, one in each of the five Arctic bioclimate subzones and the forest-tundra transition. Trends in soil factors, cover of plant growth forms (PGFs) and species diversity were examined along the summer warmth index (SWI) gradient and on loamy and sandy soils. Classification and ordination were used to group the plots and to test relationships between vegetation and environmental factors. RESULTS: Clear, mostly non-linear, trends occurred for soil factors, vegetation structure and species diversity along the climate gradient. Cluster analysis revealed seven groups with clear relationships to subzone and soil texture. Clusters at the ends of the bioclimate gradient (forest-tundra and polar desert) had many highly diagnostic taxa, whereas clusters from the Yamal Peninsula had only a few. Axis 1 of a DCA was strongly correlated with latitude and summer warmth; Axis 2 was strongly correlated with soil moisture, percentage sand and landscape age. CONCLUSIONS: Summer temperature and soil texture have clear effects on tundra canopy structure and species composition, with consequences for ecosystem properties. Each layer of the plant canopy has a distinct region of peak abundance along the bioclimate gradient. The major vegetation types are weakly aligned with described classes of the European Vegetation Checklist, indicating a continuous floristic gradient rather than distinct subzone regions. The study provides ground-based vegetation data for satellite-based interpretations of the western maritime Eurasian Arctic, and the first vegetation data from Hayes Island, Franz Josef Land, which is strongly separated geographically and floristically from the rest of the gradient and most susceptible to on-going climate change.

Changes in richness and community composition of ectomycorrhizal fungi among altitudinal vegetation types on Mount Kinabalu in Borneo.

The distribution patterns of tropical ectomycorrhizal (ECM) fungi along altitudinal gradients remain largely unknown. Furthermore, despite being an iconic site for biodiversity research, virtually nothing is known about the diversity and spatial patterns of fungi on Mt Kinabalu and neighbouring mountain ranges. We carried out deep DNA sequencing of soil samples collected between 425 and 4000 m above sea level to compare richness and community composition of ECM fungi among altitudinal forest types in Borneo. In addition, we tested whether the observed patterns are driven by habitat or by geometric effect of overlapping ranges of species (mid-domain effect). Community composition of ECM fungi was strongly correlated with elevation. In most genera, richness peaked in the mid-elevation montane forest zone, with the exception of tomentelloid fungi, which showed monotonal decrease in richness with increasing altitude. Richness in lower-mid- and mid-elevations was significantly greater than predicted under the mid-domain effect model. We provide the first insight into the composition of ECM fungal communities and their strong altitudinal turnover in Borneo. The high richness and restricted distribution of many ECM fungi in the montane forests suggest that mid-elevation peak richness is primarily driven by environmental characteristics of this habitat and not by the mid-domain effect.

Abrupt changes in the composition and function of fungal communities along an environmental gradient in the high Arctic.

Fungi play a key role in soil-plant interactions, nutrient cycling and carbon flow and are essential for the functioning of arctic terrestrial ecosystems. Some studies have shown that the composition of fungal communities is highly sensitive to variations in environmental conditions, but little is known about how the conditions control the role of fungal communities (i.e., their ecosystem function). We used DNA metabarcoding to compare taxonomic and functional composition of fungal communities along a gradient of environmental severity in Northeast Greenland. We analysed soil samples from fell fields, heaths and snowbeds, three habitats with very contrasting abiotic conditions. We also assessed within-habitat differences by comparing three widespread microhabitats (patches with high cover of Dryas, Salix, or bare soil). The data suggest that, along the sampled mesotopographic gradient, the greatest differences in both fungal richness and community composition are observed amongst habitats, while the effect of microhabitat is weaker, although still significant. Furthermore, we found that richness and community composition of fungi are shaped primarily by abiotic factors and to a lesser, though still significant extent, by floristic composition. Along this mesotopographic gradient, environmental severity is strongly correlated with richness in all fungal functional groups: positively in saprotrophic, pathogenic and lichenised fungi, and negatively in ectomycorrhizal and root endophytic fungi. Our results suggest complex interactions amongst functional groups, possibly due to nutrient limitation or competitive exclusion, with potential implications on soil carbon stocks. These findings are important in the light of the environmental changes predicted for the Arctic.

Long-term increase in snow depth leads to compositional changes in arctic ectomycorrhizal fungal communities.

Many arctic ecological processes are regulated by soil temperature that is tightly interconnected with snow cover distribution and persistence. Recently, various climate-induced changes have been observed in arctic tundra ecosystems, e.g. shrub expansion, resulting in reduction in albedo and greater C fixation in aboveground vegetation as well as increased rates of soil C mobilization by microbes. Importantly, the net effects of these shifts are unknown, in part because our understanding of belowground processes is limited. Here, we focus on the effects of increased snow depth, and as a consequence, increased winter soil temperature on ectomycorrhizal (ECM) fungal communities in dry and moist tundra. We analyzed deep DNA sequence data from soil samples taken at a long-term snow fence experiment in Northern Alaska. Our results indicate that, in contrast with previously observed responses of plants to increased snow depth at the same experimental site, the ECM fungal community of the dry tundra was more affected by deeper snow than the moist tundra community. In the dry tundra, both community richness and composition were significantly altered while in the moist tundra, only community composition changed significantly while richness did not. We observed a decrease in richness of Tomentella, Inocybe and other taxa adapted to scavenge the soil for labile N forms. On the other hand, richness of Cortinarius, and species with the ability to scavenge the soil for recalcitrant N forms, did not change. We further link ECM fungal traits with C soil pools. If future warmer atmospheric conditions lead to greater winter snow fall, changes in the ECM fungal community will likely influence C emissions and C fixation through altering N plant availability, fungal biomass and soil-plant C-N dynamics, ultimately determining important future interactions between the tundra biosphere and atmosphere.

Changes in composition and abundance of functional groups of arctic fungi in response to long-term summer warming.

We characterized fungal communities in dry and moist tundra and investigated the effect of long-term experimental summer warming on three aspects of functional groups of arctic fungi: richness, community composition and species abundance. Warming had profound effects on community composition, abundance, and, to a lesser extent, on richness of fungal functional groups. In addition, our data show that even within functional groups, the direction and extent of response to warming tend to be species-specific and we recommend that studies on fungal communities and their roles in nutrient cycling take into account species-level responses.

Long-term experimental warming alters community composition of ascomycetes in Alaskan moist and dry arctic tundra.

Arctic tundra regions have been responding to global warming with visible changes in plant community composition, including expansion of shrubs and declines in lichens and bryophytes. Even though it is well known that the majority of arctic plants are associated with their symbiotic fungi, how fungal community composition will be different with climate warming remains largely unknown. In this study, we addressed the effects of long-term (18 years) experimental warming on the community composition and taxonomic richness of soil ascomycetes in dry and moist tundra types. Using deep Ion Torrent sequencing, we quantified how OTU assemblage and richness of different orders of Ascomycota changed in response to summer warming. Experimental warming significantly altered ascomycete communities with stronger responses observed in the moist tundra compared with dry tundra. The proportion of several lichenized and moss-associated fungi decreased with warming, while the proportion of several plant and insect pathogens and saprotrophic species was higher in the warming treatment. The observed alterations in both taxonomic and ecological groups of ascomycetes are discussed in relation to previously reported warming-induced shifts in arctic plant communities, including decline in lichens and bryophytes and increase in coverage and biomass of shrubs.

Summer temperature increase has distinct effects on the ectomycorrhizal fungal communities of moist tussock and dry tundra in Arctic Alaska.

Arctic regions are experiencing the greatest rates of climate warming on the planet and marked changes have already been observed in terrestrial arctic ecosystems. While most studies have focused on the effects of warming on arctic vegetation and nutrient cycling, little is known about how belowground communities, such as fungi root-associated, respond to warming. Here, we investigate how long-term summer warming affects ectomycorrhizal (ECM) fungal communities. We used Ion Torrent sequencing of the rDNA internal transcribed spacer 2 (ITS2) region to compare ECM fungal communities in plots with and without long-term experimental warming in both dry and moist tussock tundra. Cortinarius was the most OTU-rich genus in the moist tundra, while the most diverse genus in the dry tundra was Tomentella. On the diversity level, in the moist tundra we found significant differences in community composition, and a sharp decrease in the richness of ECM fungi due to warming. On the functional level, our results indicate that warming induces shifts in the extramatrical properties of the communities, where the species with medium-distance exploration type seem to be favored with potential implications for the mobilization of different nutrient pools in the soil. In the dry tundra, neither community richness nor community composition was significantly altered by warming, similar to what had been observed in ECM host plants. There was, however, a marginally significant increase in OTUs identified as ECM fungi with the medium-distance exploration type in the warmed plots. Linking our findings of decreasing richness with previous results of increasing ECM fungal biomass suggests that certain ECM species are favored by warming and may become more abundant, while many other species may go locally extinct due to direct or indirect effects of warming. Such compositional shifts in the community might affect nutrient cycling and soil organic C storage.

Greenhouse seedlings of Alnus showed low host intrageneric specificity and a strong preference for some Tomentella ectomycorrhizal associates.

Ectomycorrhizal (ECM) fungal associates of Alnus are relatively few in comparison with those associated with other tree hosts. The composition of ECM assemblages associated with Alnus seems to change very little across the Northern Hemisphere. However, Alnus-associated ECM assemblages from the Western United States, Mexico, and Argentina tend to differ from those in eastern North America and Europe, presumably due to their different biogeographic histories. Alnus glutinosa is a northern European species subjected to diverse environmental conditions. To address intrageneric host preference within two distantly related Alnus species (Alnus acuminata and A. glutinosa), we tested the ECM colonization on seedlings of both species inoculated with natural soil from A. acuminata forests. Two tomentelloid ECM fungi from A. acuminata natural soils were determined from the anatomotyping and molecular analysis. Both species colonized A. glutinosa seedlings and presented similar relative abundances. Additional soil sequence data from A. acuminata sites suggest that a variety of tomentelloid taxa occur, including several unidentified Tomentella lineages. Maximum-likelihood and Bayesian inference analyses based on internal transcribed spacer (ITS) sequences from various locations do not reflect associations of taxa based on their biogeographic origin, and clades are in general constituted by sequences from diverse regions, including South America, Mexico, USA, and Europe. Results illustrate the probable role of specific tomentelloid fungi in the early colonization of seedlings in A. acuminata forests as well as their importance in the structure of the ECM propagule community at the sites.

Large-scale fungal diversity assessment in the Andean Yungas forests reveals strong community turnover among forest types along an altitudinal gradient.

The Yungas, a system of tropical and subtropical montane forests on the eastern slopes of the Andes, are extremely diverse and severely threatened by anthropogenic pressure and climate change. Previous mycological works focused on macrofungi (e.g. agarics, polypores) and mycorrhizae in Alnus acuminata forests, while fungal diversity in other parts of the Yungas has remained mostly unexplored. We carried out Ion Torrent sequencing of ITS2 rDNA from soil samples taken at 24 sites along the entire latitudinal extent of the Yungas in Argentina. The sampled sites represent the three altitudinal forest types: the piedmont (400-700 m a.s.l.), montane (700-1500 m a.s.l.) and montane cloud (1500-3000 m a.s.l.) forests. The deep sequence data presented here (i.e. 4 108 126 quality-filtered sequences) indicate that fungal community composition correlates most strongly with elevation, with many fungi showing preference for a certain altitudinal forest type. For example, ectomycorrhizal and root endophytic fungi were most diverse in the montane cloud forests, particularly at sites dominated by Alnus acuminata, while the diversity values of various saprobic groups were highest at lower elevations. Despite the strong altitudinal community turnover, fungal diversity was comparable across the different zonal forest types. Besides elevation, soil pH, N, P, and organic matter contents correlated with fungal community structure as well, although most of these variables were co-correlated with elevation. Our data provide an unprecedented insight into the high diversity and spatial distribution of fungi in the Yungas forests.

Phylogeny and phylogeography of the Tuber brumale aggr.

The Tuber brumale (winter truffle) is a black truffle reported from most European countries, belonging to the Melanosporum group. Its significance in the economy is ambivalent as the winter truffle has been shown to be a frequent contaminant species in the orchards of the Perigord truffle and occasionally in those of the summer truffle, yet owing to its delicate fragrance, its trade is worthy of note. The phylogeny and phylogeography of economically important truffles are relatively well-explored; however, no thorough research has been published on these aspects of the winter truffle. Therefore, here, we report the first phylogeographic analyses based on samples representing the entire distribution of the species. ITS sequences were used in this survey for haplotype and coalescent analyses, while phylogenetic analyses were based on the ITS, LSU and PKC loci. According to all loci, the samples clustered into two big clades imply the existence of two phylogenetic species. Based on our results, one of these appears to be endemic to the Carpathian Basin. In the other more widespread species, two main phylogeographic groups can be distinguished that show east-west separation with a zone of overlap in the Carpathian Basin, suggesting that they survived the latest glacial period in separate refugia.

Botrytis cinerea causes different plant responses in grape (Vitis vinifera) berries during noble and grey rot: diverse metabolism versus simple defence.

The complexity of the interaction between the necrotrophic pathogen Botrytis cinerea and grape berries (Vitis vinifera spp.) can result in the formation of either the preferred noble rot (NR) or the loss-making grey rot (GR), depending on the prevailing climatic conditions. In this study, we focus on the functional gene set of V. vinifera by performing multidimensional scaling followed by differential expression and enrichment analyses. The aim of this study is to identify the differences in gene expression between grape berries in the phases of grey rot, noble rot, and developing rot (DR, in its early stages) phases. The grapevine transcriptome at the NR phase was found to exhibit significant differences from that at the DR and GR stages, which displayed strong similarities. Similarly, several plant defence-related pathways, including plant-pathogen interactions as hypersensitive plant responses were found to be enriched. The results of the analyses identified a potential plant stress response pathway (SGT1 activated hypersensitive response) that was found to be upregulated in the GR berry but downregulated in the NR berry. The study revealed a decrease in defence-related in V. vinifera genes during the NR stages, with a high degree of variability in functions, particularly in enriched pathways. This indicates that the plant is not actively defending itself against Botrytis cinerea, which is otherwise present on its surface with high biomass. This discrepancy underscores the notion that during the NR phase, the grapevine and the pathogenic fungi interact in a state of equilibrium. Conversely the initial stages of botrytis infection manifest as a virulent fungus-plant interaction, irrespective of whether the outcome is grey or noble rot.

A modified CTAB method for the extraction of high-quality RNA from mono-and dicotyledonous plants rich in secondary metabolites.

BACKGROUND: High-quality RNA extraction from woody plants is difficult because of the presence of polysaccharides and polyphenolics that bind or co-precipitate with the RNA. The CTAB (cetyl trimethylammonium bromide) based method is widely used for the isolation of nucleic acids from polysaccharide-rich plants. Despite the widespread use of the CTAB method, it is necessary to adapt it to particular plant species, tissues and organs. Here we described a simple and generalized method for RNA isolation from mature leaf tissues of several economically important woody (17) and herbaceous plants (2) rich in secondary metabolites. High yields were achieved from small amount (up to 50 mg) of plant material. Two main modifications were applied to the basic protocol: an increase in ß-mercaptoethanol concentration (to 10%v/v) and the use of an effective DNase treatment. As opposed to similar studies, we tried to describe a more detailed protocol for isolating RNA, including the exact quantity and concentration of the reagents were used. RESULTS: Our modified CTAB method is proved to be efficient in extracting the total RNA from a broad range of woody and herbaceous species. The RNA yield was ranged from 2.37 to 91.33 µg/µl. The A260:A280 and A260:A230 absorbance ratios were measured from 1.77 to 2.13 and from 1.81 to 2.22. The RIN value (RNA Integrity Number) of the samples fell between 7.1 and 8.1, which indicated that a small degree of RNA degradation occurred during extraction. The presence of a single peak in the melt curve analyses and low standard errors of the Ct values of replicated measurements indicated the specificity of the primers to bind to the cDNA. CONCLUSIONS: Our RNA isolation method, with fine-tuned and detailed instructions, can produce high quality RNA from a small amount of starting plant material that is suitable for use in downstream transcriptional analyses. The use of an increased concentration of the reducing agent ß-mercaptoethanol in the extraction buffer, as well as the application of DNaseI-treatment resulted in a method suitable for a wide range of plants without the need of further optimalization, especially in Rhus typhina (Staghorn sumac), for which molecular-genetic studies have not yet been sufficiently explored.

Fungal palaeodiversity revealed using high-throughput metabarcoding of ancient DNA from arctic permafrost.

The taxonomic and ecological diversity of ancient fungal communities was assessed by combining next generation sequencing and metabarcoding of DNA preserved in permafrost. Twenty-six sediment samples dated 16 000-32 000 radiocarbon years old from two localities in Siberia were analysed for fungal ITS. We detected 75 fungal OTUs from 21 orders representing three phyla, although rarefaction analyses suggested that the full diversity was not recovered despite generating an average of 6677 ± 3811 (mean ± SD) sequences per sample and that preservation bias likely has considerable effect on the recovered DNA. Most OTUs (75.4%) represented ascomycetes. Due to insufficient sequencing depth, DNA degradation and putative preservation biases in our samples, the recovered taxa probably do not represent the complete historic fungal community, and it is difficult to determine whether the fungal communities varied geographically or experienced a composition shift within the period of 16 000-32 000 bp. However, annotation of OTUs to functional ecological groups provided a wealth of information on the historic communities. About one-third of the OTUs are presumed plant-associates (pathogens, saprotrophs and endophytes) typical of graminoid- and forb-rich habitats. We also detected putative insect pathogens, coprophiles and keratinophiles likely associated with ancient insect and herbivore faunas. The detection of putative insect pathogens, mycoparasites, aquatic fungi and endophytes broadens our previous knowledge of the diversity of fungi present in Beringian palaeoecosystems. A large group of putatively psychrophilic/psychrotolerant fungi was also detected, most likely representing a modern, metabolically active fungal community.