Spatial Differentiation of Heavy Metals/Metalloids, Microbial Risk Genes and Soil Microbiota in a Sulfur-Contaminated Landscape

Abstract

Legacy sulfur smelting has left behind complex contamination landscapes, yet the spatial structuring of microbial risks and adaptation strategies across soil profiles remains insufficiently understood. Microbial risk genes, including those conferring resistance to antibiotic resistance (ARGs), biocide and metal resistance (BRGs/MRGs), and virulence (VFGs), are increasingly recognized as co-selected under heavy metal stress, posing both ecological and public health concerns. In this study, we integrated geochemical analyses with metagenomic sequencing and functional annotation to jointly characterize the vertical (0–7 m) and horizontal (~2 km) distribution of heavy metals/metalloids, microbial communities, and functional risk genes at a historic smelting site in Zhenxiong, Yunnan. Heavy metals and metalloids such as arsenic (As), chromium (Cr), copper (Cu), and lead (Pb) showed clear accumulation with depth, while significantly lower concentrations were observed in both upstream and downstream locations, revealing persistent vertical and horizontal pollution gradients. Correspondingly, resistance and virulence genes were co-enriched at contaminated sites, suggesting potential co-selection under prolonged stress. LEfSe analysis revealed distinct ecological patterns: vertically, upper layers were dominated by nutrient-cycling and mildly stress-tolerant taxa, while deeper layers favored metal-resistant, oligotrophic, and potentially pathogenic microorganisms; horizontally, beneficial and diverse microbes characterized low-contamination zones, whereas heavily polluted areas were dominated by resistant and stress-adapted genera. These findings provide new insights into microbial resilience and ecological risk under long-term smelting stress.