Abstract
Soil pollution is an important stressor affecting biodiversity and ecosystem functioning. However, we lack a holistic understanding of how soil microbial communities respond to heavy metal pollution in agricultural ecosystems. Here, we explored the distribution patterns and inter-kingdom interactions of entire soil microbiome (including bacteria, fungi, and protists) in 47 paired paddy and upland fields along a gradient of legacy mercury (Hg) pollution. We found that the richness and composition of protistan community had stronger responses to Hg pollution than those of bacterial and fungal communities in both paddy and upland soils. Mercury polluted soils harbored less protistan phototrophs but more protistan consumers. We further revealed that long-term Hg pollution greatly increased network complexity of protistan community than that of bacterial and fungal communities, as well as intensified the interactions between protists and the other microorganisms. Moreover, our results consistently indicated that protistan communities had stronger responses to long-term Hg pollution than bacterial and fungal communities in agricultural soils based on structural equation models and random forest analyses. Our study highlights that soil protists can be used as bioindicators of Hg pollution, with important implications for the assessment of contaminated farmlands and the sustainable management of agricultural ecosystems.