Abundant and rare soil fungi exhibit distinct succession patterns in the forefield of Dongkemadi glacier on the central Qinghai-Tibet Plateau.

Glaciers are retreating rapidly, exposing extensive new soil habitats in glacier forefields and providing unique areas for studying primary succession. However, understanding the variation patterns and assembly mechanisms of abundant and rare fungi subcommunities along the glacier-retreating chronosequence remains a knowledge gap, especially true for the vast Qinghai-Tibet Plateau (QTP). Here, we investigated fungal communities in the glacier forefield in Dongkemadi Glaicer on the central QTP. The results showed that fungal alpha diversity exhibited a clear increasing pattern in response to increasing of distance to glacier. The percentage of abundant OTUs decreased while the percentage of rare OTUs increased, suggesting that soil development is more beneficial to the rare taxa. The distributions of both abundant and rare subcommunities exhibited a clear spatial pattern along the distance to glacier, and might be strongly controlled by multiple environmental variables, including pH, soil moisture, vegetation status, soil organic carbon, total nitrogen, and soluble reactive phosphorus. Abundant and rare fungal subcommunities were structured in different assembly regimes. Dispersal limitation processes were dominant for both abundant and rare subcommunities but with a stronger contribution to abundant subcommunity assembly. Heterogeneous selection processes contributed higher and non-dominant processes contributed lower to abundant subcommunities than to rare subcommunities. The modular structure of the fungal co-occurrence network was highly localized along the soil chronosequence. By revealing distinct diversity patterns and community assembly mechanisms of abundant and rare fungal subcommunities, our study improved our understanding of ecological succession along the glacier-retreating chronosequence.