pH-dependent distribution of soil ammonia oxidizers across a large geographical scale as revealed by high-throughput pyrosequencing

Ammonia-oxidizing archaea (AOA) and bacteria (AOB) are ubiquitous and important for nitrogen transformations in terrestrial ecosystems. However, the distribution patterns of these microorganisms as affected by the terrestrial environments across a large geographical scale are not well understood. This study was designed to gain insights into the ecological characteristics of AOA and AOB in 65 soils, collected from a wide range of soil and ecosystem types. Barcoded pyrosequencing in combination with quantitative PCR was employed to characterize the relative abundance, diversity, and community composition of archaeal 16S rRNA gene, and AOA and AOB amoA genes in 65 soil samples. The operational taxonomic unit richness and Shannon diversity of Thaumarchaeota, AOA, and AOB were highly variable among different soils, but their variations were best explained by soil pH. Soil pH was strongly correlated with the overall community composition of ammonia oxidizers, as measured by the pairwise Bray-Curtis dissimilarity across all sites. These findings were further corroborated by the evident pH-dependent distribution patterns of four thaumarchaeal groups (I.1a-associated, I.1b, I.1c, and I.1c-associated) and four AOB clusters (2, 3a.1, 10, and 12). The ratios of AOA to AOB amoA gene copy numbers significantly decreased with increasing pH, suggesting a competitive advantage of AOA over AOB in acidic soils. These results suggest that the distribution of ammonia oxidizers across large-scale biogeographical settings can be largely predicted along the soil pH gradient, thus providing important indications for the ecological characteristics of AOA and AOB in different soils.