Maomao Feng

Climate warming has led to glacier retreat worldwide. Studies on the taxonomy and functions of glacier microbiomes help us better predict their response to glacier melting. Here, we used shotgun metagenomic sequencing to study the microbial functional potential in different cryospheric habitats, i.e. surface snow, supraglacial and subglacial sediments, subglacial ice, proglacial stream water and recently deglaciated soils. The functional gene structure varied greatly among habitats, especially for snow, which differed significantly from all other habitats. Differential abundance analysis revealed that genes related to stress responses (e.g. chaperones) were enriched in ice habitat, supporting the fact that glaciers are a harsh environment for microbes. The microbial metabolic capabilities related to carbon and nitrogen cycling vary among cryospheric habitats. Genes related to auxiliary activities were overrepresented in the subglacial sediment, suggesting a higher genetic potential for the degradation of recalcitrant carbon (e.g., lignin). As for nitrogen cycling, genes related to nitrogen fixation were more abundant in barren proglacial soils, possibly due to the presence of Cyanobacteriota in this habitat. Our results deepen our understanding of microbial processes in glacial ecosystems, which are vulnerable to ongoing global warming, and they have implications for downstream ecosystems.

Glacier retreat is a visible consequence of climate change worldwide. Although taxonomic change of the soil microbiomes in glacier forefields have been widely documented, how microbial genetic potential changes along succession is little known. Here, we used shotgun metagenomics to analyse whether the soil microbial genetic potential differed between four stages of soil development (SSD) sampled along three transects in the Damma glacier forefield (Switzerland). The SSDs were characterized by an increasing vegetation cover, from barren soil, to biological soil crust, to sparsely vegetated soil and finally to vegetated soil. Results suggested that SSD significantly influenced microbial genetic potential, with the lowest functional diversity surprisingly occurring in the vegetated soils. Overall, carbohydrate metabolism and secondary metabolite biosynthesis genes overrepresented in vegetated soils, which could be partly attributed to plant–soil feedbacks. For C degradation, glycoside hydrolase genes enriched in vegetated soils, while auxiliary activity and carbohydrate esterases genes overrepresented in barren soils, suggested high labile C degradation potential in vegetated, and high recalcitrant C degradation potential in barren soils. For N-cycling, organic N degradation and synthesis genes dominated along succession, and gene families involved in nitrification were overrepresented in barren soils. Our study provides new insights into how the microbial genetic potential changes during soil formation along the Damma glacier forefield.