Ecological Succession Pattern of Fungal Community in Soil along a Retreating Glacier

Accelerated by global climate changing, retreating glaciers leave behind soil chronosequences of primary succession. Current knowledge of primary succession is mainly from studies of vegetation dynamics, whereas information about belowground microbes remains unclear. Here, we combined shifts in community assembly processes with microbial primary succession to better understand mechanisms governing the stochastic/deterministic balance. We investigated fungal succession and community assembly via high-throughput sequencing along a well-established glacier forefront chronosequence that spans 2–188 years of deglaciation. Shannon diversity and evenness peaked at a distance of 370 m and declined afterwards. The response of fungal diversity to distance varied in different phyla. Basidiomycota Shannon diversity significantly decreased with distance, while the pattern of Rozellomycota Shannon diversity was unimodal. Abundance of most frequencies OTU2 (Cryptococcus terricola) increased with successional distance, whereas that of OTU65 (Tolypocladium tundrense) decreased. Based on null deviation analyses, composition of the fungal community was initially governed by deterministic processes strongly but later less deterministic processes. Our results revealed that distance, altitude, soil microbial biomass carbon, soil microbial biomass nitrogen and [Formula: see text] –N significantly correlated with fungal community composition along the chronosequence. These results suggest that the drivers of fungal community are dynamics in a glacier chronosequence, that may relate to fungal ecophysiological traits and adaptation in an evolving ecosystem. The information will provide understanding the mechanistic underpinnings of microbial community assembly during ecosystem succession under different scales and scenario.