Statistical analyses
All statistical analyses were performed in R (v.3.4.3; R Core Team 2018). Phyloseq package (McMurdie et al. 2013). “Quantitative Insights into Microbial Ecology 2” (QIIME 2 v2018.6) (Caporaso et al. 2010). VSEARCH (Rognes et al. 2016). The rarefied data was used to calculate if the altitude gradient was significantly correlated with the fungi composition. For this purpose, Bray-Curtis distance calculation was applied, and permutation was set to 999. PERMANOVA analysis was calculated using the vegan R package (Oksanen et al. 2015).