Geospatial variation in co‐occurrence networks of nitrifying microbial guilds

Microbial communities transform nitrogen (N) compounds, thereby regulating the availability of N in soil. The N cycle is defined by interacting microbial functional groups, as inorganic N‐products formed in one process are the substrate in one or several other processes. The nitrification pathway is often a two‐step process in which bacterial or archaeal communities oxidize ammonia to nitrite, and bacterial communities further oxidize nitrite to nitrate. Little is known about the significance of interactions between ammonia‐oxidizing bacteria (AOB) and archaea (AOA) and nitrite‐oxidizing bacterial communities (NOB) in determining the spatial variation of overall nitrifier community structure. We hypothesize that nonrandom associations exist between different AO and NOB lineages that, along with edaphic factors, shape field‐scale spatial patterns of nitrifying communities. To address this, we sequenced and quantified the abundance of AOA, AOB, and Nitrospira and Nitrobacter NOB communities across a 44‐hectare site with agricultural fields. The abundance of Nitrobacter communities was significantly associated only with AOB abundance, while that of Nitrospira was correlated to AOA. Network analysis and geostatistical modelling revealed distinct modules of co‐occurring AO and NOB groups occupying disparate areas, with each module dominated by different lineages and associated with different edaphic factors. Local communities were characterized by a high proportion of module‐connecting versus module‐hub nodes, indicating that nitrifier assemblages in these soils are shaped by fluctuating conditions. Overall, our results demonstrate the utility of network analysis in accounting for potential biotic interactions that define the niche space of nitrifying communities at scales compatible to soil management.