The microbial gene diversity along an elevation gradient of the Tibetan grassland

Tibet is one of the most threatened regions by climate warming, thus understanding how its microbial communities function may be of high importance for predicting microbial responses to climate changes. Here, we report a study to profile soil microbial structural genes, which infers functional roles of microbial communities, along four sites/elevations of a Tibetan mountainous grassland, aiming to explore the potential microbial responses to climate changes via a strategy of space-for-time substitution. Using a microarray-based metagenomics tool named GeoChip 4.0, we showed that microbial communities were distinct for most but not all of the sites. Substantial variations were apparent in stress, N and C-cycling genes, but they were in line with the functional roles of these genes. Cold shock genes were more abundant at higher elevations. Also, gdh converting ammonium into urea was more abundant at higher elevations, whereas ureC converting urea into ammonium was less abundant, which was consistent with soil ammonium contents. Significant correlations were observed between N-cycling genes (ureC, gdh and amoA) and nitrous oxide flux, suggesting that they contributed to community metabolism. Lastly, we found by Canonical correspondence analysis, Mantel tests and the similarity tests that soil pH, temperature, NH₄⁺–N and vegetation diversity accounted for the majority (81.4%) of microbial community variations, suggesting that these four attributes were major factors affecting soil microbial communities. On the basis of these observations, we predict that climate changes in the Tibetan grasslands are very likely to change soil microbial community functional structure, with particular impacts on microbial N-cycling genes and consequently microbe-mediated soil N dynamics.