Distribution and Assembly Processes of Soil Fungal Communities along an Altitudinal Gradient in Tibetan Plateau

In soil ecosystems, fungi exhibit diverse biodiversity and play an essential role in soil biogeochemical cycling. Fungal diversity and assembly processes across soil strata along altitudinal gradients are still unclear. In this study, we investigated the structure and abundance of soil fungal communities among soil strata and elevational gradients on the Tibetan Plateau using Illumina MiSeq sequencing of internal transcribed spacer1 (ITS1). The contribution of neutral and niche ecological processes were quantified using a neutral community model and a null model-based methodology. Our results showed that fungal gene abundance increased along altitudinal gradients, while decreasing across soil strata. Along with altitudinal gradients, fungal α-diversity (richness) decreased from surface to deeper soil layers, while β-diversity showed weak correlations with elevations. The neutral community model showed an excellent fit for neutral processes and the lowest migration rate (R(2) = 0.75). The null model showed that stochastic processes dominate in all samples (95.55%), dispersal limitations were dominated at the surface layer and decreased significantly with soil strata, while undominated processes (ecological drift) show a contrary trend. The log-normal model and the null model (βNTI) correlation analysis also neglect the role of niche-based processes. We conclude that stochastic dispersal limitations, together with ecological drifts, drive fungal communities.