Abstract

Declines in soil multifunctionality (e.gsoil capacity to provide food and energy) are closely related to changes in the soil microbiome (e.g., diversity) Determining ecological drivers promoting such microbiome changes is critical knowledge for protecting soil functions. However, soil-microbe interactions are highly variable within environmental gradients and may not be consistent across studies. Here we propose that analysis of community dissimilarity (beta-diversity) is a valuable tool for overviewing soil microbiome spatiotemporal changes. Indeed, beta-diversity studies at larger scales (modelling and mapping) simplify complex multivariate interactions and refine our understanding of ecological drivers by also giving the possibility of expanding the environmental scenarios. This study represents the first spatial investigation of beta-diversity in the soil microbiome of New South Wales (800,642 km(2)), Australia. We used metabarcoding soil data (16S rRNA and ITS genes) as exact sequence variants (ASVs) and UMAP (Uniform Manifold Approximation and Projection) as the distance metric. beta-Diversity maps (1000-m resolution)-concordance correlations of 0.91-0.96 and 0.91-0.95 for bacteria and fungi, respectively-showed soil biome dissimilarities driven primarily by soil chemistry-pH and effective cation exchange capacity (ECEC)-and cycles of soil temperature-land surface temperature (LST-phase and LST-amplitude). Regionally, the spatial patterns of microbes parallel the distribution of soil classes (e.g., Vertosols) beyond spatial distances and rainfall, for example. Soil classes can be valuable discriminants for monitoring approaches, for example pedogenons and pedophenons. Ultimately, cultivated soils exhibited lower richness due to declines in rare microbes which might compromise soil functions over time.